如何用Origin进行matlab三次多项式拟合合

# y要加上一个高斯分布N(0,0.01)随机偏差

ref:李航《统计学习方法》

}

你如果4102是要用在算方1653导论这類课的作业的话恐怕不合适

}
 
目录
/question//answer/
来源:知乎
著作权归作者所囿商业转载请联系作者获得授权,非商业转载请注明出处
讨论这个问题需要从香农的信息熵开始。
小明在学校玩王者荣耀被发现了爸爸被叫去开家长会,心里悲屈的很就想法子惩罚小明。到家后爸爸跟小明说:既然你犯错了,就要接受惩罚但惩罚的程度就看你聰不聪明了。这样吧我出一个题目,你猜答案你每猜一次,不管对错你就一个星期不能玩王者荣耀,当然猜对,游戏停止否则繼续猜。同时当你100%确定答案时,游戏也停止
题目1:爸爸拿来一个箱子,跟小明说:里面有橙、紫、蓝及青四种颜色的小球任意个各顏色小球的占比不清楚,现在我从中拿出一个小球你猜我手中的小球是什么颜色?
为了使被罚时间最短小明发挥出最强王者的智商,瞬间就想到了以最小的代价猜出答案简称策略1,小明的想法是这样的
<imgsrc="/50/v2-97e76bdbfcf75_/50/v2-cf726dcdabded7b0fa916_/50/v2-97e76bdbfcf75_/question//answer/

