4实验第四次——修改版课件.pptVIP

4.1.2 实验指导 4.1.3 实验内容 解：方法一：用定义: syms k %定义k为符号变量 A [0.5 0.25; 0.25 0.5]; %输入矩阵 B A-k*eye length A %构造矩阵B A-KI B [ 1/2-k, 1/4] [1/4, 1/2-k] d det B ; %计算B行列式 d k^2 - k + 3/16 4.4.1 实验内容 泰勒多项式逼近 泰勒多项式逼近 4.4.2 实验内容 符号积分 符号积分由函数int来实现. 该函数的一般调用格式为: int s ：没有指定积分变量时，系统按最接近x的变量对被积函数或符号表达式s求不定积分。 int s,v ：以v为自变量，对被积函数或符号表达式s求不定积分。 int s,v,a,b ：求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷 inf 。 * 线性代数实验教学团队 * 线性代数实验教学团队 哈尔滨工程大学理学院 线性代数实验教学团队 实验四 方阵对角化与MATLAB在微积分中的应用 4.1 特征值与特征向量 4.1.1 实验目的 1.掌握MATLAB软件线性变换的几何含义； 2.掌握MATLAB软件求解特征值和特征向量的方法。 下表给出与本次实验相关的MATLAB命令 用MATLAB演示矩阵A的特征向量几何含义 eigshow A 求多项式的零点 roots 求代数方程的符号解 solve d 矩阵的特征多项式 poly 返回两个矩阵u和v，矩阵v为矩阵A的特征值所构成的对角阵，矩阵u的列为矩阵A的单位特征向量，它与v中的特征值一一对应 [u, v] eig A 以列向量形式返回方阵A的特征值 eig A 功能 指令 例4.1.1 求3阶方阵A的特征多项式，并求特征值。 解：方法一： A [11 12 13;14 15 16;17 18 19]; %输入矩阵A PA poly A %A的特征多项式 PA 1.0000 -45.0000 -18.0000 -0.0000 v roots PA %求特征多项式零点，即特征值 v 45.3965 -0.3965 0.0000 方法二： A [11 12 13;14 15 16;17 18 19]; eig A ans 45.3965 -0.3965 0.0000 图4.1线性变化几何意义 解：由附录1.2程序可以画出如下图4.1： 附录程序 1.绘制向量程序： % 绘制二维向量函数 function drawvec u plot [0;u 1 ],[0;u 2 ] ; % 画向量线段 hold on theta acos u 1 /norm u ; % 计算夹角 if u 2 0 theta 2*pi-theta; %当向量在第四象限 end fill [u 1 -0.5*cos theta+pi/12 ,u 1 ,u 1 -0.5*cos theta-pi/12 ],[u 2 -0.5*sin theta+pi/12 ,u 2 ,u 2 -0.5*sin theta- pi/12 ],black ; %用线段填充 hold off 2．线性变换程序： % 线性变换的几何含义 x [2;1];A1 [-1,0;0,1];A2 [1,0;0,-1];A3 [0.5,0;0,2];A4 [cos pi/2 ,sin pi/2 ;-sin pi/2 ,cos pi/2 ]; y1 A1*x;y2 A2*x;y3 A3*x;y4 A4*x; subplot 2,2,1 ; drawvec x ;hold on;drawvec y1 ;axis equal;axis [-3,3,-2,2] ; text x 1 ,x 2 +0.5,x ;text y1 1 ,y1 2 +0.5,y_1 ;title y_1 A_1x ;grid on; subplot 2,2,2 ; drawvec x ;hold on;drawvec y2 ;axis equal;axis [-3,3,-2,2] ; text x 1 ,x 2 +0.5,x ;text y2 1 ,y2 2 +0.5,y_2 ;title y_2 A_2x ;grid on; subplot 2,2,3 ; drawvec x ;hold on;drawvec y3 ;axis equal;axis [-3,3,-2,2] ; text x 1 ,x 2 +0.5,x ;text y3 1 -1,