把除号/改成点除./,即
t=1:1:12
Y1=3.95062*10^7-(3.92993*10^7) ./ exp(0.011122*t)
这样就是按分量计算而不是矩阵运算
解决方法如下:需要的是t和l0的元素对应计算,即矩阵乘法,那么需要在乘除运算符之前加上点号(./ 或 .*),所以改正如下:
123t = 9:1:15l0 = 3 ./ (tan(asin(0.4620 + 0.7712 * 0.6879 * cos(pi / 12 * (t - 12.24)))))plot(t,l0)