test.m

数值计算课本中第三章程序源代码

U = [ 0 1 1 0 0 0 1 1 0 0 0 0 1 1 1 1;
      0 0 1 1 0 0 0 1 1 0 1 1 1 1 0 0;
      0 0 0 0 0 1 1 1 1 1 1 0 0 1 1 0];
  
Ry = [ 0.866025  0.000000  0.500000;
       0.000000  1.000000  0.000000;
      -0.500000  0.000000  0.866025 ];  
Rz = [ 0.707107 -0.707107  0.000000;
       0.707107  0.707107  0.000000;
       0.000000  0.000000  1.000000 ];
   
subplot(2,2,1);
plot3( U(1,:), U(2,:), U(3,:) );

V = Ry*U;
subplot(2,2,2);
plot3( V(1,:), V(2,:), V(3,:) );

W = Rz*V;
subplot(2,2,3);
plot3( W(1,:), W(2,:), W(3,:) );

subplot(2,2,4);
plot3( U(1,:), U(2,:), U(3,:), 'g' );
hold on
plot3( V(1,:), V(2,:), V(3,:), 'r' );
hold on
plot3( W(1,:), W(2,:), W(3,:) );