get_phase.m

来自「该程序在matlab环境下运行可完成QPSK的调制解调过程」· M 代码 · 共 37 行

clear all;
clc;
N=1024;                      %所用的点数
Fc=4800;                     %载波频率（Hz）
fs=60000;                    %抽样频率（Hz）

fid=fopen('psk4.dat','r');   %打开gen_4psk.mdl中由gen_signal保存的数据文件
fseek(fid, 2000, 'bof');
data=fread(fid,N,'double');
fclose(fid);
data=data'; 

SA = hilbert(data);          %SA 复解析信号
t=0:1/fs:(N-1)/fs;
SB=SA.*exp(-j*2*pi*Fc*t/fs); %SB 下变频后的复基带信号
SI=real(SB);                 %同相分量
SQ=imag(SB);                 %正交分量
mag=sqrt(SI.*SI+SQ.*SQ);
x=SI./mag;                   %幅度归一化
y=SQ./mag;
figure(1)
plot(SI,SQ,'.');
grid
title('星座图')
figure(2)
plot(SI,SQ);
grid
title('轨迹图')
figure(3)
plot(x,y,'.');
title('归一化星座图')
grid
figure(4)
plot(x,y);
grid
title('归一化轨迹图')