原创力文档- 1、本文档共22页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
张旭-上机教学(改进版)
《线性代数》 —上机教学 MATLAB中基本代数运算符 运算 符号 举例 加法，a+b + 5+3 减法，a-b - 5-3 乘法，a×b * 5*3 除法，a÷b / or \ 48/4=4\48=12 乘幂，ab ^ 5^2=25 变量及数组输入 2：矩阵方式输入（请记下红色部分！） a=[1,2,3;4,5,6;7,8,0] %矩阵输入 （a为3阶方阵） b=[366;804;351] %列矩阵输入 c=[366;804;351] ‘ %行矩阵（转置）输入 det(a) %方阵行列式 inv(a) %方阵的逆 x=a\b %ax=b方程组的解 y=inv(a)*b %与x相同 r=rank(A), %矩阵的秩 disp([a,b,x]) %显示矩阵 变量及数组输入 3：矩阵的简单运算 c=inv(a) %方阵的逆阵 y=c*b %矩阵乘积 d=[a b];disp([d]) %矩阵拼接 d=a; disp([d]) %矩阵转置 g=2*a+3 %常数乘矩阵，各元素加3 p=eye(3) %3阶单位矩阵 y=a.*p %两矩阵对应元素乘积 zeros(3) %3阶零矩阵 例1 求下列矩阵列向量组的一个最大无关组. 解：在Matlab中输入： a=[1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4]; b=rref(a) 求得: 所以 是一个极大无关组,且 Warning: Could not get change notification handle for local C:\Program Files\MATLAB\R2007b\work. Performance degradation may occur due to on-disk directory change checking. a=[1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4]; b=rref(a) %不打；则计算机将显示b的值 b = 1.0000 0 0.3333 0 1.7778 0 1.0000 0.6667 0 -0.1111 0 0 0 1.0000 -0.3333 0 0 0 0 0 直接解法 利用左除运算符的直接解法 对于线性方程组Ax=b，可以利用左除运算符“\”求解： x=A\b A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; b=[13,-9,6,0]; x=A\b %不打；则计算机将显示x的值 x = -66.5556 25.6667 -18.7778 26.5556 例子三、求秩 A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; r=rank(A); r= % = 计算机不显示r的值 ??? r= | Error: Expression or statement is incomplete or incorrect. rank(A) %不打；则计算机将显示rank(A)的值 ans = 4 r r = 4 三.Help 四、上机作业 已
文档评论（0）