matlab ，已有数据，如何滤波？

思路就是设计一个低通滤波器，分界点就是50Hz，以下是一个巴特沃斯低通滤波器的程序figure; %绘制巴特沃斯滤波器频响、相位曲线及信号滤后曲线%低通滤波器技术要求，假设采样频率为fs=500， 拟定%通带截止频率为40Hz，阻带下限截止频率为60Hz%通带衰减为0.25dB，阻带衰减为30dBT=A;data=B;fs=500;Wp=2*pi*40/fs;Ws=2*pi*60/fs;Rp=0.25;Rs=30;
Omip=Wp/pi;Omis=Ws/pi; %归一化技术要求
[N，Wn]=buttord(Omip，Omis，Rp，Rs); %确定滤波器的阶数
disp(['The order of Butterworth Filtering is '，num2str(N)]);
[b，a]=butter(N，Wn); %确定Butterworth滤波器转移函数系数向量
[H，w]=freqz(b，a，512); %512点复频响应
subplot(221); %绘制幅频响应曲线
plot(w/pi，abs(H));title('幅频响应');
xlabel('w(/pi)');ylabel('|H(jw)|');axis([0，1，0，1.1]);
%set(gca，'XTickMode'，'manual'，'XTick'，[0，Omip，Omis，1]);
%set(gca，'YTickMode'，'manual'，'YTick'，[0，10^(-Rs/20)，10^(-Rp/20)，1]);
%grid;
subplot(222); %绘制相频响应曲线
plot(w/pi，angle(H)/pi);title('相频响应');
xlabel('w(/pi)');ylabel('pha(/pi)');axis([0，1，-1，1]);
%set(gca，'XTickMode'，'manual'，'XTick'，[0，Omip，Omis，1]);
%grid;
subplot(223);
plot(T，data);title('原信号波形图');axis([0，2.2，-1.1，1.1]);
subplot(224);
data1=filter(b，a，data);
plot(T，data1);title('滤波后波形图');