风向使用von Mises拟合的时候，使用最小二乘法估计初值的影响比较大，Carta在2008年的论文中提出了计算初值的方法，请问si和ci具体怎么算的，使用此方法拟合风向的matlab代码有吗