刘卫国全套配套课件MATLAB程序设计与应用第3版第3章 MATLAB矩阵处理_源程序.doc

《刘卫国全套配套课件MATLAB程序设计与应用第3版第3章 MATLAB矩阵处理_源程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第3章MATLAB矩阵处理例3-1建立随机矩阵：>>x=20+(50-20)*rand(5)>>y=0.6+sqrt(0.1)*randn(5)例3-2将[101，125]范围内的25个整数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。>>M=100+magic(5)例3-3求4阶希尔伯特矩阵及其逆矩阵。>>formatrat%以有理形式输出>>H=hilb(4)>>H=invhilb(4)>>format%恢复默认输出格式例3-4求(x+y)5的展开式。>>pascal(6)例3-5先建立5×5矩阵A，然后将A的第一行元素乘以1，第二行乘以2

2、，…，第五行乘以5。>>A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;...11,18,25,2,19];>>D=diag(1:5);>>D*A%用D左乘A，对A的每行乘以一个指定常数例3-6求方阵A的逆矩阵，且验证A与A-1是互逆的。>>A=[1,-1,1;5,-4,3;2,1,1];>>B=inv(A);>>A*B>>B*A例3-7用求逆矩阵的方法解线性方程组。>>A=[1,2,3;1,4,9;1,8,27];>>b=[5,-2,6]';>>x=inv(A)*b2例3-8用求特征值的方法解方程。

3、>>p=[3,-7,0,5,2,-18];>>A=compan(p);%A的伴随矩阵>>x1=eig(A)%求A的特征值>>x2=roots(p)%直接求多项式p的零点例3-9将X转化为稀疏存储方式。>>X=[2,0,0,0,0;0,0,0,0,0;0,0,0,5,0;0,1,0,0,-1;0,0,0,0,-5];>>A=sparse(X)例3-10根据表示稀疏矩阵的矩阵A，产生一个稀疏存储矩阵B。>>A=[2,2,1;3,1,-1;4,3,3;5,3,8;6,6,12];>>B=spconvert(A)例3-11求下列三对角线性方程组的解。>>B=[1,2,

4、0;1,4,3;2,6,1;1,6,4;0,1,2];%产生非0对角元素矩阵>>d=[-1;0;1];%产生非0对角元素位置向量>>A=spdiags(B,d,5,5)%产生稀疏存储的系数矩阵2