原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
第四章 线性代数问题求解 矩阵 线性方程组的直接解法 线性方程组的迭代法 线性方程组的符号解法 稀疏矩阵技术 特征值与特征向量 4.1 矩阵 4.1.1特殊矩阵的输入 数值矩阵的输入 零矩阵、幺矩阵及单位矩阵 生成n?n方阵： A=zeros(n), B=ones(n), C=eye(n) 生成m?n矩阵： A=zeros(m,n), B=ones(m,n), C=eye(m,n) 生成和矩阵B同样位数的矩阵： A=zeros(size(B)) 随机元素矩阵 若矩阵随机元素满足[0,1]区间上的均匀分布 生成n?m阶标准均匀分布为随机数矩阵： A=rand(n,m) 生成n?n阶标准均匀分布为随机数方阵： A=rand(n) 对角元素矩阵 已知向量生成对角矩阵： A=diag(V) 已知矩阵提取对角元素列向量： V＝diag(A) 生成主对角线上第k条对角线为V的矩阵： A=diag(V,k) 例：diag( )函数的不同调用格式 C=[1 2 3]; V=diag(C) % 生成对角矩阵 V = 1 0 0 0 2 0 0 0 3 V1=diag(V) % 将列向量通过转置变换成行向量 V1 = 1 2 3 C=[1 2 3]; V=diag(C,2) % 主对角线上第 k条对角线为C的矩阵 V = 0 0 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 生成三对角矩阵: V=diag([1 2 3 4])+diag([2 3 4],1)+diag([5 4 3],-1) V = 1 2 0 0 5 2 3 0 0 4 3 4 0 0 3 4 Hilbert矩阵及逆Hilbert矩阵 生成n阶的Hilbert矩阵： A=hilb(n) 求取逆Hilbert矩阵： B=invhilb(n) Hankel(汉克 ) 矩阵 其中：第一列的各个元素定义为C向量，最后一行各个元素定义为R。H为对称阵。 H1=hankel(C) 由 Hankel 矩阵反对角线上元素相等得出一下三角阵均为零的Hankel 矩阵 例：用Matlab语句输入Hankel矩阵 Vandermonde(范德蒙)矩阵 伴随矩阵 例：写出如下多项式的伴随矩阵 符号矩阵的输入 数值矩阵A转换成符号矩阵： B=sym(A) 例： A=hilb(3) A = 1.0000 0.5000 0.3333 0.5000 0.3333 0.2500 0.3333 0.2500 0.2000 B=sym(A) B = [ 1, 1/2, 1/3] [ 1/2, 1/3, 1/4] [ 1/3, 1/4, 1/5] 4.1.2 矩阵基本概念与性质 行列式 格式 ：d=det(A) 例：求行列式 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; det(A) ans = 0 例： tic, A=sym(hilb(20)); det(A), toc ans = 1/23774547167685345090916442434276164401754198377534864930331853312344197593106445851875857668165737734405657598672655589717656384197107933033865823241498112410235544891661547178096352577978368000000000000000000000000000
文档评论（0）