[請教]fsolve求解的結果是初值
各位大大您好
有個小小的問題想要請教
就是有個三角函數的方程式要解
使用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
→
05/11 04:09, , 1F
05/11 04:09, 1F
→
05/11 04:10, , 2F
05/11 04:10, 2F