matlab的里对一个FM信号的带通采样后的波形

我们可以使用MATLAB中的信号处理工具箱来模拟一个带通采样后的FM信号波形。下面是一段MATLAB代码，用于生成一个带有500Hz到2kHz的带通滤波器的FM信号，并对其进行采样：

% 设置采样率和时间范围
Fs = 8000; % 采样率为8kHz
T = 1; % 1秒钟的信号

% 生成FM信号
f1 = 500; % 调频频率为500Hz
f2 = 2000; % 最大调频频率为2kHz
t = linspace(0, T, T*Fs); % 生成时间序列
kf = (f2-f1)/(2*T); % 最大调频偏移
beta = kf*t.^2; % 调频偏移
fm = cos(2*pi*(f1*t + beta)); % 生成FM信号

% 设计带通滤波器
Wp = [900 1900]/(Fs/2); % 通带频率范围为900Hz到1900Hz
Ws = [800 2000]/(Fs/2); % 阻带频率范围为800Hz到2000Hz
Rp = 1; % 通带最大允许的纹波为1dB
Rs = 30; % 阻带最小衰减为30dB
[n, Wn] = buttord(Wp, Ws, Rp, Rs); % 计算巴特沃斯滤波器阶数和截止频率
[b, a] = butter(n, Wn, 'bandpass'); % 生成带通滤波器

% 对FM信号进行滤波和采样
fm_filtered = filter(b, a, fm); % 使用带通滤波器进行滤波
fm_sampled = fm_filtered(1:Fs:end); % 对滤波后的信号进行采样

% 绘制原始和采样后的信号波形
figure;
subplot(2,1,1);
plot(t, fm);
title('原始FM信号');
xlabel('时间 (s)');
ylabel('幅度');
subplot(2,1,2);
stem(0:1/Fs:T-1/Fs, fm_sampled);
title('采样后的FM信号');
xlabel('时间 (s)');
ylabel('幅度');

运行上述代码后，会生成一个包含原始FM信号和采样后的FM信号的图形。下图展示了生成的结果：