资源描述：
《信号与系统实验答案new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验三1,.利用DFT近似分析连续信号x(t)=e-2tu(t)的幅度谱并与理论值比较，将理论频谱曲线和实际计算频谱曲线绘制在一个坐标系中。（要求根据实际幅度频谱函数
2、X(jω)
3、选择合适的抽样频率，根据时域波形选择合适的窗长度，根据序列点数选择合适的DFT点数。同时，减小抽样频率，观察最终理论值与计算值间的误差变化。）fsam=50;Tp=6;N=512;T=1/fsam;t=0:T:Tp;x=exp(-2*t);X=T*fft(x,N);plot(t,x);xlabel('t');title('时域波形');w=(-N/2:N/2-1)*(2*pi/N
4、)*fsam;y=1./(j*w+2);figure;plot(w,abs(fftshift(X)),w,abs(y),'r-.');title('幅度谱');xlabel('w');legend('计算值','理论值');2.近似分析门函数信号的幅度谱，并与理论值比较，将理论频谱曲线和实际计算频谱曲线绘制在一个坐标系中，其中分别选其最高频带上限为、、时三种情况，比较结果并简单解释其区别及原因。（根据门函数的理论频谱表达式，当时值为0，并随自变量绝对值的增大呈递减趋势）fsam=16;N=512;T=1/fsam;t=-2:T:2;x=[(t>=-1)&(
5、t<=1)];X=T*fft(x,N);%消除1/T因子的影响plot(t,x);xlabel('t');title('时域波形');w=(-N/2:N/2-1)*(2*pi/N)*fsam;y=2*sin(w)./w;%理论频谱值figure;plot(w,abs(fftshift(X)),w,abs(y),'r-.');title('幅度谱');xlabel('w');legend('计算值','理论值');实验四。1.已知两有限长序x[k]=[1,2,0,-1,3,2;k=-2,-1,0,1,2,3],h[k]=[1,-1,1;k=0,1,2]。计算
6、卷积y[k]=x[k]*h[k],在同一窗口中绘出x[k]、h[k]及y[k]的波形。x=[1,2,0,-1,3,2];kx=-2:3;h=[1,-1,1];kh=0:2;y=conv(x,h);k=(kx(1)+kh(1)):(kx(end)+kh(end));subplot(3,1,1);stem(kx,x);subplot(3,1,2);stem(kh,h);subplot(3,1,3);stem(k,y);xlabel('k');ylabel('y');2.给定两个离散时间序列x[k]=0.5k{u[k]-u[k-8]},h[k]=u[k]-u[k
7、-8]，计算它们的卷积，并分别绘制x[k]、h[k]和它们的卷积y[k]的图形。k=-2:10;h=(k>=0)-(k>=8);x=(0.5.^k).*h;y=conv(x,h);ky=2*k(1):2*k(end);subplot(3,1,1);stem(k,x);xlabel(‘k’);ylabel(‘x’);subplot(3,1,2);stem(k,h);xlabel(‘k’);ylabel(‘h’);subplot(3,1,3);stem(ky,y);xlabel(‘ky’);ylabel(‘y’);3.已知连续系统的微分方程为，计算该系统的单位
8、冲激响应和单位阶跃响应，绘制响应的曲线波形图，并与理论结果比较。a=[1,4,3];b=[1];sys=tf(b,a);h=impulse(sys);g=step(sys);subplot(2,2,1);plot(h);title(‘系统的单位冲击响应h(t)’);subplot(2,2,2);plot(g);title(‘系统的单位冲击响应g(t)’);t=0:0.01:10;h=(exp(-t)-exp(-3*t))/2;g=(-exp(-t))/2+exp(-3*t)/6+1/3;subplot(2,2,3);plot(h);title(‘系统的单位
9、冲击响应的理论值h(t)’);subplot(2,2,4);plot(g);title(‘系统的单位冲击响应的理论值g(t)’);4.计算并绘制由如下微分方程表示的系统在输入信号为x(t)=(e-2t-e-3t)u(t)时的系统零状态响应并描绘曲线。a=[1,3,2];b=[8];sys=tf(b,a);t=0:0.01:10;x=exp(-2*t)-exp(-3*t);y=lsim(sys,x,t);plot(t,y);5.已知连续系统，求输入分别为和时系统的零状态响应，将响应图形分别显示在两个图形窗口，并与理论结果比较。a=[1,3,2,0];b=[4
10、,1];sys=tf(b.a);t=0:0.01:10;x1=si
