信号与系统实验报告(连续信号时域卷积).docVIP

课程实验报告 题 目： 连续时间信号的卷积 及信号的频域分析 学 院 通信与信息工程学院 学 生 姓 名 林观佑 班 级 学 号 指 导 教 师 周 波 开 课 学 院 通信与信息工程学院 日 期 2011年11月日）,转化为数值计算后则为：h（ke）=，将上式中ke用k替代，则有h（k）=，此时若令上式中的时间间隔e无限接近于零，则所得的h（k）的极限值即为两函数卷积函数值。 上机题1.已知两个信号和，试分别画出和卷积的波形（上机原程序及所画出的波形图）。 上机原程序: 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.已知两个信号和，试用数值计算法求卷积，并分别画出和卷积的波形。 （上机原程序及所画出的波形图） 上机原程序: 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：改变信号的周期，比较他们的计算结果。 （学生通过上机画出不同周期对应的计算结果并回答问题） 周期为0.45时有： a=1;tao=0.1;t=0.45; 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:10)); fn_mang=-fliplr(fn_pang(2:10)); 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