大学选修实验矩阵与处理.ppt

第4章 数值向量和数组 实验1 验证马氏链的遍历性 及求平稳分布 一、基本实验 命令： 1、矩阵输入及运算: 2、矩阵的乘方: A^n 3 、 求解平稳分布归结为求 线性方程组表示为Ax=b，A不一定是方阵。 命令格式：C=[A b]; rref(C) 题目、求解线性方程组x+2y+3z+t=3 x+4y+5z+2t=2 2x+9y+8z+3t=7 3x+7y+7z+2t=12 命令：A=[1 2 3 1;1 4 5 2;2 9 8 3;3 7 7 2];b=[3;2;7;12]; C=[A b]; rref(C) A=[1 2 3 1;1 4 5 2;2 9 8 3;3 7 7 2];b=[3;2;7;12]; C=[A b]; rref(C) ans = 1.0000 0 0 -0.5000 5.1667 0 1.0000 0 0 0.6667 0 0 1.0000 0.5000 -1.1667 0 0 0 0 0 即 x-0.5t=5.1667 y=0.6667 z+0.5t=-1.1667 命令： P=[1/3 2/3 0;1/3 0 2/3 ;0 1/3 2/3]; P^n B =[1/3 2/3 0 1;1/3 0 2/3 1;0 1/3 2/3 1]; D=[1 0 0 0;0 1 0 0;0 0 1 0]; A=B’-D’; or A=[-1+1/3 1/3 0;2/3 -1 1/3 ; 0 2/3 -1+2/3;1 1 1];b=[0;0;0;1]; C=[A b]; rref(C) 命令： p=[336/480 48/480 96/480;32/320 224/320 64/320;64/800 32/800 704/800]; p0=[0.3 0.2 0.5]; p1=p0*p p4= p0*p^4 B =[336/480 48/480 96/480 1;32/320 224/320 64/320 1;64/800 32/800 704/800 1]; D=[1 0 0 0;0 1 0 0;0 0 1 0]; A=B-D; b=[0;0;0;1]; C=[A b]; rref(C) 下面数据为1915到1973年的最大年径流量（按行从左到右）： 15600 8960 10400 10600 10800 9880 9850 10900 8810 9960 12200 7510 8640 6380 6810 8820 14400 7440 7240 6430 11000 7340 9260 5290 9130 7480 6980 9650 7260 8750 9900 7310 9040 7310 8850 7840 10700 6190 9610 7580 9990 6150 8250 6030 8980 6180 9630 9490 2340 11100 5090 10900 6490 12600 6640 7430 6760 10000 9300 实验2 时间序列的线形模型的识别 一、基本实验 1、平稳时间序列的标准化 z=[15600 8960 10400 10600 10800 9880 9850 10900 8810 9960 12200 7510 +... 8640 6380 6810 8820 14400 7440 7240 6430 11000 7340 9260 5290 +... 9130 7480 6980 9650 7260 8750 9900 7310 9040 7310 8850 7840 +... 10700 6190 9610 7580 9990 6150 8250 6030 8980 6180 9630 9490 +... 2340 11100 5090 10900 6490 12600 6640 7430 6760 10000 9300]; 命令：（1）