DSP仿真实验仿真实验截图.doc

4.7main() L=200; N=10;u=0.00002; /* 初始化数据空间 */ for(i=0;i500;i+) xi=0; di=0; ei=0; yi=0; /* 设定输入信号为正弦信号 期望输出设为比实际输入超前两个单位时间的信号 */ for(i=0;iL;i+) xi=(float)10*sin(pi*i/20); di=xi-2; /* 初始化滤波器系数 */ for(i=0;iN;i+) wi=0; /* 梯度算法 */ for(i=N;iL;i+) /*sum=0;*/ for(j=0;jN;j+) yi=yi+wj*xi-j; /*yi=sum;*/ ei=di-yi;/* 计算误差 */ for(j=0;jN;j+)/* 根据误差修正滤波器系数 */ wj=wj+2*u*ei*xi-j; 4.3 main() N=256;for(i=0; i300; i+)x_rei=0;x_imi=0;for(i=0;i=N;i+) /x_rei=exp(-i); /x_rei=1;/x_rei=sin(2*pi*i/64);if(i/32)%2=0)x_rei=(float)(i%32-16)/16;elseif(i/32)%2=1)x_rei=(float)(16-i%32)/16;/*if(i%32=0)x_rei=0;elseif(i/32)%2=0)x_rei=-1;elsex_rei=1;*/ x_imi=0;for(i=0; i300; i+)y_rei=x_rei;y_imi=x_imi;4.5 wi=2.0-2.0*i/(n-1); else wi=2.0*i/(n-1); break;/* m=3时用汉宁窗 */ case 3: for(i=0;in;i+) wi=0.5*(1.0-cos(2.0*pi*i/(n-1); break;/* m=4时用汉明窗 */ case 4: for(i=0;in;i+) wi=0.54-0.46*cos(2.0*pi*(float)i/(n-1); break;/* m=5时用布莱克曼窗 */ case 5: for(i=0;in;i+) wi=0.42-0.5*cos(2.0*pi*i/(n-1)+0.08*cos(4.0*pi*i/(n-1); break; for(i=0;in;i+) hi=hdi*wi; /* 对理想冲击响应加窗（时域相乘） */ p=pi/l; for(k=0;k=l-1;k+) wf=(pi*k)/l; re=0.0; im=0.0; for(i=0;in;i+) re=re+hi*cos(float)i*wf); im=im+hi*sin(float)i*wf); d=sqrt(pow(re,2)+pow(im,2); /* dbk=20.0*log10(d); main() /* 参数初始化 */ for(i=0;i51;i+) hdi=0; hi=0; wi=0; for(i=0;i300;i+) dbi=0; m=2; n=21; /* n=51, 为奇数 */ wc=0.2; /* wc=0.10-0.90 */ l=300; a=(n-1)/2; pi=4.0*atan(1.0); for(i=0;in;i+) /* 理想矩形函数（频域）的时域形式 */ if(i=a) hdi=wc; else b=i-a; hdi=sin(pi*b*wc)/(pi*b); switch(m) /* 根据m的值选择不同的窗函数 */ /* m=1时用矩形窗 */ c