第六讲 线性方程组和数据分析 Matlab语言程序设计 教学课件.pptVIP

二、数据分析 例10： 使用cov函数求解协方差。 A=rand(5,1); B=randn(5,1); C1=cov(A) C2=cov(A,B) X=magic(4) C3=cov(X) 二、数据分析 7. 排序 s o r t ( x ) 返回一个向量x的元素按递增排序的向量。如果元素是复数，则使用绝对值进行排序，即s o r t ( a b s ( x ) )。 [ y , i n d ] = s o r t ( x ) 返回下标向量i n d。就是y = x ( i n d )。另外向量y是x中元素按递增排序得到的。 sort(A) 对A中各列按递增排序，注意矩阵的行已被改变。 二、数据分析 例12： A=[17 3 4;3 1 12;4 12 8] sort(A) sort(A,2) [Y,I]=sort(A) 二、数据分析 按行进行数据排序 使用sortrows函数按行进行数据的排序。 Y=sortrows(X)命令将X按由小到大排序以行的方式排序。X为矩阵，返回矩阵Y，并且Y是按X的第1列由小到大，以行方式排序后生成的矩阵。 Sortrows(X,col)命令按指定列col由小到大进行排序。 [Y,I]=sortrows(X,col)命令中，Y为排序的结果，I表示第col列元素在x中位置。 二、数据分析 例13： 使用sortrows函数按行进行排序。 A=[1 7 5;6 4 2;9 8 3] sortrows(A) sortrows(A,3) sortrows(A,[2 1]) [Y,I]sortrows(A,3) 二、数据分析 最大值与最小值之差 使用函数range来求数据中最大值与最小值之差。 Y=range(X)命令求出X中最大值与最小值的差。当X为向量时，该命令返回X向量中的最大值与最小值之差；当X为矩阵时，返回X矩阵中各列元素的最大值与最小值这差。 例14 使用rangs函数求向量A以及矩阵B中向量的最大值与最小值之差。 A=[1 5 8 24] range(A) B=randn(4) range(B) 东北电力大学 理学院 徐 屹 Northeast Dianli College Matlab语言程序设计 第六讲 线性方程组和数据分析 线性方程组求解 1 数据分析 2 一、Matlab的线性方程组求解 1、线性方程组解的结构 线性方程组的解可能出现三种情形：无解、有唯一解和有无穷多组解。 这主要取决于系数矩阵的秩与增广矩阵的秩是否相等，秩与变量个数是否相等，具体地： 若 ，则线性方程组无解； 若 （n为变量个数），则线性方程组有唯一解： 若 ，则线性方程组有无穷多组解。 一、Matlab的线性方程组求解 2、线性方程组求解的matlab命令 Matlab中主要用inv,null分别求矩阵的逆和齐次方程组的基础解系。 例1 解方程组 clear; A=[1 2;4 -3];b=[23;2]; x=A\b %左除法 （x,y）=6.6364,8.1818 如果用inv命令也可 clear; A=[1 2;4 -3];b=[23;2]; x=inv(A)*b （x,y）=6.6364,8.1818 exno61 一、Matlab的线性方程组求解 例2 求不定方程组 的一个特解。 clear; A=[1 2 1;3 -2 1];b=[2;2]; x=A\b （x,y,z）=(1,0.5,0) 一、Matlab的线性方程组求解 例3 求奇异方程组 的一个特解。 clear; A=[1 2; -2 -4];b=[1;-2]; x=A\b x=inf %说明不能直接求解 如果用 clear; A=[1 2; -2 -4;0 0];b=[1;-2;0]; x=A\b 仍可求出一个特解 (x,y)=(0,0.5000) 一、Matlab的线性方程组求解 例4 求方程组 的通解。 我们首先看方程组解的结构，计算系数矩阵和增广矩阵的秩。 clear; a=[1 -1 1 -1; -1 1 1 -1;2 -2 -1 1];b=[1;1;-1]; rank(a) %系数矩阵的秩 rank([a,b]) %增广矩阵的秩 exno62 一、Matlab的线性方程组求解 计算表明，系数矩阵和