script4

Contents

Zernike polynomials

r = linspace(0,1,201);
z = zern(r);

for k=2:7
    subplot(3,2,k-1);
    plot(r,z(:,k),'k','LineWidth',2);
    axis([0 1 -1 1]);
    grid;
    ylabel(sprintf('order %d',2*k-2));
    hold on;
    plot([0 1],[0 0],'k');
    hold off;
end

Chebyshev polynomials

r = linspace(0,1,201);
z = cheb(r);

for k=1:6
    subplot(3,2,k);
    plot(r,z(:,2*k+1),'k','LineWidth',2);
    grid;
    ylabel(sprintf('order %d',2*k));
    hold on;
    plot([0 1],[0 0],'k');
    hold off;
end

Published with MATLAB® 7.6