day5 格式化代码与判断数字

在 上张贴 由 孔 祥博day5 格式化代码与判断数字已关闭评论

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
name=input("Name:")
age=input("age:")
job=input("job:")
salary=input("Salary:")
 
# print(name,age,job,salary)
msg='''
------------------info of %s----------
Name:%s
Age:%s
Job:%s
Salary:%s
----------------end------------------
'''%(name,name,age,job,salary)
 
 
if salary.isdigit():

  %s  s=string %d…

K线修炼2——其他形态

在 上张贴 由 孔 祥博K线修炼2——其他形态已关闭评论

一、上吊线 二、出水芙蓉和断头铡刀 三、覆盖线与乌云盖顶 四、包线（吞没线） 五、孕线 六、待入线、切入线、插…

MACD——八神式

在 上张贴 由 孔 祥博MACD——八神式已关闭评论

一、佛手向上 二、小鸭出水 三、漫步青云 四、天鹅展翅 五、空中揽绳 六、空中缆车 七、海底电缆 八、海底捞月

K线修炼1——剑21式

在 上张贴 由 孔 祥博K线修炼1——剑21式已关闭评论

一、上涨两颗星 实战： 二、红三兵 例子： 三、三个白武士 例子： 四、三阳开泰 例子： 五、倒锤 例子： 六…

混沌优化算法简单Python实现

在 上张贴 由 孔 祥博一条评论

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
import numpy as np
import random
import matplotlib.pyplot as plt
 
 
def myjudge(x1,x2,x3):
    myjudge=0
    a=-x1-x2-2*x3+3
    if x1>0 and x2>0 and x3>0 and a>=0:
        myjudge=1
    else:
        myjudge=0
    return myjudge
 
 
def myfunction(x1,x2,x3):
    myfunction=1-(2*x1**2+2*x2**2+x3**2+2*x1*x2+2*x1*x3-8*x1-6*x2-4*x3+9)
    return myfunction
 
 
x=np.array([[np.float(0)]*5000]*3)
TempX=np.array([np.float(0)]*3)
MaxX=np.array([np.float(0)]*3)
 
Max=np.array([np.float(0)]*100)
MaxX1=np.array([np.float(0)]*100)
MaxX2=np.array([np.float(0)]*100)
MaxX3=np.array([np.float(0)]*100)
 
 
for k in range(100):
    for z in range(100):
        for a in range(3):
            x[a][0]=random.random()
            TempX[a]=2*x[a][0]
        if myjudge(TempX[0],TempX[1],TempX[2])==1:
            break
 
    for g in range(3):
        MaxX[g]=TempX[g]
 
    MaxF=myfunction(MaxX[0],MaxX[1],MaxX[2])
    TempF=0
    for i in range(1,5000):
        for j in range(3):
            x[j][i]=4*x[j][i-1]*(1-x[j][i-1])
            TempX[j]=2*x[j][i]
 
        if myjudge(TempX[0],TempX[1],TempX[2])==1:
            TempF=myfunction(TempX[0],TempX[1],TempX[2])
        if TempF>MaxF:
            MaxX[j]=TempX[j]
            MaxF=TempF
 
    #二次载波
    for i in range(3):
        x[i][0]=random.random()
    for i in range(1,5000):
        for j in range(3):
            x[j][i]=4*x[j][i-1]*(1-x[j][i-1])
 
    for i in range(5000):
        for j in range(3):
            TempX[j]=MaxX[j]+0.0001*x[j][i]
        if myjudge(TempX[0],TempX[1],TempX[2])==1:
            TempF=myfunction(TempX[0],TempX[1],TempX[2])
        if TempF>MaxF:
            MaxX[j]=TempX[j]
            MaxF=TempF  
    
    MaxF=round(MaxF,4)
    for i in range(3):
        MaxX[i]=round(MaxX[i],4)
 
    Max[k]=MaxF
    MaxX1[k]=MaxX[0]
    MaxX2[k]=MaxX[1]
    MaxX3[k]=MaxX[2]
 
plt.subplot(2,2,1)
plt.plot(MaxX1)
plt.subplot(2,2,2)
plt.plot(MaxX2)
plt.subplot(2,2,3)
plt.plot(MaxX3)
plt.subplot(2,2,4)
plt.plot(Max)
plt.show()
 
b=Max[-1]
b1=MaxX1[-1]
b2=MaxX2[-1]
b3=MaxX3[-1]
 
print('x1:',b1,'X2:',b2,'X3:',b3,'Min:',1-b)