线性代数实验演示课课件.ppt
线性代数实验演示课第1页，此课件共47页哦画图画图矩阵与行列式运算矩阵与行列式运算简单方程组求解简单方程组求解123本次课程你将学会本次课程你将学会第2页，此课件共47页哦一、演示一、演示MATLAB的基本操作的基本操作 打开模块介绍帮助文件关闭第3页，此课件共47页哦二、数字简单计算二、数字简单计算加加+减减-乘乘*除除/右除右除 正常除法正常除法 左除左除幂幂第4页，此课件共47页哦三、常用数学函数三、常用数学函数abs(x)求求x的的绝对值绝对值sqrt(x)求求x的平方根的平方根exp(x)求求x的指数函数的指数函数sin(x),cos(x)求求x的正的正、余、余函数，函数，x为为弧度弧度tan(x),cot(x)求求x的的正、余切正、余切，x为为弧度弧度Note注意括号的配对；注意括号的配对；pi为系统变量，圆周率为系统变量，圆周率第5页，此课件共47页哦ans预设结果预设结果eps小误差小误差=2.2204e-016，可以理解为微积分中的，可以理解为微积分中的，16位位inf无穷大无穷大NaNnot a number,类似未定型（类似未定型（0/0）i，j默认虚数单位，默认虚数单位，i=j=0+1.0000i表表1 MATLAB预设变量预设变量MATLABMATLAB变量命名规则变量命名规则变量命名规则变量命名规则1.1.区分大小写区分大小写区分大小写区分大小写 2.2.第一个字符为字母第一个字符为字母第一个字符为字母第一个字符为字母3.3.可以包含下划线和数字，不能为标点和空格可以包含下划线和数字，不能为标点和空格可以包含下划线和数字，不能为标点和空格可以包含下划线和数字，不能为标点和空格第6页，此课件共47页哦2分钟完成第一次预习分钟完成第一次预习报告第二题报告第二题2第7页，此课件共47页哦四、如何赋值四、如何赋值 x=18 x=18 y=3*x2-3 y=3*x2-3 u=x+y;u=x+y;%分号含义分号含义分号含义分号含义 v=x-y;v=x-y;tan(2*u/3*v)tan(2*u/3*v)Note：注意括号和标点：注意括号和标点第8页，此课件共47页哦 A=1 2 3 4 5 五、如何输入矩阵五、如何输入矩阵 B=1;2;3;4 C=3 0 2 1;-1 4 5 2;3 5 8 7 who clear%;%;换行符换行符%查看当前变量查看当前变量%清除当前变量清除当前变量第9页，此课件共47页哦六、群运算六、群运算例如：一次求出例如：一次求出 A=0,0.25*pi,0.5*pi,0.75*pi,pi;sin(A)ans=0 0.7071 1.0000 0.7071 0.0000第10页，此课件共47页哦七、如何构造一个等差数列七、如何构造一个等差数列从初值从初值0，公差，公差0.25*pi，到终值，到终值pi结束结束 A=0:0.25*pi:pi%初值：增量：终值初值：增量：终值，A=linspace(0,pi,5)%起点，终点，等分份数起点，终点，等分份数 第11页，此课件共47页哦八、如何取出矩阵的某一行（列）八、如何取出矩阵的某一行（列）A=1,2,3;4,5,6 A=1,2,3;4,5,6 A(:,1)A(:,1)%A%A的第一列元素的第一列元素的第一列元素的第一列元素 A(2,:)%A%A的第二行元素的第二行元素的第二行元素的第二行元素第12页，此课件共47页哦总结标点符号在矩阵输入中总结标点符号在矩阵输入中的作用，完成第二题的作用，完成第二题3（3分分钟）钟）第13页，此课件共47页哦实验一实验一常见图形绘制常见图形绘制第14页，此课件共47页哦1.具有解析表达式的二维曲线具有解析表达式的二维曲线2.隐函数确定的二维曲线隐函数确定的二维曲线3.参数形式的二维曲线参数形式的二维曲线摆线摆线4.极坐标下的二维曲线，心脏线、星形线等极坐标下的二维曲线，心脏线、星形线等第15页，此课件共47页哦 X=linspace(0,2*pi,30);X=linspace(0,2*pi,30);Y=sin(X);Y=sin(X);plot(X,Y)plot(X,Y)%以以X X向量为横坐标，向量为横坐标，Y Y向量纵坐标绘图向量纵坐标绘图 xlabel(X)xlabel(X)%给给X X轴，轴，Y Y轴添加标题轴添加标题 ylabel(Y)ylabel(Y)title(Y=sinX)title(Y=sinX)%给图添加标题给图添加标题 grid on grid on%开启图像网格开启图像网格方法：描点作图方法：描点作图例例1.1.1 绘制绘制 的图像的图像.第16页，此课件共47页哦 ezplot(x2,-1,1)例例1.1.2 绘制绘制 的图像的图像第17页，此课件共47页哦 ezplot(x2-y4,-6,6,-4,4)例例1.1.3 绘制绘制 的图像的图像第18页，此课件共47页哦例例1.1.4 绘制绘制 的图像。的图像。ezplot(2*(t-sin(t),2*(1-cos(t),0,10*pi)ezplot(2*(t-sin(t),2*(1-cos(t),0,10*pi)axis(0 60-10 10)axis(0 60-10 10)第19页，此课件共47页哦v ezplot 可用于显示函数、隐函数、可用于显示函数、隐函数、参参数方程表达的曲线作图数方程表达的曲线作图v plot(X,Y)常用于描点作图的方法常用于描点作图的方法总总 结结第20页，此课件共47页哦 subplot(1,2,1)subplot(1,2,1)t=0:0.01:2*pi;t=0:0.01:2*pi;polar(t,1-cos(t),-r)polar(t,1-cos(t),-r)subplot(1,2,2)subplot(1,2,2)polar(t,t,-g)polar(t,t,-g)例例1.1.7 极坐标绘图极坐标绘图1)心脏线心脏线2)阿基米德螺线阿基米德螺线第21页，此课件共47页哦1.2.1 实验目的实验目的J 掌握掌握MATLAB软件绘制常见的三维曲线图形的软件绘制常见的三维曲线图形的 方法；方法；J 通过通过MATLAB绘制常见图形，加深对空间解析绘制常见图形，加深对空间解析 几何相关部分内容的理解几何相关部分内容的理解.第22页，此课件共47页哦第23页，此课件共47页哦 t=0:pi/50:10*pi;plot3(sin(t),cos(t),t);grid on axis square xlabel(X)ylabel(Y);zlabel(Z);例例1.2.1 绘制三维螺旋线绘制三维螺旋线第24页，此课件共47页哦 x,y=meshgrid(-8:.4:8);c=sqrt(x.2+y.2)+eps;z=sin(c)./c;mesh(x,y,z)axis square例例1.2.2 绘制墨西哥帽子图形绘制墨西哥帽子图形第25页，此课件共47页哦 x,y=meshgrid(-25:1:25,-25:1:25);z=x.2/9-y.2/4;mesh(x,y,z)例例1.2.3 绘制马鞍面绘制马鞍面第26页，此课件共47页哦 ezmesh(x2/2+y2/4)ezsurf(x2/2+y2/4)例例1.2.4 绘制椭圆抛物面绘制椭圆抛物面第27页，此课件共47页哦 sphere(50)例例1.2.5 绘制单位球面绘制单位球面第28页，此课件共47页哦 x,y,z=ellipsoid(1,0,-1,4,3,2);x,y,z=ellipsoid(1,0,-1,4,3,2);surfl(x,y,z)surfl(x,y,z)colormap copper colormap copper axis equal axis equal例例1.2.6 绘制椭球面绘制椭球面第29页，此课件共47页哦 t,z=meshgrid(0:pi/20:2*pi,-10:.5:10);t,z=meshgrid(0:pi/20:2*pi,-10:.5:10);x=2*sqrt(z.2/16+1)*cos(t);x=2*sqrt(z.2/16+1)*cos(t);y=3*sqrt(z.2/16+1)*sin(t);y=3*sqrt(z.2/16+1)*sin(t);mesh(x,y,z)mesh(x,y,z)例例1.2.7 绘制单叶双曲面绘制单叶双曲面第30页，此课件共47页哦 ezmesh(X+Y+1)例例1.2.8 绘制平面绘制平面第31页，此课件共47页哦例例1.2.9 双叶双曲面绘制双叶双曲面绘制参数方程：参数方程：第32页，此课件共47页哦例如：取例如：取ezmesh(2*tan(u)*cos(v),4*tan(u)*sin(v),4*sec(u),0,2*pi,0,pi)第33页，此课件共47页哦实验二实验二矩阵基本运算、行列式与方程组求解矩阵基本运算、行列式与方程组求解第34页，此课件共47页哦一、矩阵的输入一、矩阵的输入 eye(3)zeros(2,3)zeros(2,3)ones(2,3)ones(2,3)2.构造构造3阶单位阵阶单位阵1.输入输入3.构造构造 阶元素全为阶元素全为1的矩阵的矩阵4.构造构造 阶零矩阵阶零矩阵第35页，此课件共47页哦6.构造元素为构造元素为0-10的的 阶随机矩阵阶随机矩阵5.构造构造 阶随机矩阵阶随机矩阵rand(2,3)round(10*rand(2,3)7.构造范德蒙矩阵构造范德蒙矩阵 rot90(vander(3,4,5,6)第36页，此课件共47页哦二、矩阵的运算二、矩阵的运算 3.方阵方阵A的逆的逆inv(A)4.方阵方阵A的行列式的行列式det(A)1.矩阵矩阵A的秩的秩rank(A)2.矩阵矩阵A的转置的转置A第37页，此课件共47页哦 A=rand(3,3)B=round(rand(3,3)*100)det(A)det(B)例例2.1.2 生成元素为生成元素为0到到1间的随机矩阵间的随机矩阵 生成元素生成元素为为0到到100并且元素全是随机整数的矩阵并且元素全是随机整数的矩阵(1)求出矩阵求出矩阵 A 与与 B 的行列式的行列式.第38页，此课件共47页哦(2)求出矩阵求出矩阵A的秩；的秩；rank(A)(3)验证验证 第39页，此课件共47页哦三、行列式及其几何应用三、行列式及其几何应用*第40页，此课件共47页哦例例2.2.2 如图如图所示所示，求向，求向u=1,2,3,v=3,1,0,w=0,5,1,所构成的所构成的四面体的体积。四面体的体积。u=1 2 3;v=3 1 0;w=0 5 1;A=u;v;w;v=abs(det(A)/6 第41页，此课件共47页哦四、方程组简单求解四、方程组简单求解1.回顾利用初等变换求解方程组过程回顾利用初等变换求解方程组过程行阶梯型行阶梯型(有没有解有没有解)行最简型行最简型(解是什么解是什么)化矩阵化矩阵A为行最简型为行最简型 rref(A)第42页，此课件共47页哦例例2.3.1 利用初等变换求线性方程组利用初等变换求线性方程组 a1=1,3,-1,2;a1=1,3,-1,2;a2=3,4,-5,7;a3=1,2,4,1;a3=1,2,4,1;a4=2-3 1-6;A=a1 a2 a3 a4 A=a1 a2 a3 a4 b=1 2 3 4;b=1 2 3 4;rref(A,b)rref(A,b)第43页，此课件共47页哦2.可逆矩阵求解方程组可逆矩阵求解方程组化化A可逆可逆 X=inv(A)*b第44页，此课件共47页哦X=inv(A)*b例例2.3.1 利用求逆阵方法求线性方程组利用求逆阵方法求线性方程组第45页，此课件共47页哦3.超定方程组求解超定方程组求解(过过过过)超定：方程个数超定：方程个数比未知数多比未知数多欠定：方程个数比未知数欠定：方程个数比未知数欠定：方程个数比未知数欠定：方程个数比未知数少少少少第46页，此课件共47页哦最小二乘解最小二乘解XY 点虽然不是三条直线的交点，但到每条直线的距离平方点虽然不是三条直线的交点，但到每条直线的距离平方点虽然不是三条直线的交点，但到每条直线的距离平方点虽然不是三条直线的交点，但到每条直线的距离平方和最小和最小和最小和最小.第47页，此课件共47页哦
