北京工业大学---信息处理工程实验报告-FIR.docx

精品文档 ，知识共享 ！！ 北京工业大学 信息处理工程实验报告 FIR 专 业： 电 子 信 息 工 程 姓 名： 学 号： 指导教师： 完成日期： 设定指标设计按格式打印B数组（5分） double b[21],fl=200,fs=1000,fh=0;; p(0,fl,fh); firDesgin(b,21,LOWPASS,RECTANG,fs,fl); getch(); for(w=0;w21;w++)printf(\n b[%d]=%lf,w,b[w]); 画出幅频特性图（5分） window2(低通幅频特性,-fs/2,2,fs/2,-2, hz,|H|); xy2(BLUE); plotxy2(GREEN,2,f,firAbs(f,fs,b,N)); getch(); window2(低通相频特性,-fs/2,5,fs/2,-5, hz,FI); xy2(BLUE); plotxy2(RED,3,f,firPhase(f,fs,b,N)); getch(); 给定信号做采集（5分） for(i=0;i200;i++){x2[i]=COMPLEX((sin(150*2*pi*i/fs)),0);x3[i]=sin(150*2*pi*i/fs);} for(i=200;i512;i++){x2[i]=COMPLEX(0,0);x3[i]=0;} window2(采样,-1,-1,70,5, w,phase); xy2(RED); plotgri1(BLUE,BLUE,x2,512); plotxy2(GREEN,2,i,sin(150*2*pi*i/fs)); getch(); 信号为 前200位sin(150*2*pi*i/fs)，后312位为0 叠接相加法的实现与验证（10分） 实现 叠接相加法：重叠相加法是将待过滤的信号分割成长为?N?的若干段，，每一段都可以和有限时宽单位取样回应作卷积，再将过滤后的各段重叠相加。在实际应用中利用FFT来计算两个序列的圆周卷积从而实现计算其线性卷积，但是常遇到的问题是参加卷积的两个序列的长度相差较大，这样长度小的序列就需要补很多的零点，这样就需要大的存储量，运算时间也会变长。所以常用重叠相加法来解决。 示意图 根据实际数组参数对照，M-1=3，L-1=5，按这个参数设计程序，如下 程序设计： COMPLEX x4[15+3+(8-(3+15)%5)],h4[8],g1[8],p[23],h1[8]; int k=0;//数组的赋值与之前的方法相同 //x3[]的长度为15+3+(8-(3+15)%5)=23， for(i=0;i23;i++)x4[i]=COMPLEX(0,0); for(i=3;i18;i++)x4[i]=COMPLEX(i-2,0); for(i=0;i3;i++)h4[i]=COMPLEX(3-i,0); for(i=3;i8;i++)h4[i]=COMPLEX(0,0); for(i=0;i8;i++)g1[i]=COMPLEX(0,0);//结束各数组的赋值 fft(h4,8,1);//将只有三点的信号2进行fft变换 for(j=0;j=23-3;j=j+3){for(i=0;i8;i++)g1[i]=COMPLEX(0,0);//每次循环前清零 for(i=0;i3;i++){g1[i]=x4[i+j];}//向g数组里放三个数 fft(g1,8,1);//fft变换 for(i=0;i8;i++)h1[i]=g1[i]*h4[i];//频域相乘 fft(h1,8,-1);//逆变换回来相当于卷积 for(i=0;i6;i++){p[i+j]=p[i+j]+h1[i];}}//每次将重的点进行叠加，结果放入p for(i=3;i23;i++)printf(p[%d]=%f+%fj\n,i-3,p[i].r,p[i].i); 叠接相加法结果 验证 与示例叠接舍去法结果对比 叠接舍去法的实际滤波信号（10分） 程序：（x3为输入，b为转移函数，y为输出） for(n=0;n512;n++){y1[n]=COMPLEX(0,0);y[n]=0;} double r[52];int j;int o; for(j=0;j=512-52;j=j+32){ for(o=0;o32;o++)g[o]=x3[o+j]; convol(g,32,b,21,r,52); for(o=10;o42;o++)y[o+j]=r[o]; if(j==0){for(o=0;o10;o++)y[o]=r[o];} for(n=0;n512;n++){y1[n