小弟不是學光學的，所以想請各位大俠指點��！謝謝啦 >Hd!o"I
 
就是我用ansys計算出了鏡面的面型的數(shù)據(jù)，怎樣可以得到zernike多項式的系數(shù)，然后用zemax各階得到像差！謝謝啦！ qGzF@p(p8  

可以用matlab編程，用zernike多項式進行波面擬合，求出zernike多項式的系數(shù)，擬合的算法有很多種，最簡單的是最小二乘法，你可以查下相關(guān)資料，挺簡單的

澤尼克多項式的前9項對應象差的

非常感謝啊，我手上也有zernike多項式的擬合的源程序，也不知道對不對，不怎么會有 X//=OpS`  
function z = zernfun(n,m,r,theta,nflag) 4^ZbT  
%ZERNFUN Zernike functions of order N and frequency M on the unit circle. 1
CS\1[E  
%   Z = ZERNFUN(N,M,R,THETA) returns the Zernike functions of order N $WsyAUl  
%   and angular frequency M, evaluated at positions (R,THETA) on the 4@Bl 1b[<  
%   unit circle.  N is a vector of positive integers (including 0), and $/Llzpvny  
%   M is a vector with the same number of elements as N.  Each element QF$s([  
%   k of M must be a positive integer, with possible values M(k) = -N(k) _*wlK;`  
%   to +N(k) in steps of 2.  R is a vector of numbers between 0 and 1, n -x
Caq  
%   and THETA is a vector of angles.  R and THETA must have the same L!Gpk)}[i  
%   length.  The output Z is a matrix with one column for every (N,M) ziv*4  
%   pair, and one row for every (R,THETA) pair.  bDq<]h_7  
% Yd<9Y\W%?  
%   Z = ZERNFUN(N,M,R,THETA,'norm') returns the normalized Zernike ,S&p\(r.  
%   functions.  The normalization factor sqrt((2-delta(m,0))*(n+1)/pi), ,be$~7qS  
%   with delta(m,0) the Kronecker delta, is chosen so that the integral  $kxu-  
%   of (r * [Znm(r,theta)]^2) over the unit circle (from r=0 to r=1, g2YE^EKU~  
%   and theta=0 to theta=2*pi) is unity.  For the non-normalized