線性卷積在DSP上實(shí)現(xiàn)

1、個(gè)人資料整理僅限學(xué)習(xí)使用線性卷積在DSP上的實(shí)現(xiàn)摘要 :卷積和 <簡(jiǎn)稱卷積）是信號(hào)處理中常用的算法之一。數(shù)字卷積運(yùn)算通常采用兩種方法：線性卷積和圓卷積。DSP<DigitalSignal Processor）是一種獨(dú)特的微處理器，是以數(shù)字信號(hào)來處理大量信息的器件。其工作原理是接收模擬信號(hào)，轉(zhuǎn)換為0 或 1 的數(shù)字信號(hào)。再對(duì)數(shù)字信號(hào)進(jìn)行修改、刪除、強(qiáng)化，并在其他系統(tǒng)芯片中把數(shù)字?jǐn)?shù)據(jù)解譯回模擬數(shù)據(jù)或?qū)嶋H環(huán)境格式。它不僅具有可編程性，而且其實(shí)時(shí)運(yùn)行速度可達(dá)每秒數(shù)以千萬(wàn)條復(fù)雜指令程序，遠(yuǎn)遠(yuǎn)超過通用微處理器，是數(shù)字化電子世界中日益重要的電腦芯片。它的強(qiáng)大數(shù)據(jù)處理能力和高運(yùn)行速度，是最值得稱道

2、的兩大特色。b5E2RGbCAP關(guān)鍵詞： 線性卷積； DSP設(shè)計(jì)方案及原理：一.目的：1）加深對(duì)離散線性時(shí)不變系統(tǒng)的理解；2）設(shè)計(jì)線性卷積算法原理；3）熟悉用 TI 公司 DSP芯片進(jìn)行數(shù)字信號(hào)處理和DSP相關(guān)知識(shí)。二.要求：1）給出算法原理；2）寫出主程序；三. 原理1.1/12個(gè)人資料整理僅限學(xué)習(xí)使用2. 線性卷積2/12個(gè)人資料整理僅限學(xué)習(xí)使用3.DSP3/12個(gè)人資料整理僅限學(xué)習(xí)使用4/12個(gè)人資料整理僅限學(xué)習(xí)使用DSP 典型應(yīng)用框圖5/12個(gè)人資料整理僅限學(xué)習(xí)使用主程序MATLAB6/12個(gè)人資料整理僅限學(xué)習(xí)使用C#include<math.h>#include<s

3、tdio.h>#defineNxa/序列xk的長(zhǎng)度#defineNhb/序列hk的長(zhǎng)度#defineLengthNx+Nh-1/卷積或相關(guān)結(jié)果序列的長(zhǎng)度intxNx=。inthNh=。intyLength。voidconvolution(int*,int*,int*>。voidmain(>intk。/調(diào)用卷積程序convolution(x,h,y>。for(k=0。k<=Length-1。k+>fprintf(stdout,"%i",yk> 。/打印卷積計(jì)算結(jié)果7/12個(gè)人資料整理僅限學(xué)習(xí)使用fprintf(stdout,"

4、;n">。voidconvolution(intxNx,inthNh,intyLength>inti,k。for(k=0。k<=Length-1。k+>inttemp=0。for(i=0。i<Nx。i+>if(k-i>>=0&&(k-i><=Nh-1>temp=temp+xi*hk-i。yk=temp。 p1EanqFDPw線性卷積 matlab 實(shí)現(xiàn) <按照理論公式和conv 函數(shù)兩種方法）% 按照理論公式法一：a=4 5 6 7 8 9 。b= 7 5 8 。k=length(a> 。b

