第3章 MATLAB矩阵分析与处理.

重点
y矩阵中每一列最大的值
y向量中最大的值
最大值的位置
最大值的位置
注意：输入矩阵类型不同， 则执行的操作不同。
2.1.4 函数
因为matlab函数太多，所以要养成使用help
命令，得到有关函数的具体用法：
例：help max
2.1表达式
表达式
（即语句）：将变量、数值、函 数用操作符连接起来，就构成了表达式 。
应用：可以和其它语言程序进行数据通信。 举例：
通过MATLAB提供的函数产生矩阵
用内部函数可生成一些特殊矩阵 （函数见书上P50）
重点
通过MATLAB提供的函数产生矩阵
1、单位矩阵（
E方阵）和广义单位矩阵的
产生
重点
通过MATLAB提供的函数产生矩阵
2、随机数矩阵的产生
随机数的产生常常用在控制系统仿真以 及信号分析，是一个非常重要的手段。 MATLAB提供了很好的随机数产生函数： rand（） randn（）
A/ B A*B
1
A\B A
重点
1
*B
Matlab右除法表示形式：
C=A/B 或 C=A * i n v ( B )
Matlab左除法表示形式： C=A\B 或 C=i n v ( A ) * B
注意：只有行列式不为0的方阵才存在逆阵！！！
矩阵元素的右除、左除
a1 A a3 a2 a4
2）变量名由字母、数字和下划线构成。第一个 字母必须是英文字母。 3）有字符个数限制（版本5.0 ：最多31个字符）
2.1.2 变量

MAT
重点
（注意大小写！）
i或j: 错误：5＋j7

Matlab 仿真及其应用
主讲：陈孝敬 E-mail:chenxj9@
第3章
数学运算
主要内容：
①矩阵运算； ②矩阵元素运算；
3.1 矩阵运算
3.1.1 矩阵分析
1．向量范式定义：
x x x
1

n
k 1
xk
2 k
2

k 1 n
x
n

1/ 2


k 1
xk
向量的3种常用范数及其计算函数 在MATLAB中，求向量范数的函数为： (1) norm(V)或norm(V,2)：计算向量V的2—范数。 (2) norm(V,1)：计算向量V的1—范数。 (3) norm(V,inf)：计算向量V的∞—范数。
3.1.2 矩阵分解
矩阵分解：把矩阵分解成比较简单或对它性质比较熟悉的若干 矩阵的乘积的形式；
1．Cholesky分解： Cholesky分解是把对称正定矩阵表示成上三角矩阵的转 置与其本身的乘积，即：A=RTR，在Matlab中用函数chol 来计算Cholesky分解 例3-13 求矩阵A=pascal（4）的Cholesky分解， A=pascal(4) R=chol(A) R’*R
例3-18.求解方程组
x1 x2 3 x3 x4 1 3 x1 x2 3 x3 4 x4 4 x 5x 9 x 8x 0 2 3 4 1
解 先用Matlab函数null求出对应的齐次线性方程组的基础解 系,再利用其系数矩阵的上、下三角阵求出方程组的一个特解, 这样即可得到该方程组的通解,程序如下: >> >> >> >> >> >> A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; b=[1 4 0] ′; format rat C=null(A , ′r′); %求基础解系 [L,U]=lu(A); %A=LU,L为上三角阵,U为下三角阵 X0= U\(L\b) %用LU求出一个齐次方程的特解

2014-12-31

32/62
微分方程求解的步骤


将微分方程变换成标准型 用MATLAB描述微分方程
M-函数 入口：function dx=funmane(t,x) 匿名函数 >> f=@(t,x)[...]



求解
验证：odeset()函数
控制系统仿真与CAD 国家级精品课程
2014-12-31
2014-12-31
9/62

数值解
ห้องสมุดไป่ตู้
解析解

解析解
控制系统仿真与CAD 国家级精品课程
2014-12-31
10/62

演示：自编 funm() 求矩阵的任意函数
求
结果：左上角元素
控制系统仿真与CAD 国家级精品课程
2014-12-31
11/62
3.1 线性代数问题求解小结


线性代数很多问题可以用MATLAB语句直 接求解，和数学表示差不多一样直观 很多方法可以同时得出解析解和数值解


本节主要介绍和这门课程相关的问题 线性代数问题的MATLAB求解 代数方程求解、微分方程求解 最优化问题的求解 Laplace 变换与 z 变换问题的求解
控制系统仿真与CAD 国家级精品课程
2014-12-31
2/62
进一步学习这方面内容建议阅读


薛定宇、陈阳泉《高等应用数学问题的 MATLAB求解》（第二版），清华大学出 版社，2008。英文版：Solving Applied Mathematical Problems with MATLAB，CRC Press，2008 遇到某个MATLAB问题找不到合适的工具 箱，试在下面网址搜索

第3章数值数组及其运算习题3及解答1 要求在闭区间]2,0[ 上产生具有10个等距采样点的一维数组。

试用两种不同的指令实现。

〖目的〗●数值计算中产生自变量采样点的两个常用指令的异同。

〖解答〗%方法一t1=linspace(0,2*pi,10)%方法二t2=0:2*pi/9:2*pi %要注意采样间距的选择，如这里的2*pi/9.t1 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.2832t2 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.28322 由指令rng('default'),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。

〖目的〗●数组下标的不同描述：全下标和单下标。

●sub2ind, int2str, disp的使用。

●随机发生器的状态控制：保证随机数的可复现性。

〖解答〗rng('default')A=rand(3,5)[ri,cj]=find(A>0.5);id=sub2ind(size(A),ri,cj);ri=ri';cj=cj';disp(' ')disp('大于0.5的元素的全下标')disp(['行号 ',int2str(ri)])disp(['列号 ',int2str(cj)])disp(' ')disp('大于0.5的元素的单下标')disp(id')A =0.8147 0.9134 0.2785 0.9649 0.95720.9058 0.6324 0.5469 0.1576 0.48540.1270 0.0975 0.9575 0.9706 0.8003大于0.5的元素的全下标行号 1 2 1 2 2 3 1 3 1 3列号 1 1 2 2 3 3 4 4 5 5大于0.5的元素的单下标1 2 4 5 8 9 10 12 13 153 采用默认全局随机流，写出产生长度为1000的“等概率双位（即取-1，+1）取值的随机码”程序指令，并给出 -1码的数目。