Visual Basic 6.0: Codes of Optimization Laboratory

Visual Basic
6.0
Codes of Optimization
Laboratory
‫االعداد‬
‫علي هحود هاشن‬.‫م‬

‫ العل محع هاش‬: ‫داالعدا‬
Non-linear
equations
2
1- Bi-section Method :
xl = Val(Text1.Text)
xu = Val(Text2.Text)
f1 = (xl ^ 3) - 20
f2 = (xu ^ 3) - 20
xm = (xl + xu) / 2
10 f3 = (xm ^ 3) - 20
If (f3 * f1) > 0 Then
xl = xm
Else
xu = xm
End If
xm1 = (xl + xu) / 2
If Abs(xm1 - xm) < 0.00001 Then
Print "xm="; xm
Else
xm = xm1
GoTo 10
End If
f4 = (xm ^ 3) - 20
Ea = Abs((xm - xm1) / xm)
Print "f4="; f4
Print "Ea="; Ea
3
2- False position Method :
Private Sub Command1_Click()
xl = Val(Text1.Text)
xu = Val(Text2.Text)
f1 = (xl ^ 3) - 20
f2 = (xu ^ 3) - 20
xr = ((xu*f1)-(xl*f2))/(f1-f2)
10 f3 = (xr ^ 3) - 20
If (f3 * f1) > 0 Then
xl = xr
Else
xu = xr
End If
xr1 =((xu*f1)-(xl*f2))/(f1-f2)
If Abs(xr1 - xr) < 0.00001 Then
Print "xr="; xr
Else
xr = xr1
GoTo 10
End If
f4 = (xr ^ 3) - 20
Ea = Abs((xr- xr1) / xr)
Print "f4="; f4
Print "Ea="; Ea
4
3- Newton Raphson Method :
x = Val(Text1.Text)
10 f1 = (x ^ 3) - 20
f2 = 3 * (x ^ 2)
X1 = x - (f1 / f2)
If Abs(X1 - x) < 0.001 Then
Print " x1="; X1
Else
x = X1
GoTo 10
End If
f2 = (X1 ^ 3) - 20
Print "f2="; f2
End sub
5
4- Secant method Method :
x0 = Val(Text1.Text)
10 f1 = (x0 ^ 3) - 20
f2 = (x1 ^ 3) - 20
X2 = x1 - ((f2 * (x1 - x0)) / (f2 - f1))
If Abs(X2 - x1) < 0.000000001 Then
Print " x2="; X2
Else
x1 = X2
GoTo 10
End If
f3 = (X2 ^ 3) - 20
Print "f3="; f3
End Sub
6
Simultanous
linear
equations
7
A-Derict Methods :
1- Naẗve Guass elimination Method :
a11 = Val(Text1.Text)
b1 = Val(Text4.Text)
k = -a21 / a11
a21 = a21 + (k * a11)
a22 = a22 + (k * a12)
a23 = a23 + (k * a13)
b2 = b2 + (k * b1)
k1 = -a31 / a11
a31 = a31 + (k1 * a11)
a32 = a32 + (k1 * a12)
8
a33 = a33 + (k1 * a13)
b3 = b3 + (k1 * b1)
k2 = -a32 / a22
a31 = a31 + (k2 * a21)
a32 = a32 + (k2 * a22)
a33 = a33 + (k2 * a23)
b3 = b3 + (k2 * b2)
x3 = b3 / a33
X2 = (b2 - (a23 * x3)) / a22
X1 = (b1 - (a12 * X2) - (a13 * x3)) / a11
Print "x3="; x3
Print "x2="; X2
Print "x1="; X1
End Sub
9
2- Guass - Jordon Method :
k = -a21 / a11
a21 = a21 + (k * a11)
a22 = a22 + (k * a12)
a23 = a23 + (k * a13)
b2 = b2 + (k * b1)
k1 = -a31 / a11
a31 = a31 + (k1 * a11)
a32 = a32 + (k1 * a12)
a33 = a33 + (k1 * a13)
b3 = b3 + (k1 * b1)
10
k2 = -a32 / a22
a31 = a31 + (k2 * a21)
a32 = a32 + (k2 * a22)
a33 = a33 + (k2 * a23)
b3 = b3 + (k2 * b2)
k3 = -a12 / a22
a11 = a11 + (k3 * a21)
a12 = a12 + (k3 * a22)
a13 = a13 + (k3 * a23)
b1 = b1 + (k3 * b2)
k4 = -a13 / a33
a11 = a11 + (k4 * a31)
a12 = a12 + (k4 * a32)
a13 = a13 + (k4 * a33)
b1 = b1 + (k4 * b3)
k5 = -a23 / a33
a21 = a21 + (k5 * a31)
a22 = a22 + (k5 * a32)
a23 = a23 + (k5 * a33)
b2 = b2 + (k5 * b3)
x3 = b3 / a33
X2 = (b2 - (a23 * x3)) / a22
X1 = (b1 - (a12 * X2) - (a13 * x3)) / a11
Print "x3="; x3
11
Print "x2="; X2
Print "x1="; X1
End Sub
12
3- Grammer’s rule
D = (a11 *((a22 * a33) - (a23 * a32))) -(a12*((a21 * a33) - (a23 * a31))) +(a13 *((a21 * a32) -(a22 * a31)))
a1 = (b1 * ((a22 * a33) - (a23 * a32))) - (a12 * ((b2 * a33) - (a23 * b3))) + (a13 * ((b2 * a32) - (a22 * b3)))
a2 = (a11 * ((b2 * a33) - (a23 * b3))) - (b1 * ((a21 * a33) - (a23 * a31))) + (a13 * ((a21 * b3) - (b2 * a31)))
a3 = (a11 * ((a22 * b3) - (b2 * a32))) - (a12 * ((a21 * b3) - (b2 * a31))) + (b1 * ((a21 * a32) - (a22 * a31)))
X1 = a1 / D
X2 = a2 / D
x3 = a3 / D
Print "x3="; x3
Print "x2="; X2
Print "x1="; X1
13
B- Iterative Methods
1- Jacobi method :
X1 = Val(Text12.Text)
X1 = 1
X2 = 2
x3 = 5
10 x4 = (b1 - (a12 * X2) - (a13 * x3)) / a11
x5 = (b2 - (a21 * X1) - (a23 * x3)) / a22
14
x6 = (b3 - (a31 * X1) - (a32 * X2)) / a33
If (x4 - X1) > 0.00001 And (x5 - X2) > 0.00001 And (x6 - x3) > 0.00001 Then
X1 = x4
X2 = x5
x3 = x6
GoTo 10
Else
Print "x4="; x4
Print "x5="; x5
Print "x6="; x6
End If
End Sub
15
2- Gauss- seidel Method :
X1 = 1
X2 = 2
x3 = 5
10 x4 = (b1 - (a12 * X2) - (a13 * x3)) / a11
x5 = (b2 - (a21 * x4) - (a23 * x3)) / a22
x6 = (b3 - (a31 * X4) - (a32 * X5)) / a33
If (x4 - X1) > 0.00001 And (x5 - X2) > 0.00001 And (x6 - x3) > 0.00001 Then
16
X1 = x4
X2 = x5
x3 = x6
GoTo 10
Else
Print "x4="; x4
Print "x5="; x5
Print "x6="; x6
End If
End Sub
17
3- Sor Method :
w = 1.2
X1 = 1
X2 = 2
x3 = 5
10 x4 = X1 + ((1.2 / 12) * (1 - (a11 * X1) - (a12* X2) + (a13 * x3)))
x5 = X2 + ((1.2 / 28) * (a21 - x4 - (a22 * X2) - (a23* x3)))
x6 = x3 + ((1.2 / 76) * (76 - ( a31* 4) - (a32 * x5) - (33 * x3)))