5、k=length(b> 。n=length(a>+length(b>-1。 DXDiTa9E3da=a zeros(1,n-k>。b=b zeros(1,n-bk>。y=zeros(1,n> 。for i=1:ntemp=0 。for j=1:nif i>=jtemp=temp+a(j>.*b(i+1-j>elsebreak 。endendy(i>=tempend RTCrpUDGiT8/12個(gè)人資料整理僅限學(xué)習(xí)使用ysubplot(3,1,1>stem(y>% 按照理論公式法二：for i=1:nj=max(1,i+1-

6、bk>:min(i,k>。y2(i>=sum(a(j>.*b(i+1-j>>。endy2subplot(3,1,2>stem(y2> 5PCzVD7HxA% 利用 conv 函數(shù)a=4 5 6 7 8 9 。b= 7 5 8 。y3=conv(a,b> jLBHrnAILgsubplot(3,1,3>stem(y3>用 C 語(yǔ)言實(shí)現(xiàn)線性卷積，/* y(n>=x(n>*h(n>m-length of x(n>。n-length of h(n>。l=m+n-1 length of y(n>*/線性

7、卷積源程序(2018-12-21 17:09:57>轉(zhuǎn)載 標(biāo)簽：分類： matlab雜談%自己做出的線性卷積算法，很難得，想從網(wǎng)上下一個(gè)的，誰(shuí)知道百度文庫(kù)需要20 個(gè)金幣，真是氣憤，9/12個(gè)人資料整理僅限學(xué)習(xí)使用%然后就自己摸索著寫出了這個(gè)線性卷積源碼，解釋相當(dāng)清楚。function yn=xianxingjuanji(x1,x2>%定義函數(shù) mycomvx1=input('x1='>。x2=input('x2='>。N1=length(x1> 。%計(jì)算每個(gè)序列的長(zhǎng)度M=length(x2> 。L=N1+M-1。%卷積后的序

8、列的長(zhǎng)度f(wàn)or(n=1:L>%利用算法實(shí)現(xiàn)卷積，外層循環(huán)次數(shù)Ly(n>=0 。for(m=1:M>%內(nèi)層循環(huán)，用于計(jì)算每個(gè)y(n> 的值k=n-m+1。%因?yàn)?matlab 中序列的下標(biāo)是從1 開始的，所以所有下標(biāo)都 +1if(k>=1&k<=N1>%只用計(jì)算 >1 時(shí)的值， <1 是 x2(m>為零，不用計(jì)算y(n>=y(n>+x2(m>*x1(k>。endendendy1=conv(x1,x2> 。%用公式法計(jì)算卷積subplot(231> 。%輸出原序列 x(1>10/12個(gè)人資

9、料整理僅限學(xué)習(xí)使用stem(x1> 。grid on 。 %輸出加上網(wǎng)絡(luò)表格title('序列 x1'>subplot(232> 。%輸出原序列 x(2>stem(x2> 。grid on 。title('序列 x2'>subplot(233> 。%輸出用算法計(jì)算的結(jié)果stem(y> 。grid on 。title('線性卷積 '>subplot(234> 。%輸出用公式轉(zhuǎn)換的結(jié)果stem(y1> 。grid on 。title('conv直接卷積 '>線性卷

10、積 C 程序#include/* y(n>=x(n>*h(n>m-lengthofx(n>。n-lengthofh(n>。l=m+n-1lengthofy(n>*/voidconv(floatx,intm,floath,intn,floaty,intl>inti,j,k。for(i=0。iy=0.0。for(j=0。jk=i-j。if(k>=0&&k<N>y+=x*h。11/12個(gè)人資料整理僅限學(xué)習(xí)使用voidmain(>/參數(shù)m,n,a1,a2,b1,b2,x,h可任意改動(dòng)constintm=3。/length

11、ofx(n>constintn=3。/lengthofh(n>constintl=m+n-1。/lengthofy(n>const inta1=0。 /x(n>左 下標(biāo), 表示x(n>=x(0>,x(1>,x(2>constinta2=2。/x(n>右下標(biāo)constintb1=4。/h(n>左下標(biāo),h(n>=h(4>,h(5>,h(6>constintb2=6。/h(n>右下標(biāo)intc1。/卷積結(jié)果序列的左下標(biāo)intc2。/卷積結(jié)果序列的右下標(biāo)inti。floatx=1,1,1。floath=1,1,1。floa