[研究生入学考试]矩阵理论线性代数与计算1.docVIP

线性代数与相关计算的Matlab应用 行列式（determinant）的定义与理论计算 行列式相应具有若干性质 行列式的理论计算 利用行列式的性质将行列式化成上（下）三角形行列式 行列式在解析几何中的两个基本应用 求平面三角形的面积 ，其中为三角形三个顶点坐标，且按逆时针排列。 求空间四面体的体积 ，其中为四面体四个顶点坐标，且按右手系法则排列。 行列式的Matlab计算 在Matlab命令窗口输入矩阵A，然后用det(A) 例如： a=[2 3 4;6 7 8;9 10 11]; d=det(a) 再如 d=det(magic(5)) Matlab中如何计算行列式呢？ 可在命令窗口键入doc det 。 解线性方程组 ， 其矩阵形式为， 其中 2.1）按照克莱姆法则，可知当系数行列式时，方程组有解且是唯一的。解的形式 其中克莱姆法则的一个应用为 问题：过平面上的个不同观测点的次数不超过的多项 式存在且唯一。 假设过这些观测点的多项式为，于是有 由于 （范得蒙行列式） 所以得到问题的结论。但是如何来获得这个多项式呢？用Matlab做实验。 x=linspace(0,2*pi,5); %均匀生成[0，2*pi]区间上100个点 y=cos(x); p2=polyfit(x,y,4); %用4阶多项式拟合cos(x) xx=linspace(0,2*pi,100); y2=polyval(p2,xx); %求4阶拟合多项式在生成点x的值 yy=cos(xx); plot(x,y,go,xx,y2,b-,xx,yy,r-) title(用4阶多项式拟合cos(x)) 2.2）用求逆的办法，求解线性方程组（理论求法） （行初等变换），然后有， 但是对于高阶的线性方程组，用理论求法无法实现，为此需要借助计算机来实现。 2.3）用 Matlab求解线性方程组以及方阵求逆 输入A，b，x=A\b即可， 例1 A=magic(5),b=[1,2,3,4,5];x=A\b Ainv=inv(A),如果det(A)=0,可以求逆吗？例如A=magic(4),a=inv(A)会得到什么结果呢？ 例2 求的解。 a=[1 1 0;-3 0 3;1 -3 4];b=[1 4 2]; x=a\b 此例也可以用符号运算求解： syms x1 x2 x3; %定义符号变量 [x1,x2,x3]=solve(x1+x2=1,-3*x1+3*x3=4,x1-3*x2+4*x3=2)%solve函数还可以用于非线性解。 vap([x1,x2,x3],5) % 将解用5位有效数字表示 或 a=sym([1 1 0;-3 0 3;1,-3 4]) %定义符号变量 b=sym([1; 4 ;2]);% 定义符号变量 x=linsolve(a,b);x1=a\b; %线性方程组求解 vpa([x,x1],5) 向量、矩阵数据的生成，矩阵的代数运算以及与方阵的相关数值计算 3.1）向量数据生成 a)手工输入，例a=[1 3 7 8 10 -1]; b)命令输入，例x=linspace(-2,3,10);，x=-2:0.01:3;，x=logspace(-2,3,20);， [x,y]=meshgrid(-1:0.01:,1,-2:0.01:2); 此种数据生成在作图时非常有用,例如： x=linspace(-pi,pi,100); y=sin(x); plot(x,y,r-) title(sin曲线) 另外， [x,y]=meshgrid(-1:0.1:1); z=x.^2+y.^2; figure(1) mesh(x,y,z); figure(2) surf(x,y,z) c) 特殊矩阵生成 单位矩阵eye(n),零矩阵zeros(m,n),全一矩阵ones(m,n),魔方矩阵magic(n),pascal矩阵pascal(n),随机矩阵rand(m,n)和randn(m,n)。 d) 子矩阵的提取与整合 Matlab实现，例给定 a=magic(5);b=pascal(5); al=a(2,:);%提取a的第二行； ac=a(:,4);%提取a的第四列； alc=a(2:3,1:4); %提取a的第二行到第三行以及第一列到第四列的子矩阵。 ab=[a,b];ac=[a;b];%组装两个矩阵，注意两者的区别。 ab=horzcat(a,b);ac=vertcat(a,b)%和上面两个运算相同。 A = magic(3); B = [-5 -6 -9; -4 -4 -2]; C = eye(2