MATLAB学习必备资料.ppt

1 ）设 2 ）设 3 ）设 4 ）设 5 ）设 t ? π/6 ， 1 0 0 1 1 0 , 0 1 0 0 1 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? A1 y1 则 2 0 0 2.0 2.0 0 , 0 0.5 0 0 0.5 0.5 ? ? ? ? ? ? ? ? ? ? ? ? ? ? A2 y2 则 1 2 0 1.0 1.0 0 , 2 2 0 0 0.2 0.2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? A3 y3 则 1.0 0.5 0 1 3 2 , 0 1.0 0 2 4 2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? A4 y4 则 cos t sin t 0 0.8660 0.3660 0.5000 , sin t cos t 0 0.5000 1.3660 0.8660 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? A5 y5 则 这些矩阵的每一列代表平面上一点的横纵 坐标，四列就代表四个点。把第一数据增 补到第五列作为第五点，把这些点顺序连 以直线，就得到相应的四边形。用六个图 分别绘制出数据矩阵 , , , , 和 的 图形，就可以得到图 4.14 的图形。 图 4.14 对单元方格进行几种线性变换后生成的图形 1 y 2 y 3 y 4 y 5 y x 可以看出，矩阵 A1 使原图对纵轴生成 镜像，矩阵 A2 使原图在横轴方向膨胀并在 纵轴方向压缩；矩阵 A3 使原图在左下至右 上伸长而在其垂直方向压缩；矩阵 A4 使原 图向右方剪切变形，矩阵 A5 使原图沿反时 针方向旋转 t ? pi/6 。总之，线性变换把原来 空间的图形在各个方向以不同的比例放大 或缩小。要注意，在 x 平面中的直线，变换 到 y 平面中后，仍然是一根直线，但是其范 数和方向都可能发生变化。 例 4.5 所用的 MATLAB 程序 ea405 的部分语句 x=[0,1,1,0;0,0,1,1]; subplot(2,3,1),fill([x(1,:),0],[x(2,:),0],'r') axis equal,axis([-1.5,1.5,-1,2]),grid on A1=[-1,0;0,1],y1 A1*x subplot(2,3,2),fill([y1(1,:),0],[y1(2,:),0],'g') ... ... MATLAB 还提供了一个平面线性变换的演示函数 eigshow 。分别键入 eigshow( A2 )eigshow( A3 ) ， 它显示不同 x 的经线性变换 Ax 后生成的 y 。 x 取单 位向量，用鼠标可以拖动它转动，形成的轨迹是 单位园。同时画出 Ax 的向量，形成的轨迹如图 4.15 。 图 4.15 矩阵 A 的特征向量和特征值演示（左图 A ? A2 ，右图 A ? A3 ） ? 4.5.2 二维矩阵特征值的计算方法 在特征点处，有 ，即 （ 4-1 ） 要使这个矩阵方程有非零解，它的系数行 列式必须为零，即 。 若 则有