複數卷積c語言,c語言計算兩個複數之積

本文目錄一覽：

• 1、matlab中conv（卷積）函數的C語言實現的源代碼？請求大蝦幫助。
• 2、C語言用子函數實現卷積
• 3、如何用C語言實現數組的卷積過程~~~

matlab中conv（卷積）函數的C語言實現的源代碼？請求大蝦幫助。

function a=myconv(b,c)

bs=size(b);

cs=size(c);

i=any(bs-cs);

if i

error(‘error’)

end

i=any(~(bs-1));

if ~i

error(‘error’)

end

ko=0;

if bs(1)bs(2)

b=b’;

c=c’;

ko=1;

end

bs=size(b);

cs=size(c);

ss=2*bs(2)-1;

a=zeros(1,ss);

for i=1:cs(2)

q=zeros(1,i-1);

p=zeros(1,ss-cs(2)+1-i);

ba=[q,c,p];

ma=b(i)*ba;

a=a+ma;

end

if ko

a=a’;

end

end

測試了一下,跟conv計算的結果偏差很小(我測試的結果是10^-15左右),執行效率略低(用cputime 測試rand(1,99),差了0.2964) 商業軟體就是牛啊 真想知matlab中這個函數的源代碼

C語言用子函數實現卷積

conv(int u[],int v[],int w[], int m, int n)

{

   int i, j;

   int k = m+n-1;

   for(i=0; ik; i++)

      for(j=max(0,i+1-n); j=min(i,m-1); j++)

      {   

       w[i] += u[j]*v[i-j];

      } 

}

u[],v[]為原始數組，m，n分別為數組長度，w[]為卷積結果(w[]需初始化為0)，其長度為m+n-1

如何用C語言實現數組的卷積過程~~~

積分為線性卷積，和圓形卷積。而題目是線性卷積，然後是所求的結果個數是上面兩個數組 個數的和減去1

比如上面h數組裡面單元是5 而x數組 是4

所以肯定一點是結果是等於8個數的

result[(sizeof(h) + sizeof(x)) / sizeof(double) – 1];這個就可以說明了

第二個知識點是卷積是怎麼求的。第一步肯定是判斷兩個數組 那個長度長

conv(x, h, sizeof(x) / sizeof(x[0]), sizeof(h) / sizeof(h[0]), result); 就是實現這個目標的。

然後是長度長的放前面

好吧 我換個 數字來就把

x【】=

h【】=

然後卷積 一個是 x0*h0=1；實現語句 是第一個

for (int i = 0; i lenH; i++)

{

for (int j = 0; j = i; j++)

result[i] += x[j] * h[i – j];

}

此時 已經要轉入第二步驟了：

for (int m = lenH; m lenX; m++){

for (int j = 0; j lenH; j++)

result[m] += x[m – j] * h[j];

}

第二部 應該是 x*h+x1*h（1-1）= 這裡得h1 用0代替 但程序里 不是這樣 而是 用x*h=

好吧 我可能設置的h數組不夠長 加入 h有兩個。x有

那麼 結果 應該是x2*y1+x1*y0；

然後是第三部

是說 在要求的 結果 最後幾個數字時候 比如原題裡面 應該是有8個的。但到第二個循環才求到X得長度5個。

所以 後面應該是resual記住 數組下標 比實際小1. 所以

是這樣的

用 for (int n = lenX; n lenX + lenH – 1; n++){

for (int j = i – lenX + 1; j lenH; j++)

result[n] += x[n – j] * h[j];

}裡面的i 要改成n

for (int n = lenX; n lenX + lenH – 1; n++){

for (int j = n – lenX + 1; j lenH; j++)

result[n] += x[n – j] * h[j];

}

然後 是這樣分析的

結果等於=x(0)h(5-0)+x(1)h(5-1)+x(2)h(5-2)+x(3)h(5-3)=x(0)h(5)+x(1)h(4)+x(2)h(3)+x(3)h(2) 記住 數組不夠的地方 用0代替

copy(result, result[8], ostream_iteratordouble(cout, ” “)); 這個函數 就不想說了 自己去看stl 演算法吧

另外,虛機團上產品團購,超級便宜