ray_single.m

给出比较少见的16PSK的仿真算法

function r=ray_single
clear;
clc;
x=0;y=0;
t=0.001;
fc=10.^6; 						%Carrier frequency
wc=2*pi*fc;								
v=10000/3600;			        	 	%Receiver speed[km/h]
c=300*10.^6;						%Light speed
wm=wc*(v/c);						%Maximum shift
fm=wm/(2*pi);						%Doppler shift
c=sqrt(2);
alpha=(2*pi-pi+(2*pi*rand-pi))/4;
ph1=2*pi*rand-pi;
ph2=2*pi*rand-pi;
x=x + c*cos(wm*t*cos(alpha) + ph1);
y=y + c*cos(wm*t*sin(alpha) + ph2);
r=sqrt(x.^2 + y.^2)/sqrt(2);