數(shù)字信號處理實驗答案

1、實驗一熟悉Matlab環(huán)境一、實驗?zāi)康? .熟悉MATLAB的主要操作命令。2 .學(xué)會簡單的矩陣輸入和數(shù)據(jù)讀寫。3 .掌握簡單的繪圖命令。4 .用MATLAB編程并學(xué)會創(chuàng)建函數(shù)。5.觀察離散系統(tǒng)的頻率響應(yīng)。二、實驗內(nèi)容認真閱讀本章附錄，在MATLAB環(huán)境下重新做一遍附錄中的例子，體會各條命令的含義。在熟悉了 MATLAB基本命令的基礎(chǔ)上，完成以下實驗。上機實驗內(nèi)容：(1)數(shù)組的加、減、乘、除和乘方運算。輸入 A=1 2 3 4 , B=3 4 5 6,求C=A+B , D=A-B , E=A.*B , F=A./B, G=AAB 并用 stem語句畫出 A、B、C、D、E、F、G。clear

2、all;a=1 2 3 4;b=3 4 5 6;c=a+b;d=a-b;e=a.*b;f=a./b;g=a-b；n=1:4;subplot(4,2,1);stem(n,a);xlabel('n');xlim(0 5);ylabel('A');subplot(4,2,2);stem(n,b);xlabel('n');xlim(0 5);ylabel('B');subplot(4,2,3);stem(n,c);xlabel('n');xlim(0 5);ylabel('C');subplot(4,2,4

3、);stem(n,d);xlabel('n');xlim(0 5);ylabel('D');subplot(4,2,5);stem(n,e);xlabel('n');xlim(0 5);ylabel('E');subplot(4,2,6);stem(n,f);xlabel('n');xlim(0 5);ylabel('F');subplot(4,2,7);stem(n,g);xlabel('n');xlim(0 5);ylabel('G');(2)用MATLAB實現(xiàn)下

4、列序列：a) x(n)=0.8n0<n< 15b x(n)=e(0.2+3j)n0< n< 15c) x(n)=3cos(0.125 兀 n+0.2 兀)+2sin(0.25 兀 n+0.1 兀)0<n<15d)將c)中的x(n)擴展為以16為周期的函數(shù)x16(n)=x(n+16)，繪出四個周期。e)將c)中的x(n)擴展為以10為周期的函數(shù)x10(n)=x(n+10),繪出四個周期。clear all;N=0:15;% a) x(n)=0.8n0< n<15xa=0.8，N;040 rLJJ20QL 0 5000 r0L0- 6 -figure

5、;subplot(2,1,1);stem(N,xa); xlabel('n');xlim(0 16);ylabel('xa');% b) x(n尸e(0.2+3j)n0<n< 15xb=exp(0.2+3*j)*N);subplot(2,1,2);stem(N,xb);xlabel('n');xlim(0 16);ylabel('xb');figure;% c) x(n)=3cos(0.125 兀 n+0.2 兀)+2sin(0.25 兀 n+0.1 兀)0< n< 15xc=3*cos(0.125*pi*

6、N+0.2*pi)+2*sin(0.25*pi*N+0.1*pi);subplot(3,1,1);stem(N,xc);xlabel('n');xlim(0 16);ylabel('xc');% d)將c)中的x(n)擴展為以16為周期的函數(shù)x16(n)=x(n+16),繪出四個周期 k=0:3;m=0;for i=1:4for j=1:16 m=m+1;n(m)=N(j)+16*k(i);x16(m)=3*cos(0.125*pi*n(m)+0.2*pi)+2*sin(0.25*pi*n(m)+0.1*pi); endendsubplot(3,1,2);ste

7、m(n,x16);xlabel('n');ylabel('x16');% e)將c)中的x(n)擴展為以10為周期的函數(shù)x10(n)=x(n+10),繪出四個周期 for j=1:10x10(j)=x16(j);endfor i=1:3for m=1:10x10(i*10+m)=x10(m); endendn=1:40;subplot(3,1,3);stem(n,x10); xlabel('n');ylabel('x10');1510IIII1II1.c.Q>99?G 八).066646A "IIIIIIIIII1

