来源：小树洞谈光伏支架
各位好久不见！之前小树洞都是写一些结构或者流体等方面的计算，这次换一个话题，和大家聊聊“跟踪算法”。
最近看到一篇使用CAD方式得到跟踪器最优角度的文章，这个方法来得到一次性的结果尚且可以，但是我们知道跟踪器是一种实时追踪太阳的装置，对于各种角度如何使用理论的方式一劳永逸的得到最优值呢？
1
带倾角平单轴有何不同
熟悉平单轴跟踪算法的朋友一定对计算理论最优角度非常熟悉了，无非是使用不同的经验公式得到太阳赤纬角，时角。再推导出太阳高度角和方位角，利用这两个参数得到平单轴在不同地点不同时间的最优角度。
对于普通平单轴来说，因为组件与转动轴是相同的平面，计算算法往往使用2维的坐标转换即可完成，也就是使用平面向量。换句话说，主轴45度，组件就是45度。
但是当组件带有一定的角度，也就是带倾角平单轴时。计算最优角度则需要使用3维坐标转换公式，也就是空间向量。那么问题来了，如果主轴45度，组件还是45度吗？
2
几何原理的解释
为了方便电机控制平单轴旋转，跟踪器倾角传感器往往都是紧贴旋转轴。主轴旋转角度变化，对于组件平面来说角度也必然会跟着改变。
不管是双轴也好，平单轴也好，最优角度其实就是求组件法向量与太阳和组件中心连线的向量之间夹角最小。假设地面为参考坐标系，那么问题就转化为求解组件与地面的夹角。
只不过我们很容易就能发现，带倾角平单轴的组件角度与旋转主轴的角度变化并不是线性一致的。以带15度倾角的平单轴为例。
旋转主轴0度时，组件与地面夹角为15度。
旋转主轴45度时，组件与地面夹角为46.92度。
3
三维跟踪角度算法的原理
小树洞在这里介绍的3维坐标转换借鉴了目前可视化行业的计算机3D图形学，说白了就是将构建好的3维的物体展示到2维的显示器上。接下来说一下如何用理论的方法来计算组件与地面的夹角。
假设如下图所示：
-
组件的原法向量为OP1(x1,y1,z1)
-
地面坐标系(X,Y,Z)
-
X轴为跟踪器的旋转主轴
-
点Q，N，M分别为向量OP1到平面XZ，YZ和XY上的投影点
接下来的问题就是求当主轴转动一定的角度时，向量OP2与向量OZ的夹角。
从上面的图可以看出来，主轴绕X轴旋转时，组件法向量OP也是绕着X轴旋转，旋转360度则扫出一个圆锥面。而OP投影到平面YZ上的向量ON则长度不变，以O点为圆心旋转。
假设原先ON与Y轴夹角为Alpha，向量OP1(x1,y1,z1)，当ON转动Theta角后，OP1坐标变为OP2(x2,y2,z2)。
据此我们可以列出如下左右两对方程组：
将右边的那对(y2,z2)方程组展开得到：
再将左边那对(y1,z1)方程组带入上面的方程组得到与初始角度Alpha无关的(y1,z1,y2,z2)方程组：
再将一直没有露面的x1=x2带入，这样就得到了一个6元非齐次线性方程组：
利用行列式“OP2=常数项xOP1”整理下这个方程组得到一个旋转矩阵：
至此，关键的计算已经完成了。
4
举个栗子
下面举一个具体的例子，当跟踪器旋转主轴从放平0度转动到45度时，求带倾角15度的组件相对于地面角度。
原组件单位法向量OP1为：
-
主轴转动45度，则旋转矩阵为：
-
利用行列式OP2=旋转矩阵xOP1得到最终的组件法向量：
-
地面的单位法向量为：
-
得到两个向量（地面和OP2），可求的其夹角为46.92度。
全过程如下：
至此，带倾角平单轴真实角度的计算结束。
