10.7MHz的FM中频广播信号用ADC进行带通采样的matlab代码

以下是可能的解决方案：

% 假设输入信号已经存储在变量x中，采样率为Fs
% 带通滤波器参数
fc = 10.7e6; % 中频频率
bw = 200e3; % 带宽
f1 = fc - bw/2; % 通带左边界频率
f2 = fc + bw/2; % 通带右边界频率
n = 6; % 阶数
% 设计带通滤波器
[b,a] = butter(n,[f1 f2]/(Fs/2),'bandpass');
% 采样
y = filter(b,a,x);
% ADC采样
bits = 12; % 采样位数
Vref = 2.5; % 参考电压
q = Vref/(2^bits-1); % 量化步长
y_adc = round(y/q)*q; % 量化
% 绘制采样前后的信号频谱
N = length(x);
X = fft(x,N);
Y = fft(y,N);
Y_adc = fft(y_adc,N);
f = (0:N-1)*(Fs/N);
subplot(3,1,1);
plot(f,abs(X));
title('原始信号频谱');
subplot(3,1,2);
plot(f,abs(Y));
title('带通滤波后信号频谱');
subplot(3,1,3);
plot(f,abs(Y_adc));
title('ADC采样后信号频谱');

相关问题

10.7MHz的FM中频信号用ADC进行带通采样的matlab代码

以下是一个示例代码，可以用于将10.7MHz的FM中频信号进行带通采样：

% 设置采样率和采样时间
fs = 22.5e6;
t = 0:1/fs:0.01;

% 生成10.7MHz的中频信号
fc = 10.7e6;
fm = 1e3;
m = 2; % 调制指数
Ac = 1;
Am = 0.5;
xc = Ac*cos(2*pi*fc*t);
xm = Am*cos(2*pi*fm*t);
xif = xc.*cos(2*pi*fc*t) - xm.*sin(2*pi*fc*t);

% 带通滤波器参数
fpass = 150e3;
fstop = 200e3;
wpass = fpass/(fs/2);
wstop = fstop/(fs/2);
n = 6; % 阶数

% 设计带通滤波器
[b,a] = butter(n, [wpass, wstop], 'bandpass');

% 用带通滤波器进行采样
y = filter(b, a, xif);

% 绘制时域和频域图像
figure;
subplot(2,1,1);
plot(t, y);
xlabel('Time (s)');
ylabel('Amplitude');
title('Bandpass Sampling of 10.7MHz FM IF Signal');
subplot(2,1,2);
f = -fs/2:fs/length(y):fs/2-fs/length(y);
Y = fftshift(fft(y));
plot(f, abs(Y));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Frequency Domain of Bandpass Sampled Signal');

10.7MHzFM中频信号的AD带通采样的MATLAB

代码如下：

% 10.7MHz FM中频信号的AD带通采样

% 设置采样率和采样时间
fs = 20e6;
T = 1/fs;

% 设置带宽和中心频率
B = 200e3;
fc = 10.7e6;

% 生成FM信号
t = 0:T:1-T;
kf = B/(2*pi);
m = 5*cos(2*pi*1e6*t);
s_FM = cos(2*pi*(fc + kf.*m).*t);

% 显示FM信号时域波形
figure(1)
plot(t, s_FM);
xlabel('Time (s)');
ylabel('Amplitude');
title('FM signal (time domain)');

% 计算AD带通滤波器参数
f1 = fc - B/2;
f2 = fc + B/2;
f3dB = B/2;
N = ceil(3.3*f3dB/fs);
h = fir1(N, [f1 f2]/(fs/2));

% 进行AD带通采样
s_ADC = s_FM(1:N:end);
s_filtered = filter(h, 1, s_ADC);

% 显示AD带通采样后的信号时域波形
figure(2)
plot(t(1:N:end), s_filtered);
xlabel('Time (s)');
ylabel('Amplitude');
title('AD sampled FM signal (time domain)');

运行结果如下图所示：