单层感知器实践

import numpy as np
import matplotlib.pyplot as plt

#输入数据
X = np.array([[3,3],
              [4,3],
              [1,1],
              [0,2]])

# 添加偏置项
X = np.concatenate((np.ones((4,1)),X),axis=1)
print(X)

# 标签
Y = np.array([[1],
              [1],
              [-1],
              [-1]])
# 权值初始化，3行1列，取值范围-1到1
W = (np.random.random([3,1])-0.5) *2

print(W)

# 学习率设置
lr = 0.11
# 神经网络输出
O = 0

def update():
    global X,Y,W,lr
    O = np.sign(np.dot(X,W))
    W_C = lr * (X.T.dot(Y-O))/int(X.shape[0])
    W = W + W_C

for i in range(100):
    update() # 更新当前权值
    print(W) # 打印当前权值
    print(i) # 打印迭代次数
    O = np.sign(np.dot(X,W))
    if (O == Y).all():
        print('Finished')
        print('epoch',i)
        break

#正样本
x1 = [3,4]
y1 = [3,3]

#负样本
x2 = [1,0]
y2 = [1,2]

#计算分界线的斜率以及截距
# w0 + w1x1 + w2x2 = 0
# x1 看作x x2看作y
# y = -w0/w2 - w1/w2(x)
k = -W[1] / W[2]
d = -W[0] / W[2]

xdata = (0,5)

plt.figure()
plt.plot(xdata,xdata * k +d,'r')
plt.scatter(x1,y1,c='b')
plt.scatter(x2,y2,c='y')
plt.show()