数字信号处理（第2版）MATLAB源程序.docx

《数字信号处理（第2版）MATLAB源程序.docx》由会员分享，可在线阅读，更多相关《数字信号处理（第2版）MATLAB源程序.docx（16页珍藏版）》请在优知文库上搜索。

1、第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,x2n2);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.

2、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=min(n2)(n=1.)=x2:3三y1.y2;3、建立一个一个乘法函数，命名为sigebeng

3、.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=x1.;y2(find=min(n2)n(j)-0.4y(j-1.)+O.3(n-2)=x(n)-v(-1.)解：列断系统是否为战性、移不变系统可分别从线性和移不变性两个方面来判断,编写两个程序分别实现。(1)判断系统是否为线性系统设3个不同的输入序入X1.O0、S)和00,并且XOt)=Or1.+加式制,对I处的系统输出序列

4、为y1.D、与Ot)和Nm-若,v()=+加式/0,则系统为雄性系统。MAT1.AB程序如下n=1.:40;a=2;b=3;x1.=c*0.rn)zx2=co2xpi*0.4nx=a*x1.-H*x2;nun=(1.2x1.e=1.-1.=hcr(wn.dcn.x1.);2=fihCT2;MIbPkIK2JJMem,2的振WS):(2)判断系统是否为核不变系统MATI-AB程序如下.n=1.:40;D=9;a=3b=5;x=a*co2*piH).1.4n)4b*cM(2,p0.4*n):xd=(zcrs(I.D),x):num1.2:fcn=I.-04.0.3;)=11I(efnun.dcn.

5、x)yd=A1.1.eri1.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(1yp1.ane(b.a);%.调用函数/p1.ane求明极点urc(2Xfrcqz(b);第2章【例2-7JIJMAT1.AB求以E序列的Z变换：-V1(/：)=(2*3*)/(?)：x2（三）=-1：Xj5)=ES)=CMHI解：拜序如下，symsznw

6、XI=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)(zI)a2:x3=iztrans(X3【例2巾】已知像函数X（三）二0.976zT-2,608-7.015-0.

7、154尸I-1.400,+0j(X)z-0.072zJz0.600川MAT1.AB求具部分分式展开式并求其反变换.解：程序如下。先部分分式展开b=1.,0.;76.-2.608.1.,01.5,4154;a=1.-1.4(N(6,-4).(172J:(r4=rcsiducZ(b.a):dispCJ6tt,WispOMXM)OKXM)3.7269需蹋物器吸f岸2C-VYR-)VSinWe*H-k)S吗。用MAT1.AB可方便推计算出序列的傅甲.叶变摸，问时可画出其福僚性及相频特性.程序如F1N=6:xn=oncs(I.N):nO:N-hw=1.inpcc(0.2*pi.2kX=xntcx(-jt

8、n,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)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

9、.pi.512)H=frcqz(b,a1w)p1.oWabe1(幅度)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)%序

10、列入S)和Mn卷枳后备到序列打4的长度%对序列NI1.)作1.jSDFT%对序列h(n)作1.点DFrn1.=kNx-kn2=(kNh-kxn=%in=X)&nINx);hn=O.9n2=0n2Nhj;1.=pu2w2(NxNh-1);Xk=fTt(xn.1.);HknI1.Xhn.1.i;n=ifuXkHk);MIbPkH(3,1,1);rtcm(x11.,.,ti1.1.HS力：MibpIoK3J.2);McnMhn.VkHuCeb(I)，);WbPIot3.1.3);s4enUyn.,ti1.H义n);【例4-4已知非周期序列S)=0Tmm用MAT1.1.B闰出其时域波形并用FFT对其进

11、行频谱分析,采样点数A=128.程序如下c1.eara1.1.N=I2X;n=O:N-1.:xn=(1.9.a11;XK=m(xn.N);EagXK=ahUXK):phaXK=ang1.c;SubpIoi(1.XI);p1.ot(n.xn):x1.abck,n,KyIabchx(n);i(1.c,x(n)N=I28*subp!u1.;【例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点采怦X1.k.abX6k16.-.);box%绘刖8点DFr的幅箱特性图UUe16点IDFT(X(nT)1.abe1.fHz)Xybbd(幡度