Add files via upload

Author: Evolutionary-Intelligence

tutorials/plot_pure_nash_equilibria_for_ill_conditioned_functions.py

@@ -0,0 +1,43 @@
+"""This script has been used in Qiqi Duan's Ph.D. Dissertation (HIT&SUSTech).
+
+    Chinese: 该绘图脚本被段琦琦的博士论文（哈工大与南科大联合培养）所使用。
+"""
+import matplotlib
+import matplotlib.pyplot as plt
+
+from pypop7.benchmarks.utils import generate_xyz
+from pypop7.benchmarks.base_functions import rosenbrock
+
+
+matplotlib.rcParams['font.family'] = 'sans-serif'
+matplotlib.rcParams['font.sans-serif'] = 'SimSun'
+matplotlib.rcParams['axes.unicode_minus'] = False
+matplotlib.rcParams['font.size'] = 10  # 对应5号字体
+
+
+def cd(x):  # from https://arxiv.org/pdf/1610.00040v1.pdf
+    return 7.0 * (x[0] ** 2) + 6.0 * x[0] * x[1] + 8.0 * (x[1] ** 2)
+
+
+# helper function for 2D-plotting
+def plot_contour(func, x, y):
+    x, y, z = generate_xyz(func, x, y, 1000)
+    levels = [0, 1, 10, 100, 1000, 10000, 100000, 1000000]
+    plt.contourf(x, y, z, cmap='bone', levels=levels)
+    plt.contour(x, y, z, colors='white')
+
+
+if __name__ == '__main__':
+    ndim_problem = 2
+    bound=[-10.0, 10.0]
+    plt.figure(figsize=(2.5, 2.5))
+    plt.xlim(bound)
+    plt.ylim(bound)
+    plt.xticks(fontsize=10)
+    plt.yticks(fontsize=10)
+    plot_contour(rosenbrock, bound, bound)
+    plt.xlabel('维度1', fontsize=10)
+    plt.ylabel('维度2', fontsize=10, labelpad=-1)
+    plt.scatter([0.0], [0.0], c='green', s=24)
+    plt.savefig(str(rosenbrock.__name__) + '.png', dpi=700, bbox_inches='tight')
+    plt.show()

tutorials/plot_pure_nash_equilibria_for_nonseparable_functions.py

@@ -0,0 +1,41 @@
+"""This script has been used in Qiqi Duan's Ph.D. Dissertation (HIT&SUSTech).
+
+    Chinese: 该绘图脚本被段琦琦的博士论文（哈工大与南科大联合培养）所使用。
+"""
+import matplotlib
+import matplotlib.pyplot as plt
+
+from pypop7.benchmarks.utils import generate_xyz
+
+
+matplotlib.rcParams['font.family'] = 'sans-serif'
+matplotlib.rcParams['font.sans-serif'] = 'SimSun'
+matplotlib.rcParams['axes.unicode_minus'] = False
+matplotlib.rcParams['font.size'] = 10  # 对应5号字体
+
+
+def cd(x):  # from https://arxiv.org/pdf/1610.00040v1.pdf
+    return 7.0 * (x[0] ** 2) + 6.0 * x[0] * x[1] + 8.0 * (x[1] ** 2)
+
+
+# helper function for 2D-plotting
+def plot_contour(func, x, y):
+    x, y, z = generate_xyz(func, x, y, 500)
+    plt.contourf(x, y, z, cmap='bone')
+    plt.contour(x, y, z, colors='white')
+
+
+if __name__ == '__main__':
+    ndim_problem = 2
+    bound=[-10.0, 10.0]
+    plt.figure(figsize=(2.5, 2.5))
+    plt.xlim(bound)
+    plt.ylim(bound)
+    plt.xticks(fontsize=10)
+    plt.yticks(fontsize=10)
+    plot_contour(cd, bound, bound)
+    plt.xlabel('维度1', fontsize=10)
+    plt.ylabel('维度2', fontsize=10, labelpad=-1)
+    plt.scatter([0.0], [0.0], c='green', s=24)
+    plt.savefig(str(cd.__name__) + '.png', dpi=700, bbox_inches='tight')
+    plt.show()

