程序很简单的,有点基础的都能编。但是有一个点必须提醒你,Jacobe矩阵的计算。
两种方法:
解析Jacobe。就是需要你具体的研究非线性方程组的解析表达式,通过求导的方法,得到解析表示的Jacobe矩阵。
数值Jacobe。不需要具体研究方程组的解析式,只需要一个知道输入输出关系的黑盒子。先输入一个x0,得到y0,再输入x0+dx,得到y1,dy/dx=(y1-y0)/dx,dx很小,取10E-6即可,用割线斜率代替导数。
两种给你方法都是常用的办法,各有利弊。
你的问题解决了吗