使用MATLAB计算球谐函数Y(theta，fai，L，M)并将其保存为矩阵Y。矩阵Y是一个3维矩阵，其尺寸为180，360，（Ln_max+1）^2。此外，矩阵NNN保存了Y(L，M)在Y矩阵中第三维的位置。