[经济学]实验6.pptVIP

实验一 矩 阵 实验目的 1、学会用MATLAB软件生成矩阵. 3、熟悉MATLAB软件中关于矩阵 运算的各种命令. 2、熟悉MATLAB软件中关于矩阵 操作的各种命令. 一、矩阵的生成 （1） 输入矩阵时要以“[ ]”为其标识，即矩阵的元素应在“[ ]”内部； （2） 同行元素之间可由空格或“，”分隔，行与行间用“；”或回车符分隔； （3） 矩阵元素可为运算表达式； （4） 如不想获得中间结果，可以“ ；”结束。 1、 一般矩阵的生成: 例1： 输入矩阵 解： B=[1+9i 2+8i 3+7i;4+6j 5+5j 6+4i;7+3i 8+2i 0+j] 输入命令 在Matlab环境中定义了两个记号 i 和 j ,可以用来直接输入复数矩阵 2、特殊矩阵的生成： （1）生成单位矩阵 （２）生成全１矩阵 （3）生成全 0 矩阵 例2： 输入命令：eye(4) 输入命令：eye(3,4) 解 输入命令：zeros(4) 输入命令：zeros(3,4) 输入命令：ones(4) 输入命令：ones(3,4) 例3： 已知矩阵 试生成与矩阵A大小相同的单位阵，全0阵及全1阵； 解： 输入命令 A=[1 2 3;4 5 6;7 8 9];B1=eye(size(A)) B2=zeros(size(A)) B3=ones(size(A)) （4）生成随机矩阵 （5）生成对角矩阵 生成第 k 个对角线由向量 v 组成的对角阵， k 可以是正数，零或负数。当 k 是零时指主对角线(可简记为A= diag(v) )， k 是负数时 v 在对角线下， k 是正数时 v 在对角线上。 生成的对角矩阵的Matlab调用格式为: A= diag(v,k) 例4： （1）输入命令 A=diag([1 2 3 4]) A = 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 （2）输入命令 A=diag([1 2 3 4],1) A = 0 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 （3）输入命令 A=diag([2 3 4],-1) A = 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 0 （6）幻方矩阵（魔阵） 即每行每列和对角线上元素之和相等 命令： magic(n) 例5：magic(3) [8 1 6 3 5 7 4 9 2] （7）Vandermonde 矩阵 生成的Vandermonde矩阵的Matlat调用格式为: V=vander(C) 例6： 若向量 C=[1,2,3,4,5]，试写出该向量对应的 Vandermonde 矩阵 解： 输入命令 C=[1 2 3 4 5]; V=vander(C) V = 1 1 1 1 1 16 8 4 2 1 81 27 9 3 1 256 64 16 4 1 625 125 25 5 1 （8）希尔伯特(Hilbert) 矩阵 希尔伯特矩阵其元素为1/(i+j-1)，i,j分别为其行标和列标。 [1，1/2，1/3，……，1/n; 1/2，1/3，1/4，……，1/(n+1);1/3，1/4，1/5，……，1/(n+2);　　……1/n，1/(n+1)，1/(n+2)，……，1/(2n-1)] 例7：hilb(4) 希尔伯特矩阵是一种数学变换矩阵，正定，且高度病态（即，任何一个元素发生一点变动，整个矩阵的值和逆矩阵都会发生巨大变化），病态程度和阶数相关。 命令：hilb(n) 生成n阶希尔伯特矩阵 3、 符号矩阵的生成: 如果已经建立起数值矩阵 A，则可以由 B=sym(