[請教]fsolve求解的結果是初值

[blackman1122 (小黑鬼)]

各位大大您好

有個小小的問題想要請教

就是有個三角函數的方程式要解

使用fsolve的解 都是設定的初值

有試過把初值設定的很接近解 (從圖形上看出來的)

fsolve怎麼會解不到??

百思不得其解   所以請教各位大大

謝謝

以下是我的code
-------------------------------------------------------------

function [ks]=mymain
Yoe = 1/430 ;
Yoo = 1/34.8;
Yo = 1/92;
ke = 0.506;
ko = 0.873;

ks0=[0.2;0.2];
options=optimset('Display','off');
ks=fsolve(@myfun2,ks0,options,Yoe,Yoo,Yo,ke,ko);


function F=myfun2(ks,Yoe,Yoo,Yo,ke,ko)

F=[Yoe.*Yoo.*(cot(ke).*csc(ko)-cot(ko).*csc(ke))+
   Yo.*(cot(ks)-csc(ks)).*(Yoo.*csc(ko)-Yoe.*csc(ke))];

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

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.112.19.201
→ masspoet:試試看fmincon或fminsearch多變數函數求解
→ masspoet:不過還是有可能解不出來 主要跟起始值有關