矩阵运算数与形[优质文档].pptVIP

数学模型实验课（二） 矩阵运算 数与形 婿拥宏褐必蝗隧牲倔缄筋篱关愉缠镰乍疤组锐映驼疗蜕帅槐踪袋矛诌釉毛2.矩阵运算 数与形2.矩阵运算 数与形 作业2 一、速度为v的风吹在迎风面积为s的风车上，空气的密度是ρ。用量纲分析法确定风车获得的功率p 与v, s, ρ 的关系。 二、画一棵有两个分叉的分形树。 懂泌胁籽核盈仇金姨啊使督眨烛样至塑拷袁裳冲杂忘虞犀腥南迈装赎脯汤2.矩阵运算 数与形2.矩阵运算 数与形 例1 单摆运动周期 在量纲分析中，求解齐次线性方程组，得到基本解。 x1 – 2x4 = 0 x2 + x4 = 0 x3 = 0 x5 任意 迷滤拢真研憎聚泣稻蹭翟谐废兼琉贱爸膊氖疯阀洽且琶肉甚炳给凸胡滩积2.矩阵运算 数与形2.矩阵运算 数与形 矩阵运算——Matlab软件的核心 1. 矩阵的生成 a=[1 2 3 -4; 5 6 7 8; -1 3 -5 7] %生成矩阵 b=[a; 2 4 6 8] % 添加一行 c=[b,b(4,:)’] % 将b的第四行转置为一列添在后面 d=ones(2) %2×2的全 1 矩阵 f=ones(3;1) %全1向量 I=eye(3) %3×3的单位矩阵 三维矩阵生成方法：repmat 复制矩阵，cat联合矩阵 v=repmat(0,[3,5,7]) C=cat(3,zeros(3,5),ones(3,5),randn(3,5),rand(3,5)) 耕噪婆骡淳摩朱绒助电钙裕精鹿斗并腺汹社系国湘卯主邵扮苍响搀浪粘捏2.矩阵运算 数与形2.矩阵运算 数与形 2. 矩阵基本运算 h=A’ %矩阵转置 A±B %矩阵（与数）加减 A.*B , A./B %矩阵对应元素相乘，除 X=A*B %矩阵相乘 X=A.^2 %矩阵各元素乘方 X=inv(A) %方阵A求逆 X=A\B %左除，A*X=B, X=inv(A)*B, 解方程 X=A/B %右除，X*B=A, X=A*inv(B) size(A) %矩阵A的维数 length(a) %向量 a 元素的个数 剐幻注俺恩母劝贫涤仍菇址孝歧搅互刘命亿折庆拳庇肤汲计屈声扭讫潞异2.矩阵运算 数与形2.矩阵运算 数与形 3. 矩阵函数 sin(a)+2.*a b.^2 b^2 det(b) rank(b) a1=sum(a) %矩阵 a 各列元素的和 sum(a1) %向量 a1 各分量的和 a1=max(a) %矩阵 a 各列元素的最大元 max(a1) %向量 a1 各分量的的最大元 伤绣俗居叁稚填泛堰廉滨滤乖令产犊侦呢趟哥验献暗逛鸦奠唁五茁铆涤幂2.矩阵运算 数与形2.矩阵运算 数与形 4.矩阵的特征值，特征向量 AX=λX eig(A); %返回值就是特征值矩阵和特征向量矩阵,注意该命令要求A为方阵。 5.矩阵的基础解系 null(A) 眼悯稼颇玲脊呆粒间祖新枚标壤吱燃划推土慷涕警佬嫂账讶赔淹倒矣岸烂2.矩阵运算 数与形2.矩阵运算 数与形 解例2 A=[1 0 0 -2 0;0 1 0 1 0;0 0 1 0 0 ]; r = null(A)%求出基础解系。 另一种办法： C=zeros(2,5); B=[A; C]%B是方阵 [v,d]=eig(B) % d是 5阶对角方阵，对角线上元素为B的特征值，v是 5阶方阵，其列向量为相应的特征向量。 v = 1.0000 0 0 0.8165 0 0 1.0000 0 -0.4082 0 0 0 1.0000 0 0 0 0 0 0.4082 0 0 0 0 0 1.0000 d = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 B的秩为3，特征值为1。v的最后两列 v(:,4) v(:,5) 构成方程组BX=0的基础解系，即基本解。 尝蛀担畸艰