线性代数应用举例幻灯片.pptVIP

该方程组有三个未知数，但有五个方程，进一步分析 可以得到该方程组无解，即不存在一个二次多项式曲 线刚好能过已知的五个点。MATLAB 软件提供了一个 利用最小二乘法解决超定方程组解的方法。求系数的 公式也是 a = A\y，以找到一条二次曲线来近似地描述 已知 5 个点的变化情况。 对比插值和拟合的曲线如下图 用平面坐标系中的一个闭合图形来描述刚体，用一个矩 阵 X 来表示它。X 的一列表示刚体一个顶点的坐标。为 了使图形闭合，X 的最后一列和第一列相同；为了实现 刚体的平移运算，给矩阵 X 添加元素值都为 1 的一行， 使 X 为 3×n 矩阵。 若有矩阵： 则可以证明，矩阵 Y1 是刚体 X 沿 x 轴正方向平移 c1， 沿 y 轴正方向平移 c2 后的结果；矩阵 Y2 是刚体 X 以坐 标原点为中心逆时针转动 t 弧度的结果。 例 8 刚体的平面运动 x 0 4 6 10 8 5 3.5 6.1 6.5 3.2 2 0 y 0 14 14 0 0 11 6 6 4.5 4.5 0 0 实 例 用下列数据表示字母 A: 对 A 进行以下平面运动, 并绘制移动前后的图形。 (1) 向上移动 15, 向左移动 30; (2) 逆时针转动π /3; (3) 先逆时针转动3π /4, 然后向上平移 30, 向右平移 20。 解 构造刚体矩阵 X，平移矩阵及转动矩阵。 M1= 1 0 0 0 1 0 -30 20 1 M2= 1 0 0 0 1 0 20 30 1 R1= 0 0 1 0 0 R2= 0 0 1 0 0 MATLAB 程序 ? X=[0,4,6,10,8,5,3.5,6.1,6.5,3.2,2,0; 0,14,14,0,0,11,6,6,4.5,4.5,0,0; ones(1,12)]; % 构造刚体矩阵 X M1=[1,0,-30; 0,1,15; 0,0,1]; % 构造平移矩阵 M1 Y1=M1*X; % 计算平移结果 ? fill(Y1(1, :), Y1(2, :), ‘red’); % 绘制平移后刚体 Y1 ? plot(X(1, :), X(2, :)); % 绘制原来刚体 X hold on axis equal ? R1=[cos(pi/3), -sin(pi/3), 0; sin(pi/3), cos(pi/3), 0; 0,0,1]; % 构造转动矩阵 R1 Y2=R1*X; % 计算旋转结果 fill(Y2(1, :), Y2(2, :), ‘blue’); % 绘制转动后刚体 Y2 ? fill(Y3(1, :), Y3(2, :), ‘black’); % 绘制转动及平移后刚体 grid on hold off ? M2=[1,0,20; 0,1,30; 0,0,1]; % 构造转动矩阵 M2 R2=[cos(3*pi/4), -sin(3*pi/4), 0; % 构造转动矩阵 R2 sin(3*pi/4), cos(3*pi/4), 0; 0,0,1]; Y3=M2*R2*X; % 旋转平移后的结果 Y3 运行的结果如下 线性代数应用实例 取自《线性代数机算与应用指导（MATLAB）版》 2011.12 例 1 平板稳态温度的计算 为了计算平板形导热体的温度分布，将平板划分为许 多方格，每一个节点上的稳态温度将等于其周围四个 节点温度的平均值。由此可得出阶数与节点数相同的 线性方程组，方程的解将取决于平板的边界条件。 这