频率波数滤波matlab代码

以下是一个基于频率域的低通滤波器的MATLAB代码示例：

% 定义输入信号
Fs = 1000; % 采样率
t = 0:(1/Fs):1; % 时间向量
f1 = 10; % 信号频率
f2 = 50; % 噪声频率
x = sin(2*pi*f1*t) + sin(2*pi*f2*t); % 输入信号

% 计算频谱
N = length(x); % 信号长度
X = fft(x); % 傅里叶变换
f = (0:N-1)*(Fs/N); % 频率向量
X_mag = abs(X); % 频率幅度谱
X_phase = angle(X); % 频率相位谱

% 设计低通滤波器
Fc = 30; % 截止频率
Hd = zeros(1,N); % 滤波器系数
Hd(f<=Fc) = 1; % 截止频率以下的系数为1

% 应用滤波器
Y = ifft(X .* Hd); % 反傅里叶变换
y = real(Y); % 取实部

% 绘制结果
subplot(3,1,1);
plot(t,x);
xlabel('时间 (s)');
ylabel('幅度');
title('输入信号');
subplot(3,1,2);
plot(f,X_mag);
xlabel('频率 (Hz)');
ylabel('幅度');
title('频率幅度谱');
subplot(3,1,3);
plot(t,y);
xlabel('时间 (s)');
ylabel('幅度');
title('滤波后信号');

这段代码定义了一个包含两个正弦波的信号，其中一个是我们感兴趣的信号，另一个是噪声。然后，它使用FFT来计算信号的频谱，设计一个低通滤波器，将滤波器应用于信号，并绘制结果。请注意，此代码仅进行了低通滤波，如果需要高通滤波或其他类型的滤波器，请相应地更改滤波器系数。