MATLAB牛顿法解x^5 3x^4-5x^3=0-搜答案-神马作文网

symsxx0=2;f=x^3-3*x-1;eps=1e-6;maxcnt=1000;fx=diff(f,x);x1=x0;cnt=1;whilecnt

[x,y,z]=solve('x+y=2','y+z=1','y=0')x=2y=0z=1>>x+y+zans=3>>[x,y,z]=solve('x+y=2','y+z=1','y=1')x=1y=

你问了好多次啊?小写的问题还是?

esp表达错了吧,e的指数幂应该用exp表示.另外你这个n怎么定义的?我不大懂啊.给你分享个牛顿迭代法的程序吧,你参考下看能不能解决你的问题：

x=x0-(x0.^3+2*x0.^2+10*x0-20)/(3*x0*x0+4*x0+10);最后应该是10,不是x0

源程序如下：clearclcN=100;x=2;forii=0:Nxl=x;num=ii;x=x-(x^3-3*x+1)/(3*x^2-3);ifabs(x-xl)

你这个出什么问题了?再问：再答：Jacobian这个函数似乎不存在，你有名字为这个的.m文件么？？再问：我试了一下，可以查到。你可以运行吗？再答：那就是大小写的问题了，使用jacobian再问：膜拜。

我是用C得到结果:2.1155229/*牛顿迭代法解方程组的解x0为迭代的初值,n为迭代次数,jingdu为精度function为求根代数式,d2functoin为其导数返回最终符合一定精度的根*/d

matlab中的x=x(:)

利用“:”将矩阵元素按列取出排成一列

functionx=NewRaph(x0,f,fprime,tol)x(1)=x0;i=1;gx=x(i)-f(x(i))/fprime(x(i));while(abs(x(i)-gx)>tol)gx

兄弟,西理工的?

采用第一个.首先你的两个代码的计算过程和方法以及步骤是一致的.只不过第二个将k==N放在循环内部判断是没有必要的.放在while外面,可以节省点计算量.如果你要求结果精度高一些的话,你调用：x=nan

哈啊啊啊啊啊啊啊啊啊啊啊再问：bjijib

function x = newton_method(fun,x0,e,N)if nargin<5 N=50

Newton-Raphson求解非线性方程组matlab源程序matlab程序如下：functionhom[P,iter,err]=newton('f','JF',[7.8e-001;4.9e-001

这样吧,你不要直接通过solve求解啊,可以画一个图像啊,以x的值为横坐标,行列式值为纵坐标,或者纵坐标取一下对数（如果变化范围太大）,然后可以得出大致得到解的范围；最后可以通过其他方法计算更加精确的

x1=x-func1_1(x)/func1_1_1(x);是点除再问：Error:File:func1_1.mLine:1Column:22TheinputcharacterisnotvalidinM

定义函数functiony=nd(x)y=0.036-((x/2090.7).^(1/0.1585))-x/182000functiony=nd0(x)y=-(1/0.1585)*(x/2090.7)

%牛顿迭代法解方程y=x.^3-x-1　　x=1.5;　　formatlong;　　x1=x-func1_1(x)/func1_1_1(x);　　if(abs(x1)1e-6|abs(func1_1(

程序：clcclearx=zeros(1,25);x(1)=2;formatlongfori=1:24;x(i+1)=x(i)-(x(i)^3-3*x(i)+2)/(3*x(i)^2-3);enddi