数字信号处理技术总结.docVIP

西安理工大学 研究生课程论文/研究报告 课程名称： 数字信号分析 课程代号： 任课教师： 论文/研究报告题目： 数字信号处理 课程作业 完成日期： 2015 年 6 月 25 日 学 科： 机械电子工程 学 号： 姓 名： 成 绩： 一、作业内容 已知 求：1、该系统的阶跃响应，并画出图形； 2、该系统的冲激响应，并画出图形； 3、该系统的频率响应，画出其幅频响应和相频响应。 2.Matlab程序如下： x=ones(100); t=1:100; b=[0.001836,0.007344,0.011016,0.007374,0.001836]; %形成向量b； a=[1,-3.0544,3.829,-2.2925,0.55075]; %形成向量a； y=filter(b,a,x); %求系统阶跃响应； figure(1); subplot(2,2,1); %将窗口分成2行*2列个子图，并在一个中绘制图形； plot(t,x,g.,t,y,k-); %将x(n)和y(n)分别用绿色和黑色画在同一张图上 grid on; title(系统的阶跃响应); xlabel(n); ylabel(x(n),y(n)); [h,t]=impz(b,a,40); %求系统冲激响应； subplot(2,2,2); stem(t,h,.); grid on; title(系统的冲激响应); xlabel(n); ylabel(h(n)); [h,w]=freqz(b,a,256,whole,1); hr=abs(h); %求系统幅频响应； subplot(2,2,3); plot(w,hr); title(系统的幅频响应); xlabel(\omega/2\pi); ylabel(|{\itH}(e^{j\omega})|); grid on; hphase=angle(h); hphase=unwrap(hphase); %求系统相频响应 subplot(2,2,4); plot(w,hphase); xlabel(\omega/2\pi); ylabel(\phi(\omega)/rad); title(系统的相频响应); grid on; 输出图像为 二、设计巴特沃斯模拟低通滤波器： 通带：，； 阻带：，； ，。 要求： 1、设计出； 2、两种方法（冲激响应不变法，双线性映射法）转化成； 3、两种方法（窗函数法，频率抽样法）设计出FIR、LP、DF; 4、绘图比较幅频特性。 解： 1、参数归一化， ，， ，，； 求幅平方函数，； 求C、N，当时，C=1， 确定，得，令分母为0，， 得，， 则， 则。 2（1）、冲激响应不变法 令，，，， 则。 程序： c=[40000*pi^2]; %形成向量c； d=[1,sqrt(2)*200*pi,40000*pi^2]; %形成向量d； [h,w]=freqs(c,d,1000); b=[0,0.2449]; %形成向量b； a=[1,-1.1580,0.4112]; %形成向量a； [h1,w1]=freqz(b,a,1000); plot(w/(2*pi),20*log10(abs(h)),b);hold on; plot(w1*1000/(2*pi),20*log10(abs(h1)),r:); axis([0,500,-30,0]);title(冲激响应不变法IIR的系统的对数幅频响应); xlabel(频率/Hz);ylabel(幅值/dB); h=legend(AF,DF);grid on; （2）双线性映射法 ，，仍有，； 令，，，得，得，。 则。 程序： c=[0.3249^2*2000^2]; %形成向量c； d=[1,0.4595*2000,0.3249^2*2000^2]; %形成向量d； [h,w]=freqs(c,d,1000); b=[0.06745,0.1349,0.06745]; %形成向量b； a=[1,-1.143,0.4128]; %形成向量a； [h1,w1]=freqz(b,a,1000); plot(w/(2*pi),20*log10(abs(h)),b);hold on; plot(w1*1000/(2*pi),20*log10(abs(h1)),r:); axis([0,350,-30,0