实验三
一、实验目的
1.练习利用matlab6.5中的命令或程序分析信号频谱。
2.进一步练习M-File的建立、保存、调用。
3.加深对几种傅里叶变换的理解。
二、实验原理与方法
连续时间信号傅里叶变换(FT):
2.离散时间信号傅里叶变换(DTFT):
3.离散傅里叶变换(DFT):
4.快速傅里叶变换(FFT):
5.分析信号x(n)=(-0.5)n,-10
程序:
n=-10:10;
k=-2:0.001:2;w=pik;
x=[(-0.5).^n];
X=xexp(-jn''w);
magX=abs(X);phaX=angle(X);
figure(1);plot(w/pi,magX);
xlabel(''Frequencyinpiunits'');ylabel(''|X(ejw)|'');
title(''D?o?μ?x(n)·ù?èì?D?'');
figure(2);plot(w/pi,phaX/pi);
xlabel(''Frequencyinpiunits'');ylabel(''radians/pi'');
title(''D?o?x(n)μ??à?μì?D?'');
图像:
6.查资料弄清函数fft()的调用格式及其含义。
调用方法
?A=fft(X,N,DIM)?其中,X表示输入图像;N表示采样间隔点,如果X小于该数值,那么Matlab将会对X进行零填充,否则将进行截取,使之长度为N;DIM表示要进行离散傅立叶变换。?x=IFFT(X);
x=IFFT(X,N)
用MATLAB进行谱分析时注意:
(1)函数FFT返回值的数据结构具有对称性。
(2)做FFT分析时,幅值大小与FFT选择的点数有关,但不影响分析结果。在IFFT时已经做了处理。要得到真实的振幅值的大小,只要将得到的变换后结果乘以2除以N即可。
评分人:陈奇分数:9
|
|