Matlab写一段3✖️3相位解调代码

以下是一段3x3相位解调的Matlab代码：

% 假设输入信号为x，载波频率为fc，采样频率为fs
% 生成3个正交载波
t = :1/fs:2*pi-1/fs;
c1 = cos(2*pi*fc*t);
s1 = sin(2*pi*fc*t);
c2 = cos(2*pi*fc*t + 2*pi/3);
s2 = sin(2*pi*fc*t + 2*pi/3);
c3 = cos(2*pi*fc*t + 4*pi/3);
s3 = sin(2*pi*fc*t + 4*pi/3);

% 对输入信号进行3个正交解调
y1 = x .* c1;
y2 = x .* c2;
y3 = x .* c3;

% 对解调后的信号进行低通滤波
[b,a] = butter(4,2*fc/fs);
z1 = filter(b,a,y1);
z2 = filter(b,a,y2);
z3 = filter(b,a,y3);

% 计算相位
phase1 = atan2(-z1,s1);
phase2 = atan2(-z2,s2);
phase3 = atan2(-z3,s3);

% 输出相位
disp('相位1：');
disp(phase1);
disp('相位2：');
disp(phase2);
disp('相位3：');
disp(phase3);

希望对你有帮助！