※ 引述《sikao (si)》之銘言：
: 版上前輩大家好 再請教一個問題
: 我用plot指令作圖 但圖就跑不出來
: 想請問一下是甚麼原因
: 我寫的script如下
: for t=0:0.05:5
: T=170-22*t
: if T>=120
: G=(3.98*10^7)*exp(-6270/(8.314*(T-30)))*exp(-2.55*10^5/((T+273)*(200-T)))
: else
: G=(4.81*10^11)*exp(-6270/(8.314*(T-30)))*exp(-5.51*10^5/((T+273)*(200-T)))
: end
: plot(t,G);
: axis([0,5,10^-5,5]);
: xlabel('Time');
: ylabel('G (min^-1)')
: end
: 不知道是不是要調整y軸的範圍
: 麻煩大家非常謝謝
改寫如下
function activiteng
clc
format long
t=0:0.05:5;
T = 170-22.*t;
if T >= 120
G
=(3.98.*10^7).*exp(-6270./(8.314.*(T-30))).*exp(-2.55.*10^5./((T+273).*(200-T)))
else
G
=(4.81.*10^11).*exp(-6270./(8.314.*(T-30))).*exp(-5.51.*10^5./((T+273).*(200-T)))
end
plot(t,G);
axis([0,5,10^-5,5]);
xlabel('Time');
ylabel('G (min^-1)')
這樣應該可行