matlab在科学计算中的应用4[姜志鹏].pptVIP

第四章 线性代数问题求解;4.1 矩阵4.1.1特殊矩阵的输入;随机元素矩阵 若矩阵随机元素满足[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 矩阵;Vandermonde(范德蒙)矩阵 ;伴随矩阵 其中：P(s)为首项系数为1的多项式。 ; 例：考虑一个多项式2*x^4+4*x^2+5*x+6，试写出该多项式的伴随矩阵。 P=[2 0 4 5 6];A=compan(P) A = 0 -2.0000 -2.5000 -3.0000 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0;符号矩阵的输入 数值矩阵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 矩阵基本概念与性质;例： tic, A=sym(hilb(20)); det(A), toc ans = 1/2377454716768534509091644243427616440175419837753486493033185331234419759310644585187585766816573773440565759867265558971765638419710793303386582324149811241023554489166154717809635257797836800000000000000000000000000000000000 elapsed_time = 2.3140 高阶的Hilbert矩阵是接近奇异的矩阵。;矩阵的迹 格式： t=trace(A) 矩阵的秩 格式：r=rank(A) ％用默认的精度求数值秩 r=ra