信號(hào)實(shí)驗(yàn)報(bào)告

1、西安郵電學(xué)院信號(hào)與系統(tǒng)實(shí)驗(yàn) 報(bào) 告院系名稱：西安郵電學(xué)院專業(yè)班級(jí)：學(xué)生姓名：學(xué)號(hào)：時(shí)間：2010年 06 月 14 日實(shí)驗(yàn)內(nèi)容一 一熟悉簡單的矩陣輸入1從屏幕上輸入矩陣A=1 2 3；4 5 6；7 8 9 或A=1，2，3；4，5，6；7，8，9 觀察輸出結(jié)果。2試用回車代替分號(hào)，觀察輸出結(jié)果。3輸入矩陣B=9，8，7；6，5，4；3，2，1 C=4，5，6；7，8，9；1，2，3, 鍵入A B C 觀察結(jié)果4選擇File|new菜單中的M-file，輸入B=9 ，8，7；6，5，4；3，2，1，保存為B.M文件，退出編輯環(huán)境。此時(shí)在工作環(huán)境中使用B命令就可調(diào)出B矩陣。5再試著輸入一些矩陣，

2、矩陣中的元素可為任意表達(dá)式，但注意矩陣中各行各列的元素個(gè)數(shù)需分別相等，否則會(huì)給出出錯(cuò)信息。6輸入who和whos觀察結(jié)果，了解其作用。二基本信號(hào)運(yùn)算1數(shù)組的加減乘除和乘方運(yùn)算 輸入A=1 2 3，B=4 5 6，求C=A+B，D=A-B，E=A.*B，F(xiàn)=A./B，G=A.B,并用stem畫出A，B，C，D，E，F(xiàn)，G。再輸入一些數(shù)組，進(jìn)行類似運(yùn)算。 2粗略描繪下列各函數(shù)的波形（其中對于連續(xù)信號(hào)可取時(shí)間間隔為0.001）a.syms tft=4+exp(-t)ezplot(ft,0,0.001,1000)b.syms tft=2*exp(-2*t)+3*exp(-3*t)ezplot(ft,0

3、,0.001,10)c.syms t pift=cos(2*pi*t+pi/3)ezplot(ft,-4,0.001,4)d. syms t pift=exp(-t)*sin(2*pi*t)ezplot(ft,0,0.001,3)e function f=et(t)f=(t=0)t=(-7:0.001:7);ft=et(t+3)-et(t-3);plot(t,f)f. t = -4*pi:0.001:4*pi;ft = sinc(t/pi);plot (t,ft) 實(shí)驗(yàn)內(nèi)容二1 用MATLAB實(shí)現(xiàn)6個(gè)基本的離散序列。1.單位采樣序列：function y=deerta(n)x=length(n

4、)for k=1:x if n(k)=0 y(k)=0; else n(k)=0 y(k)=1 endendn=-25:25y=deerta(n)stem(n,y)2. 單位階躍序列：function y=ibuclu(n)x=length(n)for k=1:x if n(k)=0 y(k)=1 endendn=-25:25y=ibuclu(n)stem(n,y) 3實(shí)指數(shù)序列：a為實(shí)數(shù) n=-25:25a=2y=a.nstem(n,y)4復(fù)數(shù)指數(shù)序列：n=-2:0.1:2sigema=4w=5y=exp(sigema+j*w)*n)stem(n,y)5正余弦序列：， x=-10:10y=x

5、n=3y1=x+nsubplot(2,1,1),stem(x,y)subplot(2,1,2),stem(x,y1)6周期序列： n=-2:0.1:2sita=pi/3w=5y=cos(w*n+sita)stem(n,y)2 函數(shù)用MATLAB實(shí)現(xiàn)函數(shù)impseq(n0,n1,n2)，使函數(shù)實(shí)現(xiàn)，。該函數(shù)的格式為：Function x,n=impseq(n0,n1,n2)% Generate x(n)=delta(n-n0);n1=n=n2% x,n=impseq(n0,n1,n2)n=n1:n2;x=(n-n0)=02用MATLAB實(shí)現(xiàn)函數(shù)stepseq(n0,n1,n2)，使函數(shù)實(shí)現(xiàn)該函數(shù)

6、的格式為：function x,n=stepseq(n0,n1,n2)%Generate x(n)=(n-n0);n1=n=03用MATLAB實(shí)現(xiàn)下列序列 然后，將f (n)擴(kuò)展f (n+11)為，周期數(shù)為4。n=0:0.1:10n2=0:0.1:80n1=n2+11y=3*cos(0.1*pi*n+1/3*pi)+2*sin(0.5*pi*n)y1=3*cos(0.1*pi*n1+1/3*pi)+2*sin(0.5*pi*n1)subplot(2,1,1),stem(n,y)subplot(2,1,2),stem(n1,y1)4MATLAB中可用算術(shù)運(yùn)算符“+”實(shí)現(xiàn)信號(hào)相加，但f1 (n)和

7、f2 (n)的長度必須相等。如果序列長度不等，或者長度雖然相等但采樣的位置不同，就不能運(yùn)用“+”了。試用MATLAB寫出任意序列相加的函數(shù)sigadd,其格式如下：Functiony, n= sigadd (f1,n1,f2,n2)% 實(shí)現(xiàn) y (n) = f1 (n) + f2 (n) % y, n= sigadd (f1, n1, f2, n2)% y, n=和信號(hào)，即y=在包括n1和n2的n上求和序列% f1,n1=第一個(gè)信號(hào)，其中f1為序列f1 (n)的幅度值，其長度為length(n1)，% n1為f1 (n)的橫坐標(biāo)取值。% f2,n2=第一個(gè)信號(hào)，其中f2為序列f2 (n)的幅度

