cyn3942637
2010-05-25, 17:19
要求一个平面Z=AX+BY+C, 给出了一些三维点，要求使点到这个平面的 最大值与最小值之差 为最小 时的A,B,C的值？？:(:(
我编的算法如下，不行
x=[564.350 564.704 564.563 554.944 553.457 555.864 547.304 546.524 546.161 529.135 534.381 535.569 ]
y=[487.186 497.095 510.358 489.682 499.659 508.552 489.794 499.659 512.507 493.103 502.950 511.496 ]
z=[219.127 214.034 207.437 221.187 216.696 211.390 223.820 219.135 212.828 228.552 221.764 217.055 ]
z=X(1)*x+X(2)*y+X(3);
for j=1:length(x)
r=(X(1)*x(j)+X(2)*y(j)-z(j)+X(3))/(sqrt(X(1)^2+X(2)^2+1^2));
end
function F=myfun(X)
F=max(r)-min(r);
x0=[1,1,1] %起始点
[X,fval]=fminunc(@myfun,x0)
set(handles.edit2,'string',num2str(fval))
我编的算法如下，不行
x=[564.350 564.704 564.563 554.944 553.457 555.864 547.304 546.524 546.161 529.135 534.381 535.569 ]
y=[487.186 497.095 510.358 489.682 499.659 508.552 489.794 499.659 512.507 493.103 502.950 511.496 ]
z=[219.127 214.034 207.437 221.187 216.696 211.390 223.820 219.135 212.828 228.552 221.764 217.055 ]
z=X(1)*x+X(2)*y+X(3);
for j=1:length(x)
r=(X(1)*x(j)+X(2)*y(j)-z(j)+X(3))/(sqrt(X(1)^2+X(2)^2+1^2));
end
function F=myfun(X)
F=max(r)-min(r);
x0=[1,1,1] %起始点
[X,fval]=fminunc(@myfun,x0)
set(handles.edit2,'string',num2str(fval))