wblxx.txt

某大学开发的一个用于对微波传输线状态进行模拟的软件

VERSION 5.00
Private Sub cmdOK_Click()
  Dim pi As Double
  pi = 4 * Atn(1)
  
 If Frame8.Visible = True Then
   Text33.Text = "?"
   Text39.Text = "?"
   kw = Text35.Text
   kh = Text36.Text
   ker = Text37.Text
   kwh = kw / kh
   ka = 0.5 * (ker + 1)
   If kwh > 1 Then
   kb = 0.5 * (ker - 1) / Sqr(1 + 12 / kwh)
   kef = ka + kb
   kc = kwh + 1.393 + 0.667 * Log(kwh * 1.444)
   kzo = 120 * pi / (Sqr(kef) * kc)
   Else
   kb = 0.5 * (ker - 1) * (1 / Sqr(1 + 12 / kwh) + 0.04 * (1 - kwh) ^ 2)
   kef = ka + kb
   kzo = 60 * Log(8 / kwh + 0.25 * kwh) / Sqr(kef)
   End If
   ka1 = (kef + 0.3) / (kef - 0.258)
   ka2 = (kwh + 0.264) / (kwh + 0.8)
   kdl = 0.412 * ka1 * ka2 * kh
   kcoc = 120 * pi * kdl * Sqr(kef) / (1000 * kzo)
   Text33.Text = Format(kdl, "##########.0000000000")
   Text39.Text = Format(kcoc, "##########.0000000000")
 End If
 
 If Frame7.Visible = True Then
  Text25.Text = "?"
  Text55.Text = "?"
  cw = Text54.Text
  ch = Text53.Text
  cg = Text52.Text
  cer = Text26.Text
  cwh = cw / ch
  cgw = cg / cw
  cano = 0.619 * cwh * (Log(cwh) - 0.3853)
  cako = 4.26 - 1.453 * Log(cwh)
  If cgw > 0.3 Then
  cane = 1.565 / Exp(0.16 * Log(cwh)) - 1
  cake = 1.97 - 0.03 / cwh
  Else
  cane = 0.8675
  cake = 2.043 * Exp(0.12 * Log(cwh))
  End If
  ccoe = cw * Exp(cake) * Exp(cane * Log(cgw)) / 1000
  ccoo = cw * Exp(cako) * Exp(cano * Log(cgw)) / 1000
  cp = cer / 9.6
  cce = ccoe * Exp(0.9 * Log(cp))
  cco = ccoo * Exp(0.8 * Log(cp))
  cc1 = 0.5 * cce
  cc2 = 0.5 * (cco - 0.5 * cce)
  Text25.Text = Format(cc1, "##########.0000000000")
  Text55.Text = Format(cc2, "##########.0000000000")
  
 End If
 
 If Frame6.Visible = True Then
   Text42.Text = "?"
   Text41.Text = "?"
   Text43.Text = "?"
   Text44.Text = "?"
   jw1 = Text6.Text
   jw2 = Text7.Text
   jh = Text8.Text
   jer = Text9.Text
   jw21 = jw2 / jw1
   jwh1 = jw1 / jh
   jwh2 = jw2 / jh
   jn = Log(jw21)
   If jer = 9.6 And jw21 > 3.5 Then
   jcs = Sqr(jw1 * jw2) * (56.46 * jn - 44) / 1000
   Else
   ju = Log(jer)
   jcs = Sqr(jw1 * jw2) * ((4.386 * ju + 2.33) * jw21 - 5.472 * ju - 3.17) / 1000
   End If
   ja1 = 0.5 * (jer + 1)
   If jwh1 > 1 Then
   jb1 = 0.5 * (jer - 1) / Sqr(1 + 12 / jwh1)
   jef1 = ja1 + jb1
   jc1 = jwh1 + 1.393 + 0.667 * Log(jwh1 * 1.444)
   jzo1 = 120 * pi / (Sqr(jef1) * jc1)
   Else
   jb1 = 0.5 * (jer - 1) * (1 / Sqr(1 + 12 / jwh1) + 0.04 * (1 - jwh1) ^ 2)
   jef1 = ja1 + jb1
   jzo1 = 60 * Log(8 / jwh1 + 0.25 * jwh1) / Sqr(jef1)
   End If
    ja2 = 0.5 * (jer + 1)
   If jwh2 > 1 Then
   jb2 = 0.5 * (jer - 1) / Sqr(1 + 12 / jwh2)
   jef2 = ja2 + jb2
   jc2 = jwh2 + 1.393 + 0.667 * Log(jwh2 * 1.444)
   jzo2 = 120 * pi / (Sqr(jef2) * jc2)
   Else
   jb2 = 0.5 * (jer - 1) * (1 / Sqr(1 + 12 / jwh2) + 0.04 * (1 - jwh2) ^ 2)
   jef2 = ja2 + jb2
   jzo2 = 60 * Log(8 / jwh2 + 0.25 * jwh2) / Sqr(jef2)
   End If
   jd = jzo1 * Sqr(jef1) + jzo2 * Sqr(jef2)
   jv1 = jzo1 * Sqr(jef1) / jd
   jv2 = jzo2 * Sqr(jef2) / jd
   jalsh = 40.5 * (jw21 - 1) - 32.57 * jn + 0.2 * (jw21 - 1) ^ 2
   jl1 = jv1 * jalsh * jh / 1000
   jl2 = jv2 * jalsh * jh / 1000
   jdl = 0.3 * jalsh * jh / jd
    Text41.Text = Format(jcs, "##########.0000000000")
  Text42.Text = Format(jdl, "##########.0000000000")
    Text43.Text = Format(jl1, "##########.0000000000")
  Text44.Text = Format(jl2, "##########.0000000000")
  End If
  