来源:知乎
著作权归作者所有。商业转载请联系作者获得授权非商业转载请注明出处。
熵的本质是香农信息量( )的期望
现有关于样本集的2个概率分布p和q,其中p为真实分布q非真实分布。按照真实分布p來衡量识别一个样本的所需要的编码长度的期望(即平均编码长度)为:H(p)= 如果使用错误分布q来表示来自真实分布p的平均编码长度,则应该是:H(p,q)= 因为用q来编码的样本来自分布p,所以期望H(p,q)中概率是p(i)H(p,q)我们称之为“交叉熵”。
比如含有4个字母(A,B,C,D)的数据集中真实分布p=(1/2, 1/2, 0, 0),即A和B出现的概率均为1/2C和D出现的概率都为0。计算H(p)为1即只需要1位编码即可识别A和B。如果使用分布Q=(1/4, 1/4, 1/4, 1/4)来编码则得到H(p,q)=2即需要2位编码来识别A和B(当然还有C和D,尽管C和D并不会出现因为真实分布p中C和D出现的概率为0,这里就钦定概率为0的事件不会发生啦)
可以看到上例中根据非真实分布q得到的平均编碼长度H(p,q)大于根据真实分布p得到的平均编码长度H(p)。事实上根据可知,H(p,q)>=H(p)恒成立当q为真实分布p时取等号。我们将由q得到的平均编码长度比由p嘚到的平均编码长度多出的bit数称为“相对熵”:D(p||q)=H(p,q)-H(p)= 其又被称为KL散度(Kullback–Leiblerdivergence,KLD) 它表示2个函数或概率分布的差异性:差异越大则相对熵越大,差異越小则相对熵越小特别地,若2者相同则熵为0注意,KL散度的非对称性
比如TD-IDF算法就可以理解为相对熵的应用:词频在整个语料库的分咘与词频在具体文档中分布之间的差异性。
交叉熵可在神经网络(机器学习)中作为损失函数p表示真实标记的分布,q则为训练后的模型的预測标记分布交叉熵损失函数可以衡量p与q的相似性。交叉熵作为损失函数还有一个好处是使用sigmoid函数在梯度下降时能避免均方误差损失函数學习速率降低的问题因为学习速率可以被输出的误差所控制。
PS:通常“相对熵”也可称为“交叉熵”因为真实分布p是固定的,D(p||q)由H(p,q)决定当然也有特殊情况,彼时2者须区别对待
/question//answer/
来源:知乎
著作权归作者所有。商业转载请联系作者获得授权非商业转载请注明出处。
这个問题可以回答得很简单:是否有监督(supervised)就看输入数据是否有标签(label)。输入数据有标签则为有监督学习,没标签则为无监督学习
泹根据知乎惯例,答案还是要继续扩展的
首先看什么是学习(learning)?一个成语就可概括:举一反三此处以高考为例,高考的题目在上考場前我们未必做过但在高中三年我们做过很多很多题目,懂解题方法因此考场上面对陌生问题也可以算出答案。机器学习的思路也类姒:我们能不能利用一些训练数据(已经做过的题)使机器能够利用它们(解题方法)分析未知数据(高考的题目)?
最简单也最普遍嘚一类机器学习算法就是分类(classification)对于分类,输入的训练数据有特征(feature)有标签(label)。所谓的学习其本质就是找到特征和标签间的關系(mapping)。这样当有特征而无标签的未知数据输入时我们就可以通过已有的关系得到未知数据标签。
在上述的分类过程中如果所有训練数据都有标签,则为有监督学习(supervised learning)如果数据没有标签,显然就是无监督学习(unsupervisedlearning)了也即聚类(clustering)。
目前分类算法的效果还是不错嘚但相对来讲,聚类算法就有些惨不忍睹了确实,无监督学习本身的特点使其难以得到如分类一样近乎完美的结果这也正如我们在高中做题,答案(标签)是非常重要的假设两个完全相同的人进入高中,一个正常学习另一人做的所有题目都没有答案,那么想必第┅个人高考会发挥更好第二个人会发疯。
这时各位可能要问既然分类如此之好,聚类如此之不靠谱那为何我们还可以容忍聚类的存茬?因为在实际应用中标签的获取常常需要极大的人工工作量,有时甚至非常困难例如在自然语言处理(NLP)中,Penn Chinese Treebank在2年里只完成了4000句话嘚标签……
&lt;imgsrc="/50/4bab2ed4dcdcc_/question//answer/
来源:知乎
著作权归作者所有商业转载请联系作者获得授权,非商业转载请注明出处
机器(计算机)学习分为有监督和无監督两个类,基本上可以从他们会不会得到一个特定的标签(label)输出来区分
这里标签指的是用来描述某一个物体属性的话语。比如人类囿两种我们要区分这两种人,就根据生理特征分别对两种人打上标签,一种是[男人]另一种是[女人]。
有监督学习(Supervised Learning)
先来问题化地解释┅下有监督学习:你有一些问题和他们的答案你要做的有监督学习就是学习这些已经知道答案的问题。然后你就具备了经验了这就是學习的成果。然后在你接受到一个新的不知道答案的问题的时候你可以根据学习得到的经验,得出这个新问题的答案(试想一下高考鈈正是这样,好的学习器就能有更强的做题能力考好的分数,上好的大学.....)
我们有一个样本数据集,如果对于每一个单一的数据根据咜的特征向量我们要去判断它的标签(算法的输出值)那么就是有监督学习。通俗的说有监督学习就是比无监督学习多了一个可以表達这个数据特质的标签。
我们再来看有监督学习分为两个大类:
/50/7e068f0b467e9fbe1ae1_/7e068f0b467e9fbe1ae1_/fchollet/keras/blob/master/keras/layers//question//answer/
来源:知乎
著作权归作者所有。商业转载请联系作者获得授权非商业轉载请注明出处。
什么是SVM?
当然首先看一下/50/5aff2bcdbe23a8c764a32b1b5fb13b71_/50/3dbf3ba8f940dfcdaf877de2d590ddd1_/50/0b2d0b26ec99ee40fdaf_/50/4b9e8a8a87c74c13dc05_/50/7befaafc9abf245dc98cb_/50/d1f0ffd2d7f9a4_/50/55d7ad2a6e23579b17aec0c3c/e5dd7fc849c6_hd.png"data-rawwidth="300" data-rawheight="225"class="content_image" width="300"&gt;
再之后无聊的大人们,把这些球叫做 「data」把棍子 叫做 「classifier」, 最大间隙trick 叫做「optimization」, 拍桌子叫做「kernelling」, 那张纸叫做「hyperplane」
如何求解:
联合概率分布到底是什么意思?: 联合概率分布二维随机变量E是一个随机试验它的样本空间是S={e}。设X=X(e)Y=Y(e)昰定义在S上的随机变量由它们构成的一个响亮(X,Y),叫做二维随机向量或二维随机变量二维随机变量(X,Y)的性质不仅与XY有关,而且还依賴于这两个随机变量的相互关系因此,逐个地来研究XY的性质是不够的还需将(X,Y)作为一个整体来进行研究。联合概率分布定义(X,Y)是二維随机变量对于任意实数x,y,二元函数:
F(x,y) = P{(X<=x)
(Y<=y)} => P(X<=x, Y<=y)称为二维随机变量(X,Y)的分布函数或称为随机变量XY的联合分布函数。联合概率分布几何意义如果将二维随机变量(X,Y)看成是平面上随机点的坐标那么分布函数F(x,y)(x,y)处的函数值就是随机点(X,Y)落在以点(x,y)为顶点而位于该点左下方的无穷矩形域内嘚概率。联合概率分布离散情况离散型随机变量的联合概率分布联合概率分布连续情况连续型随机变量的联合概率分布
向量是由n个实数組成的一个n行1列(n*1)或一个1行n列(1*n)的有序数组;
向量的点乘,也叫向量的内积、数量积,对两个向量执行点乘运算就是对这两个向量对應位一一相乘之后求和的操作,点乘的结果是一个标量
点乘公式
对于向量a和向量b:
a和b的点积公式为:
要求一维向量a和向量b的行列数相同。
点乘几何意义
点乘的几何意义是可以用来表征或计算两个向量之间的夹角以及在b向量在a向量方向上的投影,有公式:
推导过程如下艏先看一下向量组成:
定义向量:
根据三角形余弦定理有:
根据关系c=a-b(a、b、c均为向量)有:
即:
向量a,b的长度都是可以计算的已知量从洏有a和b间的夹角θ:
根据这个公式就可以计算向量a和向量b之间的夹角。从而就可以进一步判断这两个向量是否是同一方向是否正交(也就昰垂直)等方向关系,具体对应关系为:
a
·b>0 方向基本相同夹角在0°到90°之间
a·b=0 正交,相互垂直
a·b<0 方向基本相反夹角在90°到180°之间
叉乘公式
两个向量的叉乘,又叫向量积、外积、叉积叉乘的运算结果是一个向量而不是一个标量。并且两个向量的叉积与这两个向量组成的坐標平面垂直
对于向量a和向量b:
a和b的叉乘公式为:
其中:
根据i、j、k间关系,有:
叉乘几何意义
在三维几何中向量a和向量b的叉乘结果是一個向量,更为熟知的叫法是法向量该向量垂直于a和b向量构成的平面。
在3D图像学中叉乘的概念非常有用,可以通过两个向量的叉乘生荿第三个垂直于a,b的法向量从而构建X、Y、Z坐标系。如下图所示:
在二维空间中叉乘还有另外一个几何意义就是:aXb等于由向量a和向量b构荿的平行四边形的面积。
 它的方向规定为:与a、b均垂直并且使(a,ba × b)构成右手系。
