配色: 字号:
2010136102张秋红实验三
2012-06-29 | 阅:  转:  |  分享 
  
实验三

一、实验目的

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

















献花(0)
+1
(本文系依米荷阳首藏)