内容发布更新时间 : 2024/11/19 6:35:25星期一 下面是文章的全部内容请认真阅读。

线性代数

案例

Cayler-Hamilton定理

【实验目的】

1．理解特征多项式的概念

2．掌握Cayler-Hamilton定理 【实验要求】掌握生成Vandermonde矩阵的vander命令、求矩阵特征多项式系数的poly()命令、求矩阵范数的norm命令及矩阵多项式运算的polyvalm命令 【实验内容】

Cayler-Hamilton定理是矩阵理论中的一个比较重要的定理，其内容为：若矩阵A的特征多项式为

f(s)?det(sI?A)?a1sn?a2sn?1???ansn?an?1

则有f(A)?0,亦即

a1An?a2An?1???anA?an?1E?0

假设矩阵A为Vandermonde矩阵，试验证其满足Cayler-Hamilton定理。 【实验方案】

Matlab提供了求取矩阵特征多项式系数的函数poly()，但是poly()函数会产生一定的误差，而该误差在矩阵多项式求解中可能导致了巨大的误差，从而得出错误的结论。

在实际应用中还有其他简单的数值方法可以精确地求出矩阵的特征多项式系数。例如，下面给出的Fadeev-Fadeeva递推算法也可以求出矩阵的特征多项式。

1?c??tr?ARk?,k?1,2,...,n?kk???R1?I,Rk?ARk?1?ck?1I,k?2,...,n

该算法首先给出一个单位矩阵I，并将之赋给R1，然后对每个k的值分别求出特

征多项式参数，并更新Rk矩阵，最终得出矩阵的特征多项式的系数ck。该算法可以直接由下面的Matlab语句编写一个poly1??函数实现：

Function c=poly1(A) [nr,nc]=size(A);

if nc==nr % 给出若为方阵，则用Fadeev-Fadeeva算法求特征多项式 I=eye(nc); R=I; c=[1 zeros(1,nc)];

for k=1:nc,c(k+1)=-1/k*trace(A*R);r=A*R+c(k+1)*I; end

1

elseif (nr==1 \\ nc==1) % 给出为向量时，构造矩阵

A=A(isfinite(A));n=length(A) ; % 出去非数或无界的特征根 c=[1 zeros(1,n)]; for j=1:n

c(2:(j+1))=c(2:(j+1))-A(j).*c(1:j);

end

else % 参数有误则给出错误信息

error (’Argument must be a vector or a square matrix.’) end.

【实验过程】

>> A = vander([1 2 3 4 5 6 7]); 运行结果： A =

1 1 1 1 1 1 1 64 32 16 8 4 2 1 729 243 81 27 9 3 1 4096 1024 256 64 16 4 1 15625 3125 625 125 25 5 1 46656 7776 1296 216 36 6 1 117649 16807 2401 343 49 7 1

>> A 运行结果： aa1 =

1.0e+009 *

0.0000 -0.0000 -0.0002 0.0287 1.1589 -6.2505 -2.4223 0.0249

如调用新的poly1()函数，则可以得出如下的精确结果。

>> aa1=poly1(A);b1=polyvalm(aa1,A);norm(B1) 运行结果： ans =

0

可见，由此得出的B矩阵就会完全等于0，故该矩阵满足Cayley-Hamilton定理。

2

小行星轨道问题

【实验目的】

1. 掌握线性方程组求解 2. 加深对正交变换的理解

3. 掌握Matlab软件中的ezplot、zplot命令的区别和适用范围 【实验要求】

掌握绘制隐函数曲线ezplot命令和彗星状轨迹图comet命令 【实验内容】

天文学家要确定一颗小行星绕太阳运行的轨道，在轨道平面内建立以太阳为原点的直角坐标系，在两坐标轴上取天文测量单位（一天文单位为地球到太阳的平均距离：9300万里）。在五个不同的时间点对小行星作了观察，测得轨道上五个点的坐标数据如下:

表 2-1 小行星观测数据 x 4.5596 5.0816 5.5546 5.9636 6.2756 y 0.8145 1.3685 1.9895 2.6925 3.5265 由开普勒第一定律知，小行星轨道为一椭圆。设方程为 a1x2?2a2xy?a3y2?2a4x?2a5y?1?0

试确定椭圆的方程并在轨道的平面内以太阳为原点绘出椭圆曲线。并应用坐标平移变换和正交变换将上例题中的二次曲线方程化为标准方程，绘椭圆轨道图，完成小行星运行的动态模拟。

【实验方案】

（1）二次曲线方程中有五个待定系数：a1，a2，a3，a4，a5。将观察所得的五个点坐标数据(xj,yj)，(j?1,2,?,5)代入二次曲线方程得到关于a1，a2，a3，

a4，a5的线性方程组

?a1x12?2a2x1y1?a3y12?2a4x1?2a5y1??1?22?a1x2?2a2x2y2?a3y2?2a4x2?2a5y2??1?22?a1x3?2a2x3y3?a3y3?2a4x3?2a5y3??1 ?22ax?2axy?ay?2a4x4?2a5y4??11424434??ax2?2axy?ay2?2ax?2ay??1255354555?15 3