終稿 信號(hào)與系統(tǒng) MATLAB 實(shí)驗(yàn)報(bào)告 YAO

1、信號(hào)與系統(tǒng) MATLAB 仿真實(shí)驗(yàn)報(bào)告姓名：姚俊濤 學(xué)號(hào)：2021302540313班級(jí)：0910學(xué)院：電氣工程學(xué)院 專業(yè)：電氣工程與自動(dòng)化目錄實(shí)驗(yàn)1： 連續(xù)時(shí)間信號(hào)的表示及可視化4、實(shí)驗(yàn)內(nèi)容4、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形4(1)沖激信號(hào)4(2)階躍信號(hào)5(3)指數(shù)信號(hào)5(4)門函數(shù)6(5)f(t)=Sa(nwt)7(6)f(t)=Sin(2*pi*f*t)8實(shí)驗(yàn)體會(huì)：9實(shí)驗(yàn)2： 離散時(shí)間信號(hào)的表示及可視化10、實(shí)驗(yàn)內(nèi)容10、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形10(1)離散抽樣序列10(2) 階躍序列11(3)指數(shù)序列11(4) f(n)=R(n) (分別取不同的N值)12(5) f(n)=Sa(nw)14(6) f

2、(n)=Sin(nw) (分別取不同的w值)15實(shí)驗(yàn)心得：17實(shí)驗(yàn)3： 系統(tǒng)的時(shí)域求解17、實(shí)驗(yàn)內(nèi)容17、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形173.2.1 畫(huà)出x(n)、h(n)、y(n)的波形173.2.2 繪出y(n)的單位抽樣響應(yīng)，并繪出H的幅頻及相頻特性曲線18實(shí)驗(yàn)心得20實(shí)驗(yàn)4： 信號(hào)的DFT分析20、實(shí)驗(yàn)內(nèi)容20、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形204.2.1 N=10時(shí)214.2.2 N=16時(shí)224.2.3 N=22時(shí)234.3 擴(kuò)展實(shí)驗(yàn) N=32()和N=64()時(shí)244.3.1 N=32時(shí)244.3.2 N=64時(shí)254.4 結(jié)果比照與差異原因分析26進(jìn)行理論值與計(jì)算值的比擬26查閱相關(guān)資料討論信號(hào)頻譜

3、分析過(guò)程中誤差原因及改善方法26實(shí)驗(yàn)心得27實(shí)驗(yàn)5： 系統(tǒng)時(shí)域解的快速卷積求法27、實(shí)驗(yàn)內(nèi)容27、問(wèn)題分析27、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形285.3.1 序列、的圖形285.3.2 N=25時(shí)295.3.3 N=34時(shí)305.3.4 N=45時(shí)31正確卷積結(jié)果325.5 結(jié)果比照與差異原因分析33、實(shí)驗(yàn)心得體會(huì)34六、信號(hào)與系統(tǒng)MATLAB仿真實(shí)驗(yàn)心得體會(huì)34實(shí)驗(yàn)1： 連續(xù)時(shí)間信號(hào)的表示及可視化1.1、實(shí)驗(yàn)內(nèi)容；分別??；分別畫(huà)出不同周期個(gè)數(shù)的波形；1.2、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形(1)沖激信號(hào)syms tf=dirac(t)ezplot(f,-5:5);xlabel('(t)');ylabe

4、l('(f)');沖激函數(shù)示意圖如下(2)階躍信號(hào)創(chuàng)立Heaviside 的M 文件，該文件如下：function f=Heaviside(t)正確定義出該函數(shù)并保存運(yùn)行后，就可調(diào)用該函數(shù)了。如先定義向量：t=-5:0.01:5然后調(diào)用 Heaviside 函數(shù)表示出該信號(hào)并繪出波形f=heaviside(t)plot(t,f)axis(-5,5,-0.2,1.2)(3)指數(shù)信號(hào)將a=1和a=2時(shí)的圖形分別畫(huà)于同一界面中，便于形成比照，syms tsubplot 121f=exp(t) % a=1ezplot(f,-5:5);xlabel('(t)');ylab

