《数字信号处理》课后上机题..docVIP

实验报告 第一章:时域离散信号和时域离散系统 *16.已知两个系统的差分方程分别为 (1) y(n)=0.6y(n-1)-0.08y(n-2)+x(n) (2) y(n)=0.7y(n-1)-0.1y(n-2)+2x(n)-x(n-2) 分别求出所描述的系统的单位脉冲响应和单位阶跃响应. 解:(可附程序) 系统差分方程的系数向量为 B1=1,A1=[1,-0.6,0.08] 系统差分方程的系数向量为 B2=[2,0,-1],A2=[1,-0.7,0.1] 调用MATLAB函数filter计算两个系统的单位脉冲响应和单位阶跃响应的程序%B1=1;A1=[1,-0.6,0.08]; %设差分方程(1)系数向量 B2=[2,0,-1];A2=[1,-0.7,0.1]; %设差分方程(2)系数向量 %系统1 xn=[1,zeros(1,30)]; ys=0; %xn=单位脉冲序列,长度N=31 xi=filtic(B1,A1,ys); %由初始条件计算等效初始条件输入序列xi hn1=filter(B1,A1,xn,xi); %调用filter解差分方程,求系统输出信号hn1 n=0:length(hn1)-1; subplot(2,2,1);stem(n,hn1,.) title((a)系统1的系统单位脉冲响应); xlabel(n);ylabel(h(n)) xn=ones(1,30); %xn=单位阶跃序列,长度N=31 sn1=filter(B1,A1,xn,xi); %调用filter解差分方程,求系统输出信号sn1 n=0:length(sn1)-1; subplot(2,2,2);stem(n,sn1,.) title((b)系统1的单位阶跃响应); xlabel(n);ylabel(s(n)) %系统2 xn=[1,zeros(1,30)]; %xn=单位脉冲序列,长度N=31 xi=filtic(B2,A2,ys); %由初始条件计算等效初始条件输入序列xi hn2=filter(B2,A2,xn,xi); %调用filter解差分方程,求系统输出信号hn2 n=0:length(hn2)-1; subplot(2,2,3);stem(n,hn2,.) title((a)系统2的系统单位脉冲响应); xlabel(n);ylabel(h(n)) xn=ones(1,30); %xn=单位阶跃序列,长度N=31 sn2=filter(B2,A2,xn,xi); %调用filter解差分方程,求系统输出信号sn2 n=0:length(sn2)-1; subplot(2,2,4);stem(n,sn2,.) title((b)系统2的单位阶跃响应); xlabel(n);ylabel(s(n)) (结果展示): 第二章:时域离散信号和系统的频域分析 *30.假设系统函数如下式: 试用MATLAB语言判断系统是否稳定. 解:(可附程序): %调用roots函数求极点,并判断系统的稳定性 A=[3,-3.98,1.17,2.3418,-1.5147]; %H(z)的分母多项式系数 p=roots(A) %求H(z)的极点 pm=abs(p); %求H(z)的极点的模 if max(pm)1 disp(系统因果稳定) else disp(系统因果不稳定) end (结果展示): p = -0.7486 0.6996 + 0.7129i 0.6996 - 0.7129i 0.6760 系统因果稳定 第三章:离散傅里叶变换(DFT) *24.给定两个序列:={2,1,1,2},={1,-1,-1,1}. (1)直接在时域计算与的卷积; (2)用DFT计算与的卷积,总结出DFT的时域卷积定理. 解:(可附程序) x1n=[2 1 1 2];x2n=[1 -1 -1 1]; %时域直接计算卷积yn: yn=conv(x1n,x2n) %用DFT计算卷积ycn: M1=length(x1n); M2=length(x2n); N=M1+M2-1; X1k=fft(x1n,N); %计算x1n的N点DFT X2k=fft(x2n,N); %计算x2n的N点DFT Yck=X1k.*X2k; ycn=ifft(Yck,N) (结果展示): yn = 2 -1 -2 2 -2 -1 2 ycn = 2.0000 -1.0000 -2.0000 2.0000 -2.0