matlab带通采样10.7MHzFM信号

以下是matlab实现10.7MHz FM信号带通采样的示例代码：

% 定义采样频率和采样时间
fs = 50e6; % 采样频率为50MHz
Ts = 1/fs; % 采样时间间隔

% 生成10.7MHz的FM信号
fc = 10.7e6; % 载波频率为10.7MHz
fm = 100e3; % 调制信号频率为100kHz
kf = 10e3; % 调频系数为10kHz/V
t = 0:Ts:1/fm; % 采样时间序列
m = cos(2*pi*fm*t); % 产生调制信号
x = cos(2*pi*(fc+kf*m).*t); % 产生FM信号

% 带通采样
fs1 = 20e6; % 采样频率1为20MHz
fs2 = 30e6; % 采样频率2为30MHz
fs3 = 40e6; % 采样频率3为40MHz
Ts1 = 1/fs1; % 采样时间间隔1
Ts2 = 1/fs2; % 采样时间间隔2
Ts3 = 1/fs3; % 采样时间间隔3
n = 0:Ts:1/fc; % 采样时间序列
y1 = x .* cos(2*pi*fs1*n); % 采样1
y2 = x .* cos(2*pi*fs2*n); % 采样2
y3 = x .* cos(2*pi*fs3*n); % 采样3

% 绘制原始信号和三个采样信号
subplot(2,2,1);
plot(t,x);
title('原始信号');
xlabel('时间(s)');
ylabel('幅度(V)');
subplot(2,2,2);
plot(n,y1);
title('采样1信号');
xlabel('时间(s)');
ylabel('幅度(V)');
subplot(2,2,3);
plot(n,y2);
title('采样2信号');
xlabel('时间(s)');
ylabel('幅度(V)');
subplot(2,2,4);
plot(n,y3);
title('采样3信号');
xlabel('时间(s)');
ylabel('幅度(V)');