实验二傅里叶分析报告及应用.docVIP

实用文案 标准文档 实验二 傅里叶分析及应用 一、实验目的 （一）掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab分析傅里叶级数展开，深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性 （二）掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔，观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建 二、实验条件 装有MATLAB的电脑 三、实验内容 1、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图（包括幅度谱和相位谱）[注：图中时间单位为：毫秒(ms)]。 符号运算法： ft=sym((t+2)*heaviside(t+2)-(t+2)*heaviside(t+1)+heaviside(t+1)-heaviside(t-1)+(2-t)*heaviside(t-1)-(2-t)*heaviside(t-2)); FW=simplify(fourier(ft)) subplot(211) ezplot(abs(FW)),grid on title(幅度谱) phase=atan(imag(FW)/real(FW)); subplot(212) ezplot(phase),grid on title(相位谱) 数值计算法： dt=0.01; t=-4:dt:4; ft=(t+2).*uCT(t+2)-(t+2).*uCT(t+1)+uCT(t+1)-uCT(t-1)+(2-t).*uCT(t-1)-(2-t).*uCT(t-2); N=2000; k=-N:N; W=pi*k/(N*dt); F=dt*ft*exp(-1i*t*W); F=abs(F); subplot(211) plot(W,F),grid on axis([-pi pi -1 3]) xlabel(w),ylabel(F(w)) title(幅度谱) phase=atan(imag(FW)/real(FW)); subplot(212) ezplot(phase),grid on title(相位谱) 2、试用Matlab命令求的傅里叶反变换，并绘出其时域信号图。 代码： t=sym(t); Fw=sym(10/(3+i*w)-4/(5+i*w)); ft=ifourier(Fw,t); ezplot(abs(ft)),grid on axis([0 3 -1 7]); xlabel(t),ylabel(F(t)) 3、已知门函数自身卷积为三角波信号，试用Matlab命令验证FT的时域卷积定理。 代码： f1=sym(heaviside(t+1)-heaviside(t-1)); F=fourier(f1); F=simplify(F)*simplify(F); subplot(2,1,1); ezplot(F); yt=sym((heaviside(t+2)-heaviside(t))*(t+2)+(heaviside(t)-heaviside(t-2))*(-t+2) ); F2=fourier(yt); F2=simplify(F2); subplot(2,1,2); ezplot(F2); 从图中可以很明显的看出时域信号先卷积在进行傅里叶变换与时域信号先进行傅里叶变换后再相乘结果一致，因此验证卷积定理正确。 4、设有两个不同频率的余弦信号，频率分别为，；现在使用抽样频率对这三个信号进行抽样，使用MATLAB命令画出各抽样信号的波形和频谱，并分析其频率混叠现象[建议：抽样信号的频谱图横坐标范围不小于-10000Hz~10000Hz或-20000*pi~20000*pi rad/s]。 答：在低抽样率时有混叠现象而高抽样率时无混叠现象。 代码： （1）当时的抽样函数 Ts=0.00025; dt=0.00001; t=-0.006:dt:0.006; a=200*pi; ft=cos(a*t); subplot(221) plot(t,ft),grid on axis([-0.006 0.006 -1.6 1.1]); xlabel(ê±??),ylabel(f(t)) title(余弦函数图像 ) N=5000; k=-N:N; w=2*pi*k/((2*N+1)*dt);