在Matlab中,三角函数与角弧度,看到例如[描述][1]的“罪”功能:>输入角弧度,指定为一个标量、向量,矩阵,或>多维数组。所以你要么保持代码但转换10°角弧度(改变α= 0.174533),或者您可以使用“信德”功能,该节点也可以计算正弦但需要角度(度)[2]作为输入。然后,你也有一个10倍的错误在你的一美元的定义:$ 10 ^{-16}$应该注意“1 e-16”。修复第一个图(速度剖面)。第二个看起来仍然奇怪……你也有一个错误的定义q $ ^{\补}_ {geo} $:单位是千瓦,所以它应该转化为q_geo = 0.05;但它不会解决所有问题。[1]https://www.mathworks.com/help/matlab/ref/sin.html [2]: https://www.mathworks.com/help/matlab/ref/sind.html # btg_ljb-3