18
If Abs(x4 - X1) > 0.00001 And Abs(x5 - X2) > 0.00001 And Abs(x6 - x3) Then
X1 = x4
X2 = x5
x3 = x6
GoTo 10
Else
Print "x4="; x4
Print "x5="; x5
Print "x6="; x6
End If
End Sub
19
‫داالعدا ‪ :‬العل محع هاش‬
‫‪Interpolation‬‬
‫‪20‬‬
A-Derict Method :
1- First order ]y = a0 + a1x [
2- Second order ] y = a0 + a1x + a2x2 [
x = Val(Text13.Text)
n = Val(Text14.Text)
If n = 1 Then
k = -a31 / a21
a31 = a31 + (k * a21)
a32 = a32 + (k * a22)
b3 = b3 + (k * b2)
a1 = b3 / a32
a0 = (b2 - (a1 * a22)) / a21
y = a0 + (a1 * x)
Print "a0=", a0
Print "a1=", a1
Print "y=", y
21
Else
k = -a21 / a11
a21 = a21 + (k * a11)
a22 = a22 + (k * a12)
a23 = a23 + (k * a13)
b2 = b2 + (k * b1)
k1 = -a31 / a11
a31 = a31 + (k1 * a11)
a32 = a32 + (k1 * a12)
a33 = a33 + (k1 * a13)
b3 = b3 + (k1 * b1)
k2 = -a32 / a22
a31 = a31 + (k2 * a21)
a32 = a32 + (k2 * a22)
a33 = a33 + (k2 * a23)
b3 = b3 + (k2 * b2)
a2 = b3 / a33
a1 = (b2 - (a23 * a2)) / a22
a0 = (b1 - (a12 * a1) - (a13 * a2)) / a11
y = a0 + (a1 * x) + (a2 * x ^ 2)
Print "a0="; a0
Print "a1="; a1
Print "a2="; a2
Print "y=", y
End If
End Sub
22
B- Newton’s divided difference polynomial :
1- Linear Interpolation ,first order ( y = a0 + a1(x-x0) )
2- Quadradic Interpolation ,second order (y = a0 + a1(x-x0) + a2(x-x0)(x-x1))
3- Cubic Interpolation , third order (y= a0+a1(x-x0) +a2(x-x0)(x-x1)+a3(x-x0)(x-x1)(x-x2))
Y0 = Val(Text5.Text)
n = Val(Text9.Text)
If n = 1 Then
a0 = Y0
a1 = (Y1 - Y0) / (X1 - X0)
y = a0 + (a1 * (x - X0))
Print "a0=", a0
Print "a1=", a1
Print "y=", y
End If
If n = 2 Then
a0 = Y0
a1 = (Y1 - Y0) / (X1 - X0)
k = (Y2 - Y1) / (X2 - X1)
a2 = (k - a1) / (X2 - X0)
y = a0 + (a1 * (x - X0)) + (a2 * (x - X0) * (x - X1))
Print "a0=", a0
Print "a1=", a1
Print "a2=", a2
Print "y=", y
End If
If n = 3 Then
a0 = Y0
a1 = (Y1 - Y0) / (X1 - X0)
23
k = (Y2 - Y1) / (X2 - X1)
a2 = (k - a1) / (X2 - X0)
k1 = (y3 - Y2) / (x3 - X2)
a3 = ((k1 - k) - (k - a1)) / (x3 - X0)
y = a0 + (a1 *(x - X0))+(a2 * (x - X0) * (x - X1))+(a3 * (x - X0)*(x - X1) * (x - X2))
Print "a0=", a0
Print "a1=", a1
Print "a2=", a2
Print "a3=", a3
Print "y=", y
End If
End Sub
24
C- Lagrangian Interpolation
1- Linear Interpolation ,first order ( y =l0*y0 + l1*y1 )
2- Quadradic Interpolation ,second order ( y = l0*y0 + l1*y1 + l2*y2 )
3- Cubic Interpolation , third order ( y = l0*y0 + l1*y1 + l2*y2 + l3*y3 )
n = Val(Text9.Text)
If n = 1 Then
l0 = (x - X1) / (x0 - X1)
l1 = (x - x0) / (X1 - x0)
y = l0 * y0 + l1 * Y1
Print "l0=", l0
Print "l1=", l1
Print "y=", y
End If
If n = 2 Then
l0 = ((x - X1) / (x0 - X1)) * ((x - X2) / (x0 - X2))
l1 = ((x - x0) / (X1 - x0)) * ((x - X2) / (X1 - X2))
l2 = ((x - x0) / (X2 - x0)) * ((x - X1) / (X2 - X1))
y = l0 * y0 + l1 * Y1 + l2 * Y2
Print "10=", 10
Print "l1=", l1
Print "l2=", l2
Print "y=", y
End If
If n = 3 Then
l0 = ((x - X1) / (x0 - X1)) * ((x - X2) / (x0 - X2)) * ((x - X3) / (x0 - X3))
l1 = ((x - x0) / (X1 - x0)) * ((x - X2) / (X1 - X2)) * ((x - X3) / (X1 - X3))
l2 = ((x - x0) / (X2 - x0)) * ((x - X1) / (X2 - X1)) * ((x - X3) / (X2 - X3))
25
l3 = ((x - x0) / (X3 - x0)) * ((x - X1) / (X3 - X1)) * ((x - X2) / (X3 - X2))
y = l0 * y0 + l1 * Y1 + l2 * Y2 + l3 * y3
Print "l0=", l0
Print "l1=", l1
Print "l2=", l2
Print "l3=", l3
Print "y=", y
End If
End Sub
26
‫‪Regression‬‬
‫‪27‬‬
1- Linear regression Deriction :
y3 = Val(Text7.Text)
n = Val(Text9.Text)
X= Val(Text10.Text)
L = (X1 * Y1) + (X2 * Y2) + (x3 * y3) + (x4 * y4)
L1 = X1 + X2 + x3 + x4
L2 = Y1 + Y2 + y3 + y4
L3 = (X1 ^ 2) + (X2 ^ 2) + (x3 ^ 2) + (x4 ^ 2)
a1 = ((n * L) - (L1 * L2)) / ((n * L3) - ((L1) ^ 2))
a0 = (L2 / n) - ((a1 * L1) / n)
Y = a0 + ( a1*X )
Print "a1=", a1
Print "a0=", a0
Print "Y=", Y
28
2- Polynomial Model :
n = Val(Text9.Text)
l = Y1 + Y2 + y3 + y4
l1 = (X1 * Y1) + (X2 * Y2) + (x3 * y3) + (x4 * y4)
l2 = ((X1 ^ 2) * Y1) + ((X2 ^ 2) * Y2) + ((x3 ^ 2) * y3) + ((x4 ^ 2) * y4)
d = X1 + X2 + x3 + x4
d1 = (X1 ^ 2) + (X2 ^ 2) + (x3 ^ 2) + (x4 ^ 2)
d2 = (X1 ^ 3) + (X2 ^ 3) + (x3 ^ 3) + (x4 ^ 3)
d3 = (X1 ^ 4) + (X2 ^ 4) + (x3 ^ 4) + (x4 ^ 4)
k = -d / n
d = d + (k * n)
d1 = d1 + (k * d)
d2 = d2 + (k * d1)
l1 = l1 + (k * l)
k1 = -d1 / n
d1 = d1 + (k1 * n)
d2 = d2 + (k1 * d)
d3 = d3 + (k1 * d1)
l2 = l2 + (k1 * l)
k2 = -d2 / d1
d1 = d1 + (k2 * d)
d2 = d2 + (k2 * d1)
d3 = d3 + (k2 * d2)
l2 = l2 + (k2 * l1)
a2 = d3 / l2
a1 = (l1 - (a2 * d2)) / d1
a0 = (l - (a1 * d) - (a2 * d1)) / n
Print " a0 =", a0 ; " a1 =", a1 ; " a2 =", a2
End Sub
29
‫أ رتقد‬
‫فهاتفها‬
‫ُزاحش‬
‫‪30‬‬

Visual Basic 6.0: Codes of Optimization Laboratory

Uploaded by

Copyright:

Available Formats

Visual Basic 6.0: Codes of Optimization Laboratory

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Visual Basic 6.0: Codes of Optimization Laboratory

Uploaded by

Copyright:

Available Formats

Visual Basic

If (f3 * f1) > 0 Then

xm1 = (xl + xu) / 2

If Abs(xm1 - xm) < 0.00001 Then

Ea = Abs((xm - xm1) / xm)

If (f3 * f1) > 0 Then

If Abs(xr1 - xr) < 0.00001 Then

Ea = Abs((xr- xr1) / xr)

If Abs(X1 - x) < 0.001 Then

Print " x1="; X1

X2 = x1 - ((f2 * (x1 - x0)) / (f2 - f1))

If Abs(X2 - x1) < 0.000000001 Then

Print " x2="; X2

a21 = a21 + (k * a11)

a22 = a22 + (k * a12)

a23 = a23 + (k * a13)

a31 = a31 + (k1 * a11)

a32 = a32 + (k1 * a12)

a31 = a31 + (k2 * a21)

a32 = a32 + (k2 * a22)

a33 = a33 + (k2 * a23)

X2 = (b2 - (a23 * x3)) / a22

X1 = (b1 - (a12 * X2) - (a13 * x3)) / a11

a21 = a21 + (k * a11)

a22 = a22 + (k * a12)

a23 = a23 + (k * a13)

a31 = a31 + (k1 * a11)

a32 = a32 + (k1 * a12)

a33 = a33 + (k1 * a13)

a31 = a31 + (k2 * a21)

a32 = a32 + (k2 * a22)

a33 = a33 + (k2 * a23)

a11 = a11 + (k3 * a21)

a12 = a12 + (k3 * a22)

a13 = a13 + (k3 * a23)

a11 = a11 + (k4 * a31)

a12 = a12 + (k4 * a32)

a13 = a13 + (k4 * a33)

a21 = a21 + (k5 * a31)

a22 = a22 + (k5 * a32)

a23 = a23 + (k5 * a33)

X2 = (b2 - (a23 * x3)) / a22

X1 = (b1 - (a12 * X2) - (a13 * x3)) / a11

10 x4 = (b1 - (a12 * X2) - (a13 * x3)) / a11

x5 = (b2 - (a21 * X1) - (a23 * x3)) / a22

10 x4 = (b1 - (a12 * X2) - (a13 * x3)) / a11

x5 = (b2 - (a21 * x4) - (a23 * x3)) / a22

x6 = (b3 - (a31 * X4) - (a32 * X5)) / a33

10 x4 = X1 + ((1.2 / 12) * (1 - (a11 * X1) - (a12* X2) + (a13 * x3)))

x5 = X2 + ((1.2 / 28) * (a21 - x4 - (a22 * X2) - (a23* x3)))

x6 = x3 + ((1.2 / 76) * (76 - ( a31* 4) - (a32 * x5) - (33 * x3)))

You might also like

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.