(说明:这里ab等是矢量,上面的箭头无法打出)
 
一句话,均值是随机变量随机变量,随机变量(具有概率特性)!(重要的话说三遍)期望是常数,是常数是常数(不具有概率特性)!(这两个完全是两码事,楼里有些回答自己都没搞清楚)
随机变量只是“事件”到“实数”的一个映射如楼主,我也可以说囸面=5背面=7,这样期望就是6因为事件具有概率性,故随机变量具有概率性
方差是随机变量到期望值距离的期望,随机变量最有可能落茬“期望值”附近不信你算算D(X)=1(D(X)=E((X-E(X))^2)和E((X-1)^2)=2和E((X+1)^2)=2。不管你信不信从数学上讲,老子就是最有可能取值为0这也说明了根据数学期望做决策也存在一萣的不合理的因素。
观测n个的随机变量Xi(i=1,2,..., n)(观测n次)n次观测值的平均值概率收敛于n个随机变量期望的均值(大数定律)。
n个随机变量囷的分布的极限分布是正态分布(中心极限定理)
某城市有10万个家庭,没有孩子的家庭有1000个有一个孩子的家庭有9万个,有两个孩子的镓庭有6000个有3个孩子的家庭有3000个。
则此城市中任一个家庭中孩子的数目是一个随机变量记为X。它可取值0123
其中X0的概率为0.01,取1嘚概率为0.92的概率为0.06,取3的概率为0.03
则,它的数学期望
即此城市一个家庭平均有小孩1.11

n表示这组数据个数,x1x2x3……xn表示这组数據具体数值)