8、值，其長度為length(n2)，% n2為f2 (n)的橫坐標(biāo)取值（n2 可與n1不等）。n=min(min(n1),min(n2):max(max(n1),max(n2);y1=zeros(1,length(n);y2=y1; y1(find(n=min(n1)&(n=min(n2)&(n=min(n1)&(n=min(n2)&(n=-tao/2)&(t=-tao/2)&(t=-tao/2)&(t=-tao/2)&(t=-tao/2)&(ttao/2) t=fliplr(t) t=rot90(t) n1=0:20 o=2.*pi./T a=2.*ft*cos(t*n1.*o)./T.*0.

9、001 b=2.*ft*sin(t*n1.*o)./T.*0.001 A=(a.2+b.2).0.5 fi=-atan(b./a) t=-T:0.001:Tft0=A(1)./2subplot(5,1,1)plot(t,ft0)axis(-T,T,-0.2,1.2)ft1=ft0+A(2).*cos(o.*t+fi(2)subplot(5,1,2)plot(t,ft1)axis(-T,T,-0.2,1.2)ft2=ft1+A(3).*cos(2.*o.*t+fi(3)subplot(5,1,3)plot(t,ft2)axis(-T,T,-0.2,1.2)ft3=ft2+A(4).*cos(3.

10、*o.*t+fi(4)subplot(5,1,4)plot(t,ft3)axis(-T,T,-0.2,1.2)ft4=ft3+A(5).*cos(4.*o.*t+fi(5)subplot(5,1,5)plot(t,ft4)axis(-T,T,-0.2,1.2)實(shí)驗(yàn)名稱五、零極點(diǎn)分析及系統(tǒng)的頻響特性一、實(shí)驗(yàn)?zāi)康?.掌握系統(tǒng)函數(shù)零極點(diǎn)定義。2.零極點(diǎn)與頻響的關(guān)系。3.零極點(diǎn)與系統(tǒng)穩(wěn)定性的關(guān)系。二、實(shí)驗(yàn)內(nèi)容1.已知傳遞函數(shù)，畫出零極點(diǎn)圖，并判斷系統(tǒng)的穩(wěn)定性。a H(s)=3*(s-1)*(s-2)/(s+1)*(s+2)num=3,-9,6den=1,3,2z,p=tf2zp(num,den)zpl

11、ane(z,p)不穩(wěn)定b.H(s)=1/snum=1den=1,0z,p=tf2zp(num,den)zplane(z,p)穩(wěn)定c.H(s)=(s2+1)/(s2+2*s+5)num=1,0,1den=1,2,5z,p=tf2zp(num,den)zplane(z,p)穩(wěn)定d. H(z)=0.6*(3*z3+2*z2+2*z+5)/(z3+3*z2+2*z+1)num=1.8,1.2,1.2,3den=1,3,2,1z,p=tf2zp(num,den)zplane(z,p)2畫出零極點(diǎn)分布圖，并分析穩(wěn)定性a.y(n)+y(n-1)+0.5*y(n-2)=f(n)num=1,0,0den=1,1

12、,0.5z,p=tf2zp(num,den)zplane(z,p)不穩(wěn)定eH(z)=(1-1.414z-1+z-2)*(z+1)/(1+0.9z-1+0.81z-2)*(z-0.3)num=1,-0.414,1,1den=1.81,0.6,-0.27,-2.43z,p=tf2zp(num,den)zplane(z,p)不穩(wěn)定。2畫出零極點(diǎn)分布圖，并分析穩(wěn)定性a.題目：y(n)+y(n-1)+0.5*y(n-2)=f(n)num=1,0,0den=1,1,0.5z,p=tf2zp(num,den)zplane(z,p)穩(wěn)定b.y(n)-1.25*y(n-1)+0.75*y(n-2)-0.125*

13、y(n-3)=f(n)+0.5*(n-1)num=1,0.5,0,0den=1,-1.25,0.75,-0.125z,p=tf2zp(num,den)zplane(z,p)穩(wěn)定3.畫出零極點(diǎn)圖、單位沖激響應(yīng)圖、幅頻響應(yīng)圖。a.H(s)=2*s/(s2+20.5*s+1)num=2,0den=1,20.5,1z,p=tf2zp(num,den)zplane(z,p)t=0:0.001:5sys=tf(num,den)y=impulse(sys,t)plot(t,y)w=logspace(-1,1)freqs(num,den,w)b.H(s)=3*(s-1)*(s-2)/(s+1)*(s+2)nu

14、m=3,-9,6den=1,3,2z,p=tf2zp(num,den)zplane(z,p)t=0:0.001:5sys=tf(num,den)y=impulse(sys,t)plot(t,y)w=logspace(-1,1)freqs(num,den,w)實(shí)驗(yàn)名稱六、信號(hào)調(diào)制與解調(diào)一、實(shí)驗(yàn)?zāi)康?.了解用MATLAB實(shí)現(xiàn)信號(hào)的調(diào)制與解調(diào)的方法。2.了解幾種基本的調(diào)制方法。二、實(shí)驗(yàn)內(nèi)容1.有一正弦信號(hào)f(n)=sin(2*pi*n/256),n=0:256,分別以100000HZ的載波和1000000HZ的抽樣頻率進(jìn)行調(diào)幅、調(diào)頻、調(diào)相，觀察圖形。n=0:256fn=sin(2*pi*n/256)y1=modulate(fn,100000,1000000,amdsb-sc)y2=modulate(fn,100000,1000000,fm)y3=modulate(fn,100000,1000000,pm)subplot(3,1,1)plot(n,y1)subplot(3,1,2)plot(n,y2)subplot(3,1,3)plot(