线性代数试后验--方阵对角化与matlab在微积分中的应用.ppt

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;dr