我们以一简单数据组来说明什么是线性回归。假设有一组数据型态为 y=y(x),其中
x={0, 1, 2, 3, 4, 5}, y={0, 20, 60, 68, 77, 110}
如果我们要以一个最简单的方程式来近似这组数据,则非一阶的线性方程式莫属。先将这组数据绘图如下
图中的斜线是我们随意假设一阶线性方程式 y=20x,用以代表这些数据的一个方程式。以下将上述绘图的 MATLAB 指令列出,并计算这个线性方程式的 y 值与原数据 y 值间误差平方的总合。
>> x=[0 1 2 3 4 5];
>> y=[0 20 60 68 77 110];
>> y1=20*x; % 一阶线性方程式的 y1 值
>> sum_sq = sum(y-y1).^2); % 误差平方总合为 573
>> axis([-1,6,-20,120])
>> plot(x,y1,x,y,'o'), title('Linear estimate'), grid
如此任意的假设一个线性方程式并无根据,如果换成其它人来设定就可能采用不同的线性方程式;所以我们 须要有比较精确方式决定理想的线性方程式。我们可以要求误差平方的总合为最小,做为决定理想的线性方 程式的准则,这样的方法就称为最小平方误差(least squares error)或是线性回归。MATLAB的polyfit函数提供了 从一阶到高阶多项式的回归法,其语法为polyfit(x,y,n),其中x,y为输入数据组n为多项式的阶数,n=1就是一阶 的线性回归法。polyfit函数所建立的多项式可以写成
从polyfit函数得到的输出值就是上述的各项系数,以一阶线性回归为例n=1,所以只有 二个输出值。如果指令为coef=polyfit(x,y,n),则coef(1)= , coef(2)=,...,coef(n+1)= 。注意上式对n 阶的多 项式会有 n+1 项的系数。我们来看以下的线性回归的示范:
>> x=[0 1 2 3 4 5];
>> y=[0 20 60 68 77 110];
>> coef=polyfit(x,y,1); % coef 代表线性回归的二个输出值
>> a0=coef(1); a1=coef(2);
>> ybest=a1*x+a0; % 由线性回归产生的一阶方程式
>> sum_sq=sum(y-ybest).^2); % 误差平方总合为 356.82
>> axis([-1,6,-20,120])
>> plot(x,ybest,x,y,'o'), title('Linear regression estimate'), grid
-
点赞数:0回答数:1
-
点赞数:0回答数:1
-
点赞数:0回答数:2
-
点赞数:0回答数:1
-
点赞数:0回答数:1
-
点赞数:0回答数:1
-
点赞数:0回答数:1
-
点赞数:0回答数:5
-
点赞数:0回答数:1
-
点赞数:0回答数:1

最新问答: W=UIt 又 I =q/t 所以 q=It 所以W=Uq 问 以上所提到的q代表什么 一个长方体水箱,低是正方形,水箱的高是4分米,它的侧面积是40平方分米,这个水箱的表面积是多少平方分米 听录音,选出你听到的句子。 ( ) 1. A. My sister is She likes red.()()()she()? 1.湖水深10m,小鸟在距湖面6m的空中飞翔,它在湖中的像距离湖面____m. 科学 声音是由( )产生的 一片只有5到选择题的阅读,急One time every month my village has a basketba Unable to log on to zhe game server.please try again later. 给几句名言要中英翻译的 张师傅7天做了8个零件,每个零件做了几分之几天. 两个球体积之和为12π,且这两个球大圆周长之和为6π,那么这两球半径之差是( ) be adj.to sth.和be adj.about sth和be adj.over sth.填入不同的介词有什么区别 以下关于会议纪要的说法正确的是(多选) 6.4一元一次方程的应用(2)1.若银行一年定期储蓄的年利率为2.25%,小丽的妈妈取出一年到期的本金及利息时,扣除了利 “长江后浪推前浪,一代更比一代强”,这句话化用的诗句是什么? 你已不再乐于敞开心扉的英文 如果地球没有了重力,大气层还存在么? 从甲地到乙地,客车要行4小时,货车要行5小时,客车和货车的时间比是______:______,客车和货车的速度比是___ 如图,二次函数y=ax的平方+2X-3的图像与x轴有一个交点在0和1之间(不含0和1),则a的 氯气与碱反应生成什么?氯气和碱反应,生成次氯酸盐和金属氯化物,这句话对吗?书上的答案说对,但不是还要生成水吗?