如何用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 = 0i <lenHi++)

{

for (int j = 0j <= ij++)

result[i] += x[j] * h[i - j]

}

此时 已经要转入第二步骤了：

for (int m = lenHm <lenXm++){

for (int j = 0j <lenHj++)

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 = lenXn <lenX + lenH - 1n++){

for (int j = i - lenX + 1j <lenHj++)

result[n] += x[n - j] * h[j]

}里面的i 要改成n

for (int n = lenXn <lenX + lenH - 1n++){

for (int j = n - lenX + 1j <lenHj++)

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_iterator<double>(cout, " ")) 这个函数 就不想说了 自己去看stl 算法吧

另外,虚机团上产品团购,超级便宜

线返闭性卷积在matlab中编写可使用以下方法：

Matlab中的漏搏裂线性卷积函数conv(xn,hn) 其中xn、hn的点数分别为N1、N2，输出点数为 N1+N2-1 。

upfirdn(xn,hn,length) 作用为把xn中的每个值乘以序列hn，然后移位相加，length表示了移位的长度。其中xn、hn的点数分别为N1、N2，输出点数为 N2+(N1-1)×length 。

线性卷积(linear convolution) 在时域描述线性系统输入和输出之间关系的一种运算。这种运算在线性系统分析和信号处理中应用很多，通常简称卷积。中文名：数字信号银型处理。

xn 长a ，hn长b ，卷积结果yn 长度谨缺旁为祥橡a+b-1

如果扮码结果的序列yn 还要参与后续的运算，例如和xn相加，那么长度就不一致了，需要截断

一般的xn或者hn非零的序列很短，两端的值为零或者趋于零，例如高斯函数，delta函数，

卷积出来的序列两端也有很多的0值，或者很小的值，可以截断

conv(a,b,'same')返回和a序列一样长的序列取于原来结果的中间

---