imaging_光学成像仿真_

function Fi = getAberratedFi(sys, maxN, c) x = linspace(-1, 1, sys.re); y = x; [X, Y] = meshgrid(x, y); [Theta, Rho] = cart2pol(X, Y); ps = find(Rho <= 1); idx = 0; Fi = zeros(sys.re); for n = 0 : maxN for m = -n : 2 : n idx = idx + 1; Fi(ps) = Fi(ps) + c(idx) * zernfun(n, m, Rho(ps), Theta(ps),'norm'); end end end