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

发布网友 发布时间：2023-02-11 14:11

我来回答

共5个回答

热心网友 时间：2024-11-26 02:39

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

下面是修改后的源代码:

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

for(i=0;i<m;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<m;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<m+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;k<l;k++)
{ s[k]=0.0; //初始化
for(i=0;i<m;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<n;i++)
{
y[i]=0.0;
for(j=0;j<m;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);
}

热心网友 时间：2024-11-26 02:40

你就是个高手了，直接看不懂，为啥用c语言呢，c#也很强呀

热心网友 时间：2024-11-26 02:40

将出错提示也贴出来才能告诉你啊

热心网友 时间：2024-11-26 02:41

//void conv(double w[],int m,double r[],int n,double s[],int l);
//void cor(double *x,int n,double *h,int m,double *y);
//int lvs(double t[], int n,double b[],double x[]);
这几个函数的声名你干嘛要注释掉？
注释掉了就没法儿调用了~~~~

另，像这样的程序干嘛不用FORTRAN啊。。。

热心网友 时间：2024-11-26 02:41

同上面，函数声明是不能注释掉的。
还有，你函数定义时引用的参数是数组，但是调用时用的是数组名，就是相当于首地址，个人觉得可能是这个错误吧~~~可以试试在函数定义时把数组参数改成地址参数~~~