利用MATLAB編寫卷積函數(shù)myconv

《利用MATLAB編寫卷積函數(shù)myconv》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、一、實(shí)驗(yàn)?zāi)康?.了解MATLAB的基本操作。2.利用MATLAB實(shí)現(xiàn)正弦信號(hào)采樣；3.利用MATLAB編寫卷積函數(shù)myconV。二、實(shí)驗(yàn)條件PC機(jī)，MATLAB7.0三、實(shí)驗(yàn)內(nèi)容(-)函數(shù)文件與腳本文件的定義，正弦信號(hào)的釆樣函數(shù)文件代碼：functionseqs=mysampling(FszL,Fx)t=l/Fs:l/Fs:L/Fs;t1=0:0.001:L/Fs;x=sin(2*pi*Fx*t);xl=sin(2*pi*Fx*tl);plot(t1,xl),holdon;stem(t,x,1r1);holdoff;

2、腳本文件代碼：Fs=40;Fx=20;L=15;forFx=20:10:80,mysampling(Fs,L,Fx);pause,End運(yùn)行結(jié)果：1080604020-02-0406-08-100.050.1015020.250.30.350.4實(shí)驗(yàn)屮遇到的問(wèn)題及改正：functionseqs=mysampling(Fs,L,Fx)t=1/Fs:1/Fs:L/Fs;x=sin(2*pi*Fx*t);plot(tzx)zholdon;stem(tzx,1r1);holdoff;一開始并沒有添加tl=0:0.001:L/

3、Fs;xl=sin(2*pi*Fx*tl);兩行代碼，想利用t=l/Fs:l/Fs:L/Fs;x=sin(2*pi*Fx*t);生成的函數(shù)直接作為需要抽樣的原函數(shù)，但是t的采樣間隔太大，沒法生成平滑的正弦曲線，所以導(dǎo)致最后仿真岀來(lái)的圖像中岀現(xiàn)不光滑的正弦曲線。經(jīng)過(guò)思考找出問(wèn)題后，另外編寫了一個(gè)正弦函數(shù)，設(shè)定的采樣間隔為0.001后，仿真出來(lái)的圖像中就出現(xiàn)了光滑的正弦曲線。(二)編寫卷積函數(shù)myconvConv函數(shù)代碼：functiony=myconvl(hfx)n=length(x);k=length(h);%定義序

4、列長(zhǎng)度f(wàn)orm=l:1:(n-1)A(m)=0;endform=n:1:n+k-1A(m)=h(m-n+l);%給h序列添(n-l)個(gè)0endB=hankel(A);fori=l:1:nC(:,i)=B(:,i);%抽取B矩陣屮從笫1列到第n列的矩陣endD=fliplr(C);%C矩陣左右翻轉(zhuǎn)E=D';%轉(zhuǎn)置y(n+k-1)=0;fora=l:1:n+k-1forb=l:1:nnumble(b)=x(b)*E(b,a);y(a)=numble(b)+y(a);endnumble=[000];貂青零endystem(

5、y);腳木代碼：x=[l-231];h=[-32-14];%定義x,h序列myconvl(hAx);運(yùn)行結(jié)果：?UntitledGy=■38-U9-911arts=-38-149-911A?I四、實(shí)驗(yàn)結(jié)論和討論通過(guò)本次實(shí)驗(yàn)運(yùn)用MATLAB實(shí)現(xiàn)采樣以及自行編寫卷積函數(shù)，基本能熟悉運(yùn)用MATLAB來(lái)進(jìn)行實(shí)驗(yàn)，并且對(duì)于采樣以及卷積的認(rèn)識(shí)更加的深入。在進(jìn)行本次卷積實(shí)驗(yàn)之前，我一直只知道利用MATLAB本身自有的conv函數(shù)來(lái)進(jìn)行卷積的運(yùn)算，并沒有想過(guò)conv函數(shù)的卷積運(yùn)算功能是如何實(shí)現(xiàn)的。通過(guò)這次自己編寫卷積函數(shù)的代碼，不僅

6、熟知了hankel,fliplr以及轉(zhuǎn)置等MATLAB中及其有用的用法，而且對(duì)于卷積函數(shù)有了更深一步的了解。并J4通過(guò)自己所寫的myconvl函數(shù)與系統(tǒng)自帶函數(shù)conv的比較，發(fā)現(xiàn)自己所編寫的函數(shù)的運(yùn)行時(shí)間要多于corw函數(shù)運(yùn)行的時(shí)間，意識(shí)到MATLAB的深?yuàn)W之處，也發(fā)現(xiàn)自己的能力還遠(yuǎn)遠(yuǎn)不足，需要在進(jìn)一步的努力。