python与机器学实践-何宇健 源代码及过程中遇到的问题

时间:2024-09-03 14:34:33

# -*- 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()