理想低通滤波器matlab代码

下面是一个简单的理想低通滤波器的 MATLAB 代码示例：

% 定义滤波器参数
fc = 1000;  % 截止频率为1000Hz
fs = 8000;  % 采样频率为8000Hz
N = 101;    % 滤波器阶数

% 计算理想低通滤波器的频率响应
h_ideal = (2*fc/fs)*sinc(2*fc*(0:N-1)/fs);

% 使用filter函数进行滤波
x = randn(1, 10000);  % 生成随机信号
y = filter(h_ideal, 1, x);  % 进行滤波

% 绘制频率响应曲线和滤波后的信号波形
f = linspace(0, fs/2, N);  % 计算频率坐标轴
H = fft(h_ideal);  % 计算滤波器的频率响应
Y = fft(y);  % 计算滤波后的信号的频谱
figure;
subplot(2,1,1); plot(f, abs(H(1:N/2+1))); title('理想低通滤波器的频率响应');
subplot(2,1,2); plot(x(1:200), 'b'); hold on; plot(y(1:200), 'r'); title('滤波前后的信号波形');
legend('原始信号', '滤波后的信号');

在这个示例中，我们首先定义了滤波器的参数，包括截止频率、采样频率和滤波器阶数。然后，我们使用`sinc`函数计算理想低通滤波器的频率响应，这里的`sinc`函数是 MATLAB 自带的函数，用于计算`sinc`函数值。接着，我们使用`filter`函数对随机生成的信号进行滤波，并绘制了滤波器的频率响应曲线和滤波后的信号波形。注意，在绘制频率响应曲线时，我们使用了`fft`函数计算滤波器的频率响应和滤波后的信号的频谱。