2024年5月26日发(作者：左源)

clear all, clear all;

M=8;

N=4*M+2;

fd=926;%an=0的fd

% fd=input('an=0的fd:测试值926：');%an=0时的最大多普勒频移

wd=2*pi*fd;

Ns=50000;

% Ns=50000;

Ts=1e-6;%样本点数和采样频率

t=[0:Ns-1]*Ts;

%f=[-Ns/2:Ns/2-1]/(Ns*Ts*fd);

t_0=0;

%t=t';

%ph=unifrnd (-pi,pi,1,M+1);%-pi~pi均匀分布随机变量

%set beita_n,an,bn,wn

uc=zeros(1,Ns);

us=zeros(1,Ns);

ruu=zeros(1,Ns);

rucuc=zeros(1,Ns);

for n=1:M+1

if(n==M+1)

beita_n(M+1)=pi/4;

an(M+1)=sqrt(2)*cos(beita_n(M+1));

bn(M+1)=sqrt(2)*sin(beita_n(M+1));

wn(M+1)=wd;

else

2024年5月26日发(作者：左源)

clear all, clear all;

M=8;

N=4*M+2;

fd=926;%an=0的fd

% fd=input('an=0的fd:测试值926：');%an=0时的最大多普勒频移

wd=2*pi*fd;

Ns=50000;

% Ns=50000;

Ts=1e-6;%样本点数和采样频率

t=[0:Ns-1]*Ts;

%f=[-Ns/2:Ns/2-1]/(Ns*Ts*fd);

t_0=0;

%t=t';

%ph=unifrnd (-pi,pi,1,M+1);%-pi~pi均匀分布随机变量

%set beita_n,an,bn,wn

uc=zeros(1,Ns);

us=zeros(1,Ns);

ruu=zeros(1,Ns);

rucuc=zeros(1,Ns);

for n=1:M+1

if(n==M+1)

beita_n(M+1)=pi/4;

an(M+1)=sqrt(2)*cos(beita_n(M+1));

bn(M+1)=sqrt(2)*sin(beita_n(M+1));

wn(M+1)=wd;

else