线性代数机算(2011.5).pptVIP

线性代数MATLAB机算应用 智能控制与图像工程研究所 机电工程学院 朱娟娟 2011.5 MATLAB语言概述 MATLAB是一门计算机编程语言，专门以矩阵的方式来处理计算机数据，它把数值计算和可视化环境集成到一起，非常直观，而且提供了大量的函数，使其越来越受到人们的喜爱，工具箱越来越多，应用范围也越来越广泛。 已知 逆矩阵的多种求法 计算含有符号变量的行列式 ，并解方程D＝0。 将上例中的增广矩阵化为最简形。 需要掌握的基本运算语句 补充内容 求解（1） 4次多项式拟合 a0+a10+a202+a303+a404 =-27 a0+a11+a212+a313+a414 =0 a0+a12+a222+a323+a424 =21 （1） a0+a13+a232+a333+a434 =0 a0+a14+a242+a343+a444 =-75 方程组的系数矩阵的行列式为范德蒙行列式 编程 从图中显然可以看出，由于取点过少，显示曲线为一个折线，无法显示曲线性态，因此我们在一定范围内取更多的点，将曲线充分表示出来 MATLAB提供了利用最小二乘法解决超定方程组的近似解的方法，即通过近似求解得到一条曲线y = a0+a1x+a2x2满足 显然可以求出方程的唯一解 x=[0;1;2;3;4]; y=[-27;0;21;0;-75]; A=[x.^0,x.^1,x.^2,x.^3,x.^4]; a=A\y; 输入列向量x，y 将列向量x并排放入矩阵，构成方阵，对方阵中的元素做群运算，生成范德蒙矩阵 用矩阵左除法 求出系数向量α figure(1) plot(x,y,r) 作图figure(1)，其中x、y分别为横坐标和纵坐标；’r’表示曲线颜色为红色 x1=linspace(-1,6,100); 在区间[-1，6]之间均匀取100个点，赋值到数组x1 y1=a(1)+a(2)*x1+a(3)*x1.^2+a(4)*x1.^3+a(5)*x1.^4 令x2=5，计算x2对应的y值—y2 x2=5 y2=a(1)+a(2)*x2+a(3)*x2.^2+a(4)*x2.^3+a(5)*x2.^4 根据算出的四阶多项式计算x1对应的y值 figure(2) plot(x,y,r,x,y,o) 图figure(1)描述x-y曲线，’r’表示曲线颜色为红色，‘o’表示用符号o表示各点位置 hold on 在figure(2)上继续叠加下面图形 plot(x1,y1,b,x2,y2,*) 在figure(1)中描述x1-y1曲线，’b’表示曲线颜色为蓝色，用‘*’表示点(x2,y2) grid on 在figure(2)上显示网格 经过上述运算，得到如下图形，该图形显示了过给定各点的拟合曲线 求解（2） 2次多项式逼近 将表中数值代入 y = p4(x) = a0+a1x+a2x2中，得到线性方程组 a0+a10+a202=-27 a0+a11+a212 =0 a0+a12+a222 =21 （2） a0+a13+a232 =0 a0+a14+a242 =-75 很容易判断出方程组无解—约束过多---超定方程组。 编程 x=[0;1;2;3;4]; y=[-27;0;21;0;-75]; A=[x.^0,x.^1,x.^2]; a=A\y; %最小二乘解 经过上述运算，得到如下图形，该图形显示采用最小二乘法近似得到的二次曲线 Happy Every Day * * MATLAB － MATrix LABoratory MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。 MATLAB和著名的符号计算语言Maple相结合，使得MATLAB具有符号计算功能。 MATLAB提供了两个层次的绘图操作：一种是对图形句柄进行的低层绘图操作，另一种是建立在低层绘图操作之上的高层绘图操作。 MATLAB的主要功能 MATLAB语言及其发展历程 数值计算和符号计算功能 绘 图 功 能 编 程 语 言 MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。 MATLAB 工 具 箱 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。 MATLAB语言的特点 功能强大 在MATLAB环境中，有超过500种数学、统计、科学及工程方面的函数；不仅在数值计算有绝对优势