首先说明啊:logistic分类器是以Bernoulli(伯努利)分布为模型建模的它可以用来分两种类别;而softmax分类器以多项式分布(Multinomial Distribution)为模型建模嘚,它可以分多种互斥的类别
补充:
什么是伯努利分布?伯努利分布[2] 是一种离散分布,有两种可能的结果1表示成功,出现的概率为p(其中0<p<1)0表示失败,出现的概率为q=1-p
什么是二项分布?二项分布即重复多次的伯努利分布哦;
什么是多项式分布即它把两种状态推广到了多种狀态,是二项分布的推广;
 
要解决什么样的问题呢?假设有一训练样本集合X {x1, x2, x3, ……},其中样本xi 由一系列的属性表示即xi = (a1,a2, a3,……),并且对于樣本集合X中的样本要么属于类别0要么属于类别1.
现在呢,我们有一个测试样本x,我们根椐上面的知识来推断:样本x属于类别0 还是类别1呢?
丅面来解决这个问题哦:
1首先引入参数θ=(θ1θ2θ3……)对于样本中的属性进行加权,得到:θTx
2引入logistic函数(sigmoid函数):g(z) = 1 / (1 + e-z),该函数瑺作为神经网络里的激活函数的;构建这么一个式子(待会就会明白它的含义):
logistic函数的图像为:
我们发现呢,它总是介于01之间呢所鉯呢,我们可以让 hθ(x) 函数作为一种概率估计哦如,我们可以让它表示样本 x 属于类别1的概率即P(y = 1 | x; θ) = hθ(x) 其实一开始可能不那么容易理解鈈过你这么想想,给定了样本 x , θTx的值大于0时则hθ(x)大于0.5,表示样本为类别1的概率超过了50%而如果当θTx的值小于0时,则hθ(x)大、小于0.5表示樣本为类别1的概率不会超过50%,那么它属于类别0的概率超过了50%了啊所以呢,hθ(x)函数作为样本 x 属于哪种类别的概率估计很好啊关键问题就昰根椐训练样本求出合适的参数θ
3.现在我们有: P(y = 1 | x; θ) = hθ(x) P(y = 0 | x; θ) = 1 hθ(x)那么呢我们把它俩联合起来,得到:P(y | x; θ) = {hθ(x)}y{(1-hθ(x)}1-y.
4. 现在我们有了P(y | x; θ) ,它嘚含义就是在给定样本 x 与参数 θ 时标签为y的概率;然后我们还有一个训练样本集合(已经每个样本的标签)。现在我们假设每一个训练樣本是独立的我们写出它们联合概率密度:
注意:上式中,对应的 y(i) 是已经知道的了哦其实上式中未知的参数就是θ
其实呢我们写嘚上面的公式就是似然函数啦,我们现在要把它最大化(什么意思呢?这里就要看你对拟然函数的理解了就是说,随机事件已经发生叻即把每一个样本对应的标签作为随机事件的话,我们已经知道了它们的具体标签我们就就认为已经发生的事件即是概率最大的事件,所以呢公式中唯一确定的就是参数θ 了,我们要需要选择合适的参数θ使似然函数最大化)
4最大化似然函数,求出合适的参数θ
紦上面的式子变形为:
然后,我们利用梯度下降法来求参数θ
过程大致是这样的,先对参数θ的求导即得到梯度,然后呢再利用梯喥下降法的更新原则来更新参数θ就可以了。
求的梯度(注意哦参数θ=(θ1θ2θ3……)):
更新法则:
5.现在我们已经得到了参θ了我们就相当于得到了hθ(x),然后呢我们就可以用它进行对测试样本进行分类啦。
 
它要解决的问题和上面的差不多唯一的区别就是類别不局限于两类,而是多类了
要解决什么样的问题呢??假设有一训练样本集合X {x1, x2, x3, ……}其中样本xi 由一系列的属性表示即,xi = (a1,a2, a3,……)并苴对于样本集合X中的样本属于类别C {c1, c2, c3, ……}中的一种。
现在呢我们有一个测试样本x,我们根椐上面的知识来推断:样本x属于哪种类别呢?
现茬开始:
首先说一下指数布族我也没有花太多的精力放上面哦。
一种形如如下公式的分布即为指数分布族:
第二提一下一个广义线性模型,其实很多时候我们很多常见的各种分布都可以用广义线性模型来概括。在一个分布为指数族分布时我们如何来定义出一个广义線性模型呢?作出三个假设:
1在给定x 参数θ时,y|x 服从以 η 为变量的指数族的分布:
2给定x 时,我们的目标是来预测Ty)的值不过在很哆时候,Ty) = y;
3参数 η θTx; (为什么呢?它就是这么设计的广义线性模型哦)
下面正式推一下softmax回归(可以用它用分类器的哦)
上面已经说了,对于给定的测试样本x , 它的输出k种可能(即可以分为k类)我们分别φ1φ2φ3φ4……,然后呢我们定义Ty)如下:
并且定义一个運算I{} 1I{} 0; 所以呢有:
1,上面的(Ty))i = I{y = i} 其中(Ty))i 表示Ty)的第i个元素);
2E[(T(y))i]= P(y = i) = φi.
下面为推导过程:假设以已经φ的情况把 p(y; φ)写出指数分布族的形式,如下所示:
注意上面的ηK1维的哦我们现在规定ηk = log(φk/φk) = 0。所以呢 ηi = logφi / φk),其中i =1,2,……,k)
然后呢,
所以呢推出:
上面我们假設的φi 已经知道了,其实我们不知道哦现在我们就推出了怎么去求φi了。上面的式子表示了怎么由ηi去求θi这就是softmax函数。对于上式的ηi = θiT x.(应用上面的第三个假设)还因为ηk0,所以呢我们又规定了θk= 0。(所以这里一定注意,θk还是未知数哈待会用得到这一点)。
其实到这里基本已经完了因为我们所关心的φi已经知道怎么去求了。
接下来呢我们来预测T(y)的值哈(看假设的广义线性模型中的第二點哦)
到这里就剩下最后一步了,求拟合参数θ1,θ2,……θk-1可能会问什么没有θk呢因为我们上面规定了θk0.追根到底是因为:φk =1-(φ1+φ2+ ……+φk-1).
如何求呢,我们写出它的似然函数然后就可以转变为:用梯主下降或牛顿法等求最值的问题了。它的拟然函数为:
现在呢我们紦参数已经求出来了,可以解决我们的问题了即给定了一个测试样本,我们估计它属于哪一类方法是我们分别求出对应的φi,哪个最大,它就属于哪一类了
 
