西农信号与系统——频谱搬移实验

% 频谱搬移实验
%% 参数设置
Fs = 1000; % 采样频率，单位为Hz
t = 0:1/Fs:0.3-1/Fs; % 时间向量

%1. 生成频率为10 Hz的信号X(t)
f_X = 10; % X(t)的频率
X = sin(2*pi*f_X*t); % X的正弦波形

%2. 生成频率为200 Hz的噪声信号C(t)
f_C = 200; % C(t)的频率
C = sin(2*pi*f_C*t); % C的正弦波形

%3. 将X(t)和C(t)组合形成Y(t)
Y = X .* C;

%4. 对Y(t)执行FFT
Y_fft = fft(Y); %计算FFT
n = length(Y); %FFT点数
f = (0:n-1)*(Fs/n); %频率范围
magnitude = abs(Y_fft); %FFT的幅度

%% 绘制原始信号和FFT的图形
figure(1);
subplot(2,2,1);
plot(t,X);
title('低频信号X(t)');
xlabel('时间 (秒)');
ylabel('幅度');

subplot(2,2,2);
plot(t,C);
title('高频信号C(t)');
xlabel('时间 (秒)');
ylabel('幅度');


subplot(2,2,3);
plot(t,Y);
title('已调信号Y(t)');
xlabel('时间 (秒)');
ylabel('幅度');

subplot(2,2,4);
plot(f,magnitude);
title('Y(t)的FFT幅度');
xlabel('频率 (Hz)');
ylabel('幅度');
axis([0 500 0 max(magnitude)]); % 限制频率轴以便更好地观察

%% 接收端
figure(2)
title('解调信号');
xlabel('时间 (秒)');
ylabel('幅度');
Re=Y .* C;
plot(t,Re);