程序运行过后的误差精度达不到，怎么修改这个程序。
本帖最后由 cppwylb 于 2012-4-23 17:17 编辑下面是相位轮廓术的程序，模拟半球，可是最后的误差分析达不到要求的精度，请高人指导一下，单位告诉的相位差改为0.2pi要怎么弄
close all;clear all;clc
N=512;x=1:N;
y=1:N;=meshgrid(x,y);R=240;r=90;
z=zeros(N);
for i=1:N
for j=1:N
if(sqrt((i-N/2)^2+(j-N/2)^2)<sqrt(R*R-r*r))
z(i,j)=.25*((sqrt(R*R-(i-N/2).*(i-N/2)-(j-N/2).*(j-N/2))-r));
end
end
end
z=0.02*z;
%figure,mesh(z);
% =meshgrid(0:0.01:4,0:0.01:4);
% Z = PEAKS(401);
% Z = Z/max(max(Z));
% z = 1.5*Z;
%h = max(max(Z))%
figure,mesh(z);xlabel('x'),ylabel('y'),zlabel('z')
f0=0.01;
noise=0;
i01=0.2+0.5*cos(2*pi*x*f0)+noise.*rand(size(z));
i02=0.5+0.5*cos(2*pi/4+2*pi*x*f0)+noise.*rand(size(z));
i03=0.5+0.5*cos(4*pi/4+2*pi*x*f0)+noise.*rand(size(z));
i04=0.5+0.5*cos(6*pi/4+2*pi*x*f0)+noise.*rand(size(z));
figure,gray256(i01);
figure,gray256(i02);
figure,gray256(i03);
figure,gray256(i04);
q2=-atan2((i01*sin(0*pi/4)+i02*sin(2*pi/4)+i03*sin(4*pi/4)+i04*sin(6*pi/4)),(i01*cos(0*pi/4)+i02*cos(2*pi/4)+i03*cos(4*pi/4)+i04*cos(6*pi/4)));
figure,mesh(q2);xlabel('x'),ylabel('y'),zlabel('z')
temp=ones(size(q2));
q21=LingXingUnwrap2(q2,temp,200,200);
figure,mesh(q21);
i1=0.5+0.5*cos(2*pi*x*f0+2*pi*z)+noise.*rand(size(z));
figure,gray256(i1)
i2=0.5+0.5*cos(2*pi/4+2*pi*x*f0+2*pi*z)+noise.*rand(size(z));
figure,gray256(i2)
i3=0.5+0.5*cos(pi+2*pi*x*f0+2*pi*z)+noise.*rand(size(z));
figure,gray256(i3)
i4=0.5+0.5*cos(6*pi/4+2*pi*x*f0+2*pi*z)+noise.*rand(size(z));
figure,gray256(i4)
q1=-atan2((i1*sin(0*pi/4)+i2*sin(2*pi/4)+i3*sin(pi)+i4*sin(6*pi/4)),(i1*cos(0*pi/4)+i2*cos(2*pi/4)+i3*cos(pi)+i4*cos(6*pi/4)));
figure,mesh(q1);xlabel('x'),ylabel('y'),zlabel('z')
q11=LingXingUnwrap2(q1,temp,200,200);
figure,mesh(q11);
b=(q11-q21)./(2*pi);
wu=(q11-q21)./(2*pi)-z;
err2=max(max(wu))
figure,mesh(q21);xlabel('x'),ylabel('y'),zlabel('z')
figure,mesh(q11);xlabel('x'),ylabel('y'),zlabel('z')
figure,mesh(b);xlabel('x'),ylabel('y'),zlabel('z')%%%恢复的物体
figure,mesh(wu);xlabel('x'),ylabel('y'),zlabel('z')%%%测量误差
补充，把单位高度引起的相位差改成0.2pi要怎么改啊，
页:
[1]