数字信号处理编程作业-机电.docVIP

数字信号处理编程作业 姓名： 学号： 单位： 专业： 导师： 实验一 离散时间系统及离散卷积 离散时间系统的单位脉冲响应 （1）选择一个离散时间系统； y(n)=x(n)+2y(n-1) y(n)=0,n0;h(n)=0,n0; （2）用笔进行差分方程的递推计算； 当x(n)=，y(n)=h(n)，得 h(0)=1+2y(-1)=1; h(1)=0+2h(0)=2; h(2)=0+2h(1)=4; h(n)=0+2h(n-1)=2n. 从而，系统函数的单位冲击响应为 h(n)=2nu(n); （3）编制差分方程的递推计算程序； n=[0:9]; y=zeros(1,length(n)); a=2; for m=1:10 if(m==1) y(1,1)=1; else y(1,m)=a*y(1,m-1); end end subplot(1,1,1),stem(n,y);title(3??¤?ìó|);xlabel(n);ylabel(h(n)); （4）在计算机上实现递推运算； 得到的结果，如图所示 （5）将程序计算结果与笔算的计算结果进行比较，验证程序运行的正确性； y y = 1 2 4 8 16 32 64 128 256 512 与笔算结果完全符合。 1-2 参考实例验证： 离散时间系统的单位脉冲响应 %单位取样序列函数 function [x,n]=impseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)==0]; %单位阶跃序列函数 function [x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)=0]; (1) 参考实例验证 [x,n]=impseq(0,-20,120); a=[1,-1,0.9];b=1; h=filter(b,a,x); subplot(2,1,1),stem(n,h);title(冲激响应);xlabel(n);ylabel(h(n)); subplot(2,1,2),zplane(b,a);title(系统零极点图); 离散系统的幅频、相频的分析方法 参考实例验证 b=[0.0181,0.0543,0.0543,0.0181]; a=[1.000,-1.76,1.1829,-0.2781]; function work1(h,N) sum=0; w=0:0.1:pi; for n=1:N sum=sum+h(n).*exp(-j*n*w); end Hm=abs(sum); Ha=angle(sum); subplot(2,1,1); plot(w/pi,Hm); title(幅度响应); grid; ylabel(幅度);xlabel(以\pi为单位的频率); subplot(2,1,2); plot(w/pi,Ha);title(相位响应);grid; ylabel(相位/pi);xlabel(以\pi为单位的频率); b=[0.0181,0.0543,0.0543,0.0181];a=[1.000,-1.76,1.1829,-0.2781]; [x,n]=impseq(0,0,100);h=filter(b,a,x); work1(h); 离散卷积的计算 (1) 参考实例验证 %矩形序列函数 function [x,n]=rectseq(n0,n1,n2,n3) n=[n2:n3]; x=[(n-n0=0)(n-n10)]; [x,n0]=rectseq(0,10,-5,50);[y,n1]=stepseq(0,-5,50);h=y.*(0.9.^n1); c=conv(x,h); subplot(3,1,1);stem(n0,x);title(输入序列); xlabel());ylabel(x(n));axis([-5,50,0,2]); subplot(3,1,2);stem(n1,h);title(冲激响应); xlabel(n);ylabel(h(n));axis([-5,50,0,2]); subplot(3,1,3);stem(-10:100,c);title(输出); xlabel(n);ylabel(y(n));axis([-10,100,0,8]); 实验二 离散傅立叶变