for i=1:1:100 T(i)=296+i; Q(T)=a+b*T+c*T.^2+d*T.^3; S(T)=(q/Q)*(t/T)*exp(-h*c*E/k*(1/T-1/t))*(1-exp(-h*c*v/(k*T)))*(1-exp(-h*c*v/(k*t)))^-1;endplot(T,S)
从给出的代码来看,出现《错误使用 / 矩阵维度必须一致。》的错误警告,其原因是
Q(T)=a+b*T+c*T.^2+d*T.^3; 中的T 无法用 T(i) 替代。因为 T 是一组向量值,而 T(i) 是一个数值。
所以,你可以试着按下列代码改一下,即
for i=1:1:100
T(i)=296+i;
Q(i)=a+b*T(i)+c*T(i).^2+d*T(i).^3;
S(i)=(q/Q(i))*(t/T(i))*exp(-h*c*E/k*(1/T(i)-1/t))*(1-exp(-h*c*v/(k*T(i))))*(1-exp(-h*c*v/(k*t)))^-1;
end
plot(T,S)