DSP大作业FIR.docVIP

要求： 1）1人1题，按本人学号的最末位选题。学号的最末位为0、3、6的同学选做第一题，学号的最末位为1、4、7、9的同学选做第二题，学号的最末位为2、5、8的同学选做第三题。 2）单独提交综合实验报告电子版。报告内容为DSP程序、软件仿真结果（输入信号及图形、滤波器类型及系数、滤波后的输出信号图形 或 频谱图）、结果分析及总结。 题目： 设输入信号为x(t)=cos(2πft) + cos(8πft) + cos(12πft)，x(n)利用MATLAB仿真分析得到，采样频率fs根据题目要求自己选择。这里的频率f为自己学号的最末两位×20 Hz。 FIR滤波器的DSP实现(低通、高通、带通、带阻) 要求：低通功能为输入x(t)=cos(3πft)+ cos(7πft)+ cos(10πft)，输出y(t)=cos(3πft)；高通功能为输入x(t)= cos(3πft)+ cos(7πft)+ cos(10πft)，输出y(t)=cos(10πft)；带通功能为输入x(t)= cos(3πft)+ cos(7πft)+ cos(10πft)，输出y(t)= cos(7πft)；带阻功能为输入x(t)= cos(3πft)+ cos(7πft)+ cos(10πft)，输出y(t)= cos(3πft)+ cos(10πft)。 答： 首先利用MATLAB程序设计出滤波器的系数。 程序说明： 1.A/D采样率fs定义为20KHZ； 2.滤波器的阶数定为79； 3.采用窗函数设计法，利用MATLAB计算出滤波器的系数，然后将该系数应用到DSP的滤波程序中去。 以下是MATLAB程序代码: b1=fir1(79,0.16); %80 point hamming window FIR filter, Wn=0.16π=2π*1550/20000 b2=fir1(79,0.26,high); %80 point high pass FIR filter, Wn=0.26π=2π*2635/20000 b3=fir1(79,[0.16 0.26]); %band pass 0.1-0.3 b4=fir1(79,[0.16 0.26], stop); %band stop 0.16-0.26 fvtool(b1,1) %Flter Visualization Tool可心以查看幅度响应、冲激响应等 fvtool(b2,1) fvtool(b3,1) fvtool(b4,1) h=buffer(round(b1*32768),8) %参数*2^15,把小数点移最高位后面后取整，然后排成8个一行便于输出。 csvwrite(table1.txt,h) %参数表输出到文件中 h=buffer(round(b2*32768),8) %参数*2^15,把小数点移最高位后面后取整，然后排成8个一行便于输出。 csvwrite(table2.txt,h) %参数表输出到文件中 h=buffer(round(b3*32768),8) %参数*2^15,把小数点移最高位后面后取整，然后排成8个一行便于输出。 csvwrite(table3.txt,h) %参数表输出到文件中 h=buffer(round(b4*32768),8) %参数*2^15,把小数点移最高位后面后取整，然后排成8个一行便于输出。 csvwrite(table4.txt,h) %参数表输出到文件中 %然后只要把参数表从table.txt中拷到代码中，每行前面加上 .word t=0:1/20000:1/20; %采样率为20KHz,一共采集1000个点 x=cos(2*pi*930*t)+cos(2*pi*2170*t)+ cos(2*pi*3100*t); %930Hz、2170Hz、3100Hz的信号相叠加 xx=buffer(round(x/3.3*32768),1) %参数*2^15,把小数点移最高位后面后取整，然后排成1个一行便于输出。 csvwrite(table5.txt,h) %参数表输出到文件中，输入信号x的离散值保存在文件table5.txt中 y1=filter(b1,1,x); %低通滤波结果 y2=filter(b2,1,x); %高通滤波结果 y3=filter(b3,1,x); %带通滤波结果 y4=filter(b4,1,x); %带阻滤波结果 %画出输入输出波形 figure(1); subpl