8、10246610121416(3) x(n)=1,-1,3,5,產(chǎn)生并繪出下列序列的樣本 a) xi(n)=2x(n+2)-x(n-1)-2x(n)>II5b) x2(n) = " nx(n - k)k 1clear all n=1:4; T=4; x=1 -1 3 5; x(5:8)=x(1:4);1234subplot(2,1,1);stem(1:8,x);grid;ILcn廠I>for i=1:4 if i-1<0x1(i)=2*x(i+2)-x(i-1)-2*x(i); elsex1(i)=2*x(i+2)-x(i-1+T)-2*x(i); end end

9、x1(5:8)=x1(1:4);subplot(2,1,2);stem(1:8,x1);grid;(4)繪出下列時間函數(shù)的圖形，對 x軸、y軸以及圖形上方均須加上適當(dāng)?shù)臉?biāo)注: a) x(t)=sin(2 兀 t) 0<t< 10sb) x(t)=cos(100 兀 t)sin(兀 t)0<t<4sta=0:0.05:10;xa=sin(2*pi*ta);subplot(2,1,1);plot(ta,xa);xlabel('t');ylabel('幅度')； tb=0:0.01:4;xb=cos(100*pi*tb).*sin(pi*tb)

10、;subplot(2,1,2);plot(tb,xb);xlabel('t');ylabel('幅度')；(5)編寫函數(shù)stepshift(n0,n1,n2/現(xiàn)u(n-n0),n1<n0<n2,繪出該函數(shù)的圖形，起點為 n1, 終點為n20n0=5;ns=1;nf=10;%ns 為起點;nf為終點；在=n=n0處生成單位階躍序列n=ns:nf;x=(n-n0)>=0;stem(n,x);(6)給一定因果系統(tǒng)H(z) =(1 +<,2z-1 +1)/(1 -0.67z -1 +0.9z-2)求出并繪制H(z)的幅頻響應(yīng)與相頻響應(yīng)。clear

11、 all;b=1,sqrt(2),1;a=1,-0.67,0.9;h,w=freqz(b,a);am=20*log10(abs(h);subplot(2,1,1);plot(w,am);ph=angle(h);subplot(2,1,2);plot(w,ph);(7)計算序列8 -2 -1 2 3和序列2 3 -1-3的離散卷積，并作圖表示卷積結(jié)果clear all;a=8 -2 -1 2 3;b=2 3 -1-3;c=conv(a,b); % 計算卷積 M=length(c)-1;n=0:1:M;stem(n,c);xlabel('n');ylabel('幅度'

12、;);(8)求以下差分方程所描述系統(tǒng)的單位脈沖響應(yīng)h(n),0< n0 50y(n)+0.1y(n-1)-0.06y(n-2)=x(n)-2x(n-1)clear all;N=50;a=1 -2;b=1 0.1-0.06;x=1 zeros(1,N-1);k=0:1:N-1;y=filter(a,b,x);stem(k,y);xlabel('n');ylabel('幅度');實驗二 信號的采樣與重建一，實驗?zāi)康? 1)通過觀察采樣信號的混疊現(xiàn)象，進一步理解奈奎斯特采樣頻率的意義。( 2)通過實驗，了解數(shù)字信號采樣轉(zhuǎn)換過程中的頻率特征。( 3)對實際的音頻文

13、件作內(nèi)插和抽取操作，體會低通濾波器在內(nèi)插和抽取中的作用。二，實驗內(nèi)容(1)采樣混疊，對一個模擬信號 Va(t)進行等間采樣，采樣頻率為 200H乙得 到離散時間信號 V(n).Va (t)由頻率為30Hz,150Hz,170Hz,250Hz,330Hz的5個 正弦信號的加權(quán)和構(gòu)成。Va(t)=6cos(60pi*t)+3sin(300pi*t)+2cos(340pi*t)+4cos(500pi*t)+10sin(660pi*t) 觀察采樣后信號的混疊效應(yīng)。程序： clear, close all t=0:0.1:20;Ts=1/2;n=0:Ts:20;V=8*cos(0.3*pi*t)+5*c

