數(shù)字信號處理實驗報告(基于MATLAB)

1、課程名稱：數(shù)字信號處理實驗實驗地點：綜合樓C407專業(yè)班級： 2014 級生物醫(yī)學(xué)工程姓名： leifeng學(xué)號： 指導(dǎo)老師：第一次實驗第一章 離散時間信號的時域分析Q1.1運行程序P1.1,以產(chǎn)生單位樣本序列un并顯示它clf； n=-10:20;u=zeros(1,10) 1 zeros(1,20);stem(n,u);xlabel('時間序號');ylabel( '振幅)； title( '單位樣本序列)； axis(-10 20 0 1.2);10.8幅0.6振0.40.20-10-505101520時間序號單位樣本序列Q1.2 命令clf, axis,

2、 title , xlabel和ylabel的作用是什么clf：清除圖形窗口內(nèi)容；axis:規(guī)定橫縱坐標(biāo)的范圍；title :使圖像面板上方顯示相應(yīng)的題目名稱;xlable：定義橫坐標(biāo)的名字；ylable :定義縱坐標(biāo)的名字。Q1.3修改程序P1.1以產(chǎn)生帶有延時11個樣本的延遲單位樣本序列udn,運行修改的程序并且顯示產(chǎn)生的序列。clf；n=0:30;u=zeros(1,11) 1 zeros(1,19);stem(n,u);xlabel(' 時間序號);ylabel(' 振幅')；title('單位樣本序列)；axis(0 30 0 1.2);單位樣本序列1

3、 .0.8、幅 0.6 . 振0.4 -0.2 0 . F 1一 1一 1一 1. I ,一 L ，一 匚 一.一，一 L 出口 0 n - 鼻 1051015202530時間序號Q1.5修改程序P1.1,以產(chǎn)生帶有超前7個樣本的延時單位階躍序列sdn。運行修改后的程序并顯示產(chǎn)生的序列。clf;n=-10:20;sd=zeros(1,3) 1 ones(1,27);stem(n,sd);xlabel('時間序號');ylabel( '振幅); title( '單位樣本序列)； axis(-10 20 0 1.2);Q1.6運行程序P1.2,以產(chǎn)生復(fù)數(shù)值的指數(shù)序列