5、el('(f)');title('f(t)=exp(t)')subplot 122f=exp(-t) %a=-1ezplot(f,-5:5);xlabel('(t)');ylabel('(f)');title('f(t)=exp(t)')(4)門函數(shù)t=-5:0.01:5;%設(shè)定時(shí)間變量t的范圍及步長(zhǎng)y=rectpuls(t,2); %用rectpulst a)命令表示門函數(shù)，默認(rèn)以零點(diǎn)為中心，寬度為aplot(t,y);%用plot函數(shù)繪制連續(xù)函數(shù)title('門函數(shù)'); %用title函數(shù)設(shè)置

6、圖形的名稱axis(-5 5 -0.5 1.5);(5)f(t)=Sa(nwt)syms tf=Sinc(pi/8*t)subplot 121ezplot(f,-15:15);xlabel('(t)');ylabel('(f)');axis(-15 15 -0.3 1.2)title('f(t)=Sa(pi/8*t)'); f=sinc(pi/18*t)subplot 122ezplot(f,-15:15);xlabel('(t)');ylabel('(f)');title('f(t)=Sa(pi/18*t

7、)');(6)f(t)=Sin(2*pi*f*t)周期信號(hào)(分別畫(huà)出不同周期個(gè)數(shù)的波形)sym ta %頻率記為af=sin(2*pi*a*t)ezplot(f,-15:15);xlabel('(t)');ylabel('(f)');axis(-15 15 -1.2 1.2)1.3實(shí)驗(yàn)體會(huì)： 通過(guò)這次實(shí)驗(yàn)，我掌握了根本信號(hào)函數(shù)圖形的繪制方法；更加直觀地認(rèn)識(shí)了個(gè)根本信號(hào)的特性曲線。實(shí)驗(yàn)2： 離散時(shí)間信號(hào)的表示及可視化2.1、實(shí)驗(yàn)內(nèi)容；分別??；分別取不同的N值；分別取不同的值；2.2、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形(1)離散抽樣序列由單位序列的定義知，只有在k = 0時(shí)

