遗传算法中初始化种群一些重要的源码

admin

1. 不可行初始种群产生的程序源代码:
   
2. 
   
3. 
   
4. Private Sub Command1_Click()          '切换到不可行性初始种群产生的窗口
   
5. Form2.Show
   
6. Form1.Hide
   
7. End Sub
   
8. Private Sub Command2_Click()          '切换到不可行性初始种群产生的窗口
   
9. Form3.Show
   
10. Form1.Hide
    
11. End Sub
    
12. 
    
13. Private Sub Form_Load()                '设置窗口背景颜色
    
14. BackColor = QBColor(3)
    
15. End Sub
    
16. 
    
17. Private Sub Form_Unload(Cancel As Integer) '提示窗口
    
18. Dim v As Integer
    
19. v = MsgBox("你真的要关闭窗口？", vbYesNo)
    
20. If v = vbNo Then
    
21. Cancel = -1
    
22. End If
    
23. End Sub
    
24. 
    
25. Option Explicit                         '定义变量
    
26. Private Declare Function timeGetTime() Lib "winmm.dll" () As Long   '声明时间函数
    
27. Dim k(1 To 100), Y(10), g(1 To 11) As Integer
    
28. Dim a, C, X, p, f, t, h, i, j, e, u, ii As Integer
    
29. Dim sh, m As Double
    
30. Dim n As String, S(1000), r(1000001) As String
    
31. Dim l(100, 10) As Single
    
32. Dim b, bb, D, z As Long
    
33. '参数注释：a，c ：参数；b，bb，D：时间；i，j：循环变量；n：染色体字符串；u：字符串长度；sh：熵值；m：概率；z：种群大小；S(1000)：存放生成的随机数；r(1000001)存放初始种群；k(1 to 100)：某一基因位上各基因的个数；Y(t)：各基因的个数
    
34. Private Sub Command1_Click()             '返回主窗口
    
35. Form1.Show
    
36. Form2.Hide
    
37. End Sub
    
38. 
    
39. Private Sub Command3_Click()             '初始化
    
40. Combo1.ListIndex = -1
    
41. Text2.Text = ""
    
42. Text7.Text = ""
    
43. Text8.Text = ""
    
44. For j = 1 To 100
    
45. k(j) = 0
    
46. t = 0
    
47. m = 0
    
48. sh = 0
    
49. n = ""
    
50. Next
    
51. For i = 1 To 100
    
52. k(i) = 0
    
53. Next
    
54. For j = 1 To 100
    
55. For i = 1 To 10
    
56.   l(j, i) = 0
    
57. Next
    
58. Next
    
59. End Sub                               '利用For循环清空各个值及数组
    
60. Private Sub Command2_Click()            '运行
    
61. b = timeGetTime()                      '调取程序开始时间
    
62. j = Combo1.ListIndex
    
63. If j = 1 Then                           '选择随机函数Rnd法
    
64. z = Text3.Text                          '给种群大小赋值
    
65. p = Val(Text9)                          '给约束赋值
    
66. For j = 1 To z
    
67.   For i = 1 To 1000
    
68.   bb = timeGetTime()()
    
69.   Randomize (bb + i)     '利用时间加循环次数作为控制因子对Rnd函数进行初始化
    
70.   X = Int(p * Rnd)
    
71.   t = X Mod 10
    
72.   If Y(t) < p Then        '当Y(t)大于各基因的约束时跳出If语句
    
73.   Y(t) = Y(t) + 1
    
74.   n = n & t
    
75.   End If
    
76.   u = Len(n)
    
77.   If u = p * p Then       '当染色体的长度达到100时跳出For循环，j=j+1
    
78.   Exit For
    
79.   End If
    
80.   Next i
    
81. n = Replace(n, "0", "a")   '将染色体中的所有字符"0"置换为"a"
    
82. r(j) = n
    
83. i = 1
    
84. n = ""
    
85. For t = 0 To 9           '将Y(t)进行清零
    
86. Y(t) = 0
    
87. Next
    
88. Next j
    
89. Open "F:\vb1\随机函数法.txt" For Output As #1   '输出初始种群
    
90. For i = 1 To z Step 1
    
91. Print #1, r(i)
    
92. Next i
    
93. Text7.Text = r(1)                             '将第一条染色体输出到窗口中
    
94. Close #1
    
95. D = timeGetTime()                           '调取程序结束时间
    
96. Text8.Text = (D - b) / 1000 & "秒"               '统计程序运行时间
    
97. shang                                           '调用子函数
    
98. Open "F:\vb1\重复.txt" For Output As #3         '输出初始种群中的重复染色体
    
99. For i = 1 To z - 1
    
100. For j = i + 1 To z
     
101. If r(i) = r(j) Then                           '如果两条染色体相同则输出
     
102. Print #3, r(i)
     
103. End If
     
104. Next
     
105. Next
     
106. Close #3
     
107. Else                                        '选择混合同余法
     
