10.7Mhz的FM中频广播信号进行带通采样的matlab

代码如下：

% 设置采样率和采样点数
fs = 50e3;
N = 1024;

% 生成10.7MHz的FM中频信号
fc = 10.7e6;
fm = 1e3;
fsig = 2*pi*fm*t;
x = cos(2*pi*fc*t + cumsum(fsig));

% 进行带通滤波
f_low = 10e3;
f_high = 20e3;
[b,a] = butter(4, [f_low, f_high]/(fs/2), 'bandpass');
y = filter(b, a, x);

% 进行采样
t_sample = (0:N-1)/fs;
y_sample = interp1(t, y, t_sample);

% 绘制时域图和频谱图
figure;
subplot(2,1,1);
plot(t_sample, y_sample);
xlabel('Time (s)');
ylabel('Amplitude');
title('Sampled FM IF signal');
subplot(2,1,2);
f = linspace(-fs/2,fs/2,N);
Y = fftshift(fft(y_sample));
plot(f,abs(Y));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Spectrum of sampled FM IF signal');

相关问题

10.7Mhz的FM中频广播信号进行AD带通采样的matlab

代码如下：

% 设置采样频率和采样时长
fs = 20e6; % 采样频率20MHz
T = 1/fs; % 采样间隔
L = 2000; % 采样时长2000个采样点
t = (0:L-1)*T; % 时间序列

% 生成10.7Mhz的FM中频广播信号
fc = 10.7e6; % 中心频率10.7MHz
df = 20e3; % 频率偏移20kHz
kf = 2*pi*df; % 调频系数
Ac = 1; % 广播信号的振幅
fm = cos(2*pi*fc*t + kf*cumsum(randn(size(t)))); % FM调制信号
xm = Ac*cos(2*pi*fc*t); % 广播信号
x = xm.*fm; % FM中频广播信号

% 进行AD带通采样
f1 = 10.695e6; % 带通下限频率10.695MHz
f2 = 10.705e6; % 带通上限频率10.705MHz
Wn = [f1/(fs/2) f2/(fs/2)]; % 归一化截止频率
b = fir1(100, Wn); % FIR滤波器系数
y = filter(b, 1, x); % 进行滤波

% 绘制原始信号和采样后的信号的时域和频域图像
figure;
subplot(2,2,1);
plot(t, x);
xlabel('时间(s)');
ylabel('幅度(V)');
title('原始信号时域图像');
subplot(2,2,2);
plot(t, y);
xlabel('时间(s)');
ylabel('幅度(V)');
title('采样后的信号时域图像');
subplot(2,2,3);
Y = fft(x);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:(L/2))/L;
plot(f, P1);
xlabel('频率(Hz)');
ylabel('幅度(V)');
title('原始信号频域图像');
subplot(2,2,4);
Y = fft(y);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:(L/2))/L;
plot(f, P1);
xlabel('频率(Hz)');
ylabel('幅度(V)');
title('采样后的信号频域图像');

FM10.7MHz的中频广播信号用ADC进行带通采样的matlab代码

以下是一个简单的MATLAB代码样例，用于进行FM广播信号的中频带通采样：

% 设置采样参数
fs = 44100;             % 采样频率
fc = 10.7e6;            % 中频频率
bw = 200e3;             % 带宽

% 生成带通滤波器
[b, a] = butter(10, [((fc-bw/2)/(fs/2)),((fc+bw/2)/(fs/2))], 'bandpass');

% 读取FM信号文件
filename = 'fm_signal.wav';
[x, fs] = audioread(filename);

% 中频带通采样
y = filter(b, a, x);

% 绘制原始信号和采样信号的频谱
N = length(x);
X = abs(fft(x))/N;
Y = abs(fft(y))/N;
f = linspace(0, fs, N);
figure;
subplot(2,1,1);
plot(f, X);
xlim([0, fc*2]);
title('原始信号频谱');
subplot(2,1,2);
plot(f, Y);
xlim([0, fc*2]);
title('中频带通采样信号频谱');

此代码将读取一个名为“fm_signal.wav”的文件，中频带通滤波后将其采样，并绘制原始信号和采样信号的频谱。请注意，此代码仅供参考，并且可能需要根据特定应用程序进行调整。