民大信号与系统实验四.doc

信号与系统实验报告 学生姓名： xxx 班级学号： xxx 指导老师： xxx 实验报告内容 实验名称：周期信号的傅里叶级数及频谱分析 傅里叶变换及性质 二、实验学时：2学时 三、实验目的： 学会运用Matlab分析傅里叶级数展开，深入理解傅里叶级数的物理意义 学会运用Matlab分析周期信号的频谱特性 学会运用Matlab求连续时间信号的傅里叶变换 学会运用Matlab求连续时间信号的频谱图 学会运用Matlab分析连续时间信号的傅里叶变换性质 四、实验内容： 分析傅里叶级数展开，深入理解傅里叶级数的物理意义：分析周期信号的频谱特性：求连续时间信号的傅里叶变换：Matlab求连续时间信号的频谱图：分析连续时间信号的傅里叶变换性质 五、实验原理： 分析傅里叶级数展开，深入理解傅里叶级数的物理意义：分析周期信号的频谱特性：求连续时间信号的傅里叶变换：Matlab求连续时间信号的频谱图：分析连续时间信号的傅里叶变换性质 六、实验结果 1.已知周期三角信号波形如图所示，试求出该信号的傅里叶级数，利用MATLAB编程实现其各次谐波的叠加，并验证其收敛性。 源程序： t=-6:0.001:6; y=(sawtooth(pi*(t+1),0.5)+1)./2; plot(t,y),grid on xlabel(t),ylabel(周期三角信号波形); axis([-6 6 -0.5 1.5]) n_max=([1 3 5 11 47]); N=length(n_max); for k=1:N n=1:2:n_max(k); b=4./(pi*pi.*n.*n); omega=pi; x=1/2+b*cos(omega*n*t); figure; plot(t,y); hold on plot(t,x); hold off; xlabel(t),ylabel(bufengheboxing) axis([-6 6 -0.5 1.5]),grid on title([max_xieboshu=,num2str(n_max(k))]) end 2.试用MATLAB分析上题中周期三角信号的频谱。当周期三角信号的周期和三角信号的宽度发生变化时，试观察分析其频谱的变化。 源程序： n=-30:30;tao=1;T=10;w1=2*pi/T; fn=T./(pi^2*n.^2.*tao).*(1-cos(pi*n.*tao./T)); subplot(311) stem(n*w1,fn),grid on title(tao=1,T=10) tao=1;T=5;w2=2*pi/T; fn=T./(pi^2*n.^2.*tao).*(1-cos(pi*n.*tao./T)); m=round(30*w1/w2); n1=-m:m; fn=fn(30-m+1:30+m+1); subplot(312) stem(n1*w2,fn),grid on title(tao=1,T=5) tao=2;T=10;w3=2*pi/T; fn=T./(pi^2*n.^2.*tao).*(1-cos(pi*n.*tao./T)); subplot(313) stem(n*w3,fn),grid on title(tao=2,T=10) 7.3编程练习 1(2).f(t)=(sin(π*t)/ π*t)^2; syms t ft=sym((sin(pi*t)/(pi*t))^2); Fw=fourier(ft); subplot(211); ezplot(abs(Fw));grid on phase=atan(imag(Fw)/real(Fw)); subplot(212); ezplot(phase);grid on 2(2).F(w)=exp^(-4*w)^2; syms t Fw=sym(exp(-2*w^2)); ft=ifourier(Fw,t) ezpolt(Fw);grid on ft = 1/4*2^(1/2)/pi^(1/2)*exp(-1/8*t^2) 4. ft=sym(heaviside(t+0.5)-heaviside(t-0.5)); subplot(211); ezplot(ft,[-2,2]); Fw=simplify(fourier(ft)); Yw=Fw*Fw; yt=ifourier(Yw); subplot(212); ezplot(yt,[-2,2]);