比例-谐振控制器,抑制谐波，有效的将谐波滤除

Harmonics_and_z分数_GEE_code Google Earth Engine代码可根据MODIS时间序列数据计算z得分和谐波分析。

此函数生成 L 阶和 M 阶的球谐函数基函数。 语法：[Ymn,THETA,PHI,X,Y,Z]=spharm4(L,M,RES,PLOT_FLAG); 输入： L - 球谐度数，[1x1] M - 球面谐波阶次，[1x1] RES - 要使用的点数的向量 [#Theta x #Phi 点]、[1x2] 或 [2x1] PLOT_FLAG - 生成球谐曲面图形的二进制标志 (DEFAULT=1) 输出： Ymn - 球谐坐标，[RES(1) x RES(2)] THETA - 圆周坐标，[RES(1) x RES(2)] PHI - 纬度坐标，[RES(1) x RES(2)] X,Y,Z - 大小、平方、球谐曲面点的笛卡尔坐标，[RES(1) x RES(2)] 注意：保持 THETA 和 PHI 的各种用法是非常重要的。 对于此函数，THETA 是方位角/经度/圆周坐标，定

这组小而基本的例程是由于我需要在地磁应用的背景下创建和操作真实的球谐系数。 这些系数被保存为向量，因为 Matlab 中的向量操作非常发达。 该工具箱支持访问对应于特定阶数 l 和阶数 m 的系数（其中负 m 对应于 sin(m*phi)，正 m 对应于 cos(m*phi) 项）； 将向量映射到矩阵并返回； 使用施密特半归一化球谐函数和基本绘图映射到网格。 由独立系数集的几个“部分”组成的向量也在某种程度上得到支持（在我的情况下，它们对应于地球不同层的参数化）。 还包含一组函数，用于将其与 FJ Simons 公开提供的 DOTM 库集成。 为了使用 SHRotateVec（球谐系数的旋转），请安装他的库http://geoweb.princeton.edu/people/simons/software.html 受 Fortran SHTOOLS 软件启发的名称http://w