[問題] 變數儲存問題 >"<~~已回收

看板MATLAB作者winddoor1227 (白小衣)時間16年前 (2009/08/12 18:14)推噓1(1推 0噓 4→)

留言5則, 2人參與討論串1/1

各位大大您好：小妹我有一點小問題想要請教大大們，我撰寫程式求出函式的根，並利用回圈，使輸入已知的data不同，而會有不同的根值被計算出來， EX: i 根 i=1 R_1 i=2 R_2 .......... 以此類推，可是我該怎麼讓計算出來的根都一一記錄起來，而不是被一直覆蓋掉，拜託各位大大幫我看一下~~ 以下是我寫的程式碼，請大大幫幫忙!! 有分為主程式跟副程式兩個檔案~~ 主 → pttex130.m R_c=[]; clear all; [h1,fval]=fsolve(@specialeqex,0); R_c=[R_c h1]; 副→specialeqex.m function f = specialeqex(R) load d:/I_t t=1996; A_t=[2000.4;2000.4;2100.3]; A_ref=[2400;2400;2400]; h_0=[30;30;30]; P=1; d=200; H=4; lnda=0.023; gamma=1; for i=1:length(A_ref) if A_ref(i) > A_t(i) % 沖蝕 A_t0=[]; for w=1:12 A= ((I_t(w,2)*exp(-lnda*(1966-(1953+w)))... +I_t(w+1,2)*exp(-lnda*(1966-(1954+w))))/2); A_t0=[A_t0 A]; end f_c=0; for j=13:42 f_1 = ((1-P*gamma*(1-exp(-R/H)))*I_t(j,2)*exp(-(P*R/d+lnda)*(t-(1953+j)))... + (1-P*gamma*(1-exp(-R/H)))*I_t(j+1,2)*exp(-(P*R/d+lnda)*(t-(1954+j))))/2; f_c=f_c+f_1; end f = A_t(i)-(sum(A_t0))*exp(-(P*R/d+lnda)*(1996-1966))-f_c; end end ps.因I_t的檔案無法PO所以願意幫我解答問題的大大，請留言我會馬上跟您聯絡!! 感激不盡 <(__)> ~~~ T__T (泣) -- 愛，是心中輕盈的詩　是幸福，是寂寞 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.190.27