求c 语言编程高手 帮看看那里出错了

2025-03-10 05:35:56
推荐回答(1个)
回答1:

改了以下几个地方:
1、把3个声明函数的地方的注释去掉了
2、把函数定义时的参数名字由m改成m1(因为m已经被定义成常量)

下面是修改后的源代码:

#include
#include
#include
#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;i { fscanf(fp1,"%f",&a);
s[i]=a;}
}

for(i=0;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;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;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;k { s[k]=0.0; //初始化
for(i=0;i {if(k-i>=0&&k-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;i {
y[i]=0.0;
for(j=0;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