If Frame5.Visible = True Then
   Text28.Text = "?"
   Text29.Text = "?"
   zw = Text31.Text
   zh = Text32.Text
   zer = Text27.Text
   zwh = zw / zh
   zswh = Sqr(zwh)
   If zwh >= 1 Then
   zcb = zw * ((9.5 * zer + 1.25) * zwh + 5.2 * zer + 7) / 1000
   Else
   zv1 = ((14 * zer + 12.5) * zwh - (1.83 * zer - 2.25)) / zswh
   zv2 = 0.02 * zer / zwh
   zcb = zw * (zv1 + zv2) / 1000
   End If
   zalb = 50 * zh * (4 * zswh - 4.21) / 1000
   Text28.Text = Format(zcb, "##########.0000000000")
  Text29.Text = Format(zalb, "##########.0000000000")
   

End If

If Frame4.Visible = True Then
  Text11.Text = "?"
  Text15.Text = "?"
  Text16.Text = "?"
  tw1 = Text10.Text
  tw2 = Text12.Text
  th = Text13.Text
  ter = Text14.Text
  twh1 = tw1 / th
  twh2 = tw2 / th
  ta1 = 0.5 * (ter + 1)
  If twh1 > 1 Then
   tb1 = 0.5 * (ter - 1) / Sqr(1 + 12 / twh1)
   tef1 = ta1 + tb1
   tc1 = twh1 + 1.393 + 0.667 * Log(twh1 * 1.444)
   tzo1 = 120 * pi / (Sqr(tef1) * tc1)
   Else
   tb1 = 0.5 * (ter - 1) * (1 / Sqr(1 + 12 / twh1) + 0.04 * (1 - twh1) ^ 2)
   tef1 = ta1 + tb1
   tzo1 = 60 * Log(8 / twh1 + 0.25 * twh1) / Sqr(tef1)
   End If
    ta2 = 0.5 * (ter + 1)
  If twh2 > 1 Then
   tb2 = 0.5 * (ter - 1) / Sqr(1 + 12 / twh2)
   tef2 = ta2 + tb2
   tc2 = twh2 + 1.393 + 0.667 * Log(twh2 * 1.444)
   tzo2 = 120 * pi / (Sqr(tef2) * tc2)
   Else
   tb2 = 0.5 * (ter - 1) * (1 / Sqr(1 + 12 / twh2) + 0.04 * (1 - twh2) ^ 2)
   tef2 = ta2 + tb2
   tzo2 = 60 * Log(8 / twh2 + 0.25 * twh2) / Sqr(tef2)
   End If
  tq = (Exp(0.0072 * tzo2) - 1 / Exp(0.0072 * tzo2)) / (Exp(0.0072 * tzo2) + 1 / Exp(0.0072 * tzo2))
  tct = tw1 * (100 / tq + 0.64 * tzo2 - 261) / 1000
  talw1 = tzo1 * Sqr(tef1) / 0.3
  talw2 = tzo2 * Sqr(tef2) / 0.3
  tu1 = 0.016 / twh1
  tu2 = twh2 * (-0.016 * twh1 + 0.064)
  tal1 = -twh2 * th * (tu1 + tu2) * talw1 / 1000
  tv1 = (0.12 * twh1 - 0.47) * twh2
  tv2 = 0.195 * twh1 - 0.357
  tv3 = 0.0283 * Sin((twh1 - 0.75) * pi)
  tal2 = th * (tv1 + tv2 + tv3) * talw2 / 1000
  Text15.Text = Format(tct, "##########.0000000000")
  Text16.Text = Format(tal1, "##########.0000000000")
  Text11.Text = Format(tal2, "##########.0000000000")
