哈尔滨理工大学电技数字信号处理实验报告

《哈尔滨理工大学电技数字信号处理实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、哈尔滨理工大学实验报告课程名称：数字信号处理学院：自动化学院专业班级：电技12-2学生姓名：学号：指导教师：马静实验课程名称数字信号处理时间：2014.11实验名称离散系统的差分方程、卷积分析和频率响应地点：B302姓名学号班级：电技12-2成绩评定：一.实验目的(1)加深对离散系统的差分方程求解过程的理解。(2)掌握离散系统的线性卷积分析方法(3)观察离散系统的频率响应。二、实验内容：(1)编制程序求解下列两差分方程的单位冲激响应，并绘出其图形。在MATLAB屮,可以用函数y=Filter(a,b,x)求解差分方程。y[n]+0

2、.75-1]+0.125)’[/?-2]=x[n]-x[n-l]y[n]=0.25{x[n-1]+x[n-2]+x[n-3]+x[n-4]}(2)应用卷积性质继续求解两系统的单位阶跃响应。在MATLAB中，可以用函数y=Conv(x,h)计算卷和。(3)求系统幅度频率响应和相位频率响应。_0.0528+0.0797厂1+0.1295厂2+0.1295厂3+0.797Z"4+0.0528z—5~1-1.8107Z-1+2.4947z-2-1.8801z"3+0.9537Z-4-0.2336z~5—三实验程序n=15;a=[l,-1]

3、;b=[l,075,0.125];x=[1zeros(1,n-1)];k=0:1:n-1;y=filter(a、cr,X)、•stem?y)；xlabel(•rV)、•ylabel(••05Im;)7【2】n=15;a=【0,0.25,0.25,0.25,0.25】；b=〔ll;x=【1zeros《1、nil二；k=0:1:nil;y=filter(a,)；M=l

4、ength(c)—1;nHO:1:M;stem(n,c)；xlacrel(•n*rylacrel(••-)/sa=【0、0.25,0.25,0.25、0.25Hb=【lkc=conv(a,k>)、•M=length(c)—1;n=0:1:M;stem(n,c)；xlabel(-n-)、•ylabel(-•03沖-)7泞=256;num=【0.05280.07970.12950.12950.7970.0528J;den=【l-1.81072k947-1.88010.9537-0.2336Hw=o:pi/k::pi;h=HbreqN{

5、ncm,den,w)/s&plot(2、2,1)；plot(w/pi,real(h))、•gridtirtle(m>r^-)xlabel(-xomelQa/x^i-r.ylacrel(-•沖's&plot《2、2,2)；勺lot(w/^5i、imag(h))/grid/title(-sor-)xlabel(•omega/XJiJ:;ylabel(••0白办-)；subplot《2,2、3)；plot(w/pi,abs(h));grid;title(1-u^E^x1)；xlabel(1omega/pi1);ylabel(1-

6、’)；subplot(2,2,4);plot(w/pi,angle(h));gridtitle('Ial»^x')；xlabel(1omega/pi1);ylabel(1Ial»(rad)1);igure1□回®虚部FileEditViewInsertToolsDesktopWindowHelp实部00.5co/ti幅度谱las倒Ides00.510co/ti0.5150-5co/ti相位谱0.51co/ti202■b=[0.0580.07970.12950.12950.7970.0528];a=[l-1.81072.4947

7、-1.88010.9537-0.2336];zplane(b,a);JFigiigure1□回®FileEditViewInsertToolsDesktopWindowHelp10鳃50实部1500.51105O-5-1Oid—.5/71OS-5co/ti1-2-101RealPart21o1■[3]b=[0.0580.07970.12950.12950.7970.0528];a=[l-1.81072.4947-1.88010.9537-0.2336];G=tf(b,a);bode(G);FileEditViewInsertToo

8、lsDesktopWindowHelp实部虚部^5laffi-505側loS1510-4-10100.5co/ti010510.5co/ti幅度谱1510MSa5BodeDiagramoooo224■imprpnl!u6eH0.5co/tid—(O