14、os(0.5*pi*t+0.6435)-10*sin(0.7*pi*t);Vn=8*cos(0.3*pi*n)+5*cos(0.5*pi*n+0.6435)-10*sin(0.7*pi*n);subplot(221)plot(t,V),),grid on , subplot(222) stem(n,Vn,- 7 -2007 / I 1-20i'X-404005101520grid on,(2)輸入信號X(n)為歸一化頻率f1=0.043, f2=0.31的兩個正弦信號相加而成, N=100,按因子M=2作抽取：(1)不適用低通濾波器；(2)使用低通濾波器。 分別顯示輸入輸出序列在時域和

15、頻域中的特性。程序：clear;N=100;M=2;f1=0.043;f2=0.31;n=0:N-1;x=sin(2*pi*f1*n)+sin(2*pi*f2*n);y1=x(1:2:100);y2=decimate(x,M, 'fir' );figure(1);stem(n,x(1:N);title('input sequence' );xlabel( 'n' );ylabel('fudu');figure(2);n=0:N/2-1;stem(n,y1);title( 'output sequence without L

16、P');xlabel( 'n' );ylabel( 'fudu');figure(3);m=0:N/M-1;stem(m,y2(1:N/M);title( 'output sequence with LP');xlabel( 'n' );ylabel( 'fudu');figure(4);h,w=freqz(x);- 10 -plot(w(1:512),abs(h(1:512);title( 'frequency spectrum of the input sequence'xlabel(

17、'w' );ylabel( 'fudu');figure(5)；h,w=freqz(y1);plot(w(1:512),abs(h(1:512);title( 'frequency spectrum of the output sequence without LP' xlabel( 'w' );ylabel( 'fudu');figure(6);h,w=freqz(y2);plot(w(1:512),abs(h(1:512);title( 'frequency spectrum of the output

18、 sequence without LP' xlabel( 'w' );ylabel( 'fudu'););););input sequence0102030405060708090100output sequence without LP05101520253035404550noutput sequence with LP1.50.5-0.5-1-1.505101520253035404550frequency spectrum of the output sequenceu d00.511.522.533.5wu dfrequency spectr

19、um of the output sequence without LPw-11 -25frequency spectrum of the output sequence without LP(3)輸入信號X(n)為歸一化頻率f1=0.043, f2=0.31的兩個正弦信號相加而成，長度N=50,內(nèi)插因子為2. (1)不適用低通濾波器；(2)使用低通濾波器。分別 顯示輸入輸出序列在時域和頻域中的特性。程序：clear,close all , N=50;L=2;f1=0.043;f2=0.31;n=0:N-1;x=sin(2*pi*f1*n)+sin(2*pi*f2*n);figure(1);s

20、tem(n,x(1:N);title('input sequence' );xlabel( 'n' );ylabel('fudu');y1=zeros(1,N*2);y1(1:2:N*2)=x;figure(2);m=0:N*L-1;stem(m,y1(1:N*L);title( 'output sequence ');-12 -xlabel( 'n' );ylabel( 'fudu');y2=interp(x,L);figure(3);m=0:N*L-1;stem(m,y2(1:N*L);titl

21、e( 'output sequence' );xlabel( 'n' );ylabel( 'fudu');figure(4);h,w=freqz(x);plot(w(1:512),abs(h(1:512);title( 'frequency spectrum of the input sequence');xlabel( 'w' );ylabel( 'fudu');figure(5);h,w=freqz(y1);plot(w(1:512),abs(h(1:512);title( 'frequ

22、ency spectrum of the output sequence');xlabel( 'w' );ylabel( 'fudu');figure(6);h,w=freqz(y2);plot(w(1:512),abs(h(1:512);title( 'frequency spectrum of the output sequence ');xlabel( 'w' );ylabel( 'fudu');- 13 -output sequenceoutput sequence32.521.510.50-0.5

23、-1-1.5-2n102030405060708090100-15 -frequency spectrum of the input sequence302520ud 151030ud 150.511.522.533.5frequency spectrum of the output sequence252010500.512.51.52w33.5-17 -frequency spectrum of the output sequence605040ud 3020111! ；10 000.511.522.533.5w二.(3)令x(n尸cos(2*pi*f*n/fs),其中f/fs=1/16,

24、即每個周期內(nèi)有16個點。試用MATLAB編程實現(xiàn)：1)作M=4倍的抽取，使每個周期變成4點。程序：clear,close all,N=100;M=4;n=0:N-1;x=cos(2*pi*n*(1/16);stem(n,x(1:N);title('input sequence');xlabel('n');ylabel('fudu');-19 -y1=x(1:4:100);y2=decimate(x,M,'fir');figure(2);m=0:N/4-1;stem(m,y1);title('output sequence

25、');xlabel('n');ylabel('fudu');figure(3);m=0:N/M-1;stem(m,y2(1:N/M);title('output sequence');xlabel('n');ylabel('fudu');figure(4);h,w=freqz(x);plot(w(1:512),abs(h(1:512);title('frequency spectrum of the input sequence');xlabel('w');ylabel(&#

26、39;fudu');figure(5);h,w=freqz(y1);plot(w(1:512),abs(h(1:512);title('frequency spectrum of the output sequence');xlabel('w');ylabel('fudu');figure(6);h,w=freqz(y2);plot(w(1:512),abs(h(1:512);title('frequency spectrum of the output sequence ');xlabel('w');yl

27、abel('fudu');- 21 -output sequenceoutput sequencen-22 -frequency spectrum of the output sequence-24 -3J " ! I3.5300.511.522.53.55040ud 302010000.511.522.5wfrequency spectrum of the output sequence1412108u u6420w2).作L=3倍的插值，使每個周期變成48點 程序：clear,close all,N=50;L=3;n=0:N-1;x=cos(2*pi*n*(1/1

28、6);figure(1);stem(n,x(1:N);title('input sequence');xlabel('n');ylabel('fudu');y1=zeros(1,N*3);y1(1:3:N*3)=x;figure(2);m=0:N*3-1;stem(m,y1(1:N*3);title('output sequence ');xlabel('n');ylabel('fudu');y2=interp(x,L);figure(3);m=0:5:N*L-1;stem(m,y2(1:5:N*L

29、);title('output sequence');xlabel('n');ylabel('fudu');figure(4);h,w=freqz(x);plot(w(1:512),abs(h(1:512);title('frequency spectrum of the input sequence');xlabel('w');ylabel('fudu');figure(5);h,w=freqz(y1);plot(w(1:512),abs(h(1:512);title('frequency

30、 spectrum of the output sequence');xlabel('w');ylabel('fudu');figure(6);h,w=freqz(y2);plot(w(1:64),abs(h(1:64);title('frequency spectrum of the output sequence ');xlabel('w');ylabel('fudu');input sequence10.80.60.40.20-0.2-0.4-0.6-0.8-11015202530354045n50-

31、 25 -output sequence1uaMr0.80.60.40.2-0.2-0.4-0.6-0.8-1 0一.n-r-50'l-i- 'I-l-LU-J-LXJl上100150output sequencen-27 -frequency spectrum of the input sequence300.511.522.5w33.5520 5 02 115uabfrequency spectrum of the output sequence0.512.533.5520 5 02 11530uaMT1.52w-29 -frequency spectrum of the

32、output sequence201000.050.10.150.20.25w0.30.350.40-30 -(4).輸入信號x(n)為歸一化頻率分別是 f1=0.04,f2=0.3的正弦信號相加而成，N=50，內(nèi)插因子為5,抽取因子 為3,給出按有理因子5/3做采樣率轉(zhuǎn)換的輸入輸出波形。程序：clear,close all,N=50;M=3;L=5;f1=0.04;f2=0.3;n=0:N-1;x=sin(2*pi*f1*n)+sin(2*pi*f2*n);y=resample(x,L,M);figure(1);stem(n,x(1:N);title('input sequence&

33、#39;);xlabel('n');ylabel('fudu');figure(2);m=0:N-1;stem(m,y(1:N);title('output sequence ');xlabel('n');ylabel('fudu');figure(3);h,w=freqz(x);plot(w(1:512),abs(h(1:512);title('frequency spectrum of the input sequence');xlabel('w');ylabel('fu

34、du');figure(4);h,w=freqz(y);plot(w(1:512),abs(h(1:512);title('frequency spectrum of the output sequence ');xlabel('w');ylabel('fudu');input sequence- 31 -output sequencenfrequency spectrum of the input sequencew-32 -frequency spectrum of the output sequence454035302520151

35、05000.511.522.533.5w-34 -實驗三快速Fourier變換(FFT)及其應(yīng)用一、實驗?zāi)康? .在理論學(xué)習(xí)的基礎(chǔ)上，通過本實驗，加深對 FFT的理解，熟悉FFT子程序。2 .熟悉應(yīng)用FFT對典型信號進行頻譜分析的方法。3 . 了解應(yīng)用FFT進行信號頻譜分析過程中可能出現(xiàn)的問題以便在實際中正確應(yīng)用FFT4 .熟悉應(yīng)用FFT實現(xiàn)兩個序列的線性卷積的方法。5,初步了解用周期圖法作隨機信號譜分析的方法。二、實驗原理與方法在各種信號序列中，有限長序列信號處理占有很重要地位，對有限長序列，我們可以使 用離散Fouier變換(DFT)。這一變換不但可以很好的反映序列的頻譜特性，而且易于用快

36、速 算法在計算機上實現(xiàn)，當(dāng)序列x(n)的長度為N時，它的DFT定義為：21_匹陽上)=工卷)磴,乙 =圖加反變換為：1 Ml伽)=WZx伏卿?有限長序列的DFT是其Z變換在單位圓上的等距采樣，或者說是序列 Fourier變換的等距 采樣，因此可以用于序列的譜分析。FFT并不是與DFT不同的另一種變換，而是為了減少 DFT運算次數(shù)的一種快速算法。 它是對變換式進行一次次分解，使其成為若干小點數(shù)的組合，從而減少運算量。常用的 FFT是以2為基數(shù)的，其長度 N =它的效率高，程序簡單，使用非常方便，當(dāng)要變換的序列長度不等于2的整數(shù)次方時，為了使用以2為基數(shù)的FFT,可以用末位補零的方法，使其 長度延

37、長至2的整數(shù)次方。(一)、在運用DFT進行頻譜分析的過程中可能產(chǎn)生三種誤差：(1) 混登序列的頻譜時被采樣信號的周期延拓，當(dāng)采樣速率不滿足Nyquist定理時，就會發(fā)生頻譜混疊，使得采樣后的信號序列頻譜不能真實的反映原信號的頻譜。避免混疊現(xiàn)象的唯一方法是保證采樣速率足夠高，使頻譜混疊現(xiàn)象不致出現(xiàn)，即在確定采樣 頻率之前，必須對頻譜的性質(zhì)有所了解，在一般情況下，為了保證高于折疊頻率的分量不會 出現(xiàn)，在采樣前，先用低通模擬濾波器對信號進行濾波。(2)泄漏實際中我們往往用截短的序列來近似很長的甚至是無限長的序列，這樣可以使用較短的DFT來對信號進行頻譜分析，這種截短等價于給原信號序列乘以一個矩形窗函

38、數(shù)，也相當(dāng)于 在頻域?qū)⑿盘柕念l譜和矩形窗函數(shù)的頻譜卷積，所得的頻譜是原序列頻譜的擴展。泄漏不能與混疊完全分開，因為泄漏導(dǎo)致頻譜的擴展，從而造成混疊。為了減少泄漏的影響， 可以選擇適當(dāng)?shù)拇昂瘮?shù)使頻譜的擴散減至最小。(3)柵欄效應(yīng)DFT是對單位圓上Z變換的均勻采樣，所以它不可能將頻譜視為一個連續(xù)函數(shù)， 就一定意 義上看，用DFT來觀察頻譜就好像通過一個柵欄來觀看一個圖景一樣，只能在離散點上看到 真實的頻譜，這樣就有可能發(fā)生一些頻譜的峰點或谷點被“尖樁的柵欄”所攔住，不能別我 們觀察到。減小柵欄效應(yīng)的一個方法就是借助于在原序列的末端填補一些零值，從而變動DFT的點數(shù)，這一方法實際上是人為地改變了對真

39、實頻譜采樣的點數(shù)和位置， 相當(dāng)于搬動了每一根“尖 樁柵欄”的位置，從而使得頻譜的峰點或谷點暴露出來。(二卜用FFT計算線性卷積用FFT可以實現(xiàn)兩個序列的圓周卷積。在一定的條件下，可以使圓周卷積等于線性卷積。一般情況，設(shè)兩個序列的長度分別為 N1和N2,要使圓周卷積等于線性卷積的充要條件是 FFT 的長度N>N1 + N2對于長度不足N的兩個序列，分別將他們補零延長到 No當(dāng)兩個序列中有一個序列比較長的時候，我們可以采用分段卷積的方法。有兩種方法： 重疊相加法。將長序列分成與短序列相仿的片段，分別用FFT對它們作線性卷積，再將分段卷積各段重疊的部分相加構(gòu)成總的卷積輸出。重疊保留法。這種方法

40、在長序列分段時，段與段之間保留有互相重疊的部分，在構(gòu)成總的卷 積輸出時只需將各段線性卷積部分直接連接起來，省掉了輸出段的直接相加。(三卜用周期圖法(平滑周期圖的平均法)對隨機信號作譜分析實際中許多信號往往既不具有有限能量，由非周期性的。無限能量信號的基本概念是隨機 過程，也就是說無限能量信號是一隨機信號。周期圖法是隨機信號作譜分析的一種方法，它 特別適用于用FFT直接計算功率譜的估值。將長度為N的實平穩(wěn)隨機序列的樣本x(n)再次分割成K段，每段長度為L,即L=N/K。每 段序列仍可表示為：- 47 -xi(n)=x(n+(i-1)L) , 0<n<L-1 , 1<i<K

41、但是這里在計算周期圖之前，先用窗函數(shù)w(n)給每段序列xi(n)加權(quán)，K個修正的周期圖定義為小/)LU£-1I z孫)M = 0其中u表示窗口序列的能量，它等于1 G儲L閘=口在此情況下，功率譜估計量可表小為(口)MLia)、實驗內(nèi)容及步驟實驗中用到的信號序列:Gaussian 序歹!Jb)衰減正弦序列Xh(n) =e-叫iiiZTtftK。<n<150, elsec)三角波序列+ 1,0門 38 -/ 4 < n < 70,如d)反三角波序列4 -凡0 < n <3 n - 3S4 « 打 « 7 Ci, else上機實驗內(nèi)容

42、：(1)、觀察高斯序列的時域和幅頻特性，固定信號 等于2, 4, 8,觀察它們的時域和幅頻特性，了解當(dāng) 特性的影響；固定q=8,改變p,使p分別等于8,xa(n)中參數(shù)p=8,改變q的值，使q分別q取不同值時，對信號序列的時域幅頻13, 14,觀察參數(shù)p變化對信號序列的時域及幅頻特性的影響，觀察p等于多少時，會發(fā)生明顯的泄漏現(xiàn)象，混疊是否也隨之出現(xiàn)? 記錄實驗中觀察到的現(xiàn)象，繪出相應(yīng)的時域序列和幅頻特性曲線。% %定義高斯序列% function Xa,Fa =gauss(p,q)% n=0:15;% Xa(n+1)=exp(-(n+1-p)，2./q);% F=fft(Xa);% Fa=ab

43、s(F);clear all;% p=8,q=2 %Xa1,Fa1= gauss(8,2);k=0:15;subplot(5,2,1);plot(k,Xa1);xlabel('n');ylabel(' 時域特性');text(10,0.5,'p=8,q=2');subplot(5,2,2);plot(k,Fa1);xlabel('n');ylabel(' 幅頻特性');text(8,3,'p=8,q=2');% p=8,q=4 %Xa2,Fa2= gauss(8,4);subplot(5,2,3);

44、plot(k,Xa2);xlabel('n');ylabel(' 時域特性');text(10,0.5,'p=8,q=4');subplot(5,2,4);plot(k,Fa2);xlabel('n');ylabel(' 幅頻特性');text(8,3,'p=8,q=4');% p=8,q=8 %Xa3,Fa3= gauss(8,8);subplot(5,2,5);plot(k,Xa3);xlabel('n');ylabel(' 時域特性');text(10,0.5,

45、'p=8,q=8');subplot(5,2,6);plot(k,Fa3);xlabel('n');ylabel(' 幅頻特性');text(8,3,'p=8,q=8');% p=13,q=8 %Xa4,Fa4= gauss(13,8);subplot(5,2,7);plot(k,Xa4);xlabel('n');ylabel(' 時域特性');text(10,0.5,'p=13,q=8');subplot(5,2,8);plot(k,Fa4);xlabel('n')

46、;ylabel(' 幅頻特性');text(8,3,'p=13,q=8');% p=14,q=8 %Xa5,Fa5= gauss(14,8);subplot(5,2,9);plot(k,Xa5);xlabel('n');ylabel(' 時域特性');text(10,0.5,'p=14,q=8');subplot(5,2,10);plot(k,Fa5);xlabel('n');ylabel(' 幅頻特性');text(8,3,'p=14,q=8');o - 一 - o

47、 .5o 1 so D. 翅-博君 知里翼之s. n.5 D.(2)、觀察衰減正弦序列xb(n)的時域和幅頻特性，a=0.1, f=0.0625,檢查譜峰出現(xiàn)位置是否正確，注意頻譜的形狀，繪出幅頻特性曲線，改變 f,使f分別等于0.4375和0.5625,觀 察這兩種情況下，頻譜的形狀和譜峰出現(xiàn)位置，有無混疊和泄漏現(xiàn)象？說明產(chǎn)生現(xiàn)象的原因。%定義衰減正弦序列% function Xb,Fb = downsin(a,f)% n=0:15;% Xb(n+1)=exp(-a.*n).*sin(2*pi*f.*n); %自然對數(shù)的底：e=:2.71828 18284 59045 23536% F =

48、fft(Xb);% Fb=abs(F);clear all;k=0:15;% a=0.1,f=0.0.0625 %Xb,Fb=downsin(0.1,0.0625);subplot(3,2,1); plot(k,Xb);xlabel('n');ylabel('時域特性');text(8,0.5,'a=0.1,f=0.0625');subplot(3,2,2); plot(k,Fb);xlabel('n');ylabel('幅值特性');text(10,3,'a=0.1,f=0.0625');% a=

49、0.1,f=0.4375 %Xb1,Fb1=downsin(0.1,0.4375);subplot(3,2,3); plot(k,Xb1);xlabel('n');ylabel('時域特性');text(8,0.5,'a=0.1,f=0.4375');subplot(3,2,4); plot(k,Fb1);xlabel('n');ylabel('幅值特性');text(10,3,'a=0.1,f=0.4375');% a=0.1,f= 0.5625 %Xb2,Fb2=downsin(0.1,0.56

50、25);subplot(3,2,5); plot(k,Xb2);xlabel('n');ylabel('時域特性');text(8,0.5,'a=0.1,f=0.5625');subplot(3,2,6); plot(k,Fb2);xlabel('n');ylabel('幅值特性');text(10,3,'a=0.1,f=0.5625');O 石。 石 0.-0.£O .5口.5( 0.4.5 0 5 0.肛 犯*w盆5J4 3 2.(3)、觀察三角波和反三角波序列的時域和幅頻特性，用 N

51、=8點FFT分析信號序列xc(n)和 xd(n)的幅頻特性，觀察兩者的序列形狀和頻譜曲線有什么異同？繪出兩序列及其幅頻特性曲 線。在xc(n)和xd(n)末尾補零，用N=16點FFT分析這兩個信號的幅頻特性，觀察幅頻特性 發(fā)生了什么變化？兩情況的FFT頻譜還有相同之處嗎？這些變化說明了什么？clear all;n=0:3;k=1:8;%定義三角波序列Xc(n+1) = n;Xc(n+5) =4-n;%定義反三角波序列Xd(n+1) = 4-n;Xd(n+5) =n;% 三角波特性 %subplot(2,2,1);plot(k-1,Xc);xlabel('n');ylabel(&

52、#39;時域特性');text(1,3,'三角波);subplot(2,2,2);plot(k-1,abs(fft(Xc);xlabel('k');ylabel('幅頻特性');text(4,10,'三角波);% 反三角波特性 %subplot(2,2,3);plot(k-1,Xd);xlabel('n');ylabel('時域特性');text(3,3,'反三角波);subplot(2,2,4);plot(k-1,abs(fft(Xd);xlabel('k');ylabel(

53、9;幅頻特性');text(4,10,'反三角波);%末尾補0,計算32點FFTXc(9:32)=0;Xd(9:32)=0;k=1:32;figure;% 三角波特性 %subplot(2,2,1);plot(k-1,Xc);xlabel('n');ylabel('時域特性');text(1,3,'三角波);subplot(2,2,2);plot(k-1,abs(fft(Xc);xlabel('k');ylabel('幅頻特性');text(4,10,'三角波')；% fan 三角波特性 %subplot(2,2,3);plot(k-1,Xd);xlabel('n');ylabel('時域特性');text(3,3,'反三角波);subplot(2,2,4);plot(k-1,abs(fft(Xd);xlabel('k');ylabel('幅頻特性');text(4,10,'反三角波')；5 0 5 0 as®40 3r 20 o1203041:14 3 2 10擔(dān)t篇fc40302nno1(4)、一個連續(xù)信號含兩個頻率分量，經(jīng)采樣得x(n尸sin