8、，f (k )的值為1，而當(dāng) 時(shí)的值為零，用MATLAB的dirac語(yǔ)句來(lái)來(lái)表示單位序列及繪制其波形。y=0,0,0,0,0,1,0,0,0,0,0'x=-5:5;stem(x,y)axis(-5,5,-0.5,1.5)title('單位脈沖序列')xlabel('n') ylabel('y(n)')(2) 階躍序列 y=0,0,0,0,0,1,1,1,1,1,1'x=-5:5;stem(x,y)axis(-5,5,-0.5,1.5)title('階躍函數(shù)')xlabel('n')ylabel(&#

9、39;y(n)')(3)指數(shù)序列 (分別取a>0及a<0)a=1; % a=1時(shí)t=-5:1:5;f=exp(a*t);subplot 121stem(t,f,'b.','MarkerSize',10)xlabel('t')ylabel('f(t)')title('f(t)=exp(a*t) a=1') a=-1; %a=-1時(shí)t=-5:1:5;f=exp(a*t);subplot 122stem(t,f,'b.','MarkerSize',10)xlabel(&

10、#39;t')ylabel('f(t)')title('f(t)=exp(a*t) a=-1')(4) f(n)=R(n) (分別取不同的N值)N=10時(shí) n=0:1:9;f=ones(10);stem(n,f,'b')axis(-5,15,0,1.5)title('R(n)n=10')N=20時(shí) n=0:1:19;f=ones(20);stem(n,f,'b')axis(-5,25,0,1.5)title('R(n)n=20')(5) f(n)=Sa(nw)n=-45:1:45;f=sin

11、c(0.1*n);stem(n,f,'.');xlabel('n');ylabel('f');axis(-50 50 -1 1)n=-45:1:45;f=sinc(0.2*n);stem(n,f,'.');xlabel('n');ylabel('f');axis(-50 50 -1 1)(6) f(n)=Sin(nw) (分別取不同的w值)w=100時(shí)n=-15:1:15;f=sin(100*n);stem(n,f,'.');xlabel('n');ylabel(&#

12、39;f');w=200時(shí)n=-15:1:15;f=sin(200*n);stem(n,f,'.');xlabel('n');ylabel('f');2.3實(shí)驗(yàn)心得：經(jīng)查閱資料得知在用MATLAB表示離散序列并將其可視化時(shí)，我們要注意：第一，與連續(xù)時(shí)間信號(hào)不同，離散時(shí)間信號(hào)無(wú)法用符號(hào)運(yùn)算來(lái)表示；第二，由于在MATLAB中，矩陣的元素個(gè)數(shù)是有限的，因此，MATLAB無(wú)法表示無(wú)限序列；第三，在繪制離散信號(hào)波形時(shí)，要使用專門繪制離散數(shù)據(jù)的stem命令，而不是plot命令。 實(shí)驗(yàn)3： 系統(tǒng)的時(shí)域求解3.1、實(shí)驗(yàn)內(nèi)容1、設(shè)，求：，并畫(huà)出、波形。2、

13、求因果線性移不變系統(tǒng)的單位抽樣響應(yīng)，并繪出的幅頻及相頻特性曲線。3.2、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形 畫(huà)出x(n)、h(n)、y(n)的波形n=-10:30; %設(shè)定變量范圍，步長(zhǎng)默認(rèn)為1h=0.9.n.*f; %hn的表示，注意向量相乘注意用點(diǎn)subplot 221 %產(chǎn)生圖形窗口2stem(n,h); %繪制hnxlabel('h(n)');axis(-10 10 0 1);x=heaviside(n)-heaviside(n-10); %階躍函數(shù)直接相減subplot 222 %產(chǎn)生圖形窗口3stem(n,x);xlabel('x(n)')axis(-10 20 0

14、 1.5); %設(shè)定坐標(biāo)軸的范圍subplot 212 y=conv(h,x); %求h和x的卷積用命令convstem(y); %繪制函數(shù)yxlabel('y(n)=x(n)*h(n)')axis(-10 20 0 10);如下列圖所示，依次為f=h(n)、x(n)、y(n)=x(n)*h(n)的圖像：經(jīng)分析，得到的卷積和結(jié)果符合理論計(jì)算結(jié)果，理論上求解卷積和的步驟為：先將x(n)反褶，然后向右平移，相乘，求和。當(dāng)門函數(shù)序列恰好全部進(jìn)入等比序列時(shí)，卷積和最大，繼續(xù)向右平移那么卷積和隨著等比數(shù)列各數(shù)相對(duì)減小而下降。3.2.2 繪出y(n)的單位抽樣響應(yīng)，并繪出H的幅頻及相頻特性

15、曲線因果線性移不變系統(tǒng)，求其單位抽樣響應(yīng)，并繪出的幅頻及相頻特性曲線。a=1 0 -0.81; %描述系統(tǒng)的差分方程的系數(shù)b=1 0 -1;%描述系統(tǒng)的差分方程的系數(shù)figure(1);h=impz(b,a,-10:10); %利用impz函數(shù)求系統(tǒng)的單位樣值響應(yīng)stem(h); %繪制函數(shù)h 的離散序列xlabel('h(n)');title('系統(tǒng)的單位樣值響應(yīng)')figure(2);freqs(b,a,:);%對(duì)連續(xù)系統(tǒng)頻率響應(yīng)H(jw)進(jìn)行分析的函數(shù)freqs()title('系統(tǒng)幅頻及相頻特性曲線')實(shí)驗(yàn)圖形如下列圖所示：實(shí)驗(yàn)結(jié)果與理論

16、解相同，相比之下使用軟件可快速求解，方便、準(zhǔn)確，更加實(shí)用。實(shí)驗(yàn)結(jié)果顯示該系統(tǒng)為低通型，相位移一直為.3.3實(shí)驗(yàn)心得通過(guò)本次實(shí)驗(yàn)，我掌握了利用MATLAB進(jìn)行時(shí)域求解及其可視化的方法，并直觀地認(rèn)識(shí)了幅度特性、相位特性的物理意義。實(shí)驗(yàn)4： 信號(hào)的DFT分析4.1、實(shí)驗(yàn)內(nèi)容計(jì)算余弦序列的DFT。分別對(duì)N=10、16、22時(shí)計(jì)算DFT，繪出幅頻特性曲線，分析是否有差異及產(chǎn)生差異的原因。4.2、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形問(wèn)題分析：連續(xù)時(shí)間傅里葉變換主要用來(lái)描述連續(xù)時(shí)間非周期信號(hào)的頻譜。任意的非周期信號(hào)，如果滿足狄里克利條件，那么，它可以被看作是由無(wú)窮多個(gè)不同頻率的周期復(fù)指數(shù)信號(hào)的線性組合構(gòu)成的，每個(gè)頻率所對(duì)應(yīng)的

17、周期復(fù)指數(shù)信號(hào)稱為頻率分量，其相對(duì)幅度為對(duì)應(yīng)頻率的的之值，其相位為對(duì)應(yīng)頻率的的相位。 N=10時(shí)實(shí)驗(yàn)程序如下：N=10;n=0:N-1;x=cos(pi/8*n);y=fft(x,N);subplot(2,1,1);stem(n,x);title('取N=10時(shí) 離散序列'); ylabel('y(n)');xlabel('n');subplot(2,1,2);stem(n,y);ylabel('Phase');xlabel('Frequency(rad)');title('取N=10時(shí) 幅頻特性曲線

18、9;);實(shí)驗(yàn)圖形如下列圖所示： N=16時(shí)實(shí)驗(yàn)程序如下：N=16;n=0:N-1;x=cos(pi/8*n);y=fft(x,N);subplot(2,1,1);stem(n,x);title('取N=16時(shí) 離散序列'); ylabel('y(n)');xlabel('n');subplot(2,1,2);stem(n,y);ylabel('Phase');xlabel('Frequency(rad)');title('取N=10時(shí) 幅頻特性曲線');實(shí)驗(yàn)圖形如下列圖所示：4.2.3 N=22時(shí)實(shí)

19、驗(yàn)程序如下：N=22;n=0:N-1;x=cos(pi/8*n);y=fft(x,N);subplot(2,1,1);stem(n,x);title('取N=22時(shí) 離散序列'); ylabel('y(n)');xlabel('n');subplot(2,1,2);stem(n,y);ylabel('Phase');xlabel('Frequency(rad)');title('取N=10時(shí) 幅頻特性曲線');實(shí)驗(yàn)圖形如下列圖所示：與補(bǔ)充的N=32、N=64時(shí)的序列圖和幅頻特性圖作以比照，然后分析選

20、取不同的N時(shí)DFT幅頻特性曲線產(chǎn)生差異的原因。4.3 擴(kuò)展實(shí)驗(yàn) N=32()和N=64()時(shí)4.3.1 N=32時(shí)實(shí)驗(yàn)程序如下：N=32;n=0:N-1;x=cos(pi/8*n);y=fft(x,N);subplot(2,1,1);stem(n,x);title('取N=32時(shí) 離散序列'); ylabel('y(n)');xlabel('n');subplot(2,1,2);stem(n,y);ylabel('Phase');xlabel('Frequency(rad)');title('取N=32時(shí)

21、幅頻特性曲線');實(shí)驗(yàn)圖形如下列圖所示：4.3.2 N=64時(shí)實(shí)驗(yàn)程序如下：N=64;n=0:N-1;x=cos(pi/8*n);y=fft(x,N);subplot(2,1,1);stem(n,x);title('取N=64時(shí) 離散序列'); ylabel('y(n)');xlabel('n');subplot(2,1,2);stem(n,y);ylabel('Phase');xlabel('Frequency(rad)');title('取N=64時(shí) 幅頻特性曲線');實(shí)驗(yàn)圖形如下列圖所

22、示：4.4 結(jié)果比照與差異原因分析進(jìn)行理論值與計(jì)算值的比擬此題為周期信號(hào)，無(wú)直流分量，所以取樣點(diǎn)數(shù)可為。從取樣點(diǎn)數(shù)N=16和N=32可以看出，取樣點(diǎn)數(shù)的不同，會(huì)造成頻率譜和相位譜的不同。當(dāng)N=16時(shí)，n=1或15時(shí)有幅度值，而在N=32時(shí)，n=2和30時(shí)有幅度值，在N=64時(shí)，n=4和60時(shí)有幅度值，得到在N=16時(shí)，其頻譜已經(jīng)和N=32時(shí)一致剛好成2倍關(guān)系，且N=16、32、64時(shí)均不產(chǎn)生混頻現(xiàn)象。而N=10或22時(shí)由于其不是周期序列的周期值，故產(chǎn)生了明顯的頻域混疊。討論信號(hào)頻譜分析過(guò)程中誤差原因及改善方法在頻譜分析過(guò)程中由于取樣頻率過(guò)低或者由于信號(hào)的截取長(zhǎng)度不當(dāng)將會(huì)產(chǎn)生誤差。取樣頻率過(guò)低,

23、可能會(huì)產(chǎn)生混頻現(xiàn)象,可以適當(dāng)提高取樣率，增加樣點(diǎn)數(shù)，來(lái)減少混疊對(duì)頻譜分析所造成的誤差。對(duì)于連續(xù)周期信號(hào)，其時(shí)域取樣必須滿足時(shí)域取樣定理：其取樣點(diǎn)數(shù)K2*N+1其中N為最高諧波分量。截取信號(hào)長(zhǎng)度不當(dāng)，會(huì)產(chǎn)生功率泄露，對(duì)周期序列進(jìn)行頻譜分析時(shí)，為防止泄露應(yīng)做到：截取的長(zhǎng)度應(yīng)取一個(gè)根本周期或根本周期的整數(shù)倍，假設(shè)待分析的周期信號(hào)事先不知道其確切的周期，那么可截取較長(zhǎng)時(shí)間長(zhǎng)度的樣點(diǎn)進(jìn)行分析，以減少功率泄露誤差。當(dāng)然，必須在取樣頻率滿足取樣定理的條件下進(jìn)行，否那么混疊與泄露同時(shí)存在給頻譜分析造成困難。4.3實(shí)驗(yàn)心得通過(guò)這次實(shí)驗(yàn)，我掌握了連續(xù)時(shí)間周期信號(hào)的傅里葉級(jí)數(shù)的物理意義和分析方法；掌握了連續(xù)時(shí)間傅里

24、葉變換的分析方法及其物理意義；進(jìn)一步學(xué)習(xí)了系統(tǒng)頻率響應(yīng)特性的計(jì)算方法和特性曲線的繪制發(fā)法；深入學(xué)習(xí)了幅度特性、相位特性的物理意義。實(shí)驗(yàn)5： 系統(tǒng)時(shí)域解的快速卷積求法5.1、實(shí)驗(yàn)內(nèi)容用快速卷積法計(jì)算系統(tǒng)響應(yīng)，：，。要求取不同的L點(diǎn)數(shù)，并畫(huà)出、波形，分析是否有差異及產(chǎn)生差異的原因。5.2、問(wèn)題分析MATLAB的conv()函數(shù)可以幫助我們快速求出兩個(gè)離散序列的卷積和。conv函數(shù)的調(diào)用格式為: f=conv (fl, f2)，其中fl為包含序列f1(k)的非零樣值點(diǎn)的行向量，f2為包含序列九(k)的非零樣值點(diǎn)的行向量，向量f那么返回序列的所有非零樣值點(diǎn)行向量。此題可利用離散傅里葉變換的相關(guān)性質(zhì)進(jìn)行

25、卷積快速運(yùn)算，用fft和ifft函數(shù)語(yǔ)句實(shí)現(xiàn)。應(yīng)用FFT實(shí)現(xiàn)數(shù)字濾波器實(shí)際上就是用FFT來(lái)快速計(jì)算有限長(zhǎng)度序列的線性卷積。X(k)x(n)FFTIFFTH(k)H(k)X(k)y(n)這種方法就是先將輸入信號(hào)x(n)通過(guò)FFT變換為它的頻譜采樣值X(k)，然后再和FIR濾波器的頻響采樣值H(k)相乘，H(k)可事先存放在存儲(chǔ)器中，最后再將乘積H(k)X(k)通過(guò)快速傅里葉變換簡(jiǎn)稱IFFT復(fù)原為時(shí)域序列，即得到輸出y(n)。5.3、實(shí)驗(yàn)程序與實(shí)驗(yàn)圖形 5.3.1 序列、的圖形 程序如下：n1=0:14;%設(shè)定n1的值一代替門函數(shù)的作用n2=0:19;%設(shè)定n2的值以代替門函數(shù)的作用x=sin(0

26、.4.*n1);%正弦函數(shù)序列y=0.9.n2; %函數(shù)y的表達(dá)式Xk=fft(x,45);%求函數(shù)x的快速傅立葉變換Yk=fft(y,45);%求函數(shù)y的快速傅立葉變換Hk=Xk.*Yk;%時(shí)域內(nèi)的卷積對(duì)應(yīng)于頻域內(nèi)的點(diǎn)乘h=ifft(Hk);%求Hk的快速傅立葉變換的反變換，即為x和y的卷積值subplot(2,1,1),stem(x);%繪制x的圖形xlabel('n');ylabel('x(n)');title('原序列 x(n)');subplot(2,1,2),stem(y);%繪制y的圖形xlabel('n');yla

27、bel('h(n)');title('原序列 h(n)');實(shí)驗(yàn)圖形如下列圖所示：5.3.2 N=25時(shí)實(shí)驗(yàn)源程序如下：clear alln1=0:14;%設(shè)定n1的值一代替門函數(shù)的作用n2=0:19;%設(shè)定n2的值以代替門函數(shù)的作用x=sin(0.4.*n1);%正弦函數(shù)序列y=0.9.n2; %函數(shù)y的表達(dá)式Xk=fft(x,25);%求函數(shù)x的快速傅立葉變換Yk=fft(y,25);%求函數(shù)y的快速傅立葉變換Hk=Xk.*Yk;%時(shí)域內(nèi)的卷積對(duì)應(yīng)于頻域內(nèi)的點(diǎn)乘h=ifft(Hk);%求Hk的快速傅立葉變換的反變換，即為x和y的卷積值stem(h);%繪制h

28、的圖形xlabel('n');ylabel('y(n)')title('取L=25時(shí) 快速卷積法處理結(jié)果');實(shí)驗(yàn)圖形如下列圖所示：5.3.3 N=34時(shí)實(shí)驗(yàn)源程序如下：clear alln1=0:14;%設(shè)定n1的值一代替門函數(shù)的作用n2=0:19;%設(shè)定n2的值以代替門函數(shù)的作用x=sin(0.4.*n1);%正弦函數(shù)序列y=0.9.n2; %函數(shù)y的表達(dá)式Xk=fft(x,34);%求函數(shù)x的快速傅立葉變換Yk=fft(y,34);%求函數(shù)y的快速傅立葉變換Hk=Xk.*Yk;%時(shí)域內(nèi)的卷積對(duì)應(yīng)于頻域內(nèi)的點(diǎn)乘h=ifft(Hk);%求Hk的

29、快速傅立葉變換的反變換，即為x和y的卷積值stem(h);%繪制h的圖形xlabel('n');ylabel('y(n)')title('取L=34時(shí) 快速卷積法處理結(jié)果');實(shí)驗(yàn)圖形如下列圖所示：5.3.4 N=45時(shí)實(shí)驗(yàn)源程序如下：clear alln1=0:14;%設(shè)定n1的值一代替門函數(shù)的作用n2=0:19;%設(shè)定n2的值以代替門函數(shù)的作用x=sin(0.4.*n1);%正弦函數(shù)序列y=0.9.n2; %函數(shù)y的表達(dá)式Xk=fft(x,45);%求函數(shù)x的快速傅立葉變換Yk=fft(y,45);%求函數(shù)y的快速傅立葉變換Hk=Xk.*Yk

30、;%時(shí)域內(nèi)的卷積對(duì)應(yīng)于頻域內(nèi)的點(diǎn)乘h=ifft(Hk);%求Hk的快速傅立葉變換的反變換，即為x和y的卷積值stem(h);%繪制h的圖形xlabel('n');ylabel('y(n)')title('取L=45時(shí) 快速卷積法處理結(jié)果');實(shí)驗(yàn)圖形如下列圖所示：用conv()函數(shù)求出兩個(gè)離散序列的卷積和。conv函數(shù)的調(diào)用格式為: f=conv (fl, f2).實(shí)驗(yàn)程序如下：clear alln1=0:14;%設(shè)定n1的值一代替門函數(shù)的作用n2=0:19;%設(shè)定n2的值以代替門函數(shù)的作用x=sin(0.4.*n1);%正弦函數(shù)序列y=0.9.n2; %函數(shù)y的表達(dá)式z=conv(x,y); %求h和x的卷積用命令convstem(z); %繪制函數(shù)yxlabel('n');ylabel('y(n)=x(n)*h(n)');title('用conv函數(shù)計(jì)算的卷積結(jié)果')實(shí)驗(yàn)圖形如下列圖所示：5.5 結(jié)果比照與差異原因分析經(jīng)與使用conv()函數(shù)求出兩個(gè)離散序列的卷積和作比照，可以看出，當(dāng)時(shí)，快速卷積法求得卷積和是正確的。