234葛雷+数学实验周15教程.doc

淮 阴 工 学 院 数学软件专业实践周 班 级： 金融1131 姓 名： 葛雷 学 号： 1134104134 教 师： 王小才　 胡平 教 师 王小才 胡平 答辩分 总 分 得 分 数 理 学 院 数学实验周题目 一、代数与矩阵综合题 生成3阶矩阵 ， 其中，系数向量利用下述方法给定：；m为你的学号最后2位。 (1) 求的秩，行列式，的伴随矩阵。 程序代码： clc clear a=fliplr(-(randi(8,3,1)-4)) A=[0,1,0;0,0,34;a] save(filename1, a,A) disp(A的秩为:) disp(rank(A)) disp(A的行列式) disp(det(A)) disp(A的伴随矩阵) disp(inv(A)*det(A)) 运行结果：A = 0 1 0 0 0 34 2 -2 -2 A的秩为: 3 A的行列式 68 A的伴随矩阵 68 2 34 68 0 0 0 2 0 (2) 向量由生成， 求方程的解。 程序代码： b=randi(20,3,1)-4 x=A\b 运行结果： x = 7.1765 -1.0000 0.1765 (3)求的特征多项式，计算，并对结果进行必要的讨论； PA=poly(A) fx=poly2str(PA,x) FA=A^3+2*A^2+68*A-68 运行结果： PA = 1.0000 2.0000 68.0000 -68.0000 fx = x^3 + 2 x^2 + 68 x – 68 FA = 0 -68 -68 -68 0 -68 -68 -68 0 (4) 确定一个区间，使该区间包含特征方程所有的根，并画出曲线的图像（显示的效果尽可能地好）； 程序代码： c=0,d=2 x=-10:0.05:10 y=x.^3+2*x.^2+68*x-68 plot(x,y,color,r,Linewidth,3)  (5) 使用函数roots 或 solve 求方程最小的实根或唯一的实根； 程序代码： syms x x1=solve(x^3+2*x^2+68*x-68==0) x=vpa(x1) 结果：x = 0.95989 - 1.4799 + 8.2856*i - 1.4799 - 8.2856*i 结论：有唯一实根x=0.95989 (6)求的特征值和化成Jordan 标准型的变换矩阵； 运行代码 D=1/2*(A+A) disp(D的特征值 ) disp(poly(D)) [V,J]=jordan(D) Q=V 结果： D的特征值 1.0000 2.0000 -257.2500 -16.5000 Q = -0.0312 - 0.0000i 0.1014 + 0.0000i -32.2741 - 0.0000i -0.9381 + 0.0000i 1.0627 - 0.0000i 2.1381 - 0.0000i 1.0000 1.0000 1.0000 (7) 根据计算的结果，指出矩阵具有的性质，并说明是否与你所了解的结论一致； 程序：K=Q*Q 结果：K = 1.0e+03 * 0.0019 -0.0000 + 0.0000i 0.0000 - 0.0000i -0.0000 - 0.0000i 0.0021 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 - 0.0000i 1.0472 K是一个对角矩阵，一致的 (8) 求二元函数 的具体表达式，说明该二元函数是那种类型的二次曲面，编写函数型m文件定义该二元函数； 程序代码： syms x y z z=[x y 1]*D*[x;y;1] 结果：z =x + 16*y + x*(y/2 + 1) + y*(x/2 + 16) – 2 类似于马鞍面的二次曲面 function z=m