本文通过Matlab代码演示了如何使用快速傅立叶变换(FFT)进行信号处理，包括信号的频域分析、滤波及逆变换过程。首先生成了一个包含两个频率分量的复合正弦信号，并对其进行了FFT变换，随后通过对变换后的频谱进行滤波处理，最后再通过逆变换得到处理后的信号。
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)
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)
结果：
扫一扫
2266
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
