gainc example.c

自己编的数值计算方法的c语言实现源程序。

#include"stdio.h"
#include"gainc.c"
void main()
{	int i,n,ns;
	double a[2][3]={{1.0,-0.8,0.0},{1.0,-0.95,0.9}};/*H1(z)=0.2/(1-0.8z^-1);H2(z)=(0.05-0.1z^-1)/(1-0.95z^-1+0.9z^-2)*/
	double b[2][3]={{0.2,0.0,0.0},{0.05,-0.1,0.0}};
	double f,x[300],y[300];
	FILE *fp;
	n=2;
	ns=2;
	gainc(b,a,n,ns,x,y,300,1);
	if((fp=fopen("gainc.txt","w"))==NULL)
	{	printf("cannot open file 'gainc.txt'!\n");
		exit(0);
	}
	for(i=0;i<300;i++)
	{	f=i*0.5/299;
		fprintf(fp," %lf",f);
	}
	fprintf(fp,"\n");
	for(i=0;i<300;i++)
	{	fprintf(fp," %lf",x[i]); 
	}
	fclose(fp);
}