3、k)=》x(n)h(k-n)A=x(k)*h(k)A(F2-2)“=-00式中,y(k)、x(k)和h(k)分别为对y(t)>x(t)和h(t)以为△吋间间隔进行采样所得的离散序列。相应的可得出结论:若x(k)和h(k)分别心在序号区间[kl,k2]和[k3,R4]有非零的值,则离散卷积(卷积和)y(t)=x(t)*h(t)有非零值的序号区间为[kl+k3,k2+k4]。上机题1.已知两个彳百号兀1(r)=£(r-1)-£(r—2)和x2(0=£(f)—£(/—1),试分别iffli出xx(r),x2(r)和卷积y(r)
4、=x{(r)*x2(r)的波形。(上机原程序及所画出的波形图)源程序:T=0.01;tl=l;t2=2;t3=0;t4=l;t二0:T:t2+t4;xl=ones(size(t)).*((t>tl)-(t>t2));x2=ones(size(t)).*((t>t3)-(t>t4));y二conv(xl,x2)*T;subplot(3,1,1),plot(t,xl);ylabel(,xl(t),);subplot(3,1,2),plot(t,x2);ylabel(,x2(t),);subplot(3,1,3),plot(t
5、,y(l:(t2+t4)/T+l));ylabel('y(t)二xl*x2');xlabel('t/s,);波形图:5ooaCXJx/s上机题2.己知两个信号x(t)=和加)二沧一"纭⑴,试用数值计算法求卷积,并分别画出x(/),M)和卷积y(t)=□)"(/)的波形。(上机原程序及所画出的波形图)源程序:t2=3;t4=ll;T=0.01;t二0:T:t2+t4;x=exp(-t)・*((t>0)一(t>t2));h=t.*exp(~t/2)・*((t>0)一(t>t4));y二conv(x,h)*T;yt=4*exp
6、(-t)+2*t・*exp(T/2*t)-4*exp(T/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(l:(t2+t4)/T+l),t,yt,'--r);legend(Jbynumerical','Theoretical');ylabel('y二x*h,);xlabel('t/s');波形图:*XIIA实验内容:(二)信号的频域分析上机题3.求周期矩
7、形脉冲信号的频谱图,已知A=1,t=0.1$,卩=0.5$(上机原程序及所画出的波形图)源程序:a二1;t£io=0.1;t=2;n0=t/tao;n=0:2*n0;xiao_n0=2*n0+l;nl=2*n0+l;fn_p=a*tao/1*(sin(n*pi*tao/t+cps*(n=0)))・/(n*pi*tao/t+cps*(n二二0));fn_pabs=abs(fn_p);fn_pang二angle(fn_p);fn_mabs=fliplr(fn_pabs(2:xiao_n0));fn_mang=-f1ip1r(
8、fn_pang(2:xiao_n0));fnabs=[fn_mabsfn_pabs];fnang=[fn_mangfn_pang];subplot(2,1,1),stem((-2*n0:2*n0),fnabs);text(15,0.03,'amplitudespectrum,);subplot(2,1,2),stem