菜鸟之路——机器学习之非线性回归个人理解及python实现

关键词：

梯度下降：就是让数据顺着梯度最大的方向，也就是函数导数最大的放下下降，使其快速的接近结果。

Cost函数等公式太长，不在这打了。网上多得是。

这个非线性回归说白了就是缩小版的神经网络。

python实现：

 import numpy as np
 import random

 def graientDescent(x,y,theta,alpha,m,numIterations):#梯度下降算法
     xTrain =x.transpose()
     for i in range(0,numIterations):#重复多少次
         hypothesis=np.dot(x,theta)          #h函数
         loss=hypothesis-y

         cost=np.sum(loss**2) / (2*m)
         print("Iteration %d / cost:%f"%(i,cost))
         graient=np.dot(xTrain,loss)/m
         theta=theta-alpha*graient
     return theta

 def getData(numPoints,bias,variance):#自己生成待处理数据
     x=np.zeros(shape=(numPoints,2))
     y=np.zeros(shape=numPoints)
     for i in range(0,numPoints):
         x[i][0]=1
         x[i][1] = i
         y[i]=(i+bias)+random.uniform(0,1)*variance
     return x,y

 X,Y=getData(100,25,10)
 print("X:",X)
 print("Y:",Y)

 numIterations=100000
 alpha=0.0005
 theta=np.ones(X.shape[1])
 theta=graientDescent(X,Y,theta,alpha,X.shape[0],numIterations)
 print(theta)

运行结果：

......输出数据太多，只截取后面十几行

Iteration 99988 / cost:3.930135
Iteration 99989 / cost:3.930135
Iteration 99990 / cost:3.930135
Iteration 99991 / cost:3.930135
Iteration 99992 / cost:3.930135
Iteration 99993 / cost:3.930135
Iteration 99994 / cost:3.930135
Iteration 99995 / cost:3.930135
Iteration 99996 / cost:3.930135
Iteration 99997 / cost:3.930135
Iteration 99998 / cost:3.930135
Iteration 99999 / cost:3.930135
[30.54541676 0.99982553]

其中遇到一个错误。

TypeError: unsupported operand type(s) for *: 'builtin_function_or_method' and 'float'

因为我第五行xTrain =x.transpose()。刚开始没加括号。直接用的xTrain =x.transpose

打印一下xTrain是

<built-in method transpose of numpy.ndarray object at 0x00000219C1D14850>

只是创建了一个transpose方法，并没有真的给x转置。加上括号就好了。再打印xTrain就能正常显示转置后的x了