连续时间信号的卷积及信号的频域分析实验报告

课程实验报告题

目：

连续时间信号的卷积及信号的频域分析学 学生姓名班级学号指导教师开课学院通信与信息工程学院日 11,ttt3,t11,ttt3,t33,t11,k33,k33,kx(n)h(kn)，记作记作y(k)x(n)h(kn)问题1：用计算机算卷积是把连续信号进行采样，得到一个个离散数值，然后用的卷积，请推导数值计算与连续信号的卷积之间的关系。生回答问题）解：连续函数x(t)和h(t)的卷积为：y(t)x(t)h(t)x()h(t)d2-1)若x(t)和h(t)分别在时间区间 和 有非零的值，则 y(t)x(t)[(tt1)(tt2)]h(t)[(tt3)(tt4)]

x()[(t1)(t2)]h(t)[(tt3)(tt4)]d要使y(t)为非零值，必须有(t1)(t2)=1和[(tt3)(tt4)=1从而，应同时满足：t

和t

t

t

，即t

和 有非零的值， 则卷积y(t)x(t)h(t)有非零值得时间区间为t t。1 2 nn

（F2-2 ]和[k 2 有非零的值，则离散卷积（卷积和）y(t)x(t)h(t)有非零值的序号区间为[kk k]。 上机题1.已知两个信号x1(t)(t1)(t2)和x2(t)(t)(t1)，试分别画出x1(t),x2(t)和卷积y(t)x1(t)x2(t)的波形。机原程序及所画出的波形图）T=0.01;t1=1;t2=2;t3=0;t4=1;x1=ones(size(t)).*((t>t1)-(t>t2));x2=ones(size(t)).*((t>t3)-(t>t4));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)');1,3),plot(t,y(1:(t2+t4)/T+1));ylabel('y(t)=x1*x2');abel('----->t/s');上机题2.已知两个信号x(t)et(t)和h(t)tet/2(t)，试用数值计算法求卷积，并分别画出x(t),h(t)和卷积y(t)x(t)h(t)的波形。机原程序及所画出的波形图）t4=11;T=0.01;x=exp(-t).*((t>0)-(t>t2));t/2).*((t>0)-(t>t4));conv(x,h)*T=4*exp(-t)+2*t.*exp(-1/2*t)-4*exp(-1/2*t);1,1),plot(t,x);ylabel('x(t)');subplot(3,1,2),plot(t,h)abel('h(t)');1,3),plot(t,y(1:(t2+t4)/T+1),t,yt,'--r');legend('bynumerical','Theoretical');abel('y=x*h');abel('----->t/s');上机题3.求周期矩形脉冲信号的频谱图，已知A1,0.1s,T0.5s机原程序及所画出的波形图）1;tao=0.1;t=0.5;=t/tao;n=0:2*n0;/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0));_pabs=abs(fn_p);fn_pang=angle(fn_p);fn_mabs=fliplr(fn_pabs(2:11));liplr(fn_pang(2:11));fnabs=[fn_mabsfn_pabs];fnang=[fn_mangfn_pang];1,1),stem((-2*n0:2*n0),fnabs);text(4,0.11,'amplitudespectrum');1,2),stem((-2*n0:2*n0),fnang);text(-2,2,'phasespectrum');abel('n');grid问题2：改变信号的周期，比较他们的计算结果。生通过上机画出不同周期对应的计算结果并回答问题）>>m=16;>>T=2;tao=0.5;>>f0=1/T;fc=2;>>dt=T/m;>>n=0:m-1;>>t=n*dt;>>xt1=(1-t/tao).*((t>=0)-(t>=tao));>>xt2=(t/tao+1-2/tao).*(t>=2-tao);>>xt=xt1+xt2=Columns1through61.0000 0.7500 0.5000 0.2500 0 Columns7through120 0 0 0 0 Columns13through160 0.2500 0.5000 >>H=ones(size(n)).*((n>=m/2-fc/f0)-(n>m/2+fc/f0));>>xn=fft(xt);>>cn=fftshift(xn/m);>>yn=cn.*H;>>yt=ifft(yn)*m;>>subplot(3,2,1),stem((n-m/2)*f0,cn);>>text(-3,0.3,'spectralcoefficientsofx(t)');>>subplot(3,2,3),plot(n-m/2*f0,H);>>text(-1,0.5,'H(f)');>>subplot(3,2,5),stem((n-m/2)*f0,yn);>>text(-3,0.3,'outputspectralcoefficients');>>xlabel('---->f/Hz')>>subplot(3,2,2),pl