迭代方程
Private Sub Command1_Click()
Dim TB1(76), TB(76), TA1(76), TA(76) As Single
Dim pa, ca, lma, pb, cb, lmb As Single
pa = 2680#: ca = 1340#: lma = 1.93: pb = 0.00152: cb = 950: lmb = 0.00165
Dim tl, ts, L, Ta0, Tb0, Tf As Single
Ta0 = 1500: Tb0 = 20
Dim dx, dt, ma, mb, hb As Single
Dim b, q, p, a, st, r, dr, ri, ab, aa As Single
Dim fa, fb, f0 As Single
Dim i As Integer
Dim time As Single
Dim f1, f2 As Single
dr = 1: dt = 0.5: r = 75
time = Val(Text1.Text)
f1 = pa * lma * dt / (ca * dt * (dr * dr))
f2 = pb * lmb * dt / (cb * dt * (dr * dr))
st = p * dt / 60
q = 0#
For i = 1 To 60
TB(i) = Tb0
Next i
For i = 61 To 75
TA(i) = Ta0
Next i
Print TA(i); TB(i)
Print CStr(time) & "分" & "time' tempreture:"
For p = 1 To 10 Step dt
ri = dr * p
For i = 2 To 59
TB(1) = 20
ab = pb * lmb / (cb * dt)
TB1(i) = (f2 / 2) * ((2 + dr / ri) * TB(i - 1) + (2 / f2 - 4) * TB(i) + (2 - dr / ri) * TB(i + 1))
Next i
Print Tab(20); Format$(i, TB(i))
If i = 60 Then
ab = pb * lmb / (cb * dt)
TB1(60) = ab * dt * ((TB(59) - TA(61)) / (2 * ri * dr) + (TB(59) - 2 * TB(60) + TA(61)) / (dr * dr)) + TB(60)
Print TB(60)
End If
If i = 61 Then
TA(61) = 780
aa = pa * lma / (ca * dt)
TA1(61) = aa * dt * ((TB(60) - TA(62)) / (2 * ri * dr) + (TB(60) - 2 * TA(61) + TA(61)) / (dr * dr)) + TA(61)
Print Tab(20); i = 61; Format$(TA(61))
End If
For i = 62 To 74
TA(61) = 1480
aa = pa * lma / (ca * dt)
TA1(i) = f1 / 2 * ((2 + dr / ri) * TA(i - 1) + (2 / f1 - 4) * TA(i) + (2 - dr / ri) * TA(i + 1))
Print TA(i)
Next i
Print Tab(20); Format$(TA(i))
TA1(75) = 4 * (TA(74) + (1 / f1 - 4) * TA(75))
Print Format$(TA(75), "###.####")
Next p
Print TA(i)
Private Sub Command1_Click()
Dim TB1(76), TB(76), TA1(76), TA(76) As Single
Dim pa, ca, lma, pb, cb, lmb As Single
pa = 2680#: ca = 1340#: lma = 1.93: pb = 0.00152: cb = 950: lmb = 0.00165
Dim tl, ts, L, Ta0, Tb0, Tf As Single
Ta0 = 1500: Tb0 = 20
Dim dx, dt, ma, mb, hb As Single
Dim b, q, p, a, st, r, dr, ri, ab, aa As Single
Dim fa, fb, f0 As Single
Dim i As Integer
Dim time As Single
Dim f1, f2 As Single
dr = 1: dt = 0.5: r = 75
time = Val(Text1.Text)
f1 = pa * lma * dt / (ca * dt * (dr * dr))
f2 = pb * lmb * dt / (cb * dt * (dr * dr))
st = p * dt / 60
q = 0#
For i = 1 To 60
TB(i) = Tb0
Next i
For i = 61 To 75
TA(i) = Ta0
Next i
Print TA(i); TB(i)
Print CStr(time) & "分" & "time' tempreture:"
For p = 1 To 10 Step dt
ri = dr * p
For i = 2 To 59
TB(1) = 20
ab = pb * lmb / (cb * dt)
TB1(i) = (f2 / 2) * ((2 + dr / ri) * TB(i - 1) + (2 / f2 - 4) * TB(i) + (2 - dr / ri) * TB(i + 1))
Next i
Print Tab(20); Format$(i, TB(i))
If i = 60 Then
ab = pb * lmb / (cb * dt)
TB1(60) = ab * dt * ((TB(59) - TA(61)) / (2 * ri * dr) + (TB(59) - 2 * TB(60) + TA(61)) / (dr * dr)) + TB(60)
Print TB(60)
End If
If i = 61 Then
TA(61) = 780
aa = pa * lma / (ca * dt)
TA1(61) = aa * dt * ((TB(60) - TA(62)) / (2 * ri * dr) + (TB(60) - 2 * TA(61) + TA(61)) / (dr * dr)) + TA(61)
Print Tab(20); i = 61; Format$(TA(61))
End If
For i = 62 To 74
TA(61) = 1480
aa = pa * lma / (ca * dt)
TA1(i) = f1 / 2 * ((2 + dr / ri) * TA(i - 1) + (2 / f1 - 4) * TA(i) + (2 - dr / ri) * TA(i + 1))
Print TA(i)
Next i
Print Tab(20); Format$(TA(i))
TA1(75) = 4 * (TA(74) + (1 / f1 - 4) * TA(75))
Print Format$(TA(75), "###.####")
Next p
Print TA(i)
解决方案 »
免费领取超大流量手机卡,每月29元包185G流量+100分钟通话, 中国电信官方发货