voacap计算落地场强和电路可靠度

用ITS HFBC软件的voacapw.exe程序实现。

软件下载http://www.greg-hand.com/versions/,下载 itshfbc_180417a.exe 

下载后程序路径:C:\itshfbc\bin_win

voacapw命令行执行的格式要求https://www.voacap.com/voacapw.html

程序调用方法:c:\itshfbc\bin_win\voacapw.exe silent c:\itshfbc voacapx.dat voacapx.out

其中voacapx.dat放在路径c:\itshfbc\run下,voacapx.out生成后自动存到c:\itshfbc\run下。

输入文件voacapx.dat的内容设置方法

先在voacap gui设置,然后用run->Circuit计算后从结果中拷贝到voacapx.dat中,格式样例如下:

COMMENT    Any VOACAP default cards may be placed in the file: VOACAP.DEF

LINEMAX      55       number of lines-per-page

COEFFS    CCIR

TIME          1    1    1    1

MONTH      201910.00

SUNSPOT     50.

LABEL     #  26 /  6NE.4      #  64 /  8NE.1

CIRCUIT   41.00N   113.00W    41.00N    74.00W  S     0

SYSTEM       1. 160. 0.10  70. 73.0 3.00 0.10

FPROB      1.00 1.00 1.00 0.00

ANTENNA       1    1    2   30     0.000[DEFAULT\CONST17.VOA  ]  0.0  500.0000

ANTENNA       2    2    2   30     0.000[DEFAULT\SWWHIP.VOA   ]  0.0    0.0000

FREQUENCY  6.07 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00

METHOD       20    0

EXECUTE

QUIT

 

输出结果文件voacapx.out的格式样例

   1.0 18.9  6.1  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 FREQ

        1F2  2F2   -    -    -    -    -    -    -    -    -    -  MODE

        5.4 12.7   -    -    -    -    -    -    -    -    -    -  TANGLE

       11.4 11.5   -    -    -    -    -    -    -    -    -    -  DELAY

        378  243   -    -    -    -    -    -    -    -    -    -  V HITE

       0.50 1.00   -    -    -    -    -    -    -    -    -    -  MUFday

        137  124   -    -    -    -    -    -    -    -    -    -  LOSS

         56   57   -    -    -    -    -    -    -    -    -    -  DBU

        -80  -67   -    -    -    -    -    -    -    -    -    -  S DBW

       -183 -160   -    -    -    -    -    -    -    -    -    -  N DBW

        102   92   -    -    -    -    -    -    -    -    -    -  SNR

        -19  -13   -    -    -    -    -    -    -    -    -    -  RPWRG

       0.92 0.94   -    -    -    -    -    -    -    -    -    -  REL

       0.00 0.00   -    -    -    -    -    -    -    -    -    -  MPROB

       0.88 0.85   -    -    -    -    -    -    -    -    -    -  S PRB

       25.0 13.2   -    -    -    -    -    -    -    -    -    -  SIG LW

       15.1  6.1   -    -    -    -    -    -    -    -    -    -  SIG UP

       26.2 15.4   -    -    -    -    -    -    -    -    -    -  SNR LW

       15.5  9.7   -    -    -    -    -    -    -    -    -    -  SNR UP

       17.0 17.0   -    -    -    -    -    -    -    -    -    -  TGAIN

       -3.7 -0.7   -    -    -    -    -    -    -    -    -    -  RGAIN

         92   86   -    -    -    -    -    -    -    -    -    -  SNRxx

其中1.0是我们设置的时间,实际时间为UT 00:30-01:30,1.0为居中时间。

6.1为当前设置的频率,18.9为最大可用频率MUF,最佳工作频率可取OWF=MUF * 0.85

DBU为监测台理论上收到的落地场强(field strength),REL为电路可靠度(circuit reliability)

我们去DBU和REL的值计算一个监测点是否合格。合格的标准是落地场强大于55dB,电路可靠度大于70%。通过计算一个服务区内的所有监测点合格的数量来计算理论覆盖率。

批量计算的方法

参考https://www.voacap.com/voacapw.html 中“Running a set of different CIRCUITS while keeping other parameters fixed“小节改变部分参数重新执行execute

定制发射台/监测站经纬度列表

将:\itshfbc\userdb\目录下Transmit.def和Receive.def的内容分别修改为关注的所有发射台和关注的所有监测点列表即可,这样在设置发射台参数的时候就可看到自定义的台列表。

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
电偶极子场强计算可以通过以下步骤实现: 1. 定义电偶极子的位置和电荷。 假设电偶极子的位置为(x,y,z),电荷为q。这些参数需要作为输入提供给程序。 2. 计算电偶极子对某一点的电场强。 定义要计算场强的点为(x0,y0,z0)。电场强计算公式为: E = (1/(4*pi*epsilon))*((q/r^3)*(2*cos(theta)*r_hat + sin(theta)*theta_hat)) 其中,r_hat和theta_hat表示径向和切向的单位向量,r和theta分别为极坐标系下偶极子和点之间的距离和夹角。 3. 在三维空间中绘制电场强矢量场。 使用matlab中的quiver函数可以在三维空间中绘制电场强矢量场,其中每个矢量的起始点为计算场强的点,终点为该点的电场强矢量。 完整的matlab代码可以如下所示: ```matlab % 电偶极子场强计算 % 定义电偶极子位置和电荷 x = 0; y = 0; z = 0; q = 1; % 定义要计算场强的点 x0 = linspace(-1,1,20); y0 = linspace(-1,1,20); z0 = linspace(-1,1,20); % 计算场强 epsilon = 8.854187817e-12; % 真空介质常数 E = zeros(length(x0), length(y0), length(z0), 3); for i = 1:length(x0) for j = 1:length(y0) for k = 1:length(z0) r = sqrt((x-x0(i))^2 + (y-y0(j))^2 + (z-z0(k))^2); theta = atan2(sqrt((x-x0(i))^2 + (y-y0(j))^2), (z-z0(k))); r_hat = [(x-x0(i))/r, (y-y0(j))/r, (z-z0(k))/r]; theta_hat = [-sin(theta)*cos(phi), -sin(theta)*sin(phi), cos(theta)]; E(i,j,k,:) = (1/(4*pi*epsilon))*((q/r^3)*(2*cos(theta)*r_hat + sin(theta)*theta_hat)); end end end % 绘制电场强矢量场 [X0,Y0,Z0] = meshgrid(x0,y0,z0); quiver3(X0,Y0,Z0,E(:,:,:,1),E(:,:,:,2),E(:,:,:,3)); ``` 注意,以上代码仅为示例,需要根据具体情况进行调整。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值