线性代数以及空间的解析几何实验.docVIP

线性代数与空间解析几何实验报告（四） 实验名称： 方阵对角化与微积分的应用（单周） 姓 学 号： 地 点 日 期： 1.实验目的： 掌握MATLAB软件线性变换的几何含义； 掌握MATLAB软件求解特征值和特征向量的方法。 掌握MATLAB软件分析方阵可否对角化的方法。 掌握MATLAB软件方阵正交化的方法； 掌握MATLAB软件化二次型为标准型的方法； 通过MATLAB软件进一步理解方阵对角化，二次型及其标准型的意义。 2.实验指导： 表4 与本实验相关的MATLAB命令 指令 功能 eig(A) 以列向量形式返回方阵的特征值 [u, v] = eig(A) 返回两个矩阵u和v，矩阵v为矩阵A的特征值所构成的对角阵，矩阵u的列为矩阵A的单位特征向量，它与v中的特征值一一对应 poly 矩阵的特征多项式 solve(d) 求代数方程的符号解 roots( ) 求多项式的零点 eigshow(A) 用MATLAB演示矩阵A的特征向量几何含义 trace 计算矩阵的迹 norm 计算矩阵范数或模 orth 求出可将矩阵化为对角阵的正交阵 eval( ) 执行括弧中字符串的有效表达式 syms x 定义x为符号变量 format rat 定义输出格式为小整数比的近似值 taylor(f,n,x,a)f进行泰勒展开 int（s,v） 以v为自变量，对被积函数或符号表达式s求不定积分 int(s,v,a,b) 求定积分 3.实验作业： 注：请你把结果答在程序与运行结果的右侧栏中，答在其它地方无效！ 题 目 1 1、求矩阵A的特征值与特征向量，并分析其几何意义随机矩阵，每个元素在1-10之间的整数。 程 序 与 运 行 结 果： A=round(rand(2,2)*10) A = 10 6 2 5 [U,V]=eig(A) U = 0.9590 -0.6632 0.2832 0.7485 V = 11.7720 0 0 3.2280 eigshow(A) 题 目 2 2、矩阵A。 其中A为随机矩阵，每个元素在1-10之间的整数；为自己学号后两位。 程 序 与 运 行 结 果： A=round(rand(3,3)*10) A = 8 8 2 4 9 4 6 7 9 [P,Q]=eig(A) P = 0.5142 0.7784 0.7784 0.5053 -0.3142 + 0.1764i -0.3142 - 0.1764i 0.6930 -0.4084 - 0.3122i -0.4084 + 0.3122i Q = 18.5564 0 0 0 3.7218 + 1.0108i 0 0 0 3.7218 - 1.0108i Y=P^-1*Q^35*P Y = 1.0e+044 * 0.7150 + 0.0000i 1.0823 + 0.0000i 1.0823 + 0.0000i 0.5884 - 0.0238i 0.8906 - 0.0361i 0.8906 - 0.0361i 0.5884 + 0.0238i 0.8906 + 0.0361i 0.8906 + 0.0361i 题 目 3 3、利用Taylor展开函数至5项，7项，并与函数图形比较（函数输入格式log(1+x)） 程 序 与 运 行 结 果： syms x f=log(1+x) f = log(1+x) taylor(f,3) ??? taylor(f,3) | Error: Missing operator, comma, or semicolon. taylor(f,3) ans = x-1/2*x^2