需要计算一个3d矩阵中所有点到直线的距离<－程序优化，该怎么操作？

Question

%% difinitionclear;x1=1;x2=100;nx=100;y1=50;y2=50;ny=100;z1=41;z2=60;nz=100;matrx1=rand(nx,ny,nz);% 所有点的坐标生成[X,Y,Z] = meshgrid(1:1:nx, 1:ny,1:nz);% e.g.% pointApositionX=Y(1,50,40);% pointApositionY=X(1,50,40);% pointApositionZ=Z(1,50,40);% 穿过两点的所有x点xCross=linspace(x1,x2,nx);% 定义直线两段点a=[x1 y1 z1];b=[x2 y2 z2];％＝＝＝＝＝＝＝＝＝%--matrxDisMatrixmatrxDisMatrix=rand(nx,ny,nz);for i=1:nxfor j=1:nyfor k=1:nzpointijk=[Y(i,j,k), X(i,j,k) ,Z(i,j,k)];[m]=pointlinedistance(a,b,pointijk);% m is the distance matrxDisMatrix(i,j,k)=m;endendend===========这里计算一个3d 矩阵中所有点到直线的距离时，并将距离保存到另一个矩阵，用的是循环，计算起来很慢；大伙帮我看看有没有优化的方法，谢谢

Answer 1

当然可以优化，不过需要对pointlinedistance函数进行相应的调整。