第二章MATLAB语言的数值运算(沈阳航空航天大学).pptVIP

MATLAB基本运算符 系统特殊变量 常用数学函数 第二章MATLAB语言的数值运算 教学目标 介绍MATLAB的两种基本的数值运算，矩阵和表达式 学习要求 掌握矩阵和多项式的构造和运算方法，能够使用常用的几种函数进行简单的问题求解 2.1 矩阵 2.2 向量 2.3 数组 2.4 多项式 2.1 矩 阵 矩阵是线性代数的基本运算单元，通常矩阵是指含有M行N列的数值和矩阵结构，每个元素可以是实数或复数，在线性代数中我们学习了矩阵的基本性质，了解了矩阵的基本运算，如：加、减、内积、逆矩阵、转置、线性方程式等。 MATLAB支持线性代数所定义的全部矩阵运算。 矩阵的信息获取 MATLAB语言的数值运算 MATLAB = matrix（矩阵）+ laboratory（实验室） 2.1.1矩阵的构造 直接输入法 三个符号：“[]”、“，”、“；” 用中括号[ ]把所有矩阵元素括起来 同一行的不同数据元素之间用空格或逗号间隔 用分号（；）指定一行结束 可分成几行进行输入，用回车符代替分号 数据元素可以是表达式，系统将自动计算结果，可以是实数 ，也可以是复数，复数用i，j 输入 矩阵的构造 A=[1 2 3 4; 5 6 7 8; 9 10 11 12; 13 14 15 16] B=[1,sqrt(25),9,13 2,6,10,7*2 3+sin(pi),7,11,15; 4,abs(-8),12,16] 由函数创建和修改矩阵 由函数创建和修改矩阵 MATLAB提供了大量的函数用于创建一些特殊的矩阵及其派生矩阵。 ① zeros生成全部元素为0的零距阵 (Zeros Array) A=zeros(n) 生成n×n零矩阵。 A=zeros(m,n)或者A=zeros([m,n]) 生成m×n的零矩阵。 A=zeros(m,n,p,…)或者A=zeros([m n p …]) 生成m×n×p×…的零矩阵。 B=zeros(size(A)) 生成和矩阵A大小相等的全零矩阵。 ② ones生成全部元素为1的距阵 A=ones (n) 生成n×n全1矩阵。 A=ones (m,n)或者A=ones ([m,n]) 生成m×n的全1矩阵。 A=ones (m,n,p,…)或者A=ones ([m n p …]) 生成m×n×p×…的全1矩阵。 B=ones (size(A)) 生成和矩阵A大小相等的全1矩阵。 ③ eye生成单位距阵 A=eye (n) 生成n×n单位矩阵。 A=eye (m,n)或者A=eye ([m,n]) 生成m×n的单位矩阵。 B=eye (size(A)) 生成和矩阵A大小相等的单位矩阵。 ④ rand生成均匀分布的随机矩阵 A=rand (n) 生成n×n随机矩阵。 A=rand (m,n)或者A=rand ([m,n]) 生成m×n的随机矩阵。 A=rand (m,n,p,…)或者A=rand ([m n p …]) 生成m×n×p×…的随机矩阵。 B=rand (size(A)) 生成和矩阵A大小相等的随机矩阵。 A=rand 不带任何参数将产生一个随机数。 提示：rand函数产生一个矩阵元素在0和1之间均匀分布的随机数的随机矩阵。 ⑤diag 生成一个对角阵或由对角线元素组成的向量 A＝diag(V) 当V为n维向量时,产生一个以向量V的元素为对角线的n维数组 A＝diag(V) 当V为n维矩阵时，产生 一个以V矩阵的主对角线元素为元素的n 维数组。 tril: 抽取主下三角; triu:抽取主上三角 矩阵的变向 rot90:旋转; fliplr:上翻; flipud:下翻 注意：在MATLAB中，不需要事先定义矩阵的维数，MATLAB自动为矩阵分配存 储空间。但如果在程序运行过程中采用零矩阵为矩阵生成的全部元素，或某一行、某一列的元素预先分配内存空间，将会大大加快MATLAB程序的运算速度。 2.1.2 矩阵元素 用下标来表示矩阵元素，同时可以利用下标对矩阵元素进行修改 A(m, n) 提取第m行，第n列元素 A(:, n) 提取第n列元素 A(m, :) 提取第m行元素 A(m1:m2, n1:n2) 提取第m1行到第m2行和第n1列 到第n2列的所有元素 A(m:end, n) 提取从第m行到最末行和第n列的子块 例： 修改矩阵A中元素的数值 A=[1 2 3 4; 5 6 7 8; 9 10 11 12; 13 14 15 16]; A(1,1)=0;A(2,2)=A(1,2)+A(2,1);A(4,4)=cos(0); 则矩阵变为： A = 0