Matlabの再帰的ループは、反復できません

Question

皆さん。私はループに問題があります。それは一度だけ実行され、反復されません。

私は宇宙船の表面温度を時間の関数として見出そうとしています。輻射熱の流れは現在の温度に依存するため、ループの次の温度は1秒前の温度に依存するはずです。これはスクリプト全体ではありませんが、ループの何が問題なのか誰にでも教えてください。

ありがとうございます！

% Starting with given surface and interior temperatures, find each new 
% temperature with each heat flux every second through one orbit. 

timeline = 1:1:time1; 
Ts = zeros(size(timeline)); 
Ts(1) = 240; % initial surface temperature, Kelvins` 

for n = 2:length(time1) 
    Ts(n) = (solarflux + IRs - (sigma*ems_mli*surfarea*(Ts(n-1)^4)))/(m*c) + Ts(n-1); 
end 

%% Plotting 

figure(1) 
plot(timeline,Ts) 
xlabel('seconds'); 
ylabel('surface temperature (Kelvins)'); 

Answer 1

time1の値を設定した後、以下のコードを試してみてください。 いくつかのエラーが修正されました：n = 2：長さ（time1）はn = 2：length（timeline）でなければなりません.time1はスカラーで、sizeの代わりにlength（timeline）そのタイムラインと同じサイズのベクトル。

timeline = 1:1:time1;        % your time array 
Ts = zeros(1,length(timeline));     % preallocate memory for Ts 
Ts(1) = 240;          % initial surface temperature, Kelvins 

for n = 2:length(timeline)       % loop through 2 to end of timeline 
Ts(n) = (solarflux + IRs - (sigma*ems_mli*surfarea*(Ts(n-1)^4)))/(m*c) + Ts(n-1); 
end 

figure         % create a figure 
plot(timeline,Ts)      % plot timeline vs Ts 
xlabel('seconds');      % X-Axis label 
ylabel('surface temperature (Kelvins)'); % Y-Axis label 

Answer 2

おかげで、みんな。私は、以前の「タイムライン」のエラーを発見したあなたの姿をプロットする。