Zernike多项式拟合横向剪切干涉

clear all; num=1024; x=[1:num]; y=[1:num]; [X1,Y1]=meshgrid(x,y); Radius=300; CenterX=400; CenterY=500; shear_value=100; X=(X1-CenterX)/Radius; Y=(Y1-CenterY)/Radius; a1=0.1; %a2=10; a2=0.1; %a3=2; a3=0.2; a4=0.7; a5=0.9; a6=0.11; a7=5; %a7=0.5; a8=0.15; %a8=15; Z1=1; Z2=2.*X; Z3=2.*Y; Z4=3^0.5.*[2*(X.^2+Y.^2)-1]; Z5=2*6^0.5.*X.*Y; Z6=6^0.5.*(X.^2-Y.^2); Z7=8^0.5.*[3.*(X.^2+Y.^2).*Y-2.*Y]; Z8=8^0.5.*[3.*(X.^2+Y.^2).*X-2.*X]; lamda=532; opd=a1.*Z1+a2.*Z2+a3.*Z3+a4.*Z4+a5.*Z5+a6.*Z6+a7.*Z7+a8.*Z8; phi=2*pi/lamda*opd; r=X.^2+Y.^2; flag1=(r<=1); flag11=(r>1); phi1=zeros(num,num); phi1(flag1)=phi(flag1); % y3=uint8(255/2*(1+0.8*cos(phi1))); % figure,imagesc(y3); % colormap gray % axis square phi2=zeros(num,num); phi2(shear_value+1:num,:)=phi1(1:num-shear_value,:); flag2=zeros(num,num); flag2=logical(flag2); flag2(shear_value+1:num,:)=flag1(1:num-shear_value,:); flag22(shear_value+1:num,:)=flag11(1:num-shear_value,:); % r=(X).^2+(Y-shear_value-1).^2; % flag2=(r<=1); % flag22=(r>1); phi_1=phi1-phi2; flag=(flag1 & flag2); phi_11=zeros(num,num); phi_11(flag)=phi_1(flag); flag=(flag11 & flag22); phi_11(flag)=pi/2; phi2=zeros(num,num); phi2(:,shear_value+1:num)=phi1(:,1:num-shear_value); flag2=zeros(num,num); flag2=logical(flag2); flag2(:,shear_value+1:num)=flag1(:,1:num-shear_value); flag22(:,shear_value+1:num)=flag11(:,1:num-shear_value); % r=(X-shear_value-1).^2+(Y).^2; % flag2=(r<=1); % flag22=(r>1); phi_2=phi1-phi2; flag=(flag1 & flag2); phi_22=zeros(num,num); phi_22(flag)=phi_2(flag); flag=(flag11 & flag22); phi_22(flag)=pi/2; % % figure,imagesc(phi_22);colormap gray % figure,imagesc(flag);colormap gray y1=uint8(255/2*(1+0.8*sin(phi_11))); y2=uint8(255/2*(1+0.8*sin(phi_22))); y3=uint8(255/2*(1+0.8*cos(phi1))); figure,imagesc(y1); colormap gray axis square %save y1 y1;%存储下 figure,imagesc(y2); colormap gray axis square %save y2 y2; figure,imagesc(y3); colormap gray axis square %save y3 y3; figure mesh(phi1);