引用回帖:

9樓: Originally posted by 鐘聲入翠微 at 2021-07-08 14:23:57
謝謝！（30000,0.05）這個初始值是自己根據(jù)實際情況設(shè)置，還是隨意設(shè)置的呢？這次我運行出來了，但是感覺結(jié)果不是我想象中的那樣。
...

這個需要嘗試，可以用 fsolve 函數(shù)隨機賦值求解。

結(jié)果不理想就要考慮方程是否正確，或者a,b取值是否合理。

function example_0708
clear;
a=0.709133065;
b=0.08555;
x0 = [1,1];
x = fsolve(@root2d,x0);
fprintf('\tx=%.8f\n',x(1));
fprintf('\ty=%.8f\n',x(2));

function F = root2d(x)
F(1) = (1+(0.0172*x(1))^2)^(-x(2)/2)-a;
F(2) = x(2)*atan(0.0172*x(1))-b;
end

end
，

引用回帖:

10樓: Originally posted by 獨孤神宇 at 2021-07-08 16:30:14
這個需要嘗試，可以用 fsolve 函數(shù)隨機賦值求解。

結(jié)果不理想就要考慮方程是否正確，或者a,b取值是否合理。

function example_0708
clear;
a=0.709133065;
b=0.08555;
x0 = ;
x = fsolve(@root2d,x0);
...

你的解答對我?guī)椭�很大，謝謝！