[問題] 利用fsolve求得的根，怎麼代入程式後半部的計算？已回收
各位大大您好：
謝謝各位忙碌的閒暇之餘，還要解答我的問題!!
我現在有一個問題想要請教各位大大一下!!
如以下所敘：
我利用fsolve求出來的根(即 h1~h4)，
我現在要將此根運用到程式後半部的運算，
請問我該怎麼插入？
我的程式碼如下：
算完 h1~h4 之後，
想要以R_c = h1~h4 帶入計算公式中，
但似乎又變成error~~無法執行~~!!
function pttex120
clear all;
[h1,fval]=fsolve(@specialeqex,0)
[h2,fval]=fsolve(@specialeqex,10)
[h3,fval]=fsolve(@specialeqex,100)
[h4,fval]=fsolve(@specialeqex,-100)
R_c=[h1,h2,h3,h4];
function f = specialeqex(R)
-----------------------------------------↓以下就是我要做的後半部運算
Ae=[];
for p=1:(t(i)-1966)
n=(P/d)*A_ref(i)*exp(lnda*(t-(1965+p))) * (( 1- P*(R_c/d ))
^((1965+p)-1966) );
Ae=[Ae n];
end
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;
拜託大大幫我解答一下~~~~冏rz
--
愛，是心中輕盈的詩
是幸福，是寂寞
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.190.27
→
06/05 14:49, , 1F
06/05 14:49, 1F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
