楼主确定公式没问题吗?我画出的图和你给出的差别比较大:
有几点不太确定:
(1)上面的公式中arg函数理解为复数的辐角(argument of a complex number),MATLAB中使用angle函数计算。
(2)上面给出-90°≤α≤90°,但α由tg2α计算得到,由于反正切函数的值域为-90°-90°,再除以2,怎样得到超出-45°≤α≤45°范围之外的结果?根据相邻值的连续性吗?
(3)考虑一个简单的情况:取φ=0,则Exi=0,Eyi=1,ζ=0。进一步计算得Exr=0,γ=90°,所以tan2α=0,sin2β=0这时为什么不取α=0而是如图取-90°呢?
因提问已过多日,不知楼主是否还需要,故暂不贴出代码,如有需要再贴。