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
呼啸庄主
- 粉丝: 83
- 资源: 4696
评论1