返回答案

在Matlab中，三角函数工作在弧度角度，参见示例描述的罪功能:

输入角度，以弧度为单位，指定为标量、矢量、矩阵或多维数组。

所以你要么保持你的代码，但转换你的10°角在弧度(更改为Alpha = 0.174533)，或使用信德函数，它也计算正弦，但取角度在度作为输入。

那么，你的定义中也有10倍的误差一个美元：10美元^ {-16}$应该注意1 e-16。这将修复第一个图形(速度剖面):

第二个看起来还是很奇怪……的定义也有错误q $ ^{\补}_ {geo} $:单位为mW，应换算为q_geo = 0.05;。这修复了方程的第一行:如果您将其余部分放在注释中，它将很好地再现线性虚线，温度为-7.619°Cz= 0。所以问题一定在第二部分，虽然我不知道在哪里。

在Matlab中，三角函数工作在弧度角度，参见示例描述的罪功能:

输入角度，以弧度为单位，指定为标量、矢量、矩阵或多维数组。

所以你要么保持你的代码，但转换你的10°角在弧度(更改为Alpha = 0.174533)，或使用信德函数，它也计算正弦，但取角度在度作为输入。

那么，你的定义中也有10倍的误差一个美元：10美元^ {-16}$应该注意1 e-16。这将修复第一个图形(速度剖面):

第二个看起来还是很奇怪……的定义也有错误q $ ^{\补}_ {geo} $:单位为mW，应换算为q_geo = 0.05;。这修复了方程的第一行:如果您将其余部分放在注释中，它将很好地再现线性虚线，温度为-7.619°Cz美元= 0米。所以问题一定在第二部分，虽然我不知道在哪里。

在Matlab中，三角函数工作在弧度角度，参见示例描述的罪功能:

输入角度，以弧度为单位，指定为标量、矢量、矩阵或多维数组。

所以你要么保持你的代码，但转换你的10°角在弧度(更改为Alpha = 0.174533)，或使用信德函数，它也计算正弦，但取角度在度作为输入。

那么，你的定义中也有10倍的误差一个美元：10美元^ {-16}$应该注意1 e-16。这将修复第一个图形(速度剖面):

第二个看起来还是很奇怪……的定义也有错误q $ ^{\补}_ {geo} $:单位为mW，应换算为q_geo = 0.05;但它不能解决所有问题。

在Matlab中，三角函数工作在弧度角度，参见示例描述的罪功能:

输入角度，以弧度为单位，指定为标量、矢量、矩阵或多维数组。

所以你要么保持你的代码，但转换你的10°角在弧度(更改为Alpha = 0.174533)，或使用信德函数，它也计算正弦，但取角度在度作为输入。

那么，你的定义中也有10倍的误差一个美元：10美元^ {-16}$应该注意1 e-16。这将修复第一个图形(速度剖面):

第二个看起来还是很奇怪……的定义也有错误q $ ^{\补}_ {geo} $:单位为mW，应换算为q_geo = 0.05;。这修复了等式的第一行:如果你把其余的放在注释中，它在z = 0的温度为-7.619°C时，可以很好地再现线性虚线。所以问题一定在第二部分，虽然我不知道在哪里。

在Matlab中，三角函数工作在弧度角度，参见示例描述的罪功能:

输入角度，以弧度为单位，指定为标量、矢量、矩阵或多维数组。

所以你要么保持你的代码，但转换你的10°角在弧度(更改为Alpha = 0.174533)，或使用信德函数，它也计算正弦，但取角度在度作为输入。

那么，你的定义中也有10倍的误差一个美元：10美元^ {-16}$应该注意1 e-16。这就固定了第一个图形(速度剖面)。

第二个看起来还是很奇怪……的定义也有错误q $ ^{\补}_ {geo} $:单位为mW，应换算为q_geo = 0.05;但它并不能解决所有问题。

在Matlab中，三角函数工作在弧度角度，参见示例描述的罪功能:

输入角度，以弧度为单位，指定为标量、矢量、矩阵或多维数组。

所以你要么保持你的代码，但转换你的10°角在弧度(更改为Alpha = 0.174533)，或使用信德函数，它也计算正弦，但取角度在度作为输入。

那么，你的定义中也有10倍的误差一个美元：10美元^ {-16}$应该注意1 e-16。这就固定了第一个图形(速度剖面)：

第二个看起来还是很奇怪……的定义也有错误q $ ^{\补}_ {geo} $:单位为mW，应换算为q_geo = 0.05;但它并不能解决所有问题。