如何用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 算法吧

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

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、把3个声明函数的地方的注释去掉了

2、把函数定义时的参数名字由m改成m1（因为m已经被定义成常量)

下面是修改后的源代码:

#includestdio.h

#includestdlib.h

#includemath.h

#define N 1199

#define m 81

#define NN 1500

void main()

{

FILE *fp1,*fp2,*fp3;

double s[NN],b[NN],rxx[2397],y[NN],y1[NN],a;

int i,k;

void conv(double w[],int m1,double r[],int n,double s[],int l);

void cor(double *x,int n,double *h,int m1,double *y);

int lvs(double t[], int n,double b[],double x[]);

if((fp1=fopen("E:\\Seis.txt","r"))==NULL) //读地震合成记录

printf("Can't open file.\n");

else

{

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

{ fscanf(fp1,"%f",a);

s[i]=a;}

}

for(i=0;im;i++) //右端项

{

if(i==0) b[i]=1.0;

else b[i]=0.0;

}

//调用相关函数，计算地震记录的自相关rxx[l]

cor(s,N,s,N,rxx);

//调用莱文森算法，计算反滤波因子,并输出

k=lvs(rxx,m,b,y);

printf("%d\n",k);

fp2=fopen("E:\\fanzibo.txt","w");

if(fp2==NULL) {printf("can't open file\n"); exit(0);}

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

{ fprintf(fp2,"%f\n",y[i]);}

fclose(fp2);

//调用褶积函数，计算反射系数

conv(y,m,s,N,y1,m+N-1);

fp3=fopen("E:\\fanzheji.txt","w");

if(fp3==NULL) {printf("can't open file\n"); exit(0);}

for(i=0;im+N-1;i++)

{ fprintf(fp3,"%f\n",y1[i]); }

fclose(fp3);

}

void conv(double w[],int m1,double r[],int n,double s[],int l) //褶积

{ int k,i;

for(k=0;kl;k++)

{ s[k]=0.0; //初始化

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

{if(k-i=0k-i=n-1)

s[k]=s[k]+w[i]*r[k-i];}

}

}

//相关计算（计算x[n]与h[m]的互相关，结果为y[n])

void cor(double *x,int n,double *h,int m1,double *y)

{

int i,j;

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

{

y[i]=0.0;

for(j=0;jm;j++)

{

if(i=j)

y[i]+=x[j]*h[j-i];

}

}

}

int lvs(double t[], int n,double b[],double x[])

{

int i,j,k;

double a,beta,q,c,h,*y,*s;

s=(double *)malloc(n*sizeof(double));

y=(double *)malloc(n*sizeof(double));

a=t[0];

if(fabs(a)+1.0==1.0)

{

free(s);free(y);

printf("fail\n");

return(-1);

}

y[0]=1.0;x[0]=b[0]/a;

for(k=1;k=n-1;k++)

{

beta=0.0;q=0.0;

for(j=0;j=k-1;j++)

{

beta=beta+y[j]*t[j+1];

q=q+x[j]*t[k-j];

}

if(fabs(a)+1.0==1.0)

{

free(s);free(y);

printf("fail\n");

return(-1);

}

c=-beta/a;s[0]=c*y[k-1];y[k]=y[k-1];

if(k!=1)

for(i=1;i=k-1;i++) s[i]=y[i-1]+c*y[k-i-1];

a=a+c*beta;

if(fabs(a)+1.0==1.0)

{

free(s);free(y);

printf("fail\n");

return(-1);

}

h=(b[k]-q)/a;

for(i=0;i=k-1;i++)

{

x[i]=x[i]+h*s[i];

y[i]=s[i];

}

x[k]=h*y[k];

}

free(s);free(y);

return(1);

}

当前标题：c语言有褶积函数吗 c语言函数详解
链接URL：http://scgulin.cn/article/ddesjdo.html