Matlab原理与工程应用第二版第二章(矩阵).pptVIP

* 考虑下面三种运算 行向量和列向量乘积 列向量和行向量乘积 行向量和矩阵乘积 * 点乘，d=a*b 或d=dot(a,b) H为b、a中所有元素乘积的组合 * 方案1 * 方案2 * 行列式MATLAB表达式det(A) * 行列式特征值 lambda=eig(A,B) 求固有频率w 例： * 矩阵的逆 inv(A)或A^-1 方程组Ax=b求解(左除) x=A\b 例： * Key * * * * 作业： 1、 * clear all;clc; format long e z=magic(5); z(:,2)=z(:,2)/sqrt(3) z(5,:)=z(3,:)+z(5,:) z(:,1)=z(:,1).*z(:,4) q=z-diag(diag(z))+2*eye(5) m=q*q d=diag(q*q) e=q.^2 min=min(min(e)) max=max(max(e)) * 作业： 2、 * w=magic(2) a=repmat(w,2,2) b=repmat(w,3,1) d=repmat(w,3,1) e=repmat(w,3,1) m=[d,e] * 作业： 3、 4、 * clear all n=6; x=magic(n) m=x; for i=2:n m(i-1,:)=x(i,:); end m(n,:)=x(1,:); m format short x=[17,31,5;6 ,5 ,4;19 ,28 9;12 ,11 ,10] h=x*(x*x)^-1*x * 作业： 5、 format short a=[16,32,33,13;5,11,10,8;9,7,6,12;34,14,15,1]; b=[91,16,5,43]; x=a\b c=det(a) * * 第二章 矩阵与MATLAB 复习 矩阵和向量 由m行n列构成的数组A称为（m×n）矩阵。总共由m×n个元素组成，矩阵的元素记为aij,，其中i表示行，j表示列。 方阵（ m=n ）、对角阵、单位阵、行矩阵和列矩阵（向量）、矩阵的转置 * 2.1 向量的创建 向量可以表示为： f=[a x b]或f =[a,x,b,…] 其中a,x,b…可以是变量（预定义且赋值）、数值、表达式或字符串 采用冒号创建向量： x=s:d:f或x=（s:d:f）或x=[s:d:f] s-起始值或初始化值；d-增量或减量值；f-结束值或终值 因此，可产生以下行向量： x=[s,s+d,s+2d,…,s+nd]其中s+nd ≤f * * 向量x中元素的个数（即向量长度）可由length(x)给出。 d省略时，默认d=1。 例：创建取值范围从0.2到1.0，步长为0.1的行向量。 x=0.2:0.1:1 n=length（x） 求转置 x=（0.2:0.1:1）’ n=length（x） x=（0.2:0.12:1） n=length（x） x=0.5:7 n=length（x） * 创建等间隔值n x=linspace(s,f,n) x=[s,s+d,s+2d,…,f=s+(n-1)d] 增量或减量值 d=（f-s）/(n-1) 系统默认n=100。 对数等间隔 x=logspace(s,f,n) x=[10s,10s+d, 10s+2d ,…,10f]系统默认n=50。 * * 向量元素的访问 b=[b1 b2 b3 … b1] b(3)等同于b(1,3) 向量元素的运算 x=[-2,1:2:9,10] 可通过保留字end访问向量中的最后一个元素 x=[-2,1:2:9,10]； x7=x(end) 向量与标量进行加减运算时，标量与向量中每一元素相加减 x=[-2,1:2:9,10]；z=x-1 (乘、除、幂运算例外) * x=[-2,1:2:9,10]； z(2)=z(2)/2; z(3:4)=z(3:4)*3-1; z 注释：向量元素z(3)和z(4)分别乘3减1，终值代替z(3)和z(4)的原始值。这种语法非常适于编写简洁的程序代码。 访问向量元素的几种方法： y=[-1,6,15,-7,31,2,-4,-5] 创建一个y的第3个到第5个元素组成的向量x * x=y(3:5) x=[y(1), y(2), y(7) y(8)] index=[1,2,7,8](常用) x=y(index) x=y[1,2,7,8](常用) 演示:(猜猜sort、find函数的用法) y=[-1,6,15,-7,31,2,-4,-5]; z=[10,20,30,40,50,60,70,80]; [yn,in]=sort(y) znew=z(in) * y=[-1,6,15,-7,31,2,-4,-5]; indx=find(y=0) s=y(indx) 思考