MATLAB最小二乘解,最小二乘算法MATLAB代码实现
对于x=zx0:m
flt0(x，uint 16(B3 * x):uint 16(B2 a2 * x))=1；%(-1)*a*x？
结束了
对于x=1:zx0
fl1(x，uint 16(B1 a 1 * x):uint 16(B4 a 4 * x))=1；
结束了
对于y=zy0:n
FLT 2(uint 16((y-B4)/a4):uint 16((y-B2)/a2)，y)=1；
结束了
对于y=1:zy0
FLT 3(uint 16((y-B1)/a1):uint 16((y-B3)/a3)，y)=1；
结束了
百分之九十处理滤波器的重合点
f02=FLT 0 FLT 2；
f03=FLT 0 FLT 3；
F12=FLT 1 FLT 2；
f13=FLT 1 FLT 3；
f1=FLT 0 FLT 1；
对于x=1:m
对于y=1:n
如果F2(x，y)==2
fl2(x，y)=0；
结束了
如果F3(x，y)==2；
flt3(x，y)=0；
结束了
如果f12(x，y)==2
fl2(x，y)=0；
结束了
如果f13(x，y)==2；
flt3(x，y)=0；
结束了
如果f1(x，y)==2；
fl1(x，y)=0；
结束了
结束了
结束了
J=flt0 .* J1 FLT 1 .*J1 flt2 .*J2第三飞行技术队.*JBOY3乐队；
图(6),网格（双(j))；
% % % % % % % % % % % % j
% - 基本方法解包裹-
[m n]=大小(j)和：
k=double(零(m，n))；
对于j=1:n
对于h=2:m
如果(J(h，j)-J(h-1，j))=pi
k(h，j)=k(h-1，j)-1；
elseif abs(J(h，j)-J(h-1，j)) pi
k(h，j)=k(h-1，j)；
elseif (J(h，j)-J(h-1，j)) (-pi)
k(h，j)=k(h-1，j)1；
结束了
结束了
结束了
对于h=1:m
对于p=2:n
如果(J(h，p)-J(h，p-1))=pi
k(h，p)=k(h，p-1)-1；
elseif abs(J(h，p)-J(h，p-1)) pi
k(高，p)=k(高，p-1)；
elseif (J(h，p)-J(h，p-1)) (-pi)
k(h，p)=k(h，p-1)1；
结束了
结束了
结束了
x=J2 ** k；
图(7)、网格（双(x))；
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %。最小二乘解包裹百分之九十
% % % % % % % % % % % % % %。边界处理
% F=spalloc(m*n，m*n，5 * m * n)；
w=1时的百分比：m*n
% x=ceil(w/n)-1；y=rem(w，n)；% w=i*n j - i j -
%如果y==0
% y=n
%结束
% if (x==0) (y~=1 y~=n)
% F(w，w)==3；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% F(w，(x1)* n y)=1；
% % F(w，(x-1)* n y)=1；
% F(w，x * n y 1)=1；
% F(w，x * n y-1)=1；
%结束
% if (x==m-1) (y~=1 y~=n)
% F(w，w)==3；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% % F(w，(x1)* n y)=1；
% F(w，(x-1)* n y=1；
% F(w，x * n y 1)=1；
% F(w，x * n y-1)=1；
%结束
% if (y==1) (x~=0 x~=m-1)
% F(w，w)==3；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% F(w，(x1)* n y)=1；
% F(w，(x-1)* n y=1；
% F(w，x * n y 1)=1；
% % F(w，x * n y-1)=1；
%结束
% if (y==n) (x~=0 x~=m-1)
% F(w，w)==3；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% F(w，(x1)* n y)=1；
% F(w，(x-1)* n y=1；
% % F(w，x * n y 1)=1；
% F(w，x * n y-1)=1；
%结束
% if (x==0) (y==1) % y~=n
% F(w，w)==2；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% F(w，(x1)* n y)=1；
% % F(w，(x-1)* n y)=1；
% F(w，x * n y 1)=1；
% % F(w，x * n y-1)=1；
%结束
% if (x==m-1) (y==1)% y~=n
% F(w，w)==2；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% % F(w，(x1)* n y)=1；
% F(w，(x-1)* n y=1；
% F(w，x * n y 1)=1；
% % F(w，x * n y-1)=1；
%结束
% if (x==m-1) (y==n)% y~=n)
% F(w，w)==2；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% % F(w，(x1)* n y)=1；
% F(w，(x-1)* n y=1；
% % F(w，x * n y 1)=1；
% F(w，x * n y-1)=1；
%结束
% if (x==0) (y==n)% y~=n
% F(w，w)==2；% F(w，(x-1)* n y=f(w，(i 1)*n j)
% F(w，(x1)* n y)=1；
% % F(w，(x-1)* n y)=1；
% % F(w，x * n y 1)=1；
% F(w，x * n y-1)=1；
%结束
% if(0x m-1)(y1y n)
% F(w，w)==4；
% F(w，(x1)* n y)=1；
% F(w，(x-1)* n y=1；
% F(w，x * n y 1)=1；
% F(w，x * n y-1)=1；
%结束
%结束
%
% %% -建立泊松方程
% dx=零(m，n)；
% dy=零(m，n)；
% bx=(零(m，n))；
% by=(零(m，n))；
lx=1的百分比：(m-1)
%代表ly=1:n
% %dx(lx，ly)=J(lx 1，ly)-J(lx，ly)；%dx(lx:)=-什么
% if J(lx 1，ly)-J(lx，ly)=pi
% bx(lx，ly)=0-1；
% elseif J(lx 1，ly)-J(lx，ly) -pi
% bx(lx，ly)=0 1；
%否则
% bx(lx，ly)=0；
%结束
% dx(lx，ly)=J(lx 1，ly)-J(lx，ly) 2*pi*bx(lx，ly)；% %注释：bb的形式！需要改变
%结束
%结束
lx=1时的百分比：m
ly=1的百分比：(n-1)
% % dy(lx，ly)=J(lx，ly 1)-J(lx，ly)；%!通尚
% if (J(lx，ly 1)-J(lx，ly))=pi
% by(lx，ly)=0-1；
% elseif (J(lx，ly 1)-J(lx，ly)-pi
% by(lx，ly)=0 1；
%否则
% by(lx，ly)=0；
%结束
% dy(lx，ly)=J(lx，ly 1)-J(lx，ly) 2*pi*by(lx，ly)；% %注释：bb的形式！需要改变
%结束
%结束
% dx(m:)=0；
% dy(:n)=0；
% r=零(m，n)；
% % R(1，2:n)=dx(1，2:n)-0 dy(1，2:n)-dy(1，1:(n-1))；
% % R(2:m，1)=dx(2:m，1)-dx(1:(m-1)，1) dy(2:m，1)-0；
% % r(1.1)=dx(LX，ly) dy(lx，ly)；
lx=1时的百分比：m
%代表ly=1:n
% if (lx==1) (ly~=1)
% R(lx，ly)=dx(lx，ly) dy(lx，ly)-dy(lx，ly-1)；
% elseif (ly==1) (lx~=1)
% R(lx，ly)=dx(lx，ly)-dx(lx-1，ly) dy(lx，ly)；
% elseif (ly==1) (lx==1)
% R(lx，ly)=dx(lx，ly) dy(lx，ly)；
%否则
% R(lx，ly)=dx(lx，ly)-dx(lx-1，ly) dy(lx，ly)-dy(lx，ly-1)；
%结束
%结束
%结束
% % % % % % % % % % % % % % % % % % % % %。解方程百分之九十！！RR=维也纳2(R，[10 10])；
g=1时的百分比：m
t=1时的百分比：n
% GG((g-1)*n t)=R(g，t)；
%结束
%结束
% GT=GG ；
%
% TT=cgs(F，GT，1e-006，100)；00
% X=零(m，n)；
g=1时的百分比：m
t=1时的百分比：n
% X(g，t)=TT((g-1)* n t)；%/2%%%%%2倍关系
%结束
%结束
%图(8),目（双(X))；
%%%%%%%%%%%%%%%%%泽尼克拟合%%%%%%
A=imresize(X，[65 65]，双线性);
zj=16
[M0号]=尺寸(一)；
coef(1:ZJ)=0；
RX=(M0 1)/2；
ry=(n0 ^ 1)/2；
R0=min(RX，RY)；
%拟合(1:M0，1:N0)=0；
fitted _ qulijiao=zeros([M0 N0])；
% fitted _ ball=zeros([M0 N0])；
r=0；k=0；
z=0；
wfront=0；
对于i=1:M0
对于j=1:N0
拉德=sqrt((i-rx)^2(j-ry)^2)；%椭圆的问题？
如果rad=R0% A(i，j)=0
r=rad/R0；
如果r~=0
theta=atan2(RX-i，j-RY)；
目标
k=k1
wfront(k)=A(i，j)；
z(1，k)=1；
z(2，k)=r * cos()；
z(3，k)=r * sin()；
z(4,k)=2*r^2-1；
z(5,k)=r^2*cos(2*theta)；
z(6,k)=r^2*sin(2*theta)；
z(7,k)=(3*r^3-2*r)*cos(theta)；
z(8,k)=(3*r^3-2*r)*sin(theta)；
z(9,k)=6*r^4-6*r^2 1；
z(10,k)=r^3*cos(3*theta)；
z(11,k)=r^3*sin(3*theta)；
z(12,k)=(4*r^4-3*r^2)*cos(2*theta)；
z(13,k)=(4*r^4-3*r^2)*sin(2*theta)；
z(14,k)=(10*r^5-12*r^3 3 * r)* cos()；
z(15,k)=(10*r^5-12*r^3 3 * r)* sin()；
z(16,k)=20*r^6-30*r^4 12*r^2-1；
目标
目标
目标
orthop=zeros(zj，k)；
orthop(1，1:k)=z(1，)；
bb(1)=wfront*orthop(1，)/(orthop(1，)*orthop(1，))；
zterm=zj
对于n=2:zterm
orthop(n，)=z(n，)；
对于m=1:n-1
aa(n，m)=z(n，)*orthop(m，)/(orthop(m，)*orthop(m，))；
orthop(n，)=orthop(n，)-aa(n，m)*orthop(m，)；
目标
bb(n)=wfront*orthop(n，)/(orthop(n，)*orthop(n，)；
目标
coef(zterm)=bb(zterm)；
对于n=1:zterm-1
coef(z term-n)=bb(z term-n)-coef(z term-n 1:z term)* aa(z term-n 1:z term，z term-n)；
目标
coef _ qulijiao=coef
coef _ qu李娇(1:4)=0；
r=0；
ZZ=零([1z项])；
对于i=1:M0
对于j=1:N0
拉德=sqrt((i-rx)^2(j-ry)^2)；
如果rad=R0
r=rad/R0；
如果r~=0
theta=atan2(RX-i，j-RY)；
目标
ZZ(1)=1；
ZZ(2)=r * cos()；
ZZ(3)=r * sin()；
zz(4)=2*r^2-1；
zz(5)=r^2*cos(2*theta)；
zz(6)=r^2*sin(2*theta)；
zz(7)=(3*r^3-2*r)*cos(theta)；
zz(8)=(3*r^3-2*r)*sin(theta)；
zz(9)=6*r^4-6*r^2 1；
zz(10)=r^3*cos(3*theta)；
zz(11)=r^3*sin(3*theta)；
zz(12)=(4*r^4-3*r^2)*cos(2*theta)；
zz(13)=(4*r^4-3*r^2)*sin(2*theta)；
zz(14)=(10*r^5-12*r^3 3 * r)* cos()；
zz(15)=(10*r^5-12*r^3 3 * r)* sin()；
zz(16)=20*r^6-30*r^4 12*r^2-1；
% fitted(i，j)=coef * ZZ ；
拟合_曲丽娇(I，j)=coef _曲李娇* ZZ ；
目标
目标
目标
如果系数(4)为0
结果_曲丽娇=拟合_曲丽娇3.25/(4 * pi)；
其他
结果_曲丽娇=(-1)*拟合_曲丽娇3.25/(4 * pi)；
目标
图(9),网格(double(Result _ qulijiao))；3仿真结果
编辑
编辑
编辑
四参考文献[1]钱晓凡，饶帆，李兴华，等。精确最小二乘相位解包裹算法[J]的缩写.中国激光， 2012, 39(2):5.
博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的矩阵实验室仿真，相关矩阵实验室代码问题可私信交流。
编辑
郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。