博主在解决一个包含13个方程和13个未知数的非线性方程组时，使用fsolve函数得到的结果精度不理想且每次计算结果波动大。为提高计算精度，寻求对fsolve函数参数调整或解法优化的建议。文中给出了自定义函数myfun的详细代码，该函数定义了方程组，以及主程序中使用fsolve的设置，如最大迭代次数和函数评估次数。
我有一个13方程13未知数的方程组需要求解，使用的是fsolve函数，设定初值之后每次计算的结果都不同且相差比较大，精度不高。
希望各位大神能帮我修改一下计算程序以提高计算精度。
下面是自定义函数的代码：
function q=myfun(p)
t3=13;t5=58;t8=94;t9=80;h1=2666.1;h2=305.6;h3=2405.9;D=3.643;%已知系数
%下面是需要求解的方程组，一共13个方程，13个未知数。
x4=p(1); x5=p(2); x6=p(3); x8=p(4); x9=p(5); t4=p(6); t6=p(7); h4=p(8); h5=p(9); h6=p(10); h8=p(11); h9=p(12); G=p(13);
q(1)=(G-D)*(h8-h9)-G*(h6-h5);
q(2)=(G-D)*x8-G*x6;
q(3)=D*h1+(G-D)*h8-G*h6+D*(h3-h2)-20000;
q(4)=(h8-h9)/(h8-h4)-0.65;
q(5)=G*t4-D*t3-(G-D)*t9;
q(6)=G*x5-(G-D)*x9;
q(7)=x8-x9;
q(8)=x5-x6;
q(9)=-2024.33+163.309*x4-4.88161*x4^2+(6.302948e-2)*x4^3-(2.913705e-4)*x4^4+t4*(18.2829-1.1692*x4+(3.24804e-2)*x4^
最低0.47元/天 解锁文章
706
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
