【电磁波】大气波导电磁波射频追踪仿真【含Matlab源码 3559期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【电磁波】基于matlab大气波导电磁波射频追踪仿真【含Matlab源码 3559期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

%% 大气传播射线跟踪法
%clc;clear all;
function xsitax = Main_ray(M0,hmax,ht,HH1,HH2,slopeElvated1,slopeElvated2,xmaxkm,sitamax)
%% 参数说明
% hmax: 最大显示高度(m)；dh：高度步长(m)；ht：发射天线高度(m)；
% HH1：波导基础层高度(m)；HH2：波导层+基础层高度/蒸发波导高度(m)；
% slopeElvated1：基础层斜率(/km)；slopeElvated2：陷获层斜率(/km)；
% xmaxkm：最大显示距离(km);ae:等效地球半径(m)
% sita0：天线初始仰角范围(°)
%% 参数输入
hmax = round(hmax);%%四舍五入返回
%hmax= 500;
%M0=330;sitamax=4;
dh = 1;
h = 0:dh:hmax;
%ht = 20;
%HH1=0;%基础层高度(m)
%HH2=50;%波导层+基础层高度(m)
%slopeElvated1=118;%基础层斜率/km
%slopeElvated2=-120;%陷获层斜率/km
%xmaxkm = 200;
%% 折射率剖面和折射率梯度
%M = Sub_evaporation(h,HH2,M0);
M = Sub_FourParaRefPro(M0,h,HH1,HH2,slopeElvated1,slopeElvated2);
%M = Sub_StanRefPro(h);
% M = hunhebodao(M0,h,HH1,HH2,HH3,slopeElvated1,slopeElvated2);
for i = 1:1:length(h)-1
g(i) = (M(i+1)-M(i))/(h(i+1)-h(i));
end
g(i+1) = g(i);
ae = 63711000;
%% 计算
sitax = zeros(1,2);
iii = 2;
sitax(1,:) = [0,0.0000001];
sitamax = sitamax180/pi;
for sita0 =-0.6:0.05:0.6
id = find(ht == h); %%返回ht=h时的i值
gi = g(id);
gi1 = g(id-1);
dhi = dh;

sita1 = sita0*pi/180;
x1 = 0; h1 = ht;    u = []; u(1,1) = x1;    u(1,2) = h1;    
for ii = 1:1:hmax
    if gi == 0

%% 梯度为0时
dhi = dh;
sita2 = sita1;
h2 = h1+dhi;
x2 = x1+cot(sita2)dh;
id = find(round(h2) == h);
gi = g(id);
gi1 = g(id-1);
else
%% 梯度不为0时
if sita1>0 && sita1<= (3pi/180)
%% 仰角大于0小于3°时，射线向上传播
r = sita1sita1+2gidhi;
if r>=0
dhi = min(h(id+1)-h1,dh);
sita2 = sqrt(sita1sita1+2gidhi);
h2 = h1+dhi;
x2 = x1+(sita2-sita1)/gi;
else
sita2 = 0;
dhi = (sita2sita2-sita1sita1)/2/gi;
h2 = h1+dhi;
x2 = x1+(sita2-sita1)/gi;
end
id = min(find( h >= h2));
gi = g(id);
gi1 = g(id-1);
elseif sita1<0 && sita1>= -3pi/180
%% 仰角小于0大于-3°时，射线向下传播
r = sita1sita1-2gi1dhi;
if r>=0
if h1 == 0
dhi = 0;
sita2 = -sita1; %% 擦地角
h2 = h1;
x2 = x1;
id = min(find(h >= h2));
gi = g(id);
sitax(iii,:) = [x1/1000,sita2]; %%保存擦地角
iii = iii+1;
else
dhi = min(h1-h(id-1),dh);
sita2 = -sqrt(sita1sita1-2gi1dhi);
h2 = h1-dhi;
x2 = x1+(sita2-sita1)/gi1;
id = min(find(h >= h2));
if id == 1
gi = g(id); gi1 = g(id);
else
gi = g(id); gi1 = g(id-1);
end
end
else
sita2 = 0;
dhi = (sita2sita2-sita1sita1)/2/gi1;
h2 = h1+dhi;
x2 = x1+(sita2-sita1)/gi1;
id = min(find(h >= h2));
gi = g(id);
end
elseif sita1 == 0
%% 仰角为0
if gi>0
%% 射线向上传播
r = sita1sita1+2gidhi;
if r>=0
dhi = min(h(id+1)-h1,dh);
sita2 = sqrt(sita1sita1+2gidhi);
h2 = h1+dhi;
x2 = x1+(sita2-sita1)/gi;
else
sita2 = 0;
dhi = (sita2sita2-sita1sita1)/2/gi;
h2 = h1+dhi;
x2 = x1+(sita2-sita1)/gi;
end
id = min(find(h >= h2));
if id == 1
gi = g(id); gi1 = g(id);
else
gi = g(id); gi1 = g(id-1);
end
elseif gi<0
%% 射线向下传播
r = sita1sita1-2gi1dhi;
if r>=0
if h1 == 0
dhi = 0;
sita2 = -sita1;
h2 = h1;
x2 = x1;
id = min(find(h >= h2));
gi = g(id+1);
sitax(iii,:) = [x1/1000;sita2]; %%保存擦地角
iii = iii+1;
else
dhi = min(h1-h(id-1),dh);
sita2 = -sqrt(sita1sita1-2gi1dhi);
h2 = h1-dhi;
x2 = x1+(sita2-sita1)/gi1;
id = min(find(h >= h2));
if id == 1
gi = g(id); gi1 = g(id);
else
gi = g(id); gi1 = g(id-1);
end
end
else
sita2 = 0;
dhi = (sita2sita2-sita1sita1)/2/gi1;
h2 = h1+dhi;
x2 = x1+(sita2-sita1)/gi1;
id = min(find(h >= h2));
gi = g(id);
end
end
elseif sita1>(3pi/180)
%% 仰角大于3°时
dhi = min(h(id+1)-h1,dh);
h2 = h1+dhi;
sq = sqrt(M(id+1)^{2-(M(1)*cos(sita1))}2);
x2 = x1+(ae/(ae+h2))*(M(1)cos(sita1)dhi)/sq;
sita2 = atan(dhi/(x2-x1));
id = min(find(h >= h2));
gi = g(id);
gi1 = g(id-1);
elseif sita1<(-3pi/180)
%% 仰角小于3°时
if h1 == 0
dhi = 0;
sita2 = -sita1;
h2 = h1;
x2 = x1;
id = min(find(h >= h2));
gi = g(id);
else
dhi = min(h1-h(id-1),dh);
h2 = h1-dhi;
sq = sqrt(M(id+1)^{2-(M(1)*cos(sita1))}2);
x2 = x1+(ae/(ae+h2))(M(1)*cos(sita1)*dhi)/sq;
sita2 = -atan(dhi/(x2-x1));
id = min(find(h >= h2));
if id == 1
gi = g(id); gi1 = g(id);
else
gi = g(id); gi1 = g(id-1);
end
end
end
end

    if h1 >= hmax-1
        break;
    end
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]张鹏，陈国平，王宇. 基于特征线法和特征参数法的管道瞬变流计算[J]. 河海大学学报（自然科学版）, 2014, 42(6): 602-607.

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除