数模讲座3.ppt

说明 1) 输入时, 矩阵同一行元素之间用逗号或空格分隔? 2) 不同行之间用分号分隔? 3) 所有的矩阵元素必须放置在一对方括号之? 4) 命令行的最后如果以分号结束时，表示要求系统不显示所输入的内容? 如果命令行的最后没有分号, 按回车键之后, 系统将立即显示所输入的内容。 创建M文件 对于大型矩阵，为便于修改，一般创建M文件: 特殊矩阵的生成 1) 全零阵 格式: B = zeros (n)? %生成n×n全零阵 B = zeros(m,n)? %生成m×n全零阵 B = zeros (size(A))? %生成与A相同大小的全零阵 2) 单位阵 格式: Y = eye(n)? %生成n×n单位阵 Y = eye (m,n)? %生成m×n单位阵 Y = eye (size(A))? %生成与A相同大小的单位阵 特殊矩阵的生成 3) 元素全为1的矩阵 格式： Y=ones(n)? Y=ones(m,n)? Y=ones(size(A))? 特殊矩阵的生成 4)均匀分布随机矩阵 函数: rand( ) 5) 正态分布随机矩阵 函数: randn( ) 6) 随机排列 p = randperm(n)? %产生1～n之间整数的随机排列. 7)线性等分向量 y = linspace(a,b)? %在(a, b)上产生100个线性等分点. y = linspace(a,b,n)?%在(a, b)上产生n个线性等分点. 特殊矩阵的生成 矩阵的缩维: 假设A是一个10×10的矩阵： (1)提取A的1行: B=A(5,:)?%提取A的第5行,所有列? (2)提取A的1列: C=A(:,4)?%提取A的第4列,所有行? (3)提取A的子矩阵: D=A(1:5,6:9)? %提取A的第1~5行,6~9列? E=A([1,3,5],[2,4,8])?%提取A的1,3,5行,2,4,8列? F=A(:,7:-1:3)?%提取A的所有行的7~3列,-1为步长? (4)X=[ ]表示X为一个空矩阵, A(:,[2,4])=[ ]? %删去A的2,4列后所成的子阵. 说明：矩阵A是已存在的矩阵，缩维是对A进行操作。 矩阵的扩维 假设AA=[1,2?3,4]BB=[2,3?4,5]CC=[3,4?5,6?7,8], DD=[9,8?7,6?10,11]? (1)扩展: 则A=[AA,BB?CC,DD] 为5×4的矩阵： (2)若B=[1?2?3?4?5], 则 C=[A,B]为: 注意：MATLAB字母 区分大小写。 矩阵运算 加法“+” 减法“-” 乘法“*” 乘方“^” 除法(左除右除) “/” “\” 转置和求逆。 矩阵运算 矩阵运算 Matlab提供了两种除法运算:左除“\”和右除“/”。一般情况下，x=A\b是方程Ax =b的解? 而x=b/A是方程xA=b 的解,A总是在分母的位置。 说明 若A与B为同型阵列时，A+B、A-B分别对对应分量进行加减? 若A与B中至少有一个为标量, 则把标量扩大为与另外一个同型的阵列，再按对应的分量进行加减? A.*B为符号数组的乘法。按参量A与B对应的分量相乘。A与B必须为同型阵列或至少有一个为标量? A.\B为数组或向量左除法, A.\B按对应分量相除?B./A为数组右除法, B./A为按对应的分量相除? A^n为矩阵的方幂，计算矩阵A的整数n次方幂?B.^n为数组或向量的方幂,计算B的整数n次方幂? Matlab二维绘图 函数plot(y) 绘制二维图形的基本命令是plot，其函数格式主要为plot (y),如果y是向量，则以y的每个元素为纵坐标，以对应元素的下标（自然数顺序）为横坐标的数对所连接的折线图； 如果y是m×n矩阵，则按列绘制连线图，共绘制出n条连线。其中第j条连线是以第j列的每个元素对应的行下标为横坐标，第 j 列的每个元素的元素值为纵坐标绘制的。 先输入矩阵A，然后在命令输入符“”后输入plot（A），再按回车键，系统即绘制出矩阵A对应的二维图形. 比如, 绘制矩阵A的二维图形: Matlab二维绘图 函数plot (x, y) 1)如果x, y是维数相同的两个向量,该命令以x, y的对应元素为横、纵坐标绘制连线图，这种方法常常用于一元函数平面图形的绘制? 2)命令中x被看成是自变量，y是函数,当自变量在定义域内取一系列值(列向量)时,得到一系列函数值(也是列向量),再利用该命令即可绘制二维图形? 3)当x是n维列向量, y是 m×n 阶矩阵时, 系统将绘制出m条折线(曲线), x是这些曲线的横坐标? 自我介绍： 姓名：梅启培 学院：土木工程与力学学院 班级：土木0701 所属团队：数模基地 全国大学生数学建模 全国大学生数学建模竞赛创办于1992年，每年一届，目前已成为全