程序.txt

基于VB语言的灰色系统建模和应用预测的研究和实现

Private Sub Command1_Click()
Dim n As Integer
n = 10
ReDim h0(n), X1(n), z(n - 1, 2), B(2, n - 1), y(n - 1, 1), g(n), f(n)
Dim c(2, 4), d(2, 4), e(2, 1), A(2)
h0(1) = Text1.Text
h0(2) = Text2.Text
h0(3) = Text3.Text
h0(4) = Text4.Text
h0(5) = Text5.Text
h0(6) = Text6.Text
h0(7) = Text7.Text
h0(8) = Text8.Text
h0(9) = Text9.Text
h0(10) = Text10.Text
'一次累加得h1 ( i)
t = 0
For i = 1 To n
  t = t + h0(i)
  X1(i) = t
Next '构造矩阵B
For j = 1 To n - 1
  B(1, j) = -(X1(j + 1) + X1(j)) / 2        '如果采用优化模型, 将B (1 , j) 的表达式修改即可
  B(2, j) = 1
  Next '构造逆矩阵
For i = 1 To 2
  For j = 1 To n - 1
z(j, i) = B(i, j)
  Next
Next
For i = 1 To 2
  For k = 1 To 2
c(i, k) = 0
'For j = 1 To n - 1
  c(i, k) = c(i, k) + B(i, j) * z(j, k)
  Text11.Text = c(1, 1)
Next
If i = 1 Then c(i, k + 2) = 1 Else c(i, k + 2) = 0
  Next
Next
For j = 2 To 4
  Text11.Text = c(1, 1)
  Text12.Text = c(2, 1)
  c(2, j) = c(2, j) - c(1, j) / c(1, 1) * c(2, 1)
Next
For i = 1 To 2
  For j = 1 To 4
  d(i, j) = c(i, j) / c(i, i)
  Next
Next
i = 2
k = 1
For j = 4 To 2 Step -1
  d(k, j) = d(k, j) - d(i, j) / d(i, i) * d(k, i)
Next '构造矩阵Y
For i = 1 To n - 1
  y(i, 1) = h0(i + 1)
Next
For i = 1 To 2
  e(i, 1) = 0
For j = 1 To n - 1
  e(i, 1) = e(i, 1) + B(i, j) * y(j, 1)
  Next
Next
k = Val(InputBox("预测长度", 5))
'求a , u
For i = 1 To 2
  A(i) = 0
  For j = 3 To 4
    A(i) = A(i) + d(i, j) * e(j - 2, 1)
  Next
Next
For t = 1 To n - 1
  'g(t) = (X1(1) - A(2) / A(1)) / Exp(A(1) * t) + A(2) / A(1)
  f(t) = X1(t + 1) - g(t)
Next
g(0) = g(n - 1)
End Sub