fir4.m

来自「这里面有很多dsp应用程序源码」· M 代码 · 共 13 行

b=fir1(79,0.2,hamming(80));	%80 point hamming window FIR filter, Wn=0.2
fvtool(b,1)			%Flter Visualization Tool可心以查看幅度响应、冲激响应等
h=buffer(round(b*32768),8)' %参数*2^15,把小数点移最高位后面后取整，然后排成8个一行便于输出。
csvwrite('table.txt',h)    %参数表输出到文件中
                           %然后只要把参数表从table.txt中拷到代码中，每行前面加上' .word '
x=4000*[ones(1,40) zeros(1,60)]'; %测试数据，方波，4000可当作十六进制的小数0.5
csvwrite('in.txt',x)        %保存测试数据到文件中，以便CCS仿真时调用
y=filter(b,1,x);                   %滤波结果
%画出输入输出波形
subplot(2,1,1);
plot(x)
subplot(2,1,2)
plot(y)