2014线性代数实验 第2次演示实验.ppt

2.3 方程组解的几何含义 例2.4 求下列线性方程组的解，并用MATLAB绘出解的情况. （1） （2） （3） （4） 解：用MATLAB的 rref 命令可以解得： 用MATLAB绘制平面的简单方法为： 方程组（1）有唯一解 方程组（2）有无穷组解 方程组（3）和（4）无解 ezmesh(‘……’) 单引号内为平面方程 在MATLAB命令窗口中运行程序 g02.m，可以得到图形： 三、特征值与特征向量、与方阵 正交对角化 线性变换几何意义：反射、平移、旋转 解：在MATLAB中运行g07.m,可以得到下列图形 例3.2 特征值与特征向量的几何含义 已知矩阵 MATLAB分析特征向量的几何含义。 ，求它们的特征值和特征向量，并用 解：用MATLAB求矩阵特征值和特征向量的方法为： 用MATLAB演示矩阵A的特征向量几何含义的命令为： （1）r=eig(A),列向量r为矩阵A的特征值 （2）[V,D]=eig(A)，对角矩阵D的对角线元素为矩阵A的特征值，矩阵V的列向量为矩阵A的特征向量。 eigshow(A) 在MATLAB命令窗口运行程序g08.m，可以分别得到图形： 3.2 方阵正交对角化 下表给出与本次实验相关的MATLAB命令 为了能够更好的求得矩阵的n次幂，需要对方阵进行对角化处理。 命令 功能 disp(‘suhuh’) 显示单引号中的字符串 [P,D]=eig(A) 矩阵D为矩阵A的特征值所构成的对角阵，矩阵P的列为矩阵A的单位特征向量，它与D中的特征值一一对应 例3.1 化方阵 为对角阵。 解： a=[2 2 -2;2 5 -4;-2 -4 5]; %输入矩阵A [d,v]=eig(a) %求特征值与特征向量 d = -0.2981 0.8944 0.3333 -0.5963 -0.4472 0.6667 -0.7454 0 -0.6667 v = 1.0000 0 0 0 1.0000 0 0 0 10.0000 四、软件在微积分中的应用 4.1 实验指导 下表给出与本次实验相关的MATLAB命令 命令 功能 taylor(f,n,x,a) 将f进行泰勒展开 syms x 定义x为符号变量 int（s,v） 以v为自变量，对被积函数或符号表达式s求不定积分 int(s,v,a,b) 求定积分 4.1 实验内容 泰勒多项式逼近 6 4 2 2 4 6 4 2 2 4 O 泰勒多项式逼近 6 4 2 2 4 6 O 4 2 2 4 例4.4.1 将 分别展开为x和x-1的幂级数. 解：计算关于x展式的前8项 ?syms?x?f=x/sqrt(1+x^2);?taylor(f,8)? ans?= -?(5*x^7)/16?+?(3*x^5)/8?-?x^3/2?+?x 计算关于x-1展式的前3项 taylor(f,3,x,1) ans = 2^(1/2)/2 + (2^(1/2)*(x - 1))/4 - (3*2^(1/2)*(x - 1)^2)/16 例4.4.2 关于 的泰勒展开式逼近. ?syms?x?a=taylor(sin(x),3) ?a?=? x?b=taylor(sin(x),5) b?= x?-?x^3/6?c=taylor(sin(x),11) c?= x^9/362880?-?x^7/5040?+?x^5/120?-?x^3/6?+?x?d=taylor(sin(x),13) d?= -x^11+?x^9/362880?-?x^7/5040?+?x^5/120?-? x^3/6 + x 解： ?ezplot(sin(x));?hold?on?ezplot(x);?hold?on?ezplot(x?-?x^3/6);?hold?on?ezplot(x^9/362880?-?x^7/5040?+?x^5/120?- … x^3/6?+?x);?hold?on?ezplot(-?x^11+?x^9/362880?-… x^7/5040?+?x^5/120?-?x^3/6?+?x