Matlab学习教程.ppt
《Matlab学习教程.ppt》由会员分享,可在线阅读,更多相关《Matlab学习教程.ppt(229页珍藏版)》请在冰豆网上搜索。
![Matlab学习教程.ppt](https://file1.bdocx.com/fileroot1/2022-10/22/c8113c71-7c8f-489e-9770-90c35400c9bf/c8113c71-7c8f-489e-9770-90c35400c9bf1.gif)
数学实验与数学实验与Matlabhttp:
/晓晓阳阳华中科技大学数学系实实验验一一矩阵运算与矩阵运算与Matlab命令命令日常矩阵及其运算日常矩阵及其运算矩阵应用实例矩阵应用实例:
榄球防护用品的生产管理榄球防护用品的生产管理应用问题应用问题一个工厂生产三种橄榄球用品一个工厂生产三种橄榄球用品:
防护防护帽、帽、垫肩、臀垫。
垫肩、臀垫。
需要不同数量的:
硬塑料需要不同数量的:
硬塑料、泡沫塑料泡沫塑料尼龙线尼龙线、劳动力。
劳动力。
为监控生产,管理者对它们之间的关系为监控生产,管理者对它们之间的关系十分关心。
十分关心。
为把握这些量的关系,他列出下面的表为把握这些量的关系,他列出下面的表原料产品关系表原料产品关系表订单订单管理者接到四份订单如上表所示。
管理者接到四份订单如上表所示。
问应该如何计算每份订单所需的原材料,问应该如何计算每份订单所需的原材料,以便组织生产?
以便组织生产?
将将表表格格写写成成矩矩阵阵形形式式计计算算输入下面输入下面Matlab指令指令A=423;132;133;322,B=35206045;10155040;20124520C=A*B请自行计算观看结果请自行计算观看结果Matlab基本指令基本指令向量的创建和运算向量的创建和运算1.直接输入向量直接输入向量x1=124,x2=1,2,1,x3=x1运行结果运行结果x1=124x2=121x3=1242.冒号创建向量冒号创建向量x1=3.4:
6.7,x2=3.4:
2:
6.7x3=2.6:
-0.8:
0运算结果运算结果x1=3.40004.40005.40006.4000x2=3.40005.4000x3=2.60001.80001.00000.20003.生成线性等分向量生成线性等分向量指令指令x=linspace(a,b,n)在在a,b区间产区间产生生n个等分点个等分点(包括端点包括端点)x=linspace(0,1,5)结果结果x=00.25000.50000.75001.0000工作空间工作空间在在Matlab窗口创建向量后并运行后,向窗口创建向量后并运行后,向量就存在于工作空间,可以被调用。
量就存在于工作空间,可以被调用。
向量的运算向量的运算设设x=x1x2x3;y=y1y2y3;为为两两个个三三维向量,维向量,a,b为标量。
为标量。
向量的数乘:
向量的数乘:
a*x=a*x1a*x2a*x3向量的平移:
向量的平移:
x+b=x1+bx2+bx3+b向量和:
向量和:
x+y=x1+y1x2+y2x3+y3向量差:
向量差:
x-y=x1-y1x2-y2x3-y3数的乘幂:
数的乘幂:
如如a2元素群运算元素群运算(四则运算)四则运算)x.*y=x1*y1x2*y2x3*y3(元元素素群群乘乘积积)x./y=x1/y1x2/y2x3/y3(元元素素群群右右除除,右边的右边的y做分母做分母)x.y=y1/x1y2/x2y3/x3(元元素素群群左左除除,左边的左边的x做分母做分母)x.5=x15x25x35(元素群乘幂元素群乘幂)2.x=2x12x22x3(元元素素群群乘乘幂幂)x.y=x1y1x2y2x3y3(元元素素群群乘乘幂幂)元素群运算(函数计算)元素群运算(函数计算)Matlab有有许许多多内内部部函函数数,可可直直接接作作用用于于向量产生一个同维的函数向量。
向量产生一个同维的函数向量。
x=linspace(0,4*pi,100);(产产生生100维向量维向量x)y=sin(x);(y也自动为也自动为100维向量维向量)y1=sin(x).2;y2=exp(-x).*sin(x);观察结果观察结果创建矩阵(数值矩阵的创建)创建矩阵(数值矩阵的创建)直接输入法创建简单矩阵。
直接输入法创建简单矩阵。
A=1234;5678;9101112B=-1.3,sqrt(3);(1+2)*4/5,sin(5);exp
(2),6观察运行结果观察运行结果创建矩阵(符号矩阵的创建)创建矩阵(符号矩阵的创建)用指令用指令“syms”说明符号变量。
说明符号变量。
symsa11a12a13a14a21a22a23a24a31a32a33a34b11b12b13b14b21b22b23b24b31b32b33b34A1=a11a12a13a14;a21a22a23a24;a31a32a33a34,B1=b11b12b13b14;b21b22b23b24;b31b32b33b34运行运行矩阵的运算矩阵的运算(矩阵的加减、数乘、乘积矩阵的加减、数乘、乘积)C=A1+B1D=A1-B1symsc,cA=c*A1A2=A1(:
1:
3),B1G=A2*B1矩阵的运算矩阵的运算(矩阵的加减、数乘、乘积矩阵的加减、数乘、乘积)A,A_trans=AH=123;210;123,K=123;210;231h_det=det(H),k_det=det(K),H_inv=inv(H),K_inv=K-1矩阵的运算矩阵的运算(左除和右除左除和右除)左除左除“”:
求矩阵方程求矩阵方程AX=B的解;(的解;(A、B的行要保持一致)的行要保持一致)解为解为X=AB;当当A为方阵且可逆时有为方阵且可逆时有X=AB=inv(A)*B;右除右除“/”:
求矩阵方程求矩阵方程XA=B的解的解(A、B的列要保持一致)的列要保持一致)解为解为X=B/A,当当A为方阵且可逆时有为方阵且可逆时有X=B/A=B*inv(A)矩阵的运算矩阵的运算(左除和右除左除和右除)求矩阵方程:
求矩阵方程:
设设A、B满足关系式:
满足关系式:
AB2B+A,求求B。
其中其中A=301;110;014。
解:
有解:
有(A-2I)BA程序程序:
A=301;110;014;B=inv(A-2*eye(3)*A,B=(A-2*eye(3)A观察结果:
观察结果:
分块矩阵分块矩阵(矩阵的标识矩阵的标识)1.矩阵元素的标识矩阵元素的标识:
A(i,j)表示矩阵表示矩阵A的第的第i行行j列的元素;列的元素;2.向量标识方式向量标识方式A(vr,vc):
vr=i1,i2,ik、vc=j1,j2,ju分分别别是是含含有矩阵有矩阵A的行号和列号的单调向量。
的行号和列号的单调向量。
A(vr,vc)是是取取出出矩矩阵阵A的的第第i1,i2,ik行行与与j1,j2,ju列交叉处的元素所构成新矩阵。
列交叉处的元素所构成新矩阵。
分块矩阵分块矩阵(矩阵的标识矩阵的标识)取取出出A的的1、3行行和和1、3列列的的交交叉叉处处元元素素构成新矩阵构成新矩阵A1程序程序nA=10112;01-123;n30510;23121,nvr=1,3;vc=1,3;nA1=A(vr,vc)观察结果观察结果分块矩阵分块矩阵(矩阵的标识矩阵的标识)将将A分分为为四四块块,并并把把它它们们赋赋值值到到矩矩阵阵B中中,观察运行后的结果。
观察运行后的结果。
程序程序nA11=A(1:
2,1:
2),A12=A(1:
2,3:
5),nA21=A(3:
4,1:
2),A22=A(3:
4,3:
5)nB=A11A12;A21A22结果结果分块矩阵分块矩阵(矩阵的修改和提取)(矩阵的修改和提取)修改矩阵修改矩阵A,将它的第将它的第1行变为行变为0。
程序:
程序:
nA=10112;01-123;n30510;23121nA(1,:
)=00000;A删除上面矩阵删除上面矩阵A的第的第1、3行。
行。
程序:
程序:
nA(1,3,:
)=结果结果生成特殊矩阵生成特殊矩阵全全1阵阵nones(n),ones(m,n),ones(size(A)全零阵:
全零阵:
nzeros(n),zeros(m,n),zeros(size(A)n常常用于对某个矩阵或向量赋常常用于对某个矩阵或向量赋0初值初值单位阵:
单位阵:
neye(n),eye(m,n)随机阵:
随机阵:
nrand(m,n),rand(n)=rand(n,n)用用于于随随机机模模拟,常和拟,常和rand(seed,k)配合使用配合使用。
生成特殊矩阵生成特殊矩阵将将rand指令运行多次,观察结果。
指令运行多次,观察结果。
程序:
程序:
ny1=rand(1,5),y2=rand(1,5),nrand(seed,3),x1=rand(1,5),nrand(seed,3),x2=rand(1,5)结果结果常用矩阵函数常用矩阵函数det(A):
方阵的行列式;方阵的行列式;rank(A):
矩阵的秩;矩阵的秩;eig(A):
方阵的特征值和特征向量;方阵的特征值和特征向量;trace(A):
矩阵的迹;矩阵的迹;rref(A):
初等变换阶梯化矩阵初等变换阶梯化矩阵Asvd(A):
矩阵奇异值分解。
矩阵奇异值分解。
cond(A):
矩阵的条件数;矩阵的条件数;数据的简单分析数据的简单分析1.当当数数据据为为行行向向量量或或列列向向量量时时,函函数数对对整个向量进行计算整个向量进行计算.2.当当数数据据为为矩矩阵阵时时,命命令令对对列列进进行行计计算算,即即把把每每一一列列数数据据当当成成同同一一变变量量的的不不同同观观察值。
察值。
max(求求最最大大)、min(求求最最小小)、mean(求求平平均均值值)、sum(求求和和)、std(求求标标准准差差)、cumsum(求求累累积积和和)、median(求求中中值值)、diff(差差分分)、sort(升升序序排排列列)、sortrows(行升序排列行升序排列)等等等等。
数据的简单分析数据的简单分析观观察察:
生生成成一一个个36的的随随机机数数矩矩阵阵,并并将将其其各各列列排排序序、求求各各列列的的最最大大值值与与各各列列元素之和。
元素之和。
程序程序nrand(seed,1);A=rand(3,6),nAsort=sort(A),Amax=max(A),Asum=sum(A)结果结果实验二实验二函数可视化与函数可视化与Matlab作图作图函数的可视化函数的可视化nf(x),g(x)是周期函数吗?
观察它们的图象。
是周期函数吗?
观察它们的图象。
n程序程序clf,x=linspace(0,8*pi,100);F=inline(sin(x+cos(x+sin(x);y1=sin(x+cos(x+sin(x);y2=0.2*x+sin(x+cos(x+sin(x);plot(x,y1,k:
x,y2,k-)legend(sin(x+cos(x+sin(x),0.2x+sin(x+cos(x+sin(x),2)n令令绘制平面曲线绘制平面曲线(plot指令指令)plot(x,y):
n以以x为横坐标、为横坐标、y为纵坐标绘制二维图形为纵坐标绘制二维图形nx,y是同维数的向量;是同维数的向量;plot(y):
n相当于相当于x=1,2,length(y)时情形。
时情形。
绘制平面曲线绘制平面曲线(绘制多个图形(绘制多个图形)1.plot(x,y1;y2;),nx是是横横坐坐标标向向量量,y1;y2;是是由由若若干干函函数数的纵坐标拼成的矩阵的纵坐标拼成的矩阵2.plot(x,y1),holdon,plot(x,y2),holdoff3.plot(x,y1,x,y2,)4.plotyyn两个坐标系,用于绘制不同尺度的函数。
两个坐标系,用于绘制不同尺度的函数。
绘制平面曲线绘制平面曲线(线型、点形和颜色的控制)(线型、点形和颜色的控制)(线型、点形和颜色的控制)(线型、点形和颜色的控制)plot(x,y,颜色线型点形颜色线型点形)plot(x,y,颜颜色色线线型型点点形形,x,y,颜色线型点形颜色线型点形,)句句柄柄图图形形和和set命命令令改改变变