第2章 矩阵和数组运算.pptVIP

第2章 矩阵和数组运算 2.1 矩阵函数和运算 矩阵创建： 直接输入法 矩阵编辑器法 函数创建法 矩阵保存和提取：save/load 矩阵运算:+、-、*、\和/ 矩阵处理函数 直接输入法 A=[2 3 4 5; 3 4 5 6; 4 5 6 7] A=[2+3i 3+4i 4 1+2i 2+3i 5 7 2+3i 9] x=2,y=pi/3,A=[2,5,cos(y); x,2*x,3*x;x,x/2,x/4] 注意：变量名大小写有区别！！ 矩阵编辑器法 函数创建法 zeros(m,n):m行n列的零矩阵； ones (m,n):m行n列的1矩阵； eye(m,n):m行n列矩阵，对角线元素为1，其余0； rand(m,n): m行n列随机矩阵; diag:生成对角阵或由对角线元素组成的向量； …….. load和save提取和保存矩阵法 load 数据文件名 在数据文件中，数据以ASCII代码方式存放，数据应排成矩阵形式，数据与数据之间用空格分开，数据每行仅包含矩阵一行，且各行元素个数必须相等。 矩阵变量的名称与数据文件名相同。 Save 文件名 变量名 [-ascii] 不加-ascii 则将变量以二进制方式存放到文件中（mat文件），否则以文本方式将变量对应的数据存放到文件中。 Load和save的一个例子 load d:\data1.txt data1 data1 = 1 2 3 4 5 6 7 8 9 矩阵运算 加法：A+B 减法：A-B 乘法：A*B 除法：A\B=(B’/A’)’,如果A和B均是方阵，且A非奇异，则A\B=inv(A)*B,B/A=B*inv(A).inv表示求逆矩阵。 乘方：A^p 转置：A’ 矩阵运算例子 矩阵处理函数 det(A):求矩阵A的行列式的值； dot(A,B):求矩阵A和B的点积； eig(A):求矩阵A的特征值和特征向量； norm(A,1):求矩阵A的1-范数； rank(A):求矩阵A的秩； trace(A):求矩阵A的迹； expm(A):求矩阵A的指数eA logm(A):求矩阵A的对数； sqrtm(A):求矩阵A的平方根 lu(A):对矩阵作LU分解； inv(A):求矩阵A的逆 矩阵运算：解方程组 解下列方程组: X+Y=5; 2X+Y=7; 应用Cramer法则求解: A=[1 1;2 1]; C=[5; 7]; 解法 [X Y] =（inv(A)*C）’ 结果：[X Y]=[2 3] 2.2 数组函数与运算 数组构造 1:0.5:3 ans = 1.0000 1.5000 2.0000 2.5000 3.0000 linspace(1,10,5) ans = 1.0000 3.2500 5.5000 7.7500 10.0000 logspace(2,5,4) ans = 100 1000 10000 100000 数组运算 数组函数 abs(x):数x的绝对值（向量x的长度、复数x的模）； angle(x):复数x的相角； conj(x):复数x的共轭复数； gcd(x,y):整数x和y的最大公约数； lcm(x,y):整数x和y的最小公倍数； round(x):四舍五入致最近整数； sin(x):求正弦函数； …. 数组元素标识（引用） 字符数组 字符串数组－基本元素是字符串。 字符串数组可用于计算过程和结果的可视化。 一维字符数组 country=‘China’ 二维字符数组 names=[‘Li Yi’;’Hu Xu’] 或names=char(‘Li Yi’;’Hu Ku’) strcat:将两个或多个字符数组组合在一起； str2num:字符串转换成数字； mat2str:数组转换成字符串; … 数据格式化写入字符串sprintf s=sprintf(格式说明，数据，…) [s,errmsg]=sprintf(格式说明，数据，…) 其中“格式说明”与C语言的格式输出语句printf类似。 格式控制读字符串sscanf A=sscanf(s,format) A=sscanf(s,format,si