本文目錄一覽:
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語言的問題(卷積)
裡面的你可以到Turbo C裡面 用F7單步執行,就可以看到相應的步驟了。
其他的編譯軟體也有相應的調試方法
卷積積分的C語言程序
代碼並不長,建議單步運行去看看,看走到哪一步時結果與你想像的不一致,再分析那一段代碼…
如何用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 演算法吧
另外,虛機團上產品團購,超級便宜
信號與系統初學,實在不明白卷積的時候對誰積分的問題….
你覺得卷積積分實在不好理解,那我們另外尋找一條理解途徑。卷積積分說白了就是為求電路響應,為了理解抽象的卷積積分公式:
e(t)※h(t)=∫ e(a)h(t-a)da
我們可選一條簡單題目,從一題多解殊途同歸的角度來理解卷積方法之正確性。例如表達式為e^(-αt)的電流源接到R//C的電路,求C上電壓響應。我們不妨用三種方法求解。①用微分方程求解電容電壓響應。②用拉氏變換求電容電壓響應。③用卷積積分求解電容電壓響應。結果發現①②③答案完全一致: Uc(t)=( R/1-αRC)[ e^(-αt)-e^(-t/RC) ]。這樣就從殊途同歸角度理解了卷積積分求解電路響應確實可行,且卷積思想方法正確無疑。
原創文章,作者:GEXOE,如若轉載,請註明出處:https://www.506064.com/zh-tw/n/331176.html