数字信号处理概要.doc

西安理工大学 研究生课程论文/研究报告 课程名称： 数字信号分析 课程代号： 任课教师： 论文/研究报告题目： 数字信号处理 课程作业 完成日期： 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(22,2); stem(t,h,.); grid on; title(系统的冲激响应); xlabel(n); ylabel(h(n)); [h,w]=freqz(b,a,256,whole,1); hr=abs(h); %求系统幅频响应； subplot(22,3); plot(w,hr); title(系统的幅频响应); xlabel(\omega/2\pi); ylabel(|{\itH}(e^{j\omega})|); grid on; hphase=angle(h); hphase=unwrap(hphase); %求系统相频响应 subplot(22,4); plot(w,hphase); xlabel(\omega/2\pi); ylabel(\phi(\omega)/rad); title(系统的相频响应); grid on; 输出图像为 二、设计巴特沃斯模拟低通滤波器： 通带：，； 阻带：，； ，。 要求： 1、； 2、； 3、 （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]);title(双线性映射法的系统的对数幅频响应); xlabel(频率/Hz);ylabel(幅值/dB); h=legend(AF,DF);grid on; 3（1）、窗函数法 设定，； 求单位抽样响应：； 截断：由，查表知矩形窗满足要求，其表达式为： ，取M=10，得； 程序： Wp=0.2; %归一化截止频率，相当于除以π Ws=0.6; %归一化截止频率 Wc=(Wp+Ws)/2; %-6dB截止频率 delta_w=Ws-Wp; %过渡带宽 N=10; b=fir1(N-1,Wc,low,rectwin(N)); omega=linspace(0,pi,512); %频率抽样512个点 mag=freqz(b,1,omega); %计算频率响应 magdb=20*log10(abs(mag)); figure; stem([0:N-1],b,.); grid on; axis([0,N,-0.2,0.5]); xlabel(n);ylabel(h(n));title(单位抽样响应); figure; plot(omega/pi,magdb