python牛顿法解非线性方程组_用牛顿迭代法解非线性方程组

题目：

用牛顿迭代法解非线性方程组

有两个非线性方程,未知数是x1,x2：

(15x1+10x2)/[(40-30x1-10x2)^2×(15-15x1)]=5e-4;

(15x1+10x2)/[(40-30x1-10x2)×(10-10x2)]=4e-2.

大侠们帮帮忙,实在整不出来

解答：

function [r,n]=mulNewton(x0,eps)

if nargin==1

eps=1.0e-4;

end

r=x0-myf(x0)*inv(dmyf(x0));

n=1;

tol=1;

while tol>eps

x0=r;

r=x0-myf(x0)*inv(dmyf(x0));

tol=norm(r-x0);

n=n+1;

if(n>100000)

disp('迭代步数太多,方程可能不收');

return;

end

function f=myf(x)

x1=x(1);

x2=x(2);

f1=(15*x1+10*x2)-((40-30*x1-10*x2)^2*(15-15*x1))*5e-4;

f2=(15*x1+10*x2)-((40-30*x1-10*x2)*(10-10*x2))*4e-2;

f=[f1 f2];

function df=dmyf(x)

x1=x(1);

x2=x(2);

df=[ (3*(27*x1^2 + 12*x1*x2 - 66*x1 + x2^2 - 14*x2 + 60))/4,...

(3*x1*x2 - 3*x2 - 21*x1 + 9*x1^2 + 32)/2;(-3)*(4*x2 - 9),(-2)*(6*x1 + 4*x2 - 15)]

----------------------------------------------------------------------------------------------------------

[r,n]=mulNewton([0.2 0.6],0.001)

df =

30.9600 13.3600

19.8000 22.8000

df =

33.1741 13.7408

22.0443 24.5260

df =

35.4490 14.2583

21.4349 24.9400

df =

35.2223 14.2212

21.1843 24.7600

df =

34.9859 14.1678

21.2334 24.7120

df =

34.9828 14.1660

21.2601 24.7231

r =

0.1203 0.4785

n =

python牛顿法解非线性方程组_用牛顿迭代法解非线性方程组相关推荐

python牛顿法解非线性方程组_matlab实现牛顿迭代法求解非线性方程组.pdf
matlab实现牛顿迭代法求解非线性方程组.pdf matlab 实现牛顿迭代法求解非线性方程组实现牛顿迭代法求解非线性方程组已知非线性方程组如下 3*x1-cosx2*x3-1/20 x12-81 ...
matlab牛顿法解非线性方程组,matlab实现牛顿迭代法求解非线性方程组.pdf
资源描述 matlab 实现牛顿迭代法求解非线性方程组实现牛顿迭代法求解非线性方程组已知非线性方程组如下 3*x1-cos(x2*x3)-1/2=0 x1^2-81*(x2+0.1)^2+sin(x ...
c语言牛顿迭代法求解非线性方程组,利用牛顿迭代法求解非线性方程组
最近一个哥们,是用牛顿迭代法求解一个四变量方程组的最优解问题,从网上找了代码去改进,但是总会有点不如意的地方,迭代的次数过多,但是却没有提高精度,真是最近一个哥们,是用牛顿迭代法求解一个四变量方程组 ...
非线性方程组牛顿迭代法matlab,matlab实现牛顿迭代法求解非线性方程组
<matlab实现牛顿迭代法求解非线性方程组>由会员分享,可在线阅读,更多相关<matlab实现牛顿迭代法求解非线性方程组(5页珍藏版)>请在人人文库网上搜索. 1.matla ...
python牛顿法解非线性方程组_牛顿迭代法解非线性方程组（MATLAB版）
牛顿迭代法,又名切线法,这里不详细介绍,简单说明每一次牛顿迭代的运算:首先将各个方程式在一个根的估计值处线性化(泰勒展开式忽略高阶余项),然后求解线性化后的方程组,最后再更新根的估计值.下面以求解最简 ...
python牛顿法解非线性方程组_科学网—求解多元非线性方程组F(x)=0的Newton-Raphson方法及其MATLAB实现 - 王福昌的博文...
科学网对公式支持不太好,在博客园有相同博文牛顿迭代法可以推广到多元非线性方程组 $boldsymbol{F}(boldsymbol{x})=boldsymbol{0}$的情况,称为牛顿-- 拉夫逊方 ...
Matlab利用牛顿迭代法求解非线性方程组
我们得首先了解牛顿迭代法的原理是什么: 在这里,我以二阶非线性方程组为例: f1(x,y)=0 f2(x,y)=0,求解x,y 原理假设方程组的一组近似解为(x0,y0),将方程f1(x,y)=0与 ...
二元牛顿迭代法matlab,牛顿迭代法解二元方程组以及误差分析 matlab实现
电子科大数值分析第一次作业类似一元函数的泰勒展开: f(x) f(x0) f'(x0)(x x0) 二元函数的展开为: f(x,y) f(x0,y0) [(x x0) (y y0)]f(x,y)x ...
python牛顿迭代法推导求根的迭代公式_推导牛顿迭代法
推导牛顿迭代法牛顿法是方程求根的一个有力方法,常常能快速求出其他方法求不出或者难以求出的解. 假定有一个函数 y=f ( x ) ,方程 f(x)=0 在 x = r 处有一个根,对于此根,我们先 ...

python牛顿法解非线性方程组_用牛顿迭代法解非线性方程组

python牛顿法解非线性方程组_用牛顿迭代法解非线性方程组相关推荐

最新文章

热门文章