P0=[4 5 6 7 8 9 10 11 12];%定义固定点座标
%xs1 ys1 zs1 xs2 ys2 zs2 xs3 ys3 zs3
v1=[1 1 1];%定义直线向量
L=-10*(rand(15,3)*2-1)%产生15*3的随机位移矩阵,定义不同动点到固定点P的距离,负号表示方向
%行向量表示每次实验中三个位移传感器的位移,列向量表示每个位移传感器15次实验中的位移
P1=zeros(15,9)%预设内存给动点集,每一行分别为新的xs1 ys1 zs1 xs2 ys2 zs2 xs3 ys3 zs3
for i=1:length(L)
P1(i,1)=P0(1)+L(i,1)*v1(1)/norm(v1);P1(i,2)=P0(2)+L(i,1)*v1(2)/norm(v1);P1(i,3)=P0(3)+L(i,1)*v1(3)/norm(v1)%计算固定点1对应的动点各分量
P1(i,4)=P0(1)+L(i,2)*v1(1)/norm(v1);P1(i,5)=P0(5)+L(i,2)*v1(2)/norm(v1);P1(i,6)=P0(6)+L(i,2)*v1(3)/norm(v1)%计算固定点2对应的动点各分量
P1(i,7)=P0(1)+L(i,3)*v1(1)/norm(v1);P1(i,8)=P0(8)+L(i,3)*v1(2)/norm(v1);P1(i,9)=P0(6)+L(i,3)*v1(3)/norm(v1)%计算固定点3对应的动点各分量
end
计算结果如下:
P1 =
8.8976 9.8976 10.8976 2.8224 6.8224 7.8224 8.5420 15.5420 13.5420
9.1505 10.1505 11.1505 6.7370 10.7370 11.7370 -1.3335 5.6665 3.6665
3.6444 4.6444 5.6444 2.2208 6.2208 7.2208 9.7200 16.7200 14.7200
0.7765 1.7765 2.7765 1.8151 5.8151 6.8151 0.8256 7.8256 5.8256
-1.0115 -0.0115 0.9885 1.1346 5.1346 6.1346 0.3361 7.3361 5.3361
8.2735 9.2735 10.2735 4.5711 8.5711 9.5711 -0.2573 6.7427 4.7427
3.2053 4.2053 5.2053 8.8056 12.8056 13.8056 8.7985 15.7985 13.7985
4.3534 5.3534 6.3534 7.1295 11.1295 12.1295 5.1572 12.1572 10.1572
9.6361 10.6361 11.6361 -0.7728 3.2272 4.2272 6.7728 13.7728 11.7728
5.8807 6.8807 7.8807 8.0140 12.0140 13.0140 0.5351 7.5351 5.5351
7.9008 8.9008 9.9008 0.2378 4.2378 5.2378 4.7920 11.7920 9.7920
0.6019 1.6019 2.6019 3.5573 7.5573 8.5573 -0.7417 6.2583 4.2583
6.1799 7.1799 8.1799 -1.7289 2.2711 3.2711 7.6737 14.6737 12.6737
3.6705 4.6705 5.6705 8.8708 12.8708 13.8708 6.7274 13.7274 11.7274
7.8608 8.8608 9.8608 4.6619 8.6619 9.6619 8.0930 15.0930 13.0930
P1的行向量分别为xs1 ys1 zs1 xs2 ys2 zs2 xs3 ys3 zs3,列向量分别为TEST 1 TEST 2...TEST 15