【正文】 HPASS){line2(0,db1,db1)。line2(0,db2,db2)。line2(f1,0,f1,)。line2(f2,0,f2,)。}else{line2(0,db1,db1)。line2(0,db2,db2)。line2(f1,0,f1,)。line2(f2,0,f2,)。line2(f3,0,f3,)。line2(f4,0,f4,)。}getch()。frame2()。for(i=0。iL。i++)//用printf()格式化把系統(tǒng)函數的系數H數組和級連數L顯示出來；{ for(j=0。j2。j++) { for(k=0。k5。k++) { printf(H[%d][%d][%d]=%.2f\n,i,j,k,H[i][j][k])。 } }}getch()。window2(相頻,600,20,600.,20.,hz,db,BLUE,BLUE)。xy2(BLUE)。plotxy2(RED,2,f,btwphHs(f,fs,H,L))。//畫幅頻特性曲線getchar()。frame2()。}void IIRLv(int bandType,double db1,double db2,double fs,double f1,double f2,double f3,double f4){ order=btwOrder(bandType,db1,db2,fs,f1,f2,f3,f4)。//返回值為N L=btwAf(b,order)。//級連數 btwC23(c,bandType,order,db1,fs,f1,f2,f3,f4)。//其中oroder為N,即階數。函數返回的是Wc btwAf2Df(H,L,b,c)。 chouyang(100.)。 JILian_()。 getchar()。}void IIRLv2(int bandType,double db1,double db2,double fs,double f1,double f2,double f3,double f4){ order=btwOrder(bandType,db1,db2,fs,f1,f2,f3,f4)。//返回值為N L=btwAf(b,order)。//級連數 btwC23(c,bandType,order,db1,fs,f1,f2,f3,f4)。//其中oroder為N,即階數。函數返回的是Wc btwAf2Df(H,L,b,c)。 chouyang2(fs)。 JILian_2()。}void main(){IIR(1,2000.,200.,300.,400.,500.)。IIR(2,2000.,200.,300.,400.,500.)。IIR(3,2000.,200.,300.,400.,500.)。IIR(4,2000.,200.,300.,400.,500.)。IIRLv(1,2000.,200.,300.,400.,500.)。//對方波進行低通濾波。IIRLv2(2,2000.,300.,400.,500.,600.)。//對多次正弦諧波進行高通濾波。}void JILian_(void)//實現L級級聯(lián){ int i,j,n。 window2(函數圖形顯示,1,2,200,2,i,x[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,x,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。dft(X,x,N,1)。for(i=0。iN。i++)Y[i]=abs(X[i])。 window2(函數圖形顯示,2,25,200,2,i,X[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,Y,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。 for(i=0。i=L1。i++) { a[0]=0。 for(j=1。j5。j++)//歸一化 {a[j]=H[i][1][j]/H[i][1][0]。}//y[n]系數歸一 for(j=0。j5。j++) {b0[j]=H[i][0][j]/H[i][1][0]。}//x[n]系數歸一 IIRT(x,y,a,b0,N,5)。 for(n=0。nN。n++) x[n]=y[n]。 } for(n=0。nN。n++) y[n]=x[n]。 dft(X,y,N,1)。for(i=0。iN。i++)Y[i]=abs(X[i])。 window2(函數圖形顯示,1,10,200,2,i,y[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,y,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。 window2(函數圖形顯示,10,200,200,1,i,X[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,Y,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。}void JILian_2(void)//實現L級級聯(lián),第二個濾波實現{ int i,j,n。 window2(函數圖形顯示,1,10,200,10,i,x[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,x,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。 dft(X,x,N,1)。for(i=0。iN。i++)Y[i]=abs(X[i])。 window2(函數圖形顯示,1,400,200,1,i,X[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,Y,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。 for(i=0。i=L1。i++) { a[0]=0。 for(j=1。j5。j++)//歸一化 {a[j]=H[i][1][j]/H[i][1][0]。}//y[n]系數歸一 for(j=0。j5。j++) {b0[j]=H[i][0][j]/H[i][1][0]。}//x[n]系數歸一 IIRT(x,y,a,b0,N,5)。 for(n=0。nN。n++) x[n]=y[n]。 } for(n=0。nN。n++) y[n]=x[n]。 dft(X,y,N,1)。for(i=0。iN。i++)Y[i]=abs(X[i])。 window2(函數圖形顯示,1,20,200,20,i,y[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,y,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。 window2(函數圖形顯示,1,600,200,1,i,X[i])。 xy2(RED)。//畫xy軸。 plotgri2(BLUE,RED,Y,N)。//宏函數工具方式,注意x是形式上的宏參數，它指明后面函數的自變量。 getch()。 frame2()。}void chouyang(double fc) //對方波信號進行抽樣{ int i,T=2。 double dt=(*fc)。 double t=0。 N=fc*T。 for(i=0。iN。i++) { y[i]=0。 }// for(i=0。iN。i++) { t=i*dt。 x[i]=COMPLEX(fangbo(t),0)。 }}void chouyang2(double fs){ int i,f1=100,f2=200,f3=300。 for(i=0。iN。i++) { y[i]=0。 }// for(i=0。iN。i++) {////輸入信號x由f1=100,f2=200,f3=300,3個頻率的正弦波合成，幅度分別為1,2,3 x[i]=COMPLEX(1*sin(2*M_PI*f1*i/fs)+2*sin(2*M_PI*f2*i/fs)+3*sin(2*M_PI*f3*i/fs),0)。 }}void plotgri2(COLORREF gridcolor,COLORREF linecolor,COMPLEX p[],int N){int i。HPEN pen1=CreatePen(PS_SOLID,1,gridcolor),oldpen=(HPEN)SelectObject(,pen1)。HPEN pen2=CreatePen(PS_SOLID,1,linecolor)。for(i=0。iN。i++){ moveto2(i,0)。 lineto2(i,p[i].r)。}}double fangbo(double t)//一個周期方波{ int h=1。 if(t=0amp。amp。t=2)return t=h?:。 else if(t2||t0) return 。 else return 。}//通用IIR濾波器 void IIRT(COMPLEX input[],COMPLEX output[],double a[],double b[],int N,int Ne)//N為數據點數，Ne為0狀態(tài)最大負輸入{ int i,n。 for(n=0。nN。n++) { for(i=1。iNe。i++) { if(ni) output[n]+=a[i]*0。 else output[n]=output[n]+a[i]*output[(ni)]。 } for(i=0。iNe。i++) { if(ni)output[n]+=b[i]*0。 else output[n]=output[n]+b[i]*input[(ni)]。 } }}訓練十五 線性卷積計算自定數組：信號復數組 x[32]={1,0,2,0,3,0,4,0,5,0,6,0,7,0,8,0,9,0,10,0,11,0,12,0,13,0,14,0,15,0}。系統(tǒng)復數組 h[32]={3,0,2,0,1,0}。信號數據實際長度為15，系統(tǒng)數據實際長度為3，15+31=17，所以要做32的FFT，故兩數組長度都聲明為32。for(i=0。i15。i++)x2[i]=COMPLEX(i+1,0)。for(i=15。i32。i++)x2[i]=COMPLEX(0,0)。for(i=0。i3。i++)h2[i]=COMPLEX(3i,0)。for(i=3。i32。i++)h2[i]=COMPLEX(0,0)。直接計算程序：//。printf((3) 15個信號數據代表長序列，用一次FFT實現數據序列的卷積\n)。//信號復數組 x[32]={1,0,2,0,3,0,4,0,5,0,6,0,7,0,8,0,9,0,10,0,11,0,12,0,13,0,14,0,15,0}。//系統(tǒng)復數組 h[32]={3,0,2,0,1,0}。//信號數據實際長度為15，系統(tǒng)數據實際長度為3，//15+31=17，所以要做32的FFT，故兩數組長度都聲明為32。COMPLEX x2[32],h2[32]。for(i=0。i15。i++)x2[i]=COMPLEX(i+1,0)。for(i=15。i32。i++)x2[i]=COMPLEX(0,0)。for(i=0。i3。i++)h2[i]=COMPLEX(3i,0)。for(i=3。i32。i++)h2[i]=COMPLEX(0,0)。fft(x2,32,1)。fft(h2,32,1)。for(i=0。i32。i++)x2[i]=x2[i]*h2[i]。fft(x2,32,1)。for(i=0。i20。i++)printf(x2[%d]=%f+j%f,\n,i,x2[i].r,x2[i].i)。getch()。直接做FFT求出的卷積結果疊接舍去實現：//。//，檢驗下面疊接舍去法做出的結果。printf((4) 15個數據用疊接舍去法做出的結果，檢查與第3步結果是否一致\n)。//長序列信號：x2[15+3+(8(3+15)%5)]={[0,0,0,0,0,0,(1,0,2,0],3,0,4,0,5,0),6,0,7,0,//