《实验二傅里叶分析及应用.doc》由会员分享，可在线阅读，更多相关《实验二傅里叶分析及应用.doc（8页珍藏版）》请在三一办公上搜索。
1、实验二 傅里叶分析及应用一、实验目的一掌握使用Matlab进展周期信号傅里叶级数展开和频谱分析1、学会使用Matlab分析傅里叶级数展开，深入理解傅里叶级数的物理含义2、学会使用Matlab分析周期信号的频谱特性二掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质1、学会运用Matlab求连续时间信号的傅里叶变换2、学会运用Matlab求连续时间信号的频谱图3、学会运用Matlab分析连续时间信号的傅里叶变换的性质（三） 掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进展分析 2、学会运用MATLAB改变抽样时间间隔，观察抽
2、样后信号的频谱变化 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)
3、,grid ontitle(幅度谱)phase=atan(imag(FW)/real(FW);subplot(212)ezplot(phase),grid ontitle(相位谱)数值计算法：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*e*p(-1i*t*W);F=abs(F);subplot(211)plot(W,F),grid ona*is(-pi pi -1
4、 3)*label(w),ylabel(F(w)title(幅度谱)phase=atan(imag(FW)/real(FW);subplot(212)ezplot(phase),grid ontitle(相位谱)2、试用Matlab命令求的傅里叶反变换，并绘出其时域信号图。代码：t=sym(t);Fw=sym(10/(3+i*w)-4/(5+i*w);ft=ifourier(Fw,t);ezplot(abs(ft),grid ona*is(0 3 -1 7);*label(t),ylabel(F(t)3、门函数自身卷积为三角波信号，试用Matlab命令验证FT的时域卷积定理。代码：f1=sym
5、(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、设有两个不同频
6、率的余弦信号，频率分别为，；现在使用抽样频率对这三个信号进展抽样，使用MATLAB命令画出各抽样信号的波形和频谱，并分析其频率混叠现象建议：抽样信号的频谱图横坐标围不小于-10000Hz10000Hz或-20000*pi20000*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 ona*is(-0.006 0.006 -1.6 1.1);*label(),ylabel(f(t)
7、title(余弦函数图像 )N=5000;k=-N:N;w=2*pi*k/(2*N+1)*dt);Fw=dt*ft*e*p(-1i*t*w);subplot(222)plot(w,abs(Fw)grid ona*is(-10000 10000 0 0.008);*label(w),ylabel(F(w)title(余弦函数的频谱 )t2=-0.006:Ts:0.006; fst=cos(2.*pi.*100*t2); subplot(223) plot(t,ft,:),hold onstem(t2,fst),grid ona*is(-0.006 0.006 -1.5 1.5)*label(),
8、ylabel(Fs(t)title(抽样后的信号),hold offFsw=Ts*fst*e*p(-1i*t2*w);subplot(224) plot(w,abs(Fsw),grid ona*is(-10000 10000 0 0.008) *label(w),ylabel(Fs(w) title(抽样信号的频谱)2当时的抽样函数Ts=0.00025;dt=0.00001;t=-0.006:dt:0.006;a=3800*pi;ft=cos(a*t);subplot(221)plot(t,ft),grid ona*is(-0.006 0.006 -1.6 1.1);*label(),ylab
9、el(f(t)title(余弦函数图像 )N=5000;k=-N:N;w=2*pi*k/(2*N+1)*dt);Fw=dt*ft*e*p(-1i*t*w);subplot(222)plot(w,abs(Fw)grid ona*is(-10000 10000 0 0.008);*label(w),ylabel(F(w)title(余弦函数的频谱 )t2=-0.006:Ts:0.006; fst=cos(2.*pi.*100*t2); subplot(223) plot(t,ft,:),hold onstem(t2,fst),grid ona*is(-0.006 0.006 -1.5 1.5)*l
10、abel(),ylabel(Fs(t)title(抽样后的信号),hold offFsw=Ts*fst*e*p(-1i*t2*w);subplot(224) plot(w,abs(Fsw),grid ona*is(-10000 10000 0 0.008) *label(w),ylabel(Fs(w) title(抽样信号的频谱)5、结合抽样定理，利用MATLAB编程实现信号经过冲激脉冲抽样后得到的抽样信号及其频谱建议：冲激脉冲的周期分别取4*pi/3 s、pi s、2*pi/3 s三种情况比照，并利用构建信号。1当冲激脉冲的周期取4*pi/3 s时代码：wm=2;wc=1.2*wm;Ts=4
11、*pi/3;dt=0.1;ft=sinc(t1/pi).*(uCT(t1+10)-uCT(t1-10);N=500;k=-N:N;w=2*pi*k/(2*N+1)*dt);n=-100:100;nTs=n*Ts;fst=sinc(nTs/pi).*(uCT(nTs+10)-uCT(nTs-10);subplot(221);plot(t1,ft,:),hold onstem(nTs,fst),grid ona*is(-4 4 -1 1);*label(),ylabel(fs(t)title(Sat抽样后的信号),hold off;Fsw=Ts*fst*e*p(-1i*nTs*w);subplot
12、(222)plot(w,abs(Fsw),c),grid ona*is(-10 10 -3 10)*label(w),ylabel(Fs(w)title(Sa(t)抽样信号的频谱)t=-10:dt:10;f=fst*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t);subplot(223)plot(t,f,m),grid on;a*is(-10 10 -2 9);*label(t),ylabel(f(t)title(由f(nTs)信号重建得到Sa(t)信号)2当冲激脉冲的周期取pi时wm=2;wc=1.2*wm;T
13、s=pi;dt=0.1;ft=sinc(t1/pi).*(uCT(t1+10)-uCT(t1-10);N=500;k=-N:N;w=2*pi*k/(2*N+1)*dt);n=-100:100;nTs=n*Ts;fst=sinc(nTs/pi).*(uCT(nTs+10)-uCT(nTs-10);subplot(221);plot(t1,ft,:),hold onstem(nTs,fst),grid ona*is(-4 4 -1 1);*label(),ylabel(fs(t)title(Sat抽样后的信号),hold off;Fsw=Ts*fst*e*p(-1i*nTs*w);subplot(
14、222)plot(w,abs(Fsw),c),grid ona*is(-10 10 -3 10)*label(w),ylabel(Fs(w)title(Sa(t)抽样信号的频谱 )t=-10:dt:10;f=fst*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t);subplot(223)plot(t,f,m),grid on;a*is(-10 10 -2 9);*label(t),ylabel(f(t)title(由 f(nTs)信号重建得到 Sa(t)信号 )3当冲激脉冲的周期取2*pi/3时wm=2;wc=1
15、.2*wm;Ts=2*pi/3;dt=0.1;ft=sinc(t1/pi).*(uCT(t1+10)-uCT(t1-10);N=500;k=-N:N;w=2*pi*k/(2*N+1)*dt);n=-100:100;nTs=n*Ts;fst=sinc(nTs/pi).*(uCT(nTs+10)-uCT(nTs-10);subplot(221);plot(t1,ft,:),hold onstem(nTs,fst),grid ona*is(-4 4 -1 1);*label(),ylabel(fs(t)title(Sat抽样后的信号),hold off;Fsw=Ts*fst*e*p(-1i*nTs*
16、w);subplot(222)plot(w,abs(Fsw),c),grid ona*is(-10 10 -3 10)*label(w),ylabel(Fs(w)title(Sa(t)抽样信号的频谱 )t=-10:dt:10;f=fst*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t);subplot(223)plot(t,f,m),grid on;a*is(-10 10 -2 9);*label(t),ylabel(f(t)title(由 f(nTs)信号重建得到 Sa(t)信号 )6、周期三角信号如下列图所示注
17、：图中时间单位为：毫秒(ms)：1试求出该信号的傅里叶级数自己求或参见课本P112或P394，利用Matlab编程实现其各次谐波如1、3、5、13、49的叠加，并验证其收敛性；代码： t=-10:0.01:10;omega=pi;y=abs(sawtooth(pi*0.5*t,0.5);plot(t,y),grid on;a*is(-10,10,0,1.5);n_ma*=1,3,5,13,49;N=length(n_ma*);for k=1:N n=1:2:n_ma*(k); b=4./(pi*n).2); *=b*cos(omega*n*t); figure; plot(t,y); hold
18、 on; *=*+1/2; plot(t,*); hold off; a*is(-10,10,0,1.5); title(=,num2str(n_ma*(k);end2用Matlab分析该周期三角信号的频谱三角形式或指数形式均可。当周期三角信号的周期如由2ms1ms或由2ms4ms和宽度如2ms1ms分别变化时，试观察分析其频谱的变化。答：其频谱图中间局部的幅值降低，两侧的局部幅值上升。代码：dt=0.01;t=-4:dt:4;ft=(t=-1&t0&t=-0.5&t0&t=0.5).*(1-t);%subplot(2,1,1)%plot(t,ft);grid onn=2000;k=-n:n;
19、w=pi*k/(n*dt);f=dt*ft*e*p(-1i*t*w);f=abs(f);%subplot(2,1,2)plot(w,f);a*is(-20 20 0 1.1),grid on;title(周期三角信号的周期1ms )四、实验结论和讨论本次实验总体难度较大，但是数据上根本没有问题，除了在选择坐标长度时有些问题之外，总体没什么问题。图像清晰完整，结果也比拟明显。五、实验思考本次实验让我更加熟悉了MATLAB的根本用法和一些常用的数学计算函数，在此根底上也让我更加深入的对傅里叶变换，傅里叶级数，抽样函数，抽样定理等加深了认识，在一些细节用法的方面加深了印象。本次实验难度比拟大，但是难度大的同时收获也非常丰富。今后，我要更加熟悉常用函数以及根本规则，争取快速又保质保量完成任务。
