clc;
clear all;
close all;
fs=1000;
t=0:1/fs:1000/fs;
signal=sin(2*pi*50*t);
figure(1);
plot(t,signal,'-*');
N=30; %%滤波器结束为N=30;
hold on;
Hd = my_lowpass;
output_1=filter(Hd,signal);
output_2=filter(Hd,[signal zeros(1,N/2)]);
plot(t,output_2(1,(N/2+1):end),'-*r');
title('滤波后的波形');
figure(2);
plot(t,signal,'-*b');
hold on;
plot(t,output_1,'-*r');
![](https://img-blog.csdn.net/20170908151510644?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvZGFtb25fbWE=/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center)