线性代数实验_行列式矩阵.pptVIP

第二篇 线性代数实践 实验在线性代数中的重要性 利用软件工具进行实验有以下的一些好处： 好处一：对于低价（三阶及以下）的线性代数问题，MATLAB能提供图形帮助，有利于牢固地掌握概念。 好处二：对于高阶的问题，MATLAB能提供计算程序，方便而简捷，节省时间。 好处三：由于解题快捷，在课程中可以较多地放进线性代数的应用实例。扩展学生的视野，提高学习的目的性和积极性。 ‘线性代数实践’的预期效果 线性代数抽象吗？用了MATLAB后，你会知道它的概念都基于空间形象。 线性代数冗繁吗？学了MATLAB后，你会懂得它的计算全可有简明程序。 线性代数枯燥吗？ 精通MATLAB后，你会发现它的应用极其广泛又精彩。 首先着重于对低阶概念的理解，要在二维和三维空间内体会线性代数的定义。 结合相应的MATLAB程序，弄清低阶的算法，然后再引伸到高阶方程中去，进一步搞清其算法和程序应有的扩展。 对于应用问题，不必全看，可结合自已能理解的问题先看。 矩阵元素的引用 矩阵元素的引用 人口迁徙问题 例 假设一个城市的总人口数是固定不变，但人口的分布情况变化如下：每年都有5％的市区居民搬到郊区；而有15％的郊区居民搬到市区。若开始有700000人口居住在市区，300000人口居住在郊区。请利用分析： （1）10年后市区和郊区的人口各是多少？ （2）30年后、50年后市区和郊区的人口各是多少？ （3）分析（2）中数据相似的原因。 解：这个问题可以用矩阵乘法来描述。令人口变量 其中 为市区人口所占比例， 为郊区人口所占比例。在n+1年的人口分布状态为： 用矩阵乘法可写成： 开始市区和郊区的人口数为 可以得到n年后市区和郊区的人口分布： 因此10年后的人口可用程序计算如下： A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; X10=A^10*X0 程序运行的结果为： 市区和郊区人口数约为：744630和255370。 无限增加时间n，市区和郊区人口之比将趋向一组常数0.25/0.75。为了弄清为什么它趋向于一个稳态值，可以将A对角化。令 ，其中Λ为对角矩阵，则有 对角矩阵的幂次可以化为元素的幂次 所以，它就很容易计算。 程序la24 % 分析n年后城市人口分布 A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; [P,lamda]=eig(A); syms n % 定义符号变量n Xn=P*lamda.^n*inv(P)*X0 % .^n对矩阵lamda中所有元素进行幂运算 计算结果为： 随n增大后一项(4/5)^n趋近于零。 直线和平面的快速绘制程序 平面曲线的快速绘制程序 ezplot(‘ ‘,[a,b]) 引号中函数可以只有一个自变量，代表显函数 ezplot(‘f(x)’, [a,b]) 系统将在 a x b的范围内画出 f = f(x) 引号中的函数若有两个自变量，那就代表隐函数，其典型格式为 ezplot(‘f(x,y)’, [a,b]) 系统将在 a x b的范围内画出 f(x,y)=0。 [a,b]的默认值为[-2π, 2π] 曲线快速绘制举例 ezplot(‘x1+0.2*x2^3+1’) 画出 在x1=[-2π, -2π]的曲线 画多条曲线可按下列方法编程 s1=x1+0.2*x2^3+1 % 方程1 s2=3*x1+2*x2+3 % 方程2 ezplot(s1),hold on % 画方程1,保持 ezplot(s2),grid on % 画方程2,加网格 [x1,x2]=solve(s1,s2) % 解联立方程1,2 解的结果 方程及方程组的几何意义 二元一次方程在几何上表示的是一根直线，则两个二元一次方程组在几何上则表示两根直线的位置关系： 相交====〉有惟一解 平行====〉无解 重合====〉无穷多解 例1 求解下列三个线性方程组 （a） （b） （c） 用ezplot(s1),hold on, ezplot(s2),命令