在Matlab中,三角函数与角弧度,例如描述的罪功能:
输入角弧度,指定为一个标量、向量,矩阵,或多维数组。
所以你要么保持代码但转换10°角弧度(改变α= 0.174533),或者您可以使用信德函数,它也计算正弦但需要角度在度作为输入。
然后,你也有一个10倍误差的定义一个美元:10美元^ {-16}$应该注意的1 e-16。修复第一个图(速度剖面):
第二个看起来仍然奇怪……你也有一个错误的定义q $ ^{\补}_ {geo} $:单位是千瓦,所以它应该转换为q_geo = 0.05;。这个修正方程的第一行:如果你把其余的评论,它再现了线性,虚线,温度为-7.619°Cz美元m = 0。所以这个问题必须在第二部分,虽然我不能找出。