matlab对10.7MHzFM中频信号进行AD采样

下面是一个简单的 MATLAB 代码，用于对10.7MHz FM中频信号进行AD采样：

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

% 生成10.7MHz FM中频信号
fc = 10.7e6; % 载频频率为10.7MHz
fm = 1e3; % 调制信号频率为1kHz
kf = 2*pi*50e3; % 频偏系数为50kHz/V
t = 0:T:1/fm; % 生成调制信号的时间轴
m = cos(2*pi*fm*t); % 生成调制信号
x = cos(2*pi*fc*t + kf*cumsum(m)*T); % 生成中频信号

% 进行AD采样
n = length(x); % 信号长度
t = (0:n-1)*T; % 生成时间轴
y = x + 0.1*randn(size(x)); % 加入噪声
bits = 12; % 采样精度为12位
Vref = 2; % 参考电压为2V
step = Vref/(2^bits); % 量化步长
q = round(y/step)*step; % 进行量化

在这个代码中，我们首先设置了采样率和采样时间。然后，我们生成了一个10.7MHz FM中频信号，其中包括一个1kHz的调制信号。接下来，我们对信号进行了AD采样，使用了12位精度和2V的参考电压，并加入了一些噪声。最后，我们对采样后的信号进行了量化。

相关问题

matlab对10.7MHzFM中频信号进行AD带通采样

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

% 生成10.7 MHz FM中频信号
f_IF = 10.7e6; % 中频频率
fs = 100e6; % 采样率
t = 0:1/fs:0.01; % 时间向量
x_IF = cos(2*pi*f_IF*t + 2*pi*500*t.^2); % 10.7 MHz FM中频信号

% 设计带通滤波器
f_low = 10e6; % 带通下限频率
f_high = 12e6; % 带通上限频率
Wn = [f_low f_high]/(fs/2); % 归一化截止频率
b = fir1(100, Wn); % FIR滤波器系数

% 进行AD带通采样
x_sampled = downsample(filter(b, 1, x_IF), 10); % 采样率降低到10 MHz

% 绘制采样后的信号频谱
NFFT = 2^nextpow2(length(x_sampled)); % FFT点数
X = fft(x_sampled, NFFT)/length(x_sampled); % FFT
f = fs/2*linspace(0,1,NFFT/2+1); % 频率向量
plot(f/1e6,20*log10(abs(X(1:NFFT/2+1)))); % 绘制幅度谱
xlabel('Frequency (MHz)');
ylabel('Magnitude (dB)');
title('10.7 MHz FM中频信号的采样频谱');

这段代码首先生成了一个10.7 MHz FM中频信号，然后设计了一个带通滤波器，带通下限频率为10 MHz，带通上限频率为12 MHz。最后，使用`filter`函数对中频信号进行滤波，然后使用`downsample`函数将采样率降低到10 MHz。最后，绘制了采样后的信号频谱。

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

代码如下：

% 设置采样频率和带宽
fs = 22.05e6;   % 采样频率
bw = 200e3;     % 信号带宽

% 生成频率为10.7MHz的FM信号
t = 0:1/fs:0.01;
fdev = 75e3;
fm = sin(2*pi*10.7e6*t + 2*pi*fdev*cumsum(randn(size(t))));

% 生成抽样时钟
fclk = bw*4;    % 抽样时钟频率
tclk = 0:1/fclk:0.01;
clk = square(2*pi*fclk*tclk);

% 对FM信号进行AD带通采样
n = length(clk);
sampled_fm = zeros(1,n);
for i=1:n
    if clk(i) > 0
        sampled_fm(i) = fm(round((i-1)*fs/fclk)+1);
    end
end

% 绘制FM信号和采样后的信号的时域和频域图像
figure();
subplot(2,2,1);
plot(t,fm);
title('FM信号');
xlabel('时间(t)');
ylabel('幅值');
xlim([0 0.01]);

subplot(2,2,2);
f = linspace(-fs/2,fs/2,length(fm));
fm_fft = abs(fftshift(fft(fm)));
plot(f/1e6,fm_fft);
title('FM信号频谱');
xlabel('频率(MHz)');
ylabel('幅值');
xlim([-1 1]*20);

subplot(2,2,3);
plot(tclk,clk);
title('抽样时钟');
xlabel('时间(t)');
ylabel('幅值');
xlim([0 0.01]);

subplot(2,2,4);
f = linspace(-fs/2,fs/2,length(sampled_fm));
sampled_fm_fft = abs(fftshift(fft(sampled_fm)));
plot(f/1e6,sampled_fm_fft);
title('采样后的信号频谱');
xlabel('频率(MHz)');
ylabel('幅值');
xlim([-1 1]*100e3);

该代码生成一个频率为10.7MHz的FM信号，并采用22.05MHz的采样频率进行AD带通采样，采样带宽为200kHz。采样时钟频率为采样带宽的4倍，即800kHz。采样过程中，每当抽样时钟上升沿出现时，就对FM信号进行采样，并将采样结果存储到数组中。采样后，绘制FM信号和采样后的信号的时域图和频域图，以观察采样效果。

运行该代码后，可以得到如下图所示的结果：


从图中可以看出，采样后的信号频谱中只有10.7MHz的基带信号和其两侧的镜像频率成分，符合AD带通采样的理论。同时，采样后的信号在时域上也能够恢复出原始的FM信号的形态。