10.7MHz的FM中频带通采样的matlab

代码：

fs = 44.1e3; % 采样频率
fc = 10.7e6; % 中频频率
bw = 100e3; % 带宽
t = 0:1/fs:0.1; % 采样时间
x = cos(2*pi*fc*t + 0.5*pi*cos(2*pi*1e3*t)); % FM信号

% 中频带通滤波器
[b,a] = butter(6, [fc-bw/2, fc+bw/2]/(fs/2), 'bandpass');
y = filter(b, a, x);

% 采样
fs1 = bw; % 采样频率为带宽
ts = 1/fs1; % 采样时间间隔
n = length(y)/fs1; % 采样点数
t1 = 0:ts:(n-1)*ts; % 采样时间
z = y(1:fs/fs1:end); % 采样

% 频谱分析
N = length(z); % 采样点数
f = (0:N-1)*fs1/N; % 频率坐标
Z = fft(z)/N; % 快速傅里叶变换
Z = Z(1:N/2+1); % 取一半
Z(2:end-1) = 2*Z(2:end-1); % 归一化
plot(f/1e3,abs(Z)); % 绘制频谱图
xlabel('频率/kHz');
ylabel('幅度');
title('FM信号频谱');

说明：

1. 首先定义采样频率、中频频率、带宽和采样时间，生成一个带有调频的信号。
2. 对信号进行中频带通滤波，保留调频信号的中心频率及其周围的一定带宽。
3. 对滤波后的信号进行采样，采样频率为带宽。
4. 对采样后的信号进行频谱分析，得到频谱图。

注意事项：

1. 由于中频带通滤波器的阶数和截止频率的选择会影响滤波效果，需要根据具体情况进行调整。
2. 采样频率的选择需要满足奈奎斯特采样定理，即采样频率应该大于信号带宽的两倍。在本例中，采样频率为100kHz，信号带宽为10kHz，符合要求。
3. 在频谱分析时，需要进行归一化处理，得到正确的幅度值。