信号课程实验报告.docVIP

课程实验报告 题 目： 连续时间信号的卷积 及信号的频域分析 学 院 通达学院 学 生 姓 名 沈元强 班 级 学 号 指 导 教 师 兰吕瑞 开 课 学 院 通信与信息工程学院 日 期 2011年6月17日）和，试分别画出和卷积的波形。 （上机原程序及所画出的波形图） 解：原程序如下： T=0.01; t1=1;t2=2; t3=0;t4=1; t=0:T:t2+t4; x1=ones(size(t)).*((tt1)-(tt2)); x2=ones(size(t)).*((tt3)-(tt4)); y=conv(x1,x2)*T; subplot(3,1,1),plot(t,x1); ylabel(x1(t)); subplot(3,1,2),plot(t,x2); ylabel(x2(t)); subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1)); ylabel(y(t)=x1*x2); xlabel(-----t/s); 波形图如下： 上机题2.已知两个信号和，试用数值计算法求卷积，并分别画出和卷积的波形。 （上机原程序及所画出的波形图） 解：原程序如下： 数值法： syms tao; t=sym(t,positive); xt=exp(-t); ht=t*exp(-t/2); xh_tao=subs(xt,t,tao)*subs(ht,t,t-tao); yt=int(xh_tao,tao,0,t) yt = 4*exp(-t)+2*t*exp(-1/2*t)-4*exp(-1/2*t) 波形图： t2=3;t4=11; T=0.01; t=0:T:t2+t4; x=exp(-t).*((t0)-(tt2)); h=t.*exp(-t/2).*((t0)-(tt4)) y=conv(x,h)*T; yt=4*exp(-t)+2*t.*exp(-1/2*t)-4*exp(-1/2*t); subplot(3,1,1),plot(t,x); ylabel(x(t)); subplot(3,1,2),plot(t,h); ylabel(h(t)); subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1),t,yt,--r); legend(by numerical,Theoretical); ylabel(y=x*h); xlabel(----t/s); 波形图如下： 实验内容：（二） （上机原程序及所画出的波形图） 解：原程序如下： a=1;tao=0.1;t=0.5; n0=t/tao; n=0:2*n0; fn_p=a*tao/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0)); fn_pabs=abs(fn_p); fn_pang=angle(fn_p); fn_mabs=fliplr(fn_pabs(2:11)); fn_mang=-fliplr(fn_pang(2:11)); fnabs=[fn_mabs fn_pabs]; fnang=[fn_mang fn_pang]; subplot(2,1,1),stem((-2*n0:2*n0),fnabs); text(4,0.11,amplitude spectrum); subplot(2,1,2),stem((-2*n0:2*n0),fnang); text(-2,2,phase spectrum); xlabel(n);grid 波形图如下： 问题2：改变信号的周期，比较他们的计算结果。 直接在程序中改变T会报错。应在第4行输入：xian_no=2*no+1,而且将程序中的11改为xian_no （学生通过上机画出不同周期对应的计算结果并回答问题） 解： 1．T=0.1 原程序： a=1;tao=0.1;t=0.1; n0=t/tao; n=0:2*n0; xian_n0=2*n0+1 fn_p=a*tao/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0)); fn_pabs=abs(fn_p); fn_pang=angle(fn_p); fn_mabs=fliplr(fn_pabs(2: xian_n0)); fn_mang=-fliplr(fn_pang(