例68 Lambert函数=lam(x),x是变量,是 方程we=x的解,对不同的x,求解 然后绘制(x) ↓求解策略和过程 使用for循环 使用匿名函数描述 e生成w向量|0=1 ambert(x) 绘制函数曲线 → MATLAB求解语句 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 例 6.8 Lambert 函数 ， 是变量， 是 方程 的解， 对不同的 ，求解 然后绘制 求解策略和过程 使用for循环 使用匿名函数描述 生成w向量 绘制函数曲线 MATLAB求解语句：
MATLAB表述 ∈>>y [];xx=0:.05:10;x0=0 h=optimset; h Display=off for x=xx f=@(w)w*exp(w)-x; y1=fsolve(f, xo, h); x0=y1;y=[y,y1]; end P lot(xx, y →直接求解,使用1 ambert函数 a > yO=lambertw(xx); plot(xx, yo) 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 MATLAB表述 直接求解，使用lambertw函数
4例6.9 →数值方法求解(2 x2+ 1=0 0.75x3-y+0.9=0 ↓选择变量x1=x,x2=y 把原始常微分方程组ODEs)变为 x2+x2-1=0 变成矩阵形式 0.75x3-x2+0.9=0 f(a) x+x5-1 0.75x1-x2+0.9=0 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 例 6.9 数值方法求解 选择变量 把原始常微分方程组(ODEs)变为 变成矩阵形式
描述方程的方法 M-函数 function y=myfun(x) y=[x(1)*x(1)+x(2)*x(2)-1 0.75*x(1)^3-x(2)+0.9]; 匿名函数 >>f=Q(p)[p(1)*p(1)+p(2)*p(2)-1;, 0.75*P(1)^3-p(2)+0.9]; Inline函数 ∈>>f=in1ine('[p(1)*P(1)p(2)*p(2)-1; 0.75*P(1)^3p(2)+0.9];’,p3) 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 描述方程的方法 M-函数 匿名函数 Inline函数
当初值选为xo=1,2]T EE > OPT=optimset; OPT. LargeScale=off>; [x, Y,c,d]=fsolve(f, [1; 2], OPT) 当使用另一个搜索初始点co=[-1,0T >>[x,Y,c,d]= fsolve(f,[-1,0]’,0PT) x, norm(y), kk=d. funcCount 注意:选择不同的初值可以得出不同的结果 2/20/2021星期六, 2008-9-6,13:09:29 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:09:29 当初值选为 当使用另一个搜索初始点 注意：选择不同的初值可以得出不同的结果