gain example.c

来自「自己编的数值计算方法的c语言实现源程序。」· C语言 代码 · 共 37 行

#include"stdio.h"
#include"gain.c"
/*传递函数为：H(z)=(-0.1z^-1)/(1+0.9z^-2)*/
void main()
{	int i;
	double a[]={1.0,0.0,0.9},b[]={0.0,-0.1};
	double f,x[300],y[300];
	FILE *fp;
	gain(b,a,1,2,x,y,300,1);
	
	/*在gain,txt 中输出x[i]、y[i]数据*/
	if(( fp=fopen("gain.txt","w"))==NULL)
	{	printf("cannot open file 'gain.txt'!\n");
		exit(0);
	}
	for(i=0;i<300;i++)
	{	f=i*0.5/299;/*ATTENTION!!!!!!!!!!!!!!!!!!!! THE APRAMETER HERE (0.5) MUST BE THE SAME AS IN GAIN.C!!!!!!*/
		fprintf(fp,"%lf ",f);
	}
	fprintf(fp,"\n");
	for(i=0;i<300;i++)
	{	fprintf(fp,"%lf ",x[i]);
	}
	fprintf(fp,"\n");
	for(i=0;i<300;i++)
	{	fprintf(fp,"%lf ",y[i]);
	}
	fclose(fp);
}