tutorials/plot_pure_nash_equilibria_for_nonsmooth_functions.py

@@ -0,0 +1,41 @@
+"""This script has been used in Qiqi Duan's Ph.D. Dissertation (HIT&SUSTech).
+
+    Chinese: 该绘图脚本被段琦琦的博士论文（哈工大与南科大联合培养）所使用。
+"""
+import matplotlib
+import matplotlib.pyplot as plt
+
+from pypop7.benchmarks.utils import generate_xyz
+
+
+matplotlib.rcParams['font.family'] = 'sans-serif'
+matplotlib.rcParams['font.sans-serif'] = 'SimSun'
+matplotlib.rcParams['axes.unicode_minus'] = False
+matplotlib.rcParams['font.size'] = 10  # 对应5号字体
+
+
+def non_smooth(x):
+    return abs(x[0] - x[1]) - min(x)
+
+
+# helper function for 2D-plotting
+def plot_contour(func, x, y):
+    x, y, z = generate_xyz(func, x, y, 500)
+    plt.contourf(x, y, z, cmap='bone')
+    plt.contour(x, y, z, colors='white')
+
+
+if __name__ == '__main__':
+    ndim_problem = 2
+    bound=[-10.0, 10.0]
+    plt.figure(figsize=(2.5, 2.5))
+    plt.xlim(bound)
+    plt.ylim(bound)
+    plt.xticks(fontsize=10)
+    plt.yticks(fontsize=10)
+    plot_contour(non_smooth, bound, bound)
+    plt.xlabel('维度1', fontsize=10)
+    plt.ylabel('维度2', fontsize=10, labelpad=-1)
+    plt.plot([-10.0, 10.0], [-10.0, 10.0], c='green', linewidth=3)
+    plt.savefig(str(non_smooth.__name__) + '.png', dpi=700, bbox_inches='tight')
+    plt.show()

tutorials/plot_pure_nash_equilibria_for_rotated_functions.py

@@ -0,0 +1,43 @@
+"""This script has been used in Qiqi Duan's Ph.D. Dissertation (HIT&SUSTech).
+
+    Chinese: 该绘图脚本被段琦琦的博士论文（哈工大与南科大联合培养）所使用。
+"""
+import matplotlib
+import matplotlib.pyplot as plt
+
+from pypop7.benchmarks.rotated_functions import ellipsoid
+from pypop7.benchmarks.rotated_functions import generate_rotation_matrix
+from pypop7.benchmarks.utils import generate_xyz
+
+
+matplotlib.rcParams['font.family'] = 'sans-serif'
+matplotlib.rcParams['font.sans-serif'] = 'SimSun'
+matplotlib.rcParams['axes.unicode_minus'] = False
+matplotlib.rcParams['font.size'] = 10  # 对应5号字体
+font_size = 10
+
+
+# helper function for 2D-plotting
+def plot_contour(func, x, y):
+    levels = [0, 50000, 20000000, 50000000,
+            100000000, 150000000, 200000000]
+    generate_rotation_matrix(ellipsoid, 2, 72)
+    x, y, z = generate_xyz(func, x, y, 500)
+    plt.contourf(x, y, z, cmap='bone', levels=levels)
+    plt.contour(x, y, z, colors='white')
+
+
+if __name__ == '__main__':
+    bound=[-10.0, 10.0]
+    plt.figure(figsize=(2.5, 2.5))
+    plt.xlim(bound)
+    plt.ylim(bound)
+    plt.xticks(fontsize=font_size)
+    plt.yticks(fontsize=font_size)
+    plot_contour(ellipsoid, bound, bound)
+    plt.xlabel('维度1', fontsize=font_size)
+    plt.ylabel('维度2', fontsize=font_size, labelpad=-1)
+    plt.scatter([0.0], [0.0], c='green', s=24)
+    plt.savefig(str(ellipsoid.__name__) + '.png',
+                dpi=700, bbox_inches='tight')
+    plt.show()