小卫星多普勒频偏MATLAB仿真程序及参考文献

MALAB仿真程序，可以根据小卫星轨道高度，产生对应的多普勒频偏。

clear
clc
close all;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=6378140;%地球赤道半径；
H=524000;%卫星距离地面高度；
R=r+H;%卫星轨道半长轴；
% theta=（-20:20）;%卫星当前位置与观察者位置的地心角角距；
% d=sqrt（R^2+r^2-2*R*r*cos（theta*pi/180））;%卫星到观察者的距离；
% alpha=asin（（R*cos（theta*pi/180）-r）./d）*180/pi;%?——观察者对卫星的仰角；
alpha=0:180;%?——观察者对卫星的仰角；
d=sqrt（R^2-r^2+r^2*sin（alpha*pi/180）.^2）-r*sin（alpha*pi/180）;%卫星到观察者的距离；
theta=acos（（sin（alpha*pi/180）.*d+r）/R）*180/pi;%卫星当前位置与观察者位置的地心角角距；
theta（alpha>90）=-theta（alpha>90）;
GM=398603e9;%地球引力常数；
v=sqrt（GM/R）;%卫星绕地运动的线速度；
vd=v*cos（（alpha+theta）*pi/180）;%卫星绕地运动的线速度在d 方向的分量。
t=theta*pi/180*R/v;%时间
f_send=146000000;%卫星发送载波频率
C=3e8;%真空中光速；
f_bias=f_send*vd/C;%多普勒频偏
f_recv=f_send*（1+vd/C）;%地面接收载波频率
atten=-10*log10（d.^2/H^2）;%功率相对过定点时变化
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure（1）;
plot（alphad/1000）;%卫星到观察者的距离与仰角关系
xlabel（‘观察者对卫星的仰角（度）‘）;
ylabel（‘卫星到观察者的距离（km）‘）;
figure（2）;
plot（alphaabs（vd/1000））;%卫星到观察者的距离与仰角关系
xlabel（‘观察者对卫星的仰角（度）‘）;
ylabel（‘卫星相对观察者的速度（km/s）‘）;
figure（3）;
plot（alphaf_bias）;%多普勒频偏
xlabel（‘观察者对卫星的仰角（度）‘）;
ylabel（‘多普勒频偏（Hz）‘）;
figure（4）;
plot（alphaf_bias）;%地面接收载波频率
xlabel（‘观察者对卫星的仰角（度）‘）;
ylabel（‘地面接收载波频率（MHz）‘）;
figure（5）;
plot（alphaatten）;%功率相对过定点时变化
xlabel（‘观察者对卫星的仰角（度）‘）;
ylabel（‘功率相对过定点时变化（dB）‘）;
figure（6）;
plot（talpha）;%功率相对过定点时变化

ylabel（‘观察者对卫星的仰角（度）‘）;
xlabel（‘90度过境时时间（s）‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1839  2018-03-29 17:36  xiaoweixin\Dopler_sin.asv

     文件       1839  2018-03-29 17:37  xiaoweixin\Dopler_sin.m

     文件       1619  2018-03-29 17:55  xiaoweixin\Dopler_sin2.m

     文件    3698454  2018-03-28 15:18  xiaoweixin\紫丁香2号卫星VHF_UHF星地测控链路设计与实现.pdf

     目录          0  2018-03-29 17:38  xiaoweixin

----------- ---------  ---------- -----  ----

              3703751                    5