End If

If Frame3.Visible = True Then
  Text21.Text = "?"
  Text22.Text = "?"
  Text23.Text = "?"
  Text49.Text = "?"
  sw1 = Text17.Text
  sw2 = Text18.Text
  sh = Text19.Text
  ser = Text20.Text
  swh1 = sw1 / sh
  swh2 = sw2 / sh
  sswh1 = Sqr(swh1)
  sswh2 = Sqr(swh2)
  spwh1 = swh1 * swh1
  spwh2 = swh2 * swh2
  sq1 = Exp(-1 / 3 * Log(swh1))
  sq2 = Exp(-2 / 3 * Log(swh1))
  sq3 = Exp(-2 / 3 * Log(swh2))
  sv1 = (37.61 * swh2 - 13.42 * sswh2 + 159.38) * Log(swh1)
  sv2 = spwh2 * swh2
  sv3 = 74 * swh2 + 130
  sv4 = -60 + 0.5 / swh2 - 0.375 * swh1 * (1 - swh2)
  sc = (0.25 * sw1 * (sv1 + sv2 + sv3) * sq1 + sv4) / 1000
  su1 = (165.6 * swh2 + 31.2 * sswh1 - 11.8 * spwh1) * swh2 - 32 * swh1 + 3
  sal1 = sh * su1 * sq2 / 1000
  su2 = (165.6 * swh1 + 31.2 * sswh1 - 11.8 * spwh1) * swh2 - 32 * swh1 + 3
  sal2 = sh * su2 * sq3 / 1000
  su3 = 337.5 + (1 + 7 / swh1) / swh2
  su4 = 5 * swh2 * Cos(0.5 * pi * (1.5 - swh1))
  sal3 = -sh * (su3 - su4) / 1000
  Text21.Text = Format(sc, "##########.0000000000")
  Text22.Text = Format(sal1, "##########.0000000000")
  Text23.Text = Format(sal2, "##########.0000000000")
  Text49.Text = Format(sal3, "##########.0000000000")
End If

If Frame2.Visible = True Then
  Text45.Text = "?"
  w = Text1.Text
  h = Text2.Text
  a = Text3.Text
  b = Text4.Text
  er = Text5.Text
  wh = w / h
  wbh = (w - b) / h
  a = 0.5 * (er + 1)
   If wh > 1 Then
   b1 = 0.5 * (er - 1) / Sqr(1 + 12 / wh)
   ef1 = a + b1
   c1 = wh + 1.393 + 0.667 * Log(wh * 1.444)
   zo1 = 120 * pi / (Sqr(ef1) * c1)
   Else
   b1 = 0.5 * (er - 1) * (1 / Sqr(1 + 12 / wh) + 0.04 * (1 - wh) ^ 2)
   ef1 = a + b1
   zo1 = 60 * Log(8 / wh + 0.25 * wh) / Sqr(ef1)
   End If
   If wbh > 1 Then
   b2 = 0.5 * (er - 1) / Sqr(1 + 12 / wbh)
   ef2 = a + b2
   c2 = wbh + 1.393 + 0.667 * Log(wbh * 1.444)
   zo2 = 120 * pi / (Sqr(ef2) * c2)
   Else
   b2 = 0.5 * (er - 1) * (1 / Sqr(1 + 12 / wbh) + 0.04 * (1 - wbh) ^ 2)
   ef2 = a + b2
   zo2 = 60 * Log(8 / wbh + 0.25 * wbh) / Sqr(ef2)
   End If
   aln = 2 * h * (1 - zo1 * Sqr(ef1 / ef2) / zo2) ^ 2
   Text45.Text = Format(aln, "##########.0000000000")
End If
End Sub