回归代码_线性回归进阶版,多项式线性回归讲解与实现(附完整代码)

每天给小编五分钟,小编用自己的代码,带你轻松学习深度学习!本文将会带你做完一个深度学习进阶版的线性回归---多项式线性回归,带你进一步掌握线性回归这一深度学习经典模型,然后在此基础上,小编将在下篇文章带你实现神经网络,并且用它实现对数据集的训练。野蛮智能,小白也能看懂的人工智能。

eea6d5a559a34f0e3196151bee772e89.png

本文代码环境:Anaconda3+pytorch1.0,python版本为3.4-3.6,如果你还没有搭建好环境,可以参考小编的:文科生也能看懂的深度学习入门宝典:pytorch+tensorflow快速上手

关于多项式线性回归的前生今世小编已经在深度学习模型速成,三分钟解决经典线性回归模型(附完整代码)进行了讲解。本文要讲的是一个进阶版的线性回归---多项式线性回归。

多项式线性回归

多项式线性回归的最大优点就是把原来的x1,x2,x3,这些数据更加多层次的描述变成了,例如x1,就变成了三个对应权重分别乘以x1的一次方,x1的二次方,x1的三次方。这样做的好处是可以更加准确的描述数据的特点。方便我们在选择其他模型前,对数据有一个更好地预估。

多项式线性回归代码实现

讲解完了多项式线性回归的基本原理,下面进入代码的实现。代码如下:

d52d319c077f3a02cb5b82daa7e3ca76.png

代码略微有点长,但还是希望大家耐心看完。首选是导入必要的包,这一步操作和小编的深度学习模型速成,三分钟解决经典线性回归模型一文中的操作是一样一样的。然后make features​函数是将原来的数据进行升维操作。将原来的数据由{x1,x2,x3,x4},变成{[x1**1,x1**2,x1**3],[x2**1,x2**2,x2**3]},然后数据就从一个维度变为了三个维度,这样可以更加具体的描述数据。画出来的曲线也更加具体,准确。

然后定义数据,定义好我们的数据集:x,和y=wx+b,这是我们希望将x作为数据集后训练得到的最拟合的结果。(注意,这里的w和x是矩阵,所以她们的乘法要用mm,mm指的是矩阵乘法)。最终通过,f(x)这一函数得到结果。

然后get_batch函数定义了一组训练集,和对应的结果y。并将其返回。poly_model是我们定义的模型。criterion指的是损失函数,optimizer选择了梯度下降法进行优化。最后在while(ture)下进行训练。

测试结果如下:

d6ff2d66bf17da49a0d364b6363bb204.png

为了方便大家对比,我把蓝点向上移了一个单位。通过对比,我们发现,结果还是很准确的。野蛮智能,小白也能看懂的人工智能。欢迎大家评论。

#ifndef FUNCTION_H_ #define FUNCTION_H_ #include #include #include "polyfit.h" #include using namespace std; dxs::dxs() { ifstream fin("多项式拟合.txt"); fin>>n; x=new float[n]; y=new float[n]; for(int i=0;i>x[i]; } for(i=0;i>y[i]; } cout<>nn; m=nn+1; u=new float*[m]; for(i=0;i<m;i++) { u[i]=new float[m+1]; }//创建m行,m+1列数组 } void dxs::dfine() { for(int i=0;i<m;i++) { for(int j=0;j<m+1;j++) { u[i][j]=0; } } for(i=0;i<m;i++) { for(int j=0;j<m;j++) { for(int k=0;k<n;k++) { u[i][j]=u[i][j]+pow(x[k],j+i); } } } for(i=0;i<m;i++) { for(int k=0;k<n;k++) { u[i][m]=u[i][m]+pow(x[k],i)*y[k]; } } } void dxs::show() { for(int i=0;i<m;i++) { for(int j=0;j<m+1;j++) { cout<<u[i][j]<<" ";//<<endl; } cout<<endl; } ////显示具有m行m+1列u数组的各元素值 } void dxs::select_main(int k,float **p,int m) { double d; d=*(*(p+k)+k); //cout<<d; int l=k; int i=k+1; for(;i fabs(d)) { d=*(*(p+i)+k); l=i; } else continue; } if(d==0) cout<<"错误"; else { if(k!=l) { for(int j=k;j<m+1;j++) { double t; t=*(*(p+l)+j); *(*(p+l)+j)=*(*(p+k)+j); *(*(p+k)+j)=t; } } } } void dxs::gaosi() { for(int k=0;k<m;k++) { select_main(k,u,m);//调用列主元函数 for(int i=1+k;i<m;i++) { // *(*(p+i)+k)=(float) *(*(p+i)+k) / *(*(p+k)+k); u[i][k]=(float) u[i][k] / u[k][k]; } for(i=k+1;i<m;i++) { for(int j=k+1;j=0;i--) { float a=0; for(int j=i+1;j<m;j++) { //a=a + (*(*(p+i)+j) * *(*(p+j)+m)); a=a+u[i][j] * u[j][m]; } //*(*(p+i)+n-1)= (*(*(p+i)+n-1) - a) / *(*(p+i)+i); u[i][m]= (u[i][m] -a) / u[i][i]; } cout<<"方程组的解为:"<<endl; for(i=0;i<m;i++) { cout<<"a"<<i+1<<"="; cout<<u[i][m]<<endl; // l[i]=*(*(p+i)+n-1); } cout<<"y="<<u[0][m]; for(i=1;i<m;i++) { cout<<showpos<<u[i][m]<<"x"; if(i!=1)cout<<"^"<<noshowpos<<i; } cout<<endl; } dxs::~dxs() { delete[]x,y; delete []*u; } #endif
拟合功能说明: 对于形如y=a1*x^n+a2*x^(n-1)+...+an*x+a(n+1)的函数,若a1,a2,...,a(n+1)未知,则可以由n+1对不同的(x,y)值求出这些a的值。更进一步的说,给你一条X-Y轴平面里的曲线,或是,已知一系列的(x,y)值,要你用一个近似的一元函数拟合它,该怎么办?这个时候你不知道该一元函数的最高次数是多少,即不知道未知数系数a1,a2,...,a(n+1)的个数。在这里显然要先假设一个n的值,再来算出这些系数值:有n+1个未知数系数要求,就得有n+1对(x,y)的值,这样可以精确算出唯一的解。但是,如果给你的(x,y)值对数多于n+1对该怎么办(因为你碰到的情况几乎都是这样的)?你算出来的近似函数得兼顾所有的(x,y)取样点,使得误差最小。你不能只取所有(x,y)取样点中的n+1对计算,而扔掉了多出来的那些取样点。该程序的拟合功能就是完成这个任务的。它会计算所有的取样点,最后得出一个误差最小的函数。 执行过程如下: 先读取右边两个编辑框内的(x,y)取样点,再弹出一个对话框让你假设一个一元函数,即最高次数n,然后程序会开始计算。注意,若n设的太大,计算过程可能会很久,会造成程序假死现象。得出的结果会直接以函数的形式显示在下面的长条形编辑框内。你可以通过点击“←画图”和“画图→”来进行比较拟合的近似程度。注意,你设置的n值不是越大,就拟合得越好,也不是越小,拟合得越好,原因你细想一下就会明白。程序所用的计算原理是最小二乘法,所用的计算方法是矩阵运算。你可以百度相关知识。 该程序的另外一个“辨识”功能,用到的也是最小二乘法,传递函数等相关知识是大学工科学生会学到的。 其他不懂的问题请联系作者:hastings1986@163.com。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值