采用matlab编程，其主函数如下，可以模拟各阶的zernike多项式： acd[rjeT
%Display the Zernike function Z(n=5,m=1) UcRP/LR%C
clc F.9} jd{
clear {{G)Ry*pb
a=5;%%%%%%%%%%Z的阶数下标 ~k"+5bHa*
b=1;%%%%%%%%%%Z的阶数的上标 0=Z[6Q@ :
x = -1:0.01:1; z0z@LA4k6@
[X,Y] = meshgrid(x,x); x~5uc $
[theta,r] = cart2pol(X,Y); As:O|!F
idx = r<=1; iq#{*:1
z = nan(size(X)); D6"=2XR4n
z(idx) = zernfun(a,b,r(idx),theta(idx)); e4z`:%vy
figure(1) >)>f~ >
pcolor(x,x,z), shading interp &f*o rM:
axis square, colorbar [Vd$FDki
xlabel('X'); k<gH*=uXY'
ylabel('Y'); b1& t k~ D
title(['Zernike function Z^a_b','(r,\theta)']) -ng1RA>
figure(2) iRQ!J1SGcG
mesh(x,x,z) l_I)d7
xlabel('X'); d "wA"*8~y
ylabel('Y'); M0V<Ay\%O
title(['Zernike function Z^a_b','(r,\theta)'])
