方程：
a1x1+b1x2+c1x3+d1x4=e1
a2x1+b2x2+c2x3+d2x4=e2
a3x1+b3x2+c3x3+d3x4=e3
x1、x2、x3、x4為未知數(shù),取值為0-9000范圍內(nèi)整數(shù),其他全部為常數(shù)
現(xiàn)在要求用vb進行求解,求出最優(yōu)的x1、x2、x3、x4值,即假設(shè)：
a1x1+b1x2+c1x3+d1x4值為s1
a2x1+b2x2+c2x3+d2x4值為s2
a3x1+b3x2+c3x3+d3x4值為s3
(s1-e1)*(s1-e1)+(s2-e2)*(s2-e2)+(s3-e3)*(s3-e3)的值最小
VB最原始的方法代碼為：
VB code
m1 = -1:m2 = -1:m3 = -1:m4 = -1
aaa = -1
For n1 = 0 To 9000
For n2 = 0 To 9000
For n3 = 0 To 9000
For n4 = 0 To 9000
s1 = a1 * n1 + b1 * n2 + c1 * n3 + d1 * n4
s2 = a2 * n1 + b2 * n2 + c2 * n3 + d2 * n4
s3 = a3 * n1 + b3 * n2 + c3 * n3 + d3 * n4
If aaa = -1 Then
m1 = n1:m2 = n2:m3 = n3:m4 = n4
aaa = (s1 - e1) * (s1 - e1) + (s2 - e2) * (s2 - e2) + (s3 - e3) * (s3 - e3)
Else
If aaa < (s1 - e1) * (s1 - e1) + (s2 - e2) * (s2 - e2) + (s3 - e3) * (s3 - e3) Then
m1 = n1:m2 = n2:m3 = n3:m4 = n4
aaa = (s1 - e1) * (s1 - e1) + (s2 - e2) * (s2 - e2) + (s3 - e3) * (s3 - e3)
End If
End If
Next
Next
Next
Next
MsgBox "最優(yōu)解為：" + Trim(m1) + "," + Trim(m2) + "," + Trim(m3) + "," + Trim(m4)