四步相移 matlab代码

四步相移是一种用于频率调制和解调的数字信号处理技术。在Matlab中，可以用以下代码实现四步相移：

% 生成信号
fs = 1000; % 采样频率
t = 0:1/fs:1-1/fs; % 时间向量
f1 = 10; % 信号频率
f2 = 20; % 载波频率
A = 1; % 信号幅度
x = A*sin(2*pi*f1*t); % 生成信号

% 相移
phi = pi/2; % 相位偏移
y = A*sin(2*pi*f2*t + phi); % 相移后的信号

% 解调
z1 = A*sin(2*pi*f2*t); % 解调信号1
z2 = A*cos(2*pi*f2*t); % 解调信号2
phi_hat = atan2(mean(z1),mean(z2)); % 估计相位偏移

% 显示结果
figure;
subplot(3,1,1);
plot(t, x);
title('原始信号');
subplot(3,1,2);
plot(t, y);
title('相移后的信号');
subplot(3,1,3);
plot(t, z1, 'r', t, z2, 'b');
title('解调结果');
legend('sin(2*pi*f2*t)', 'cos(2*pi*f2*t)');
xlabel('相位偏移估计值：');
text(0.5, 0.5, ['\hat{\phi} = ', num2str(phi_hat)]);

以上代码首先生成了一个正弦信号x，并在频率为f2的载波上对该信号进行了相移。然后使用解调技术得到了解调信号z1和z2，并计算出了估计的相位偏移phi_hat。最后，使用Matlab的绘图函数将原始信号、相移后的信号和解调结果显示在了三个子图中，以便直观地观察结果。