python与机械教育初探_python与机器学实践-何宇健 源代码及过程中遇到的问题

# -*- coding: utf-8 -*-

"""

Spyder Editor

This is a temporary script file.

"""

import numpy as np

import matplotlib.pyplot as plt

#第一步 获取与处理数据

x,y=[],[] #定义存储输入数据和目标数据的数组

for sample in open('a.txt','r'): #遍历数据集并保存

_x,_y=sample.split(",")

x.append(float(_x))

y.append(float(_y))

x,y=np.array(x),np.array(y) #转化为numpy数组

x=(x-x.mean())/x.std() #标准化

#将原始数据以散点图的形式画出

plt.figure()

plt.scatter(x,y,c="g",s=6)

plt.show

#第二步：选择与训练模型

#模型：多项式拟合   多项式拟合散点是线性回归很小的一部分

x0=np.linspace(-2,4,100) #在(-2,4)这个区间上取100个点作为画图的基础

#核心代码  仔细研究   建立回归模型

def get_model(deg):#得到模型，这一步尤其重要，要仔细分析

return lambda input_x=x0: np.polyval(np.polyfit(x,y,deg),input_x)

def get_cost(deg,input_x,input_y):#返回损失值

return 0.5*((get_model(deg)(input_x)-input_y)**2).sum()

test_set={1,4,10}

for d in test_set:

print(get_cost(d,x,y))

#第三步：评估与可视化结果

plt.scatter(x,y,c="g",s=20)#s是点的大小即size

for d in test_set:

plt.plot(x0,get_model(d)(),label="degree={}".format(d))

plt.xlim(-2,4)

plt.ylim(1e5,8e5)

plt.legend()

plt.show

过程中遇到的问题：

1.在获取与处理数据的过程中，文件老是找不到，在这里有两种常见的错误

(1)FileNotFoundError: [Errno 2] No such file or directory: 'C：/a.txt'

错误分析与解决：可能是文件路径写错了，也可能是文件名写错了，检查自己设置的文件名后缀，我的错误就是没注意文件名后缀，结果我的文件名实际上是a.txt.txt，所以说一直提示我找不到文件，后来在保存文件的地方重新设置显示文件扩展名，顺利解决问题，可以将数据文件和py源文件放在一个文件夹下，就可以只输入文件名，而不输入文件路径

(2)SyntaxError: (unicode error) 'unicodeescape' codec can't decode bytes in position 2-3: truncated \UXXXXXXXX escape

错误分析与解决：

从他的博客取经的

原因：

window 读取文件可以用\，但是在字符串中\是被当作转义字符来使用，所以’d:\a.txt’会被转义成’d:\a.txt’这是正确路径，所以不会报错。而‘C:\Users\FrankYuan\Pictures\Camera Roll\WIN_20161010_08_51_57_Pro.jpg ’中经过转义之后可能就找不到路径的资源了，例如\t可能就转义成tab键了。

解决办法

python在描述路径时可以有多种方式，现列举常见的三种

方式一:转义的方式

'd:\\a.txt'

方式二:显式声明字符串不用转义

'd:r\a.txt'

方式三:使用Linux的路径/         最推荐

'd:/a.txt'

我强烈推荐第三种写法，这在Linux和window下都是行的通的。

2.标签的正确拼写 label   这个错误犯过很多次了，每次拼写成lable，导致报错，去网上搜还出来一堆看起来特别合理的解释，说什么IDE问题，哎，实际上就是粗心大意拼写错误

3.#核心代码  仔细研究   建立回归模型，多看多思考

def get_model(deg):#得到模型，这一步尤其重要，要仔细分析

return lambda input_x=x0: np.polyval(np.polyfit(x,y,deg),input_x)

def get_cost(deg,input_x,input_y):#返回损失值

return 0.5*((get_model(deg)(input_x)-input_y)**2).sum()