#include <stdio.h>
#include "iirbcf.c"
#include "gainc.c"

int main(int argc, char *argv[])               //帶通濾波器
{
	int i,k,n,ns;                              //n是每級濾波器的級數（取4）；ns是級聯數（取2）
	double a[50],b[50],x[1000],y[1000]         //a，b是濾波器系數存放數組（二維）；x，y是輸入輸出序列的存放數組
	double f1,f2,f3,f4,fs,flc,fls,fhc,fhs,freq,db  //flc:通帶下邊界頻率；fhc:通帶上邊界頻率；flc:阻帶下邊界頻率；flc:阻帶上邊界頻率；
	                                               //fs:采樣頻率;db:阻帶衰減
    char fname[40];
	FILE *fp;
    db=50;
	ns=2;
	n=4;
	flc=10;
	fhc=500;
	f1=0;
	f2=flc;
	f3=fhc;
	f4=0;
	fs=2000;
	f1=f1/fs;
    f2=f2/fs;
	f3=f3/fs;
	f4=f4/fs;
	iirbcfpass(ns,n,f1,f2,f3,f4,db,b,a);
	for (k=0;k<ns ;k++ )
	{
		printf("\nsection %d\n\n",k+1);
		for (i=0;i<=n ;i++ )
		{
			printf("b[%d][%d]=%10.7if",k,i,b[k*(n+1)+i]);
			if (((i%2==0)&&(i!=0))
			{
				printf("\n");
			}
		}
		printf("\n");
	}
	printf("\nenter file name f magnitude response\n");        //輸入響應
	scanf（"%s,fname"）;
	if((fp=fopen(fname,"w"))==NULL)
	{
		printf("cannot open file %s\n",fname);
		exit(0);
	}
	gainc(b,a,n,ns,x,y,1000,2)                                          //2代表是幅頻響應+相頻響應
	for (i=0;i<100 ;i++ )
	{
		freq=i*0.5/1000.0;
		fprintf(fp,"%if    %if\n",freq,x[i]);
	}
	fclose(fp);
}