MATLAB设计FIR滤波器时遇到的问题？

我想用的是Hamming窗 FIR带通滤波器，下限截止频率是0.5Hz，上限截止频率是35Hz，采样频率是1000Hz.下面是我做的部分代码（不知道对不对，麻烦您帮忙看一下），还有我想测试一下，用三个30Hz、50Hz、70Hz的正弦波信号混合出一个信号测试一下这个滤波器，希望在MATLAB中以程序的方式实现，能观察到输入波形和输出波形就行。希望得到您的帮助，谢谢。不知道怎么在MATLAB中进行测试，调用的代码不会写fc1=0.5; fc2=35; fs=1000;  %定义滤波器带通下限截止频率fc1=0.5Hz，带通上限截止频率fc2=35Hz，采样频率fs=1000Hzc=31;                   %定义滤波器阶数32阶w1=2*pi*fc1/fs;w2=2*pi*fc2/fs;      %参数转换，将模拟滤波器的技术指标转换为数字滤波器的技术指标window=hamming(c+1);        %调用hamming窗函数h=fir1(c,[w1/pi w2/pi],window);  %使用标准响应的加窗设计函数fir1freqz(h,1,512);      %数字滤波器频率响应