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