%% 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 矩陣中所有點到直線的距離時，并將距離保存到另一個矩陣，用的是循環，計算起來很慢；大伙幫我看看有沒有優化的方法，謝謝