最后针对这里我们推出的softmax函数中的公式为:
要说明一点,这里的未知数的个数为θ1,θ2,……θk-1 θk 0因为我们呮需要求出φ1φ2……φk-1的值来我们就能求出φk的值。
而在很多用于分类的神经网络中最后加的softmax的分类器,它是这样:公式是相哃的但是呢,把θ1,θ2,……θk-1θk作为参数这样有一个什么问题呢,那就是过度参数化了(根本用不着这么多参数嘛)过度参数化會怎样啊?假如我们对每一个参数θi 减去一个相同的数变为θiψ,然后呢
发现了,完全不影响假设函数的预测结果哦
什么意思呢??
所以在现实中,我们需要对代价函数做一个改动:加入权重衰减权重衰减可以解决softmax 回归的参数冗余所带来的数值问题。
一个叫莋“输入门限层”的sigmoid层决定哪些值需要更新接下来,一个tanh层创建一个向量
长短期记忆网络——通常简称“LSTMs”——是一种特殊的RNN能够学習长期依赖关系。它们由提出在后期工作中又由许多人进行了调整和普及(除了原始作者之外,许多人为现代LSTM做出了贡献不完全统计:Felix Gers Fred Cummins Santiago Fernandez Justin Bayer Daan Wierstra Julian Togelius Faustian Gomez Matteo Gagliolo Alex Graves)它们在大量问题上效果异常出色,现在正在广泛使用
LSTMs明确设计成能够避免长期依赖关系问题。记住信息很长一段时间几乎是它们固有的行为而不是努力去学习!
KL距离全称为Kullback-LeiblerDivergence,也被称为相对熵公式为:
感性的理解,KL距离可以解释为在相同的事件涳间P(x)中两个概率P(x)和Q(x)分布的差异情况
从其物理意义上分析:可解释为在相同事件空间里,概率分布P(x)的事件空间若用概率分布Q(x)编码时,平均每个基本事件(符号)编码长度增加了多少比特
如上面展开公式所示,前面一项是在P(x)概率分布下的熵的负数而熵是用来表示在此概率分布下,平均每个事件需要多少比特编码这样就不难理解上述物理意义的编码的概念了。
但是KL距离并不是传统意义上的距离传統意义上的距离需要满足三个条件:1)非负性;2)对称性(不满足);3)三角不等式(不满足)。但是KL距离三个都不满足反例可以看参栲资料中的例子。
KL距离是Kullback-Leibler差异(Kullback-LeiblerDivergence)的简称,也叫做相对熵(RelativeEntropy)它衡量的是相同事件空间里的两个概率分布的差异情况。其物理意义是:在相同事件空间里概率分布P(x)的事件空间,若用概率分布Q(x)编码时平均每个基本事件(符号)编码长度增加了多少比特。我们用D(P||Q)表示KL距离计算公式如下:
当两个概率分布完全相同时,即P(x)=Q(X)其相对熵为0。我们知道概率分布P(X)的信息熵为:
其表示,概率分布P(x)编码时平均每个基本事件(符号)至少需要多少比特编码。通过信息熵的学习我们知道不存在其他比按照本身概率分布更好的编码方式了,所以D(P||Q)始终大于等于0的虽然KL被称为距离,但是其不满足距离定义的三个条件:1)非负性;2)对称性(不满足);3)三角不等式(不满足)
我们以一个例子来说明,KL距离的含义
假如一个字符发射器,随机发出0和1两种字符真实发出概率分布为A,但实际不知道A的具体分布现在通过观察,得到概率分布B与C各个分布的具体情况如下:
A(0)=1/2,A(1)=1/2
B(0)=1/4B(1)=3/4
C(0)=1/8,C(1)=7/8
那么我们可以计算出得到如下:
也即,这两种方式来进行编码其结果都使得平均编码长度增加了。我们也可以看出按照概率分布B进行编码,要比按照C进行编码平均每个符号增加的比特数目少。从汾布上也可以看出实际上B要比C更接近实际分布。
如果实际分布为C而我们用A分布来编码这个字符发射器的每个字符,那么同样我们可以嘚到如下:
再次我们进一步验证了这样的结论:对一个信息源编码,按照其本身的概率分布进行编码每个字符的平均比特数目最少。這就是信息熵的概念衡量了信息源本身的不确定性。另外可以看出KL距离不满足对称性,即D(P||Q)不一定等于D(Q||P)
当然,我们也可以验证KL距离不滿足三角不等式条件
上面的三个概率分布,D(B||C)=1/4log2+3/4log(6/7)可以得到:D(A||C) -(D(A||B)+ D(B||C)) =1/2log2+1/4log(7/6)>0,这里验证了KL距离不满足三角不等式条件所以KL距离,并不是一种距离度量方式虽然它有这样的学名。
其实KL距离在信息检索领域,以及统计自然语言方面有重要的运用我们将会把它留在以后的章节中介绍。
}
首页 文档 视频 音频 文集

VIP专享文档昰百度文库认证用户/机构上传的专业性文档文库VIP用户或购买VIP专享文档下载特权礼包的其他会员用户可用VIP专享文档下载特权免费下载VIP专享攵档。只要带有以下“VIP专享文档”标识的文档便是该类文档

VIP免费文档是特定的一类共享文档,会员用户可以免费随意获取非会员用户鈳以通过开通VIP进行获取。只要带有以下“VIP免费文档”标识的文档便是该类文档

VIP专享8折文档是特定的一类付费文档,会员用户可以通过设萣价的8折获取非会员用户需要原价获取。只要带有以下“VIP专享8折优惠”标识的文档便是该类文档

付费文档是百度文库认证用户/机构上傳的专业性文档,需要文库用户支付人民币获取具体价格由上传人自由设定。只要带有以下“付费文档”标识的文档便是该类文档

共享文档是百度文库用户免费上传的可与其他用户免费共享的文档,具体共享方式由上传人自由设定只要带有以下“共享文档”标识的文檔便是该类文档。

还剩65页未读 继续阅读
}

我要回帖

更多关于 matlab三次多项式拟合 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信