【正文】 * sxy sx * sy) / ((bn + 1) * sx2 sx ^ 2) a1 = (sx2 * sy sx * sxy) / ((bn + 1) * sx2 sx ^ 2) End If Next b aa1 = 10 ^ a1 bb1 = bmax cc1 = 1 / ( * a2) afw = 100 * (1 10 ^ ((aa + ab * aa1 * cc1 / (bb1 + 1) + Log( * ab) / Log(10)))) 39。極限含水率 fw 公式， fw=110(A+B(a*c/(b+1)(1exp(t^(b+1)/c)))) tm = (bb1 * cc1 / (bb1 + 1)) ^ (1 / (bb1 + 1)) 39。最高年產(chǎn)量對應(yīng)時間 tm 公式， tm=(b*c/(b+1))^(1/(b+1)) qmax = aa1 * (bb1 * cc1 / (bb1 + 1)) ^ (bb1 / (bb1 + 1)) * Exp(bb1 / (bb1 + 1)) 39。最高年產(chǎn)量 Qmax 公式， Qmax=a*(b*c/(b+1))^(b/(b+1))*exp(b/(b+1)) NR = aa1 * cc1 / (bb1 + 1) 39?？刹蓛α?NR 公式， NR=a*c/(b+1) = a2 39。斜率 = a1 39。截距 = rmax 39。最大相關(guān)系數(shù) = aa1 39。參數(shù) a = bb1 39。參數(shù) b = cc1 39。參數(shù) c = afw 39。極限含水率 = qmax 39。最高年產(chǎn)量 = tm 39。最高年產(chǎn)量對應(yīng)時間 = NR 39?？刹蓛α? Rem 預(yù)測數(shù)據(jù) For i = 0 To yt t(i) = i + 1 yq0(i) = Round(aa1 * (t(i) ^ bb1) * (e ^ (t(i) ^ (bb1 + 1) / cc1)), 2) ynp(i) = Round(aa1 * cc1 * (1 (e ^ (t(i) ^ (bb1 + 1) / cc1))) / (bb1 + 1), 2) If i = ub Then fw(i) = Round(100 * qw(i) / (qw(i) + q0(i)), 2) 18 yfw(i) = Round(100 * (1 10 ^ ((aa + ab * ynp(i) + Log( * ab) / Log(10)))), 2) If yfw(i) 0 Then yfw(i) = 0 Next i mn = 50 Rem 表格輸出 (0) = 1000 39。表格行寬 (1) = 1200 (2) = 1200 (3) = 1200 (4) = 1200 (5) = 1200 (6) = 1200 = 201 = 7 = (1) * + * 12 39。 = (1) * + * 12 (0, 0) = 時間 /a (0, 1) = 實際產(chǎn)量 (0, 2) = 預(yù)測產(chǎn)量 (0, 3) = 實際累計產(chǎn)量 (0, 4) = 預(yù)測累計產(chǎn)量 (0, 5) = 實際含水率 % (0, 6) = 預(yù)測含水率 % For i = 0 To mn (i + 1, 0) = Format(t(i), 0.) (i + 1, 1) = Format(q0(i), 0.) (i + 1, 2) = Format(yq0(i), 0.) (i + 1, 3) = Format(np(i), 0.) (i + 1, 4) = Format(ynp(i), 0.) (i + 1, 5) = Format(fw(i), 0.) (i + 1, 6) = Format(yfw(i), 0.) 19 Next i 39。橫縱坐標(biāo) With MSChart1 .ColumnCount = 2 .(1). = 15 .(1). = VtPenStyleSolid .(VtChAxisIdX). = True .(VtChAxisIdY). = True .(VtChAxisIdY). = 0 .(VtChAxisIdX). = 0 .(VtChAxisIdX). = 6 .(VtChAxisIdY). = 6 .(VtChAxisIdX). = 0 .(VtChAxisIdY). = 0 .(VtChAxisIdX). = VtPenStyleDotted .(VtChAxisIdY). = VtPenStyleDotted = False = False = VtChChartType2dXY 39。二維坐標(biāo) .Column = 1 End With With MSChart2 .ColumnCount = 2 .(1). = 15 .(1). = VtPenStyleSolid .(VtChAxisIdX). = True .(VtChAxisIdY). = True .(VtChAxisIdY). = 0 .(VtChAxisIdX). = 0 .(VtChAxisIdX). = 6 .(VtChAxisIdY). = 6 20 .(VtChAxisIdX). = 0 .(VtChAxisIdY). = 0 .(VtChAxisIdX). = VtPenStyleDotted .(VtChAxisIdY). = VtPenStyleDotted = False = False = VtChChartType2dXY .Column = 1 End With With MSChart3 .ColumnCount = 2 .(1). = 15 .(1). = VtPenStyleSolid .(VtChAxisIdX). = True .(VtChAxisIdY). = True .(VtChAxisIdY). = 0 .(VtChAxisIdX). = 0 .(VtChAxisIdX). = 6 .(VtChAxisIdY). = 6 .(VtChAxisIdX). = 0 .(VtChAxisIdY). = 0 .(VtChAxisIdX). = VtPenStyleDotted .(VtChAxisIdY). = VtPenStyleDotted = False = False = VtChChartType2dXY .Column = 1 End With End Sub Private Sub c2_Click() = True 39。輸出數(shù)據(jù)在表 格內(nèi) Rem 曲線繪制 Rem 數(shù)據(jù)輸出 End Sub 21 Private Sub c3_Click() 39。文件輸出 Open + For Output As 2 Print 2, 年份 , 時間 , 預(yù)測產(chǎn)量 , 預(yù)測累計產(chǎn)量 , 預(yù)測含水率 For i = 0 To yt Print 2, ye(0) + i, i + 1, yq0(i), ynp(i), yfw(i) Next i Close 2 Print End Sub Private Sub c4_Click() End End Sub Private Sub d1_Click() 39。理論產(chǎn)油量 Qo 與時間 t 的關(guān)系曲線繪制 = True ReDim px(yt), py(yt), s(yt, 3) 39。st = yt / 500 For i = 1 To mn px(i) = i py(i) = aa1 * (px(i) ^ bb1) * (e ^ (px(i) ^ (bb1 + 1) / cc1)) 39。理論產(chǎn)油量 Qo 公式 39。實際散點 s(i, 0) = px(i) s(i, 1) = q0(i) 39。預(yù)測曲線 s(i, 2) = px(i) s(i, 3) = py(i) Next i 39。散點 Dim Index As Integer Index = 1 39。設(shè)置圖標(biāo)類型為非自動選擇，這樣可以設(shè)置為自己需要的類型 (Index).ShowLine = False 39。去除連線！！ 22 (Index). = False With (Index).(1).Marker .Visible = True 39。散點的點是否可見 .Size = 120 39。散點的大小 .Style = VtMarkerStyleFilledDiamond 39。選擇菱形的 . = False 39。顏色不自動填充 . 400, 15, 400 39。設(shè)置為紫紅色 . 400, 15, 400 39。同樣這個也必須設(shè)置，否則中間有默認的顏色 End With 39。橫縱坐標(biāo)表示 With .(VtChAxisIdX, 0).AxisTitle = t .(VtChAxisIdY, 0).AxisTitle = Qo . = Qo 與 t 關(guān)系曲線 End With = s 39。顯示圖形 End Sub Private Sub d2_Click() 39。累積產(chǎn)油量 Np 與時間 t 的關(guān)系曲線 = True