[問題]怎麼把fsolve求出來的根運用到後面程式的計算？？已回收
各位大大您好：
我利用fsolve求出某函數的根(即 h1~h4)，
我現在要將此根運用到程式後半部的運算，
請問我該怎麼插入？
我的程式碼如下：
算完 h1~h4 之後，
想要以R_c = h1~h4 帶入計算公式中，
function pttex120
clear all;
[h1,fval]=fsolve(@specialeqex,0)
[h2,fval]=fsolve(@specialeqex,10)
[h3,fval]=fsolve(@specialeqex,100)
[h4,fval]=fsolve(@specialeqex,-100)
function f = specialeqex(R)
Ce=[];
for k=1:(t-1966)
n1=(I_t( 12 +k,2)/R_c)*P(i)*gamma*(1-exp(-R_c/H)) + Ae(k);
Ce=[Ce n1];
end
Cd=[];
for w=1:(t-1966)
n2=Ce(j)*exp(-lnda*(t-(1966+w)));
Cd=[Cd n2];
end
CD=P*sum(Cd);
R(i)=( A_t(i)-A_ref )/ CD;
例如說：R_c = h1 (這我試過了～～不行....冏rz)
帶到n1的公式裡，
n1=(I_t( 12 +k,2)/R_c)*P(i)*gamma*(1-exp(-R_c/H)) + Ae(k);
--
愛，是心中輕盈的詩
是幸福，是寂寞
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.190.27
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
