重庆分公司,新征程启航
为企业提供网站建设、域名注册、服务器等服务
include stdio.h
为阜城等地区用户提供了全套网页设计制作服务,及阜城网站建设行业解决方案。主营业务为网站设计制作、成都网站制作、阜城网站设计,以传统方式定制建设网站,并提供域名空间备案等一条龙服务,秉承以专业、用心的态度为用户提供真诚的服务。我们深信只要达到每一位用户的要求,就会得到认可,从而选择与我们长期合作。这样,我们也可以走得更远!
#include stdlib.h
conv(int *u,int *v, int m, int n);
void main()
{
int i = 0;
int a[4] = {3,10,13,22};
int b[5] = {6,2,15,4,1};
int m = sizeof(a)/sizeof(a[0]);
int n = sizeof(b)/sizeof(b[0]);
conv(a,b,m,n);
}
conv(int u[],int v[], int m, int n)
{
int i, j;
int k = m+n-1;
int w[10]={0};
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];
}
printf("%d/n",w[i]);
}
}
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
main(){
int a[10],b[10],c[10];
int i;
for(i=0;i10;i++){
printf("Enter a arr number %d :",i+1);
scanf("%d",a[i]);
printf("Enter b arr number %d :",i+1);
scanf("%d",b[i]);
}
for(i=0;i10;i++){
c[i]=a[i]*b[i]; //得到相乘结果存入C数组.
}
for(i=0;i10;i++)
printf("\n%d",c[i]);//输出结果
}
积分为线性卷积,和圆形卷积。而题目是线性卷积,然后是所求的结果个数是上面两个数组 个数的和减去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 算法吧
另外,虚机团上产品团购,超级便宜