matlab应用之求解非线性方程需要注意的问题

2011-11-3 15:41:47 2805 matlab

0 刚才自己把参数输入错误，用fsolve求解的时候老是提示找不到解，实际上这个问题是无解的。原问题可以简化为求解如下方程: 1899392341sin(a))/600745 - (55375870cos(a))/120149 - 47944874479/240298=0'，在0到pi/2之间的解。其实在[0 pi/2]之间是没有解的，画图可以看出： ezplot('(1899392341sin(a))/600745 - (55375870cos(a))/120149 - 47944874479/240298',[0 pi/2]) 以后一定要细心。一个用二分法求解方程根的程序如下： function root=MultiRoot(f,a,b,eps) if(nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if(f1==0) root=a; end if(f2==0) root=b; end if(f1f2>0) disp('两端点函数值乘积大于0!'); return; else tol=1; fun=diff(sym(f)); ddf=diff(sym(fun)); fa=subs(sym(f),findsym(sym(f)),a); fb=subs(sym(f),findsym(sym(f)),b); dfa=subs(sym(fun),findsym(sym(fun)),a); dfb=subs(sym(fun),findsym(sym(fun)),b); if(dfa>dfb) root=a; else root=b; end while(tol>eps) r1=root; fx=subs(sym(f),findsym(sym(f)),r1); dfx=subs(sym(fun),findsym(sym(fun)),r1); ddfx=subs(sym(ddf),findsym(sym(ddf)),r1); root=r1-fxdfx/(dfxdfx-fxddfx); tol=abs(root-r1); end end 例子如下： MultiRoot('x^2-0.5',0,1,1e-7) ans = 0.7071 0
2011-11-3 15:41:47　　评论淘帖0 举报相关推荐 • 非线性方程的求解 3671 • matlab求解非线性方程组问题 3083 • 基于并联机器人非线性方程求解 2 • labview求解非线性方程组 4995 • 用fsolve求解非线性方程组 6499 • 关于非线性方程的解法 3012 • 问题:matlab实现牛顿迭代法求解非线性方程 14378 • 线性方程组求解fortran程序 0 • C语言求解线性方程组 8295 • MATLAB矩阵运算、线性方程组求解、特征值与特征向量 3648