信号与系统测试答案.docVIP

实验三 1,. 利用DFT近似分析连续信号x(t)=e-2tu(t)的幅度谱并与理论值比较，将理论频谱曲线和实际计算频谱曲线绘制在一个坐标系中。 （要求根据实际幅度频谱函数|X(jω)|选择合适的抽样频率，根据时域波形选择合适的窗长度，根据序列点数选择合适的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)*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)(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]。计算卷积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-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.已知连续系统的微分方程为，计算该系统的单位冲激响应和单位阶跃响应，绘制响应的曲线波形图，并与理论结果比较。 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(‘系统的单位冲击响应的理论值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(