在Matlab中,三角函数与角弧度,看到例如[描述][1]的“罪”功能:>输入角弧度,指定为一个标量、向量,矩阵,或>多维数组。所以你要么保持代码但转换10°角弧度(改变α= 0.174533),或者您可以使用“信德”功能,该节点也可以计算正弦但需要角度(度)[2]作为输入。然后,你也有一个10倍的错误在你的一美元的定义:$ 10 ^{-16}$应该注意“1 e-16”。修复第一个图(速度剖面):[![在这里输入图像描述][3]][3]第二个看起来仍然奇怪……你也有一个错误的定义q $ ^{\补}_ {geo} $:单位是千瓦,所以它应该被转换成“q_geo = 0.05;”。这个修正方程的第一行:如果你把其余的评论,它再现了线性,虚线,温度为-7.619°C for z = 0。所以这个问题必须在第二部分,虽然我不能找出。[1]https://www.mathworks.com/help/matlab/ref/sin.html [2]: https://www.mathworks.com/help/matlab/ref/sind.html btg_ljb-3 [3]: https://i.stack.imgur.com/pJdLH.png