10.7Mhz的FM中频广播信号进行AD带通采样的matlab

代码如下：

% 设置采样频率和采样时长
fs = 20e6; % 采样频率20MHz
T = 1/fs; % 采样间隔
L = 2000; % 采样时长2000个采样点
t = (0:L-1)*T; % 时间序列

% 生成10.7Mhz的FM中频广播信号
fc = 10.7e6; % 中心频率10.7MHz
df = 20e3; % 频率偏移20kHz
kf = 2*pi*df; % 调频系数
Ac = 1; % 广播信号的振幅
fm = cos(2*pi*fc*t + kf*cumsum(randn(size(t)))); % FM调制信号
xm = Ac*cos(2*pi*fc*t); % 广播信号
x = xm.*fm; % FM中频广播信号

% 进行AD带通采样
f1 = 10.695e6; % 带通下限频率10.695MHz
f2 = 10.705e6; % 带通上限频率10.705MHz
Wn = [f1/(fs/2) f2/(fs/2)]; % 归一化截止频率
b = fir1(100, Wn); % FIR滤波器系数
y = filter(b, 1, x); % 进行滤波

% 绘制原始信号和采样后的信号的时域和频域图像
figure;
subplot(2,2,1);
plot(t, x);
xlabel('时间(s)');
ylabel('幅度(V)');
title('原始信号时域图像');
subplot(2,2,2);
plot(t, y);
xlabel('时间(s)');
ylabel('幅度(V)');
title('采样后的信号时域图像');
subplot(2,2,3);
Y = fft(x);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:(L/2))/L;
plot(f, P1);
xlabel('频率(Hz)');
ylabel('幅度(V)');
title('原始信号频域图像');
subplot(2,2,4);
Y = fft(y);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:(L/2))/L;
plot(f, P1);
xlabel('频率(Hz)');
ylabel('幅度(V)');
title('采样后的信号频域图像');