關于求解非線性方程組的非負解問題

是 否 +2 論壇幣

k人 參與回答

經(jīng)管之家送您一份

應屆畢業(yè)生專屬福利!

求職就業(yè)群

趙安豆老師微信：zhaoandou666

經(jīng)管之家聯(lián)合CDA

送您一個全額獎學金名額~ !

立即領取

感謝您參與論壇問題回答

經(jīng)管之家送您兩個論壇幣！

+2 論壇幣

想解一個非線性方程組，求非負解，自己試著編了程序，現(xiàn)有如下疑問，特請教高人指點。程序代碼如下：function F = myfun1(x)
F=[x(1)*cos(x(1))^2+x(2)*cos(x(2))^2-x(1);x(1)*sin(x(1))^2+x(3)*cos(x(3))^2-x(2);x(2)*sin(x(2))^2+x(3)*sin(x(3))^2-x(3);x(1)+x(2)+x(3)-1]

x0=[1;0;0];
options=optimset('display','iter');
[x,fval]=fsolve('myfun1',x0,options)
結果執(zhí)行如下：
x =
    0.6807
    0.2944
    0.0248
但是當我設初值為x0=[0;0;1]時，結果為：
x =
   -0.4704
   -0.0971
    1.5674

請問如何設置才能使得在設置可以使得在初值x0=[0;0;1]的情況下，得到非負解？

掃碼加我 拉你入群

請注明：姓名-公司-職位

以便審核進群資格，未注明則拒絕

分享0 收藏0 回帖

關鍵詞：非線性方程組 線性方程組 非線性方程 線性方程 方程組 方程組

請大家?guī)兔?/div>

youngblues 發(fā)表于 2015-5-7 09:45
請大家?guī)兔?/blockquote>

非線性方程的解可能不存在, 也可能不唯一, 你的問題很難解決.

enxizheng 發(fā)表于 2015-5-7 11:38
非線性方程的解可能不存在, 也可能不唯一, 你的問題很難解決.

非常感謝

用1stOpt，加大于零的范圍約束就行了，不需猜初值：

1. Parameter x(3)>0;
   
2. Function
   
3. x(1)*cos(x(1))^2+x(2)*cos(x(2))^2-x(1);
   
4. x(1)*sin(x(1))^2+x(3)*cos(x(3))^2-x(2);
   
5. x(2)*sin(x(2))^2+x(3)*sin(x(3))^2-x(3);x(1)+x(2)+x(3)-1;

復制代碼

x1: 0.68074188638042
x2: 0.294444475972497
x3: 0.0248136376327703