108. p = Val(Text9)                                '给初始种群赋值
     
109. a = Text4.Text                                '给参数赋值
     
110. a = 4 * a + 1
     
111. C = Text5.Text
     
112. C = 2 * C + 1
     
113. z = Text3.Text
     
114. For j = 1 To z
     
115. bb = timeGetTime()
     
116. Randomize (bb + i)                   '利用时间加循环次数作为控制因子对Rnd函数进行初始化
     
117. X = Int(Rnd * (p * 10))                       '产生0到99的随机数
     
118. S(1) = X
     
119. For i = 1 To 127
     
120. S(i + 1) = (a * S(i) + C) Mod 128             '对随机数进行运算
     
121. If S(i + 1) < 100 Then
     
122. u = S(i + 1) Mod 10
     
123. n = n & u
     
124. End If
     
125. Next i
     
126. n = n & (X Mod 10)
     
127. n = Replace(n, "0", "a")                        '将染色体中的字符"0"置换为"a"
     
128. r(j) = n
     
129. i = 1
     
130. n = ""
     
131. Next j
     
132. Open "F:\vb1\混合同余法.txt" For Output As #2  '输出不可行初始种群
     
133. For i = 1 To z Step 1
     
134.   Print #2, r(i)
     
135. Next i
     
136. Text7.Text = r(1)                              '将第一条染色体输出在窗口上
     
137. Close #2
     
138. D = timeGetTime()                            '调取程序结束时间
     
139. Text8.Text = (D - b) / 1000 & "秒"                '统计程序运行时间
     
140. shang                                       '调用子函数
     
141. Open "F:\vb1\重复.txt" For Output As #4       '输出重复的染色体
     
142. For i = 1 To z - 1
     
143. For j = i + 1 To z
     
144.   If r(i) = r(j) Then                            '如果两条染色体相同则输出
     
145.   Print #4, r(j)
     
146.   End If
     
147. Next
     
148. Next
     
149. Close #4
     
150. End If
     
151. End Sub
     
152. Private Function shang()                        '子函数计算熵值
     
153. t = 0
     
154. For j = 1 To 100
     
155.   For i = 1 To z
     
156.    Select Case Mid(r(i), j, 1)            '统计r(i)中第i基因位上的各基因个数
     
157.    Case "1"
     
158.    g(1) = g(1) + 1                             '若是字符"1",则g(1)加1，统计字符的个数
     
159.    Case "2"
     
160.    g(2) = g(2) + 1
     
161.    Case "3"
     
162.    g(3) = g(3) + 1
     
163.    Case "4"
     
164.    g(4) = g(4) + 1
     
165.    Case "5"
     
166.    g(5) = g(5) + 1
     
167.    Case "6"
     
168.    g(6) = g(6) + 1
     
169.    Case "7"
     
170.    g(7) = g(7) + 1
     
171.    Case "8"
     
172.    g(8) = g(8) + 1
     
173.    Case "9"
     
174.    g(9) = g(9) + 1
     
175.    Case "a"
     
176.    g(10) = g(10) + 1
     
177.    Case Else
     
178.    g(11) = 0
     
179.    End Select
     
180.    Next
     
181.    For i = 1 To 10                          '不同基因的个数统计
     
182.     If g(i) <> 0 Then
     
183.     k(j) = k(j) + 1
     
184.     End If
     
185.    Next
     
186.     For i = 1 To 10                        '不同基因位上不同基因的概率统计
     
187.      l(j, i) = g(i) / z
     
188.     Next
     
189.     For i = 1 To 10                        '将g(10)清零
     
190.      g(i) = 0
     
191.     Next
     
192.   If k(j) > t Then                          '统计各基因位上基因的最大值
     
193.    t = k(j)
     
194.   End If
     
195. Next
     
196. For j = 1 To 100                           '将二维数组中的各个值赋给m
     
197. For i = 1 To k(j)
     
198.   m = l(j, i)
     
199.   If  m <> 0 Then
     
200.   sh = 1 / (Log(t) * 100) * Log(m) * (-m) + sh   '计算熵值
     
201.   End If
     
202. Next
     
203. Next
     
204. Text2.Text = sh                            '输出熵值
     
205. End Function
     
206. Private Sub Command4_Click()                    '在窗口中直接查看重复的初始种群
     
207. Shell "notepad.exe " & App.Path & "\重复", vbNormalFocus
     
208. End Sub
     
209. Private Sub Command6_Click()       '在窗口中查看由混合同余法产生的初始种群                  
     
210. Shell "notepad.exe " & App.Path & "\混合同余法", vbNormalFocus
     
211. End Sub
     
212. Private Sub Command7_Click()       '在窗口中查看由随机函数法产生的初始种群
     
213. Shell "notepad.exe " & App.Path & "\随机函数法", vbNormalFocus
     
214. End Sub
     
215. Private Sub Command8_Click()       '查看每条染色体的个数
     
216. Load Form4
     
217. Form4.Show
     
218. End Sub

复制代码