4、。clf；c=-(1/12)+(pi/6)*i;K=2;n=0:40;x=K*exp(c*n);subplot(2,1,1);stem(n,real(x);xlabel('時間序號 n' );ylabel( '振幅');title( 實部)；subplot(2,1,2);stem(n,imag(x);xlabel('時間序號 n' );ylabel( '振幅');title( 虛部)；實部Q1.7哪個參數(shù)控制該序列的增長或衰減率？哪個參數(shù)控制該序列的振幅?參數(shù)C控制序列增長或衰減率參數(shù)K控制序列的振幅參數(shù)C控制該序列的增長或衰減率

5、；參數(shù)K控制該序列的振幅。Q1.8若參數(shù)C更改為(1/12)+(pi/6)*i ,將會發(fā)生什么情況？50dlh-50 L010【rrr152025時間序號n虛部30354010050-50 L0r10152025時間序號n303540Q1.9運算符real和imag的作用是什么？運算符real和imag的作用分別是提取運算數(shù)值x實部和虛部。Q1.10命令subplot的作用是什么？subplot (a, b, c),其中a代表圖像分a行顯示，b代表圖像分b列顯示，c代表第a行的第 b列圖像。Q1.11運行程序P1.3以產(chǎn)生實數(shù)值的指數(shù)序列clf；n=0:35;a=1.2;K=0.2;x=K*a

6、.An;stem（n,x）;xlabel（'時間序號 n' ）;ylabel（ '振幅）;Q1.12哪個參數(shù)控制該序列的增長或衰減率？哪個參數(shù)控制該序列的振幅？參數(shù)a控制該序列的增長或者衰減率， 參數(shù)k控制該序列的振幅。Q1.13算術(shù)運算符A和A之間的區(qū)別是什么？運算符 表示a的指數(shù)是n這個序列；運算符“人”表示a的指數(shù)分別是n這個序列中的每 一個數(shù)字，即對應(yīng)的每一個元素。Q1.14若參數(shù)a小于1,會發(fā)生什么情況？將參數(shù) a更改為0.9,將參數(shù)k更改成20,再次運行程序P1.3。clf；n=0:35;a=0.9;K=20;x=K*a.An;stem（n,x）;xlabe

7、l（'時間序號 n' ）;ylabel（ '振幅）；Q1.1620余Q1.15該序列的長度是多少？怎樣才能改變它？該序列的長度是35;通過改變“ n=0:35”中的“ 35”這個數(shù)字可以改變序列的長度。Q1.16 使用MATLAB命令sum(s.*s)可計算用向量s表示的實數(shù)序列sn的能量。試求在習(xí)題Q1.11和習(xí)題Q1.14中產(chǎn)生的實數(shù)值指數(shù)序列xn的能量。clf；n=0:35;a=0.9;K=20;x=K*a.An;stem(n,x);sum(s.*s)xlabel('時間序號 n' );ylabel( '振幅);結(jié)果：ans =2.3106

8、e+004clf；n=0:35;a=1.2;K=0.2;x=K*a.An;s=stem(n,x);sum(s.*s)xlabel('時間序號 n' );ylabel( '振幅);結(jié)果：ans =2.3719e+004Q1.17運行程序pi.4,以產(chǎn)生正弦序列并顯示它。n=0:40;f=0.1;phase=0;A=1.5;arg=2*pi*f*n-phase;x=A*cos(arg);clf;stem(n,x);axis(0 40 -2 2);grid;title( ' 正弦序列 ' );xlabel(' 時間序號 n' );ylabel(

9、' 振幅 ' );axis;Q1.18 該序列的頻率是多少？怎樣可以改變它？哪個參數(shù)控制該序列的相位？哪個參數(shù)控制該序列的振幅？該序列的周期是多少？該序列的頻率是0.1Hz ;通過改變f的值可以此正弦序列的頻率； 參數(shù)phase控制該序列的初相位；參數(shù)A控制該序列的振幅；該序列的周期是10s。Q1.19 該序列的長度是多少？怎樣可以改變它？該序列的長度是41 個時間單位，通過改變n 的參數(shù)可以改變該序列的長度。Q1.21 axis的作用axi斷口 grid命令的作用是什么是規(guī)定像顯圖示的橫縱坐標(biāo)的范圍； grid 的作用是顯示圖像上面的“網(wǎng)格” 。Q1.22修改程序P1.4,以

10、產(chǎn)生一個頻率為0.9的正弦序列并顯示它。把此序列和習(xí)題Q1.17中產(chǎn)生的序列相比較。修改程序 P1.4,以產(chǎn)生一個頻率為1.1的正弦序列并顯示它。把此序列與 Q1.17 中產(chǎn)生的序列相比較，評價你的結(jié)果。正弦序列正弦序列時正該圖彳與f=0.1時的圖像一樣，因為該正弦序列的最小周期是2pi的，而當(dāng)f=0.9與f=1.1好是f=0.1的整數(shù)倍。Q1.23修改上述程序，以產(chǎn)生長度為50,頻率為0.08,振幅為2.5,相移為90度的一個正弦序列并顯示它。該序列的周期是多少？n=0:50;f=0.08;phase=pi/2;A=2.5;arg=2*pi*f*n-phase;x=A*cos(arg);cl

11、f；stem(n,x);axis(0 40 -3 3);grid；title( '正弦序列)；xlabel('時間序號n');ylabel('振幅)； axis;周期為12.5sQ1.24在程序P1.4中用plot命令彳t替stem命令，運行新程序。新圖形與 Q1.7中產(chǎn)生的圖形 有什么區(qū)別？正弦序列原圖像是離散的，新圖像是連續(xù)的，說明 plot命令是用平滑的線將各點連接起來產(chǎn)生連續(xù)波， 而stem命令則是用各點來產(chǎn)生離散波。Q1.25在程序P1.4中用stairs命令彳弋替stem命令，運行新程序。新圖形與Q1.17和Q1.24中產(chǎn)生的圖形有什么區(qū)別？stai

12、rs函數(shù)是用直線將相應(yīng)的點連接起來，類似于階梯狀。Q1.26編寫一個MATLAB程序，以產(chǎn)生并顯示一個長度為100的隨機信號，該信號在區(qū)間-2,2中均勻分布。clf;x=4*rand(1,100)-2; stem(x);Q1.27編寫一個MATLAB程序，以產(chǎn)生并顯示一個長度為75的高斯隨機信號，該信號正太分布且均值為0,方差為3clf;x=3*randn(1,75); stem(x);8-8 L 01020304050607080Q1.36運行程序P1.7以產(chǎn)生掃頻正弦序列 xnn=0:100;a=pi/2100;b=0;arg=a*n.*n+b*n;x=cos(arg);clf;stem(

13、n,x);axis(0,100,-1.5,1.5);title('掃頻正弦信號); xlabel('時間序號n');ylabel('振幅');grid;axis;掃頻正弦信號時間序號n第二次實驗第二章離散時間系統(tǒng)的時域分析Q2.1對M =2,生成輸入xn= s1n+s2n的輸出信號，輸入 xn的那個分量被改離散時間系統(tǒng)抑制？幅振幅振00答：xn被該離散時間系統(tǒng)抑制的分量為信號#2的高頻分量Q2.2.線性時不變系統(tǒng)變?yōu)閥n = 0.5(xn - xn - 1)后，對輸入xn = s1n+s2n的影響是什么?信號#1信號#221 0-1-2050100幅振2

14、1 0-1-2050100時間序號n輸入信號時間序號n輸出信號210-1-2幅振210-1-2050100050100時間序號n時間序號n答：對于輸入的影響是-該系統(tǒng)現(xiàn)在是一個高通濾波器。它通過高頻率的輸入分量信號#2,而不是低頻分量信號#1.Q2.3 s1n和s2n頻率取其他值結(jié)果。M=4 f1=0.1 f2=0.9信號#121幅 0振-1210-1-2幅 振050時間序號n輸入信號21幅 o振-1-2050時間序號n幅 振050100時間序號n輸出信號2 r-" I Ji ndl I -2 tc050100時間序號n-2M=10 f1=0.1 f2=0.9信號#1信號#2幅振2

15、10-22 1 O-1-250005000幅振幅振2 1 O-1-25000Q2.4修改程序P2.1 ,用一個長度為101、最低頻率為0、最高頻率為0.5的掃頻正弦信號作為輸入信號，計算其輸出信號。你能用該系統(tǒng)對掃頻信號的響應(yīng)來解釋Q2.1和Q2.2的結(jié)果嗎?輸入信號21.51；10.5-0.5-1.5時間序號幅 振幅 振21.510.50-0.5-1-1.5-2050100時間序號nQ2.7運行程序P2.3,對由力口權(quán)輸入得到的 yn與在相同權(quán)系數(shù)下輸出 y1n和y2n相力口得至U的 ytn進(jìn)行比較，這兩個序列是否相等？該系統(tǒng)是線性系統(tǒng)嗎？加權(quán)輸入：a x1n+b x2n的輸出50F 0 振

16、-500510152025303540加權(quán)輸入：a y1n+b y2n 50 ,!LLLLL小甌中<?.7芋0幅 n鼻*。牛&L.*：,1©.一中q9_；£<振 0小二J3.1jg、.二 讓心 -50 LIccE-15X 100510152025303540差信號x 10-15差信號5 r 11：鼻 J I I 13” H幅 o I -119 K. II I I I f0振 a工“1I，.-5 1c二cc0510152025303540時間序號n答：可以看出這加權(quán)輸入與輸出的圖像幾乎是相差無幾的，觀察差信號的圖像，縱坐標(biāo)已達(dá)到 10A(-15),可以看出

17、是相差無幾的，所以該系統(tǒng)是線性系統(tǒng)。Q2.8用三組不同的權(quán)系數(shù) a和b的值以及三組不同的輸入頻率重做Q2.7答：改動后a和b值及輸入頻率和相應(yīng)的運行結(jié)果均在下圖：根據(jù)圖進(jìn)行分析，加權(quán)輸入與輸 出的圖像幾乎是相差無幾的，觀察差信號的圖像，縱坐標(biāo)都達(dá)到10A(-15)或10A(-14),可以看出是相差無幾的。所以系統(tǒng)為線性系統(tǒng)。1.20加權(quán)輸入：a x i n+b x 2n的輸出2.3.-201015200-2001015x 105,ILilJ.-5100-100100-100-220r25303540Ad加權(quán)輸入：a y 1 n+b y 2 n5j；I2025303540差信號E%,',

18、_,I-I.國卬/，fJI11I.，. A10r1520時間序號25303540fTflr4 I10152025303540x 2n的輸出加權(quán)輸入：a x 1 川+b. I,.1110-14x 10If加權(quán)輸入：a y 1 n+b V 2 nJr15a：- BIf2025r303540差信號中占 用 J 由10152025時間序號n303540加權(quán)輸入：a x1n+b x2n的輸出10100-100時間序號nQ2.9初始條件非零時圖像:加權(quán) 輸入：a x1n+b x2n的輸出加權(quán)輸入：a y1n+b y2n50 r10幅振差信號時間序號n答：觀察上圖，看差彳t號，振幅達(dá)到 3,可見加權(quán)輸入與輸

19、出的圖像相差較大，不是線性系統(tǒng)。y1n,y2n 和Q2.11假另一個系統(tǒng)為yn=xnxn-1修改程序P2.3,計算這個系統(tǒng)的輸出序列yn。比較yn和ytn。這兩個序列是否相等？該系統(tǒng)是線性系統(tǒng)嗎？clf;n=0:200;a=2;b=-3;s1=cos(2*pi*0.1*n);s2=cos(2*pi*0.4*n);x=a*s1+b*s2;x1=0 s1 0;x2=0 0 s1;x3=0 s2 0;x4=0 0 s2;x5=0 x 0;x6=0 0 x;y1=x1.*x2;y1=y1(2:202);y2=x3.*x4;y2=y2(2:202);y=x5.*x6;y=y(2:202);yt=a*y1

20、-b*y2;d=y-yt;subplot(3,1,1)stem(n,y);ylabel('振幅')；title(' 加權(quán)輸入： acdot x_1n+bcdot x_2n 的輸出 ');subplot(3,1,2)plot(n,yt);ylabel('振幅')；title(' 加權(quán)輸入： acdot y_1n+bcdot y_2n');subplot(3,1,3)plot(n,d);xlabel(' 時間序號 n')ylabel('振幅')；title(' 差信號 ');加權(quán)輸入：a

21、 x1n+b x2n的輸出幅 振加權(quán)輸入：a y1n+b y2n-5 1c02040608010012014016018020010幅0振020406080100120140160180200時間序號n-10答：觀察上圖，看差信號，可見加權(quán)輸入與輸出的圖像相差較大，不是線性系統(tǒng)。Q2.12運行程序P2.4并比較輸出序列yn和ydn-10。這兩個序列之間有什么關(guān)系？該系統(tǒng)是時不變系統(tǒng)嗎？clf;n=0:40;D=10;a=3.0;b=-2;x=a*cos(2*pi*0.1*n)+b*cos(2*pi*0.4*n);xd=zeros(1,D) x;num=2.2403 2.4908 2.2403;

22、den=1 -0.4 0.75;ic=0 0;y=filter(num,den,x,ic);yd=filter(num,den,xd,ic);d=y-yd(1+D:41+D);subplot(3,1,1)stem(n,y);ylabel('振幅')；title('輸出 yn');grid;subplot(3,1,2)stem(n,yd(1:41);ylabel('振幅')；title('由于延時輸入 xn',num2str(D),'的輸出');grid; subplot(3,1,3)stem(n,d);xlabel

23、('時間序號 n');ylabel('振幅')；title('差值信號');grid;20-40 0f1J . 111i；23.PT1輸出yn0-205101520253035由于延時輸入xn10的輸出50-50 0. . . 一.iL,f 1J1,16.- 1 11 11J 1400510303540152025差值信號-1 00510303540152025時間序號n答：系統(tǒng)為線性時不變系統(tǒng),yd為y時移10后的序歹U。Q2.17考慮另一個系統(tǒng)：yn=nxn+xn-1修改程序P2.4,以仿真上面的系統(tǒng)并確定該系統(tǒng)是否為時不變系統(tǒng)。clf;n=

24、0:40;D=10;a=3.0;b=-2;x=a*cos(2*pi*0.1*n)+b*cos(2*pi*0.4*n);xd=zeros(1,D) x;num=n 1;den=1 0;y=filter(num,den,x);yd=filter(num,den,xd);d=y-yd(1+D:41+D);subplot(3,1,1)stem(n,y);ylabel('振幅')；title('輸出 yn');grid;subplot(3,1,2)stem(n,yd(1:41);ylabel('振幅')；title('由于延時輸入 xn',

25、num2str(D),'的輸出');grid; subplot(3,1,3)stem(n,d);xlabel('時間序號 n');ylabel('振幅')；title('差值信號');grid;輸出yn由于延時輸入xn10的輸出差值信號時間序號n答：差值為零系統(tǒng)為線性時不變系統(tǒng)。第三次實驗Q2.19運行程序P2.5,生成式（2.15）所給離散時間系統(tǒng)的沖激響應(yīng)。運行程序：幅 振沖擊響應(yīng)時間序號Q2.20修改程序 P2.5,產(chǎn)生以下因果線性時不變系統(tǒng)的沖激響應(yīng)的前45個樣本: yn+0.71yn-1-0.46yn-2-0.62yn-

26、3=0.9xn-0.45xn-1+0.35n-2+0.002xn-3修改程序P2.5 ,產(chǎn)生如下因果線性時不變系統(tǒng)的沖擊響應(yīng)的前45個樣本:111111.11 J1小丁沖擊響應(yīng)21.51 0.5 幅 振0-0.5-1-1.5051015202530354045時間序號Q2.21利用filter命令編寫一個 MATLAB程序，生成式（2.17）給出的因果線性時不變系統(tǒng)的 沖激響應(yīng)，計算并畫出前40個樣本。把你的結(jié)果和Q2.20中得到的結(jié)果相比較。沖擊響應(yīng)前40個樣本。Q2.22編寫一個MATLA弧序，生成并畫出式（2.11）給出的因果線性時不變系統(tǒng)的階躍響應(yīng)。用該程序計算并畫出式（2.15）給出

27、的線性時不變系統(tǒng)的階躍響應(yīng)的前40個樣本。階躍響應(yīng)前40個樣本。沖擊響應(yīng)Q2.23運行程序P2.6,計算輸出序列yn和y2n以及差值信號dn , yn和y2n相等嗎?四階實現(xiàn)的輸出0-14010差值信號x 10幅 振幅 振-1 00.5幅 0振-0.51OC1L隼r11士1,4,1 1J1015202530351,1中11 . K1.1|,1.一由11.1 11 J級聯(lián)實現(xiàn)的輸出10152025303540答：可以看出這加權(quán)輸入與輸出的圖像幾乎是相差無幾的，觀察差信號的圖像，縱坐標(biāo)已達(dá)到 10A(-14),可以看出是相差無幾的，所以兩信號相等。時間序號nQ2.25用任意的非零初始向量ic, i

28、c1和ic2來重做Q2.23ic=1 1y1=filter(num1,den1,x,ic);y2=filter(num2,den2,y1,ic);-1 0rO 0<JlJ*0 .z3 - 口口.j' r.%#2ri-1.卜,66%J¥四階實現(xiàn)的輸出0t5101520253035-2 040>? o?o g,一,¥U2 口 T 11-.J上匚.丁飛二產(chǎn)> ii 也、y01JOJcrju5(J級聯(lián)實現(xiàn)的輸出0510152025303540-2 0Q 口卜-a中>g oq T T)9i-_。 q.J- 1cL .XL*'燈口 F-,差值信號

29、0510303540152025時間序號nic1=2 6四階實現(xiàn)的輸出10幅 0振-10051015202530354010幅 0振-10J-:JI-11 q.1 11J.0510152025303540時間序號n差值信號ic2=7 4-1 0I1j%9,1 1' 1 kA .16p1 1.出1J、四階實現(xiàn)的輸出051015202530354010-101,j)廣1rii h 1J lc. .1-f*KJV級聯(lián)實現(xiàn)的輸出010152025303540幅 振差值信號時間序號n不相等兩個輸出之間Q2.26修改程序P2.6,將兩個二階系統(tǒng)順序顛倒并在零初始條件下重復(fù)此過程。有差別嗎?四階實現(xiàn)

30、的輸出0 -510152025303540-1 00.5幅 0振0-0.5C1<1 n.自凸4Q八口丫 X 1了 (46 J中1-3 dL5103035152025時間序號n40Q2.28運行程序P2.7,對序列hn和xn求卷積，生成yn,并用FIR濾波器hn 對輸入xn濾波，求得y1n0 yn和y1n有差別嗎？為什么要使用對xn補零后得 到的x1n作為輸入來產(chǎn)生y1n?20-20 09用卷積得到的輸出100-102468101214時間序號n由濾波生成的輸出20-20 0r一100-102468101214時間序號n答：yn與y1n沒有差別。序列hn和xn卷積得到的yn的長度hn和xn

31、的長度之和再減去 1,即9+7-1=15,而用FIR濾波器hn對輸入信號xn濾波得到的y1n長度與xn相同，故要對 xn補零后得到的x1n(長度擴展為15),濾波后所得輸出與 y1n與y2n相同。Q2.30程序中為什么要使用命令for和end?For后常為循環(huán)語句需執(zhí)行的的循環(huán)次數(shù)，end表示循環(huán)語句執(zhí)行完畢在累加的過程中若|h(k)|的值小于10A(-6)S(K)已經(jīng)收斂，所以不必再檢查之后的值，所以用for ,end來進(jìn)行判斷Q2.3邪序P2.所定的沖激響應(yīng)所對應(yīng)的離散時間系統(tǒng)是什么？運行程序P2.8,生成該沖激響應(yīng)。該系統(tǒng)穩(wěn)定嗎？若hK不小于10A(-6)且圖形顯示了一個衰減的沖激 響應(yīng)

32、，用一個較大的N值再次運行程序P2.8。答：該系統(tǒng)穩(wěn)定。Q2.33 考慮用差分方程 yn=xn-4xn-1+3xn-2+1.7yn-1-yn-2描述的離散時間系統(tǒng)。修改程序P2.8,計算并畫出上述系統(tǒng)的沖激響應(yīng)。該系統(tǒng)穩(wěn)定嗎？沖擊響應(yīng)時間序號n4-4 -：:020406080100120140160180200時間序號n答：該系統(tǒng)不穩(wěn)定。第四次實驗第三章離散時間信號的頻域分析Q3.1在程序P3.1中，計算離散時間傅里葉變換的原始系列是什么？MATLAB命令pause的作用是什么？答；原始序歹U H(jw)=(2+eA(-jw)/(1-0.6eA(-jw) pause的作用是暫時停止直到用戶按任

33、何鍵。Q3.2運行程序P3.1,求離散時間傅里葉變換的實部，虛部以及幅度和相位譜。離散時間福利 葉變換是w的周期函數(shù)嗎？若是，周期是多少？描述這四個圖形表示的對稱性。答；離散時間傅里葉變換是3的周期函數(shù)，周期為 2兀。H(ej )的實部86幅4振I120-4-3-2-101234 omegaApiH(ej )的虛部8H(e )幅度譜6產(chǎn)4 振20i-4-3-2-101234 omega/pi 相位譜H(e j )2立目1勺立津0為姜-1瓜八-2-4-3-2-101234 omega/piQ3.3修改程序p3.i，在范圍0 w W W 口內(nèi)計算如下序列的離散時間傅里葉變換；U(ej )0.7 0

34、.5e j 0.3e j2 e j31 0.3e j 0.5e j2 0.7e j3并重做習(xí)題 Q3.2。討論你的結(jié)果。你能解釋相位譜中的跳變嗎？MATLAB命令unwrap可以移除跳變。試求移除后的相位譜。H(ej )的實部0.50-0.5-1H(ej )幅度譜 omegaApiH(ej )的虛部0.50-0.50.30.40.50.60.70.80.9 omegaApi-1 00.10.100.10.20.30.40.50.60.70.80.91 omega/pi相位譜H(ej )2 O-2 相的位單為度-4 00.10.20.30.40.50.60.70.80.91 omega/piH(

35、ej )幅度譜 omega/pi答；因為離散時間傅里葉變換是3的周期函數(shù)，周期為 2兀，當(dāng)計算的相位在頻率范圍-71,兀 之外時，就會自動加2兀發(fā)生相位譜跳變。P3.2,對程序生成的圖形中的兩個軸加Q3.6通過加入合適的注釋語句和程序語句，修改程序標(biāo)記。哪個參數(shù)控制時移量?答：D控制時移量Q3.7運行修改后的程序并討論你的結(jié)果。60原序列的幅度譜40幅 振201600-1-0.500.5w/pi原序列的相位譜4時移后的幅度譜40幅振2012 O-2-4 位相的位單為度弧以-1-0.50w/pi0.51位相的位單為度弧以0-1-0.500.5w/pi時移后的相位譜420-2-41-1-0.500

36、.5w/pi答：通過觀察知道當(dāng) D增大時，序列的幅度譜不變；相位譜變密，周期縮短。Q3.8選取不同的時移值重做 Q3.7D=560幅振D=50原 序列的幅度譜幅振402060時移后的幅度譜o O 4 2譜 位 相 的 列譜 位 相 的 后 移 時位相的位單為度弧以4 2 O-2-4位相的位單為度弧以4 260原序列的幅度譜1幅 振40幅振200 -1-0.500.5w/pi600-1時移后的幅度譜4020-0.500.51w/pi時移后的相位譜原序列的相位譜1-1-0.500.5w/pi位相的位單為度弧以-4-120-2-0.500.51w/piP3.2,對程序生成的圖形中的兩個軸Q3.10通

37、過加入合適的注釋語句和程序語句，修改程序加標(biāo)記。哪個參數(shù)控制頻移量?答:w0控制頻移量Q3.11運行修改后的程序并討論你的結(jié)果原序列的幅度譜4原序列的相位譜頻移后的相位譜2 O-2-4 譜位相位單為度弧以-1-0.500.51w/piw/pi譜位相位單為度弧以答：在參數(shù)w0的控制下，離散時間傅里葉變換的幅度譜和相位譜都隨著控制參數(shù)左移。Q3.12選取不同的頻移值，重做 Q3.11W0=0.2pi原序列的幅度譜150100幅振50頻移后的幅度譜100啊 50振譜位相位單為度弧以譜位相位單為度弧以4 2 04 2 0 2 4 - -W0=pi幅振P3.4,對程序生成的圖形中的兩個軸Q3.14通過加

38、入合適的注釋語句和程序語句，修改程序 加標(biāo)記。幅度譜的乘積200150啊100 振50卷積后序列的幅度譜200150多100 振50位相的位單為度弧以00.5/列的相位譜Q3.15運行修改后的程序并討論你的結(jié)果。幅度譜的乘積400010003000x 1 x10O 2卷積后序列的幅度譜40000-1-0.500.5w相位譜的和430002000100010-1-0.500.5w卷積后序列的相位譜4-1-0.500.5w12 O-2 度痹H1H-42 O-2度幅位相后積卷-4 :-1-0.500.51w時域信號的卷積等于時域信號在頻域上的乘積，時域信號的乘積等于時域信號在頻域上的第五次實驗Q3.

39、23編寫一個MATLAB程序，計算并畫出長度為 N的L點離散傅里葉變換 Xk的值，其中L N,然后計算并畫出L點離散傅里葉逆變換Xk。對不同長度N和不同的離散傅里葉變換長度L,運行程序。討論你的結(jié)果。序歹 U; xn = 0 1 2 3 4 5 6 7 8 9 0 0 0 0 00；XK模xn度 幅時間序列nXK相位1005-51010時間序列nxn模度 幅度 幅1015時間序列nxn相位度 幅度 幅00時間序列n時間序列n時間序列n序列:xn=1/3 1/3 1/3 xnXK模時間序列n時間序列nxn模時間序列n時間序列nQ3.26在函數(shù)circshift中，命令rem的作用是什么?答：re

40、m(x,y)是用y對x求余數(shù)函數(shù)。Q3.27解釋函數(shù)circshift怎樣實現(xiàn)圓周移位運算。答：在輸入序列x由M的位置開始被循環(huán)移位。如果 M> 0 ,則circshift刪除從矢量x最左邊開 始的M個元素和它們附加在右側(cè)的剩余元素，以獲得循環(huán)移位序列。 如果如果M<0,則circshift首先通過x的長度來彌補 M ,即序列x最右邊的長度的 M樣品從x中刪除和所附在其余的 M個 樣本的右側(cè)，以獲得循環(huán)移位序列。Q3.28在函數(shù)circshift中，運算符=的作用是什么?答：二是不等于的意思。Q3.29解釋函數(shù)circonv怎樣實現(xiàn)圓周卷積運算。答：輸入是兩個長度都為 L的向量x1

41、和x2,它是非常有用的定期延長X2的函數(shù)。讓x2p成為x2延長無限長的周期的序列。從概念上講，在定點時間上通過時序交換后的x2p的長度L交換x2p序列和x2tr等于1的元素。然后元素1至L的輸出向量y是通過取x1和獲得的長度為 L的sh矢量之間的內(nèi)積得到通過循環(huán)右移的時間反轉(zhuǎn)向量x2tr。對于輸出樣本 Yn的1WNWL時，右循環(huán)移位的量為n-1個位置上。Q3.31運行修改后的程序并驗證圓周時移運算Q3.32通過加入合適的注釋語句和程序語句，修改程序P3.8,對程序生成的圖形中的兩個軸加標(biāo)記。時移量是多少？nn答：時移量為5Q3.33運行修改后的程序并驗證離散傅里葉變換的圓周時移性質(zhì)。原序列的離散傅里葉變換幅畫周移位后