数学实验 课件10.pptVIP

实验十：寻优(Optimization) 10.1 光的折射定律 law of light refraction 光在同一种介质(medium)中走最短路线（直线），聪明！ 从一种介质进入另一种介质发生折射，傻了？ 问题： 设光在两种介质中的速度比为 u，两介质分界线(separatrix)为直线。从第一介质的A点到第二介质的B点，走怎样的路线时间最短？ T=f(x)=Sqrt[x^2+a^2] +Sqrt[(c-x)^2+b^2]/u 求 f(x) 的最小值(minimum 1 均匀搜索法 区间(interval) [a,b]内取n个点, d为步长(step) a0=a, a1=a+d, a2=a+2d ,…, an=a+nd, an+1=b 依次计算对应的yk=f(a+kd),找到最佳的yk,把对应的xk作为最佳点x*的近似值. 误差|x*-xk|d 均匀搜索法 2) 0.618 搜索法(求每一点的 f 值需花成本 确定最初的搜索区间[a,b] 记d=b-a. 插入两分点 x1= a+0.382d (0.382点）x2= a+0.618d （0.618点） 求出 f(x1),f(x2). 如果f(x1)优于f(x2),用[a,x2] 代替[a,b]。 x1已是0.618 点，只须再求0.382点的值。 否则，用[x1,b]代替[a,b]。 不断重复以上过程。 3)求导法 求出使导函数y’=f’(x)等于0的x=x*和使y’不存在的x=x*值 如果在x*两侧有 xx*,f’(x)0;xx*,f’(x)0, 则x=x*是f(x)的极小值点 可以用迭代法求f’(x)的根 如果不容易求出f(x)的二阶导数值，可用如下方法将f(x)在x*附近展成二次函数： 求出y0=f(x0),y’0=f’(x0),如果y’0≠0在x0附近选一点x0+d,计算y1=f(x0+d)的值。 确定二次函数g(x)=ax2+bx+c的系数a,b,c使满足条件：g(x0)=y0, g’(x0)=2ax0+b=y’0, g(x0+d)=y1 如果a0则用二次函数g(x)的的最小值点 X1=-b/2a=x0-y’0/2a=x0-y’0d2/(2(y1-y0-dy’0)) 作为最小值点的近似值。用x1代替x0重复上述过程。 如果a0,g(x)的最小值点在端点，我们直接用x0+d代替x0重复上述过程 4） 利用Mathematica： 画出 y=f(x) 的图象，观察最小值点的近似位置x=c. 运行语句 FindMinimum[f,{x,c}] 10.2 奔向最优点 问题：设A,B,C是平面上三点。求点P(x,y)到三点距离之和 s=f(x,y) 最小。 Mathematica语句：FindMinimum[f[x,y],{x,x0},{y,y0}] 多元搜索：最速下降法 从初始位置 P0(x0,y0)出发 求P0点及点(x0+d,y0), (x0,y0+d)的f值s0,s1,s2。 u=(s1-s0)/d,v=(s2-s0)/d. 梯度方向G=(u,v). 最快下降方向 - G=(-u,-v) 选步长 h。从 P0 到 P1(x0-hu,y0-hv)。 从 P1出发再前进。 10.3 最小二乘法 数据点集t={(xi,yi)}. 用直线 y=kx+b 去拟合，使所有 kxi+b-yi 的平方和 s(k,b)最小。 Fit[t,{1,x},x] 用下列语句求达到最小值的a,b * * * * *