Matlab复合'if'语句未给出预期结果[关闭]

Question

在为本文http://adsabs.harvard.edu/abs/1992PASJ...44..493K编写方程（5）的程序时，

我想我犯了一些错误（下）。我不知道如何在等式中包含“或”。

    if alpha_sgp<alpha<=(alpha_sgp-180)
        if p<=180-t
          P=p+t;
        else %p>180-t
          P=p+t-180;
        end
    elseif alpha<=alpha_sgp | (alpha_sgp-180)<alpha   
        if p>=t
          P=p-t;
        else %p<t
          P=p-t+180;
        end
    end

任何帮助将不胜感激。

Answer 1

我觉得你的说法有问题：

if alpha_sgp<alpha<=(alpha_sgp-180)

虽然这在 Matlab 中的语法是正确的，但它并没有按照你的想法做。

简单示例：

a=1;b=2;c=3;
disp(a<b<c)

将返回 1 (true)。但是

disp(a<c<b)

也将返回1。这与操作顺序有关，并且在将结果与c 进行比较之前，将评估a<b 类型的子表达式（到1）。因此

a<b<c

相当于

(a<b) < c

相当于

If ( (a < b) & (1 < c)) | ( (a >= b ) & (0 < c) )

你真正想要的第一行是

(alpha_sgp<alpha) & (alpha<=(alpha_sgp-180))

除了确保正确地将度数转换为弧度等。我建议您坚持使用弧度 - 使用内置的 pi - 而不是以度数工作，并且必须使用像 @ 这样的函数987654334@ 而不是 sin。

另外 - 你对“不知道如何包含for”提出了一个奇怪的评论 - 当没有for 时，等式中只有一个or。这可能是您的错字吗？

最后 - 请注意，第二组条件基本上与第一组相反 - 当(alpha_sgp < alpha) & (alpha <= (alpha_sgp-180)) 为假时，反之亦然。所以在实现你的公式时，你可以用 else 替换整个 elseif 行。

如果您需要进一步说明，请随时发表评论。