【卫星信号】模拟卫星信号传播研究（Matlab代码实现）

👨‍🎓个人主页：研学社的博客    

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

本文用于模拟信标卫星在地球电离层中传播研究。标准卫星轨道参数和地球磁场模型用于确定传播路径。使用这些输入来计算测量字段实现。信号调制是通过沿轨道参数确定的平移方向扫描场来合成的。包括文档讲解和Matlab代码实现。详细讲解见第4部分。

📚2 运行结果

 

例如，图2.1显示了300公里(默认)电离层截距点的经纬度(见图4.4)。洋红色的五角星标志着车站的位置。空间站的位置几乎在轨道上，这一事实表明该通行证几乎在头顶上，这对科学家来说是可取的。c分析。

图2.1(参见图4.5)显示了叠加在磁偏角图上的轨迹(洋红色)。要执行这段代码，IGRF11 Demo代码生成的?le Bz300必须位于提示符处。该通道穿过地磁赤道，但由于磁场几何形状的变化，会遇到一系列与磁场有关的角度。

这可以在图2.1(见图4.6)中更直接地看到，图2.1显示了传播方向与磁场方向之间夹角的余弦(布里格斯-帕金角)。

 

 

 

 

 

矢量是视速度的缩放投影。在假定轨道变化在凌日过程中不变的情况下，空间结构可以转换为时间序列，由参考网格处的接收天线测量。

图2.2和2.2的上帧分别显示了去除2?模棱两可。图2.2的下一帧显示了在0.75秒的滑动间隔内测量的SI指数。

即使驱动不规则结构在名义上是均匀的，局部闪烁结构也会发生变化，这是第4.5.3节中讨论的大规模结构的结果。在图2.2中，上框中的品红叠加是一个512点居中的箱车平均值，用于捕获大规模结构。可以清楚地看到，大尺度相结构与强度变化相吻合。剩余相位变化是相同的吗?有闪烁，虽然划分不严格。

 

 

 

 

 

 

部分代码：

function [rfn,rfnModelID,Ns] = AtmosphericCorrectionCRPL( x, z, varargin )
% USAGE: [rfn,rfnModelID,Ns] = AtmosphericCorrectionCRPL( x, z )
%        [rfn,rfnModelID,Ns] = AtmosphericCorrectionCRPL( x, z, curve, rfnModel, NSurf, href)
%   curve   ='y' => make apply curved earth correction, else flat earth
%   rfnModel= 0  => Standard model, 1=> Exponential atmosphere
%   NSurf   = refractivity at surface (303 nominal)
%   href    = duct height (<100 m)
%   NOTE:  refractivity=(refractive index-1)*1.e6
%
% Translates model radial refractive index to rectangular
% coordinates (x,z) with origin on spherical earth at source location
% Default is standard CRPL model with 303 refractivity units at surface
% This model is close to 4/3 earth to 1000 km height increasing Ns 
% corresponds to denser atmosphere (more refraction)
% Surface duct is from standard model
%------------------CRPL model--------------------------
Re = 6378166;
if isempty(varargin)
    curve='y';
    rfnModel=0;
    NSurf=303;
    href=0;
    rfnModelID='CPRL';
elseif length(varargin)==4
    curve=varargin{1};
    rfnModel=varargin{2};
    if rfnModel==1
        rfnModelID='CRPL_Exp';
    end
    NSurf=varargin{3};
    href =varargin{4};
else
    error('AtmosphericCorrectionCRPL')
end

if( lower(curve) == 'y' )
   CF  = 2*Re*(sqrt(1-(x/Re)^2)-1);
   eta = z+Re*(sqrt(1+CF*z./(z+Re).^2)-1);
else
   eta = z;
end

if rfnModel==0 | rfnModel==1
    N=natmCRPL(eta,rfnModel,NSurf,href);
    rfn = N*1.e-6+1;
else
    error('Not Available')
end
return

 

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]于超,罗昕,王周.GLONASS卫星导航信号传播模型研究[J].通讯世界,2015(20):34-35.

[2]Charles Rino (2023). Beacon Satellite Propagation Simulation 

[3]许正文. 电离层对卫星信号传播及其性能影响的研究[D].西安电子科技大学,2005.

🌈4 Matlab代码实现