10.7MHzFM中频的带通采样matlab

代码如下：

% 设置参数
fc = 10.7e6; % 中频频率
fs = 44.1e3; % 采样频率
BW = 15e3; % 中频带宽

% 设计带通滤波器
h = fir1(100, [fc-BW/2, fc+BW/2]/(fs/2));

% 读取信号并进行带通滤波
x = audioread('filename.wav');
y = filter(h, 1, x);

% 进行采样
N = length(y);
n = 0:N-1;
ts = 1/fs;
t = n*ts;
y_sampled = y.*cos(2*pi*fc*t);

% 绘制原始信号和采样信号的频谱图
Y = fft(y);
Y_sampled = fft(y_sampled);
f = (0:N-1)*(fs/N);
subplot(2,1,1);
plot(f, abs(Y));
title('原始信号频谱图');
xlabel('频率(Hz)');
ylabel('幅度');
subplot(2,1,2);
plot(f, abs(Y_sampled));
title('采样信号频谱图');
xlabel('频率(Hz)');
ylabel('幅度');