matlab/简单的滤波程序

 x_old =  1:100

y_old = sin(x_old)+sin(100*x_old)

subplot(1,4,1)

plot(x_old, y_old)

x_new = 1:200

y_old = double(y_old)

z_fft = fft(y_old,200)

z_fft = fftshift(z_fft)

subplot(1,4,2)

plot(x_new, z_fft)

z_fft(1:80) = 0

z_fft(120:200) = 0

subplot(1,4,3)

plot(x_new, z_fft)

z_fft = double(z_fft)

z_ifft = ifftshift(z_fft)

z_ifft = double(z_ifft)

y_new = ifft(z_ifft, 200)

for j=1:100

        y_2(j) = y_new(j)

 end

subplot(1,4,4)

plot(x_old, y_2) 

结果：