两序列相关c程序.cpp

来自「这里面包括了地球物理方面常见的三个程序」· C++ 代码 · 共 29 行

#include <stdio.h>
# define nx 3
# define ny 3
# define nrxy 5

void  main()
{	int i,k,n;
    float x[nx],y[ny];
    float rxy[nrxy];
	float sum;
	 x[0]=2;  x[1]=3;  x[2]=-1;
	 y[0]=-1; y[1]=3;  y[2]=1;
	for(i=0; i<nrxy; i++)
	rxy[i]=0.0;
	for(n=0; n<nrxy; n++)
	{   sum=0.0;
		for(i=0;i<nx;i++)
		{
		    k=i-ny+n+1;
		    if(k<0 || k>=ny) continue;
		    sum+=x[i]*y[k];
		}
		rxy[n]=sum;
	}

	for (i=0;i<nrxy;i++)
    printf("%9.1f",rxy[i]);
}