杭電實(shí)驗(yàn)線性卷積與圓周卷積的計(jì)算

1、題目題目：已知兩個(gè)有限長(zhǎng)序列x(n)=δ(n)2δ(n1)3δ(n2)4δ(n3)5δ(n4)h(n)=δ(n)2δ(n1)δ(n2)2δ(n3)計(jì)算以下兩個(gè)序列的線性卷積和圓周卷積（1）x(n)⑤y(n)(2)x(n)⑥y(n)(3)x(n)⑨y(n)(4)x(n)⑩y(n)學(xué)弟學(xué)妹們??！不要盲目的復(fù)制進(jìn)去啊學(xué)弟學(xué)妹們??！不要盲目的復(fù)制進(jìn)去啊~~~~~~這里應(yīng)該先創(chuàng)建一個(gè)文檔保存函數(shù)這里應(yīng)該先創(chuàng)建一個(gè)文檔保存函數(shù)circonv再創(chuàng)建一

2、個(gè)文檔保存再創(chuàng)建一個(gè)文檔保存調(diào)用函數(shù)調(diào)用函數(shù)cirshiftd然后再把主程序拉進(jìn)去然后再把主程序拉進(jìn)去才能調(diào)用以上才能調(diào)用以上2個(gè)函數(shù)喔?jìng)€(gè)函數(shù)喔~~~還是好好學(xué)習(xí)吧！還是好好學(xué)習(xí)吧！●調(diào)用函數(shù)調(diào)用函數(shù)circonvfunctionyc=circonv(x1x2N)%用直接法實(shí)現(xiàn)圓周卷積%y=circonv(x1x2N)%y:輸出序列%x1x2:輸入序列%N:圓周卷積的長(zhǎng)度iflength(x1)Nerrendiflength(x2)Ne

3、rrend%以上語(yǔ)句判斷兩個(gè)序列的長(zhǎng)度是否小于Nx1=[x1zeros(1Nlength(x1))]%填充序列x1(n)使其長(zhǎng)度為N，序列h(n)的長(zhǎng)度為N1序列x(n)的長(zhǎng)度為N2x2=[x2zeros(1Nlength(x2))]%填充序列x2(n)使其長(zhǎng)度為Nn=[0:1:N1]x2=x2(mod(nN)1)%生成序列x2((n))N，鏡像，可實(shí)現(xiàn)對(duì)x(n)以N為周期的周期延拓，加1是因?yàn)镸ATLAB向量下標(biāo)只能從1開(kāi)始。H=ze

4、ros(NN)%生成N行N列的零矩陣fn=1:1:NH(n:)=cirshifted(x2n1N)%該矩陣的k行為x2((k1n))Nendyc=x1H%計(jì)算圓周卷積●調(diào)用函數(shù)調(diào)用函數(shù)cirshiftdfunctiony=cirshiftd(xmN)%直接實(shí)現(xiàn)序列x的圓周移位%y=cirshiftd(xmN)%x:輸入序列，且它的長(zhǎng)度小于N%m:移位位數(shù)?函數(shù)（函數(shù)（2）x(n)⑥y(n)x(n)⑥y(n)clearallN1=5N2=