• 豆搜网 > 非线性方程组Matlabfsolve机械
  • 非线性方程组Matlabfsolve机械

    免费下载 下载此文档 文档格式:PDF   更新时间:2008-05-01   下载次数:0   点击次数:4
    文档基本属性
    文档语言:
    文档格式:pdf
    文档作者:Share
    关键词:
    主题:
    备注:
    点击这里显示更多文档属性
    非线性方程组—docsou.comb—fsolve—机械 CAD 论坛—docsou.com
    非线性方程组求解-docsou.comb-fsolve
    实例一: ①建立文件 fun.m: docsou.comn y=fun(x) y=[x(1)-0.5*sin(x(1))-0.3*cos(x(2)), ... x(2) - 0.5*cos(x(1))+0.3*sin(x(2))]; ②>>clear;x0=[0.1,0.1];fsolve(@fun,x0,docsou.comt('fsolve')) 注:...为续行符 m 文件必须以 docsou.comn 为文件头, 调用符为@; 文件名必须与定义的函数名相同; fsolve()主要求解复杂非线性方程和方程组,求解过程是一个逼近过程.
    实例二: ①建立文件 fun.m
    docsou.comn F=myfun(x) F=[x(1)-3*x(2)-sin(x(1));2*x(1)+x(2)-cos(x(2))];
    ②然后在命令窗口求解:
    >> x0=[0;0]; >> docsou.comns=docsou.comt('docsou.comy','iter'); >> [x,fv]=fsolve(@myfun,x0,docsou.comns) %docsou.comB 显示的优化过程 %设定求解初值 %设定优化条件 %优化求解 Norm of step
    First-order Trust-docsou.comn docsou.comn docsou.comunt f(x) docsou.comty docsou.com 0 3 1 2 1 1 6 0.000423308 0.5 0.0617 1 2 9 5.17424e-010 0.00751433 4.55e-005 1.25 3 12 9.99174e-022 1.15212e-005 9.46e-011 1.25 docsou.comn docsou.comed: first-order docsou.comty is less than docsou.comns.TolFun. x= 0.4966 0.0067 fv = 1.0e-010 * 0.3161 0.0018
    非线性方程组—docsou.comb—fsolve—机械 CAD 论坛—docsou.com
    实例三:
    求下列非线性方程组在(0.5,0.5) 附近的数值解. (1) 建立函数文件 myfun.m. docsou.comn q=myfun(p) x=p(1); y=p(2); q(1)=x-0.6*sin(x)-0.3*cos(y); q(2)=y-0.6*cos(x)+0.3*sin(y); (2) 在给定的初值 x0=0.5,y0=0.5 下,调用 fsolve 函数求方程的根. x=fsolve('myfun',[0.5,0.5]',docsou.comt('docsou.comy','off')) x= 0.6354 0.3734
    非线性方程组—docsou.comb—fsolve—机械 CAD 论坛—docsou.com
    实例四:
    【例 5.6.3-1 】 求 解 二 元 函 数 方 程 组
    2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2
    f 1 ( x, y ) sin( x y ) 0 的零点. f 2 ( x, y ) cos( x y ) 0
    (0)从三维坐标初步观察两函数图形相交 情况
    x=-2:0.05:2;y=x;[X,Y]=meshgrid(x ,y); %产生 x-y 平面上网点坐标 F1=sin(X-Y);docsou.coms(X+Y); F0=zeros(size(X)); surf(X,Y,F1), docsou.coml('x'),docsou.coml('y'), view([-31,62]),hold on, surf(X,Y,F2),surf(X,Y,F0), docsou.comterp, hold off
    -1.5
    -1
    -0.5
    0
    0.5
    1
    1.5
    2
    图 5.6.3-1
    两个二元函数 0 等位线的交点图
    (2)从图形获取零点的初始近似值 在图 5.6.3-1 中,用 ginput 获取两个函 数 0 等位线(即三线组中间那条线)交点的 坐标.

    下一页

  • 下载地址 (推荐使用迅雷下载地址,速度快,支持断点续传)
  • 免费下载 PDF格式下载