实验3_非线性方程AX=0的解法..docVIP

线性方程组AX=B的数值解法 1.实验描述 1．P93.1，2,3:通过矩阵可表示立方体的坐标位置，与另一矩阵相乘可实现立方体坐标位置进行变换 2．P108.1:不通过行变换就能解决三角线性方程。 3.p109.7：将单位矩阵表示成列矩阵，通过对目标矩阵分别求解得出列矩阵从而得到目标矩阵的逆矩阵。 4. 120.3：将单位矩阵表示成列矩阵，通过分解成上下三角矩阵对目标矩阵分别求解得出列矩阵从而得到目标矩阵的逆矩阵。 5.p120.4:应用程序3.3求解基尔霍夫电流。 6.p129.4:应用高斯-赛德尔迭代法求解带状方程。 2.实验内容 P93.1. 单位立方体位于第一卦限，一个顶点在原点。首先，以角度沿y轴旋转立方体，然后再以角度沿z轴旋转立方体。求旋转后立方体的8个顶点的坐标，并与例3.10的结果比较。它们的区别是什么？试通过矩阵一般不满足交换律的事实进行解释。使用plot3命令画出3个图形。 P93.2. 设单位立方体位于第一卦限，其中一个顶点位于坐标原点。首先以角度沿x轴旋转立方体，然后再以角度沿z轴旋转立方体。求旋转后立方体的8个顶点的坐标。使用plot3画出这3个立方体。 P93.3. 四面体的坐标为（0,0,0），（1,0,0），（0,1,0），（0,0,1）。首先以弧度0.15沿y轴旋转，然后再以弧度-1.5沿z轴旋转，最后以弧度2.7沿x轴旋转。求旋转后的顶点坐标。使用plot3画出这4个立方体。 P108.1 .许多科学应用包含的矩阵带有很多零。在实际情况中很重要的三角形线性方程组有如下形式： …… 构造一个程序求解三角形线性方程组。可假定不需要变换。而且可用第k行消去第k+1行的。 p109.7 下面的习题虽然是针对3x3维矩阵的，但其概念可用于NxN维矩阵。如果矩阵A非奇异，则存在。而且。设，，是的列，而，，是I的列。方程可表示为 则上式等价于三个线性方程组，， 这样求等价于求解三个线性方程组。 使用程序2.2或上题中的程序求解下面每个矩阵的逆。通过计算 和使用命令inv(A)检查答案，并解释可能的差异。 （a） （b） 120.3 修改程序3.3，使得它可以通过重复求解N个线性方程组 其中J=1，2，…，N 来得到， 则 、 而且 ，保证对LU分解只计算一次。 p120.4 基尔霍夫电压定律说明电路网络中任意单向闭路的电压和为零。现有如下电路线性方程： 如果 使用程序求解电流 p129.4 利用高斯-赛德尔迭代法求解下列带状方程。 3.实验结果及分析 P93.1: 算法: 令X=zeros(8,3);X([5:8,11,12,15,16,18,20,22,24])=1;d=[1 2 4 3 1 5 6 8 7 5 6 2 4 8 7 3]; i=0。 判断i100是否成立，若成立，执行步骤（3）；若不成立， r1=[cos(i*pi/600) -sin(i*pi/600) 0;0 1 0;-sin(i*pi/600) 0 cos(i*pi/600)];U=X*r1;plot3(U(d,1),U(d,2),U(d,3));drawnow，i=i+1,返回步骤（2）.\ i=0. 判断i100是否成立，若成立，执行步骤（4）；若不成立，r2=[cos(i*pi/400) -sin(i*pi/400) 0;sin(i*pi/400) cos(i*pi/400) 0;0 0 1];W=U*r2;plot3(W(d,1),W(d,2),W(d,3));drawnow，i=i+1,返回步骤（3）. subplot(2,2,1);plot3(X(d,1),X(d,2),X(d,3))subplot(2,2,2);plot3(U(d,1),U(d,2),U(d,3));subplot(2,2,3)； plot3(W(d,1),W(d,2),W(d,3)); xlabel(x)；ylabel(y)；zlabel(z)；view(3); rotate3d。 图1.1(左上) 初始立方体 图1.2（右上），沿y轴旋转 图1.3 （左下）沿z轴旋转 表1.1：第次旋转后立方体的坐标 0 0.8660 0.8660 0 0.5000 0.5000 1.3660 1.3660 Y 0 0 1.0000 1.0000 1.0000 0 0 1.0000 Z 0 -0.5000 -0.5000 0 0.8660 0.8660 0.3660 0.3660