非常感谢啊,我手上也有zernike多项式的拟合的源程序,也不知道对不对,不怎么会有 gG#M-2P
function z = zernfun(n,m,r,theta,nflag) f\Qi()
%ZERNFUN Zernike functions of order N and frequency M on the unit circle. z 6p.{M
% Z = ZERNFUN(N,M,R,THETA) returns the Zernike functions of order N rx%lL
% and angular frequency M, evaluated at positions (R,THETA) on the FN G]
% unit circle. N is a vector of positive integers (including 0), and Z/%>/
% M is a vector with the same number of elements as N. Each element &n['#7 <(!
% k of M must be a positive integer, with possible values M(k) = -N(k) lLnD%*03
% to +N(k) in steps of 2. R is a vector of numbers between 0 and 1, IF<jq\M
% and THETA is a vector of angles. R and THETA must have the same H=*;3gM,'
% length. The output Z is a matrix with one column for every (N,M) iZ&