信号与系统实验二傅里叶分析及应用实验报告分析.doc

实验二 傅里叶分析及应用 一、实验目的 （一）掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab分析傅里叶级数展开，深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性 （二）掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔，观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建 二、实验条件 装用Matlab R2015a的电脑。 程序： y=(sawtooth(pi*(t+1),0.5)+1)/2; plot(t,y),xlabel(t),ylabel(三角波信号); axis([-2,2,0,1.1]); n_max=[1,3,5,11,47]; N=length(n_max); for k=1:N n=1:2:n_max(k); b=4./((pi.^2)*(n.^2)); x=0.5+b*cos(pi*n*t); figure,plot(t,x);hold on; plot(t,y);hold off; xlabel(t),ylabel(部分和的波形); axis([-3,3,0,1]); title([最大谐波数=,num2str(n_max(k))]); End 结果： 用Matlab分析该周期三角信号的频谱[三角形式或指数形式均可]。当周期三角信号的周期（如由2ms(1ms或由2ms(4ms）和宽度（如2ms(1ms）分别变化时，试观察分析其频谱的变化。 程序 n=-20:20;T=10;w1=2*pi/T; fn=2./(n.*n.*w1.*w1); subplot(311),stem(n*w1,fn),title(T=10);%axis([-10,10,-0.1,20]); T=5;w1=2*pi/T;fn=2./(n.*n.*w1.*w1); subplot(312),stem(n*w1,fn),title(T=5);axis([-20,20,-0.1,3]); T=1;w1=2*pi/T;fn=2./(n.*n.*w1.*w1); subplot(313),stem(n*w1,fn),title(T=1);axis([-100,100,-0.01,0.1]); 结果： 2、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图（包括幅度谱和相位谱）[注：图中时间单位为：毫秒(ms)]。 程序：1. ft=sym((t+2)*(heaviside(t+2)-heaviside(t+1))+(heaviside(t+1)-heaviside(t-1))+(-t+2)*(heaviside(t-1)-heaviside(t-2))); fw=simplify(fourier(ft)); subplot(311);ezplot(ft,[-3,3]); subplot(312);ezplot(abs(fw),[-10,10,0,4]);title(·ù?è?×); phase=atan(imag(fw)/real(fw)); subplot(313);ezplot(phase);title(?à???×); 程序2： dt=0.01; t=-4:dt:4; ft=(t+2).*(heaviside(t+2)-heaviside(t+1))+(heaviside(t+1)-heaviside(t-1))+(-t+2).*(heaviside(t-1)-heaviside(t-2)); n=2000; k=-n:n; w=pi*k/(n*dt); f=dt*ft*exp(-1j*t*w); f=abs(f); plot(w,f); axis([-10 10 0 3.2]); 的傅里叶反变换，并绘出其时域信号图。 [注意：(1)写代码时j(i] 程序： syms t; fw=sym(10/(3+iw)-4/(5+iw)); ft=simplify(ifourier(fw,t)); ezplot(abs(ft)); ft=(t=-0.5)-(t=0.5); f=conv(ft,ft)*dt;n=length(f);tt=(0:n-1)*dt-10; subplo