数字信号处理（第2版）MATLAB源程序.docx
第1聿MIH1.b埸程求IW1-1和【例12】两个序列的和、末积。I、建立一个主函数命名为main.m程序如下：n1.=-k4;n2=24;x2=(1.-1.02.1.0.-1.();y3,n=sigadd(x1jI2j2):y,n=sigcheng(x1.,nI,x2>n2);S1.1.bP1.o(221)；StcnWn1.x1.Ar):grid;x1.abc1.('n);y1.abc1('x);S1.IbPI3(2.22);stc11Un2.x2,.r);grid:X1.abeKmy1.abc1.('x2');subp1.ot(2.23):sem(n.y3A,);grid;IiI1.式'序列加工X1.abeICn');y1.abe1.(,y3,);subp1.ot(2.2.4);s1.em(n.y.,.):grid;Ii1.ke序列积*)：x1.abe1.Cn,);Wabd)；2、建立个加法函数,命名为sigadd.m：程序如下functiony3.n)*sigadd(x1.n1tx2m2)n=min(min(nI),min(n2)>:nax(max(nIj.max(n2);y1.=zcros(I.IcngihCn»>2=zeros(I,1.engthen);y1(find(n>=min(nI)Wn<=maxnI)=1)=xI：y2(find(n>=min(n2)>(n<=max(n2)>=1.)>=x2:>3三y1.÷y2;3、建立一个一个乘法函数，命名为sigebeng.m：程序如下，fu1.kmy,n=sigcbeng(x1.nIax2j2)n=min(mi(n1),min(n2):max(max(n1Xmax(n2);y1=zeros(I,gth(n>);),2=zeros(IJengh(n);y1.(find(n>=min(n)Wn<=maxnI)=1.)>=x1.;y2(find<(n>=min(n2)><n<=max(n2)=1.)=x2:y三y1.*y2;【例113利用MAT1.AB计算整位脓冲序列次吁I)在-3加3他用内各点的取以W：煽写如下程序计以函数收低。ks=-3;kc=3:n=1:k=kkc;K=KkM=O:令n不等于I时的所有X值为零n等于I时工值为1McnMk,x);ibbc1.('k');【例117已切系统的养分方程如下，试判断该系统是否为移不变浅性系统,>(j)-0.4y(j-1.)+O.3>(n-2)=x(n)-v(-1.)解：列断系统是否为战性、移不变系统可分别从线性和移不变性两个方面来判断,编写两个程序分别实现。(1)判断系统是否为线性系统设3个不同的输入序入X1.O0、±S)和00,并且XOt)=Or1.+加式制,对I处的系统输出序列为y1."D、与Ot)和Nm-若,v()=+加式/0,则系统为雄性系统。MAT1.AB程序如下n=1.:40;a=2;b=3;x1.=c<w(2>>*0.rn)zx2=co2xpi*0.4n>x=a*x1.-H>*x2;nun=(1.2x1.e=1.-<).4.03;>1.=hcr(wn.dcn.x1.);>>2=fihCT<num.dcn.x2);%fihcrFftIi用来求系统的响应y=fi1.kr(nun.den.x);yt=a*y1.+b*>2;MIbPkIK2JJ>Mem<ty):Iiuer轮出y的一幡)：SUbP1.OK2J.2);Stan(IytKIiUcC输出*y1.+b*>,2的振WS'):(2)判断系统是否为核不变系统MATI-AB程序如下.n=1.:40;D=9;a=3b=5;x=a*co2*piH).1.4n)4b*c<M(2,p0.4*n):xd=(zcrs(I.D),x):num1.2:<fcn=I.-04.0.3;)=11I(efnun.dcn.x)yd=A1.1.er<num.den.Ad):HibpIoK2.1.1Xstcimy>i1.c('ff1.入为Wn)时的输出')；MibpIoKX1,2);Mem(W检入为双nD耐的输而)：【例1-19如果系统传递函数为/;)=½-71.+0.8r'+0.5;试用MAT1.AB绘制系统的领率响应特性曲战和手、极点图，并判断系统的稳定性.解:MAT1.AB程序如下，b=1.2=I.O.8.O.5;ftgurt(1>yp1.ane(b.a);%.调用函数/p1.ane求明极点urc(2Xfrcqz(b);第2章【例2-7JIJMAT1.AB求以E序列的Z变换：-V1(/：)=(2*÷3*)/(?)：x2（三）=-1：Xj5)="ES)=CMHI解：拜序如下，symsznwXI=2n+3nzXI=ZtranS(K1):rvt1.y(XI);x2=n1.X2=z11ms(x2);rctty(X2):x3=M2:XAZgnMX3);prctty(X3):x4=cxp(j4w*n)4=ztrans(x4);PrCHy(X4)【例210】未下列函数的Z反变换；X1.j鼠=，_X1.闫12-1(-ay5(z+1.)2解：程序如下XsymsnzaX1.=Z/(24);x1=izrans(X1)X2=z(z-a)2;x2=iztrans(X2)X3三z*(x-I)(z÷I)a2:x3=iztrans(X3>【例2巾】已知像函数X（三）二"0.976zT-2,608-7.015-0.154尸I-1.400,+0j(X)z-0.072zJz>0.600川MAT1.AB求具部分分式展开式并求其反变换.解：程序如下。先部分分式展开b=1.,0.<76.-2.608.1.,01.5,4154;a=1.-1.4(N(6<X>,-4).(172J:(r4>=rcsiducZ(b.a):dispCJ6tt,Wisp<ndi即C极点%d即(p,)dispf4数Misp(k)运行的结果为系数6.1074-1.6605-7.173«极点0.60(X>OMXM)OKXM)3.7269需蹋物器吸f岸2C-VYR-)VSinWe*H-k)S吗。用MAT1.AB可方便推计算出序列的傅甲.叶变摸，问时可画出其福僚性及相频特性.程序如F1N=6:xn=oncs(I.N):n»O:N-hw=1.inpcc(0.2*pi.2<X>kX=xntcx(-jtn,wK%ubp1.ot(311).Mcmn.xn.,b');x1.abdCn,)r1.abdx(n)();subpk1.<31.2).pkH(w.ab5XXT)KIabd(WWabdr福频特性)：subp1.u1.<313).pkM(w.ang1.e<X>')Jdabd(WKy1.abdr相频特性工【例2-18】某离故线性移不变(1.SI)系统.当激励N")=M")时.系统的零状态响应为Nm=(2-0.5-+(-1.5)"1).求系统函数M:)并确定系统的养分方程，并利用MAT1.AB求频率呼应的幅较特性曲戏。用MAT1BA求藤翔特性曲线的程序如F：b=(2OO.5)a=(1.I-O.75W=IinsPaCCW.pi.512)H=frcqz(b,a1w)p1.o<(w/pi.abMH)>Wabe1('幅度)X1.abeICNormaIizedfrequency')【例2-21】已知HR数字池波法的系统函数，所出该泄波2i的必联型结构图借助MAT1.AB的tf2sos函数可以将直接型结构转换为级联型结构，上述例题的程序如下:b=8.-1.11.-2a=11.-5/4.3/4.-1/8)sos,g=t,2sos(b.a)第4章，【例4-3已知序列Ar(Io=SM(OAOK1式)()=0.9*(),试利用快速卷积法计Mr这两个序列的卷积yS>=Mn)”(n).解：程序如下，Nx=15:Nh=20;%产生序列x(11)%产生序列Mn)%序列入S)和Mn>卷枳后备到序列打4的长度%对序列NI1.)作1.jSDFT%对序列h(n)作1.点DFrn1.=kNx-kn2=(kNh-kxn=%in<0.44n1.)n1.>=X)&nI<Nx);hn=O.9n2<n2>=0n2<Nhj;1.=p<>u'2<xtp(>w2(Nx÷Nh-1);Xk=fTt(xn.1.);HknI1.Xhn.1.i;>n=ifuXkHk);MIbPkH(3,1,1);rtcm(x11.,.,ti1.1.H'S力：MibpIoK3J.2);McnMhn.VkHuCeb(I>)，);WbPIot3.1.3);s4enUyn.,ti1.H'义n)');【例4-4已知非周期序列'S)=0<Tmm用MAT1.1.B闰出其时域波形并用FFT对其进行频谱分析,采样点数A'=128.程序如下c1.eara1.1.N=I2X;n=O:N-1.:xn=(1.9.a11;XK=m(xn.N);EagXK=ahUXK):phaXK=ang1.c<XK>SubpIoi(1.XI);p1.ot(n.xn):x1.abck,n,KyIabch'x(n)');i(1.c<,x(n)N=I28*subp!u1.<I.X2);k'D:Iengt1.XnugXKH1.;MenHkjnagXK."/):X1.abcHVXyhbcHTX1.k):1.idd,X(k)N=12R->【例4-5】已知连续时间15号，/=8s8w+8s1.6Rcos2011/选样采样频率,=64Hz变换区间N=16,32.64三种情况进行谱分析.分别打印具幡频特性,弁进行分析和讨论.解r程序如下。Fs=64iT=14;N=1.6;n=0:N-1.;%FFT的变换区间N=16x6nT=cos(8*pi*n*T)Hx*(164pi*n*THws(20*pi4n*T):x6(c)1.6点采怦X<*1.6=1.T«x(mT);Xa6=tTshif(X6k1.6):Tp=N*T:F=1./Tp:k=N2:Na1.k=k*F:%计算x6nT的16点DFT%将岑顿率带到拓m中心%频率分辨率F%产牛164DkT对拔的果样点知率(以零担率为中心)MibpiotOJJ)en>1.k.ab>X6k16>.-.);box<x>%绘刖8点DFr的幅箱特性图UUe<16点IDFT(X(nT)1.abe1.fHz)'Xybbd('幡度
