10.7Mhz的FM中频广播信号进行带通采样的matlab
时间: 2023-11-03 07:32:01 浏览: 35
代码如下:
% 10.7 MHz FM intermediate frequency signal
fc = 10.7e6; % carrier frequency
fm = 15e3; % modulation frequency
kf = 75e3; % frequency deviation constant
fs = 500e3; % sampling frequency
t = 0:1/fs:0.01; % time axis
% Generate FM signal
x = cos(2*pi*fc*t + 2*pi*kf*cumsum(sin(2*pi*fm*t)));
% Bandpass sampling
wc = 2*pi*fm;
Ts = 1/(2*fm);
n = 0:1/fs:0.01;
m = sin(wc*n).*x;
y = downsample(m,fs/(2*fm));
% Plot results
subplot(2,1,1)
plot(t,x)
title('FM signal')
xlabel('Time (s)')
ylabel('Amplitude')
subplot(2,1,2)
plot(n,y)
title('Bandpass sampled signal')
xlabel('Time (s)')
ylabel('Amplitude')
相关问题
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 10.7 MHz的中频广播信号进行带通采样:
%% 设置采样参数
fs = 22.05e6; % 采样率
fc = 10.7e6; % 中频频率
f_cutoff = 200e3; % 信号带宽
f_pass = 50e3; % 通带截止频率
f_stop = 250e3; % 阻带截止频率
%% 设计滤波器
n = 100; % 滤波器阶数
ripple = 0.1; % 通带波动
attenuation = 60; % 阻带衰减
Wp = f_pass/(fs/2); % 归一化通带截止频率
Ws = f_stop/(fs/2); % 归一化阻带截止频率
[b,a] = ellip(n,ripple,attenuation,Wp,'high'); % 高通Elliptic滤波器设计
%% 读取采样数据
data = load('samples.mat');
x = data.samples;
%% 过滤信号
y = filter(b,a,x);
%% 显示结果
figure;
subplot(211);
plot(x);
title('原始信号');
subplot(212);
plot(y);
title('过滤后信号');
%% 保存结果
save('filtered.mat','y');