信号与系统实验三连续时间周期信号的傅里叶级数.docxVIP

实验三连续时间周期信号的傅里叶级数 一、实验目的 掌握连续时间周期信号的傅里叶级数的展开和合成，理解吉布斯现象，掌握 周期矩形脉冲信号的频谱及脉冲宽度、周期对周期信号频谱的影响。 实验内容 实验内容 1、周期信号的傅里叶级数的展开和合成 画出如下图对称方波（取E=l、T=l），并采用有限项傅里叶级数对原函数进 行逼近，画出对称方波的1、3、5、7、9、11次谐波的傅里叶级数合成波形，观 察吉布斯现象。 sum=O; t=-3:0.01:3; E 二 1 ;T= 1 血二 T/2;w=2*3」4159/T; for n=l:l fn=(2*E*ta/T)*sin(w *ta*n/2)/(w *ta*n/2); f=(E*ta/T)+cos(n*w*t)*fn-E/2; sum=sum+f; end plot(t,sum) sum=O; t=-3:0.01:3; E=1 ；T=1 ;ta=T/2;w=2*3」4159/T; for n=l:3 fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2); f=(E*ta/T)-i-cos(n*w*t)*fn-E/2; sum=sum+f; end plot(t,sum) sum=O; t=-3:0.01:3; E=1 ；T=1 ;ta=T/2;w=2*3」4159/T; for n=l:5 fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2); f=(E*ta/T)-i-cos(n*w*t)*fn-E/2; sum=sum+f; end plot(t,sum) sum=O; t 二 3:0.01:3; E= 1 ;T=1 ;ta 二 T/2; w 二2*3.14159/T; for n=l:7 fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2); 仁(E*ta/T)+cos(n*w*t)*fn?E/2; sum二sum+f; end plot(t,sum) sum=0; t=-3:0.01:3; E= 1 ;T= 1 ;ta=T/2;w=2*3.14159/T; for n=l:9 fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2); f=(E*ta/T)+cos(n*w*t)*fn-E/2; sum=sum+f; end plot(t,sum) sum=O; t=-3:0.01:3; E= 1 ;T= l;ta=T/2;w=2*3」4159/T; for n= 1:11 fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2); f=(E*ta/T)+cos(n*w*t)*fn-E/2; sum=sum+f; end plot(t,sum) 2、周期矩形脉冲信号的频谱 取E=l, t=1,画出周期矩形脉冲(教材P83图3?6)的傅里叶级数的频谱(教 材P83图3?7)； n二 12:12; E= 1 ;t= 1 ;T=5*t;w=2/T; fn=(E*t/T)*sinc(w*t*n/2); stem(n,fn,filled); hold on k二12:0.01:12; f= abs(E*t/T)*sinc(w*t*k/2); plot(k,f,J); 取E=l, t=1,画出教材P85图3-8(a)； 012:0.01:12; y=u(t+l/4)-u(t-l/4)+u(t-l 9/4)-u(t-2 l/4)-u(t+19/4)+u(t+2 l/4)+u(t-39/4)-u(t-41/4)-u( t+39/4)+u(t+41/4); subplot(2,1,1 );plot(t,y);axis([-12 12 -0.1 1.1 ]);xlabel(t);ylabel(^f(t)^); n=-12:12; E= 1 ;t= 1 ;T= 10*t; w=2/T; fn=(E*t/T)*sinc(w*t*n/2); subplot(2,1,2);stem(n,fn,,filled,);hold on; k二 12:0.01:12; f= abs(E* t/T) * sinc( w * t* k/2); xlabel(,w,);ylabel(Fn,); 取E=], T=l,画出教材P85图3-8(c)o 匸? 12:0.01:12; y=u(t+1 /4)-u(t- l/4)+u(t-39/4)-u(t-4 l/4)-u(t+39/4)+u(t+41/4); subplot(2,1,1 );plot(t,y);axis([-12 12 -0.1 l?l]);xlabel(t);ylabel(f(t)); n二 12:12; E= 1 ;t= 1 ;T=5*t;w=2