基于MATLAB的线性代数应用案例
用MATLAB解决线性代数问题实验报告
实验三使用MATLAB解决线性代数问题学院：数计学院班级：1003班姓名：黄晓丹学号：1051020144实验目的：学习MATLAB有关线性代数运算的指令，主要学习运用MATLAB解决矩阵除法，线性方程组的通解，矩阵相似对角化问题，以及解决投入产出分析等应用问题。
实验内容：矩阵转置：A=[1 2;3 4];B=[4 3;2 1];>> A',B'ans =1 32 4ans =4 33 1矩阵加减：A-Bans=-3 -11 3矩阵乘法：A*B,A.*B(数组乘法)||比较矩阵乘法与数组乘法的区别ans=8 520 13ans=4 66 4矩阵除法：A\B,B./Aans=-6 -55 4ans=4 1.50.6667 0.25特殊矩阵生成：zeros（m，n）||生成m行n列的矩阵ones（m，n）||生成m行n列的元素全为一的矩阵eye（n）||生成n阶单位矩阵rand（m，n）||生成m行n列[0 ，1]上均匀分布随机数矩阵zeros(2,3)ans =0 0 00 0 0>> ones(3,3)ans =1 1 11 1 11 1 1>> eye(3)ans =1 0 00 1 00 0 1>> rand(2,4)ans =Columns 1 through 30.9501 0.6068 0.89130.2311 0.4860 0.7621Column 40.45650.0185矩阵处理：trace（A）||返回矩阵的迹diag（A）||返回矩阵对角线元素构成的向量tril（A）||提取矩阵的下三角部分triu（A）||提取矩阵的上三角部分flipud（A）||矩阵上下翻转fliplr（A）||矩阵左右翻转reshape(A,m,n)||将矩阵的元素重排成m行n列矩阵A=[1 2 3;4 5 6;7 8 9];>> t=trace(A),d=diag(A),u=triu(A)t =15d =159u =1 2 30 5 60 0 9 flipud(A),fliplr(A)ans =7 8 94 5 61 2 3 ans =3 2 16 5 49 8 7矩阵特征值与标准型：[V,D]=eig（A）||返回矩阵特征值与特征向量[V J]=Jordan(A)||返回矩阵的相似变换矩阵和若尔当标准型A=[1 2;3 4];>> [V,D]=eig(A)V =-0.8246 -0.41600.5658 -0.9094D =-0.3723 00 5.3723>> [V,J]=jordan(A)V =0.2389 0.76110.5222 -0.5222J =5.3723 00 -0.3723线性方程组求解A=[1 2 1;3 -2 1];B=[1;4];x=A\B x =1.2500 ||求一特解-0.1250>> A=[1 2;3 -2;1 -1];B=[1;4;2];x=A\Bx = ||求得一最小二乘近似解1.2838-0.1757：方阵的相似对角化及应用：A=[1 1/4 0;0 1/2 0;0 1/4 1];[P,T]=eig(A) P =1.0000 0 -0.40820 0 0.81650 1.0000 -0.4082T =1.0000 0 00 1.0000 00 0 0.5000求得三个特征值1,1,0.5，对应特征向量（1,0,0），（0,0,1），（-0.4028,0.8165，-0.4082），由于三个特征向量线性无关，从而A 可相似对角化，即p-1AP=T.那么A∧n=p[1 0 0;0 1 0;0 0 0]p-1，计算的P*diag（[1,1,0]）*inv（P）ans =1.0000 0.50000 00 0 00 0.5000 1.0000所以得到近似解。
Matlab在线性代数中的应用
利用Matlab的控制设计方法，如PID控制、状态反馈控制等，可以 设计出有效的控制系统。
THANKS
感谢观看
利用Matlab的图像处理函数，可以从图像中提取 特征，如边缘、角点等，用于目标检测和识别。
在控制系统中的应用
系统建模
使用Matlab的控制系统工具箱，可以对系统进行建模，如线性时 不变系统、非线性系统等。
系统分析和仿真
通过Matlab的控制系统函数，可以对系统进行稳定性分析、控制 性能分析和仿真测试。
向量运算
向量的基本运算
包括向量的加法、减法、数乘、向量的模等。
向量的内积和外积
内积和外积是描述向量之间关系的运算，用于计算向量的长度、角 度等。
向量运算的实际应用
向量运算在物理、工程等领域有广泛应用，如描述物体运动轨迹、计 算力的合成等。
特征值与特征向量
01
特征值和特征向量 的定义
特征值和特征向量是描述矩阵特 性的重要概念，用于描述矩阵变 换的性质。
04
Matlab在线性代数中的优势与 局限性
优势
高效计算能力
Matlab提供了强大的矩阵运算 和数值计算功能，使得线性代
数问题的求解更加高效。
可视化工具
Matlab内置了丰富的可视化工 具，可以直观地展示线性代数 中的向量、矩阵和线性变换等 概念。
易于学习和使用
Matlab的语法相对简单，使得 线性代数运算变得容易理解和 实现。
解的精度和稳定性
Matlab在线性方程组求解过程中考虑了精 度和稳定性问题，能够提供可靠的解。
向量运算和特征值问题
向量运算
Matlab支持向量的基本运算 ，如加法、减法、数乘、点 积等。
#实验1 MATLAB续：MATLAB在线性代数中的应用
6 x3 5 x3
6
x4
0 0
x3 5x4 6x5 0
x4 5x5 1
A=[5 6 0 0 0 15600 01560 00156 0 0 0 1 5];
B=[1 0 0 0 1]'; R_A=rank(A) %求秩 X=A\B %求解
例2
求解方程组
3xx112xx2253xx333xx44
>>B=null(A,'r') %求解空间的有理基
得到：B =
2
5/3
-2 -4/3
1
0
0
1
MATLAB初步
于是，我们得到原线性方程组的解：
syms k1 k2 % 定义两个符号 X=k1*B(:,1)+k2*B(:,2) %写出方程组的通解 pretty(X) %让通解表达式更加精美
求解的完整代码如下：
%V已经被归一化为单位向量了
例 4： 求矩阵
MATLAB初步
2 1 1
A
0 4
2 1
0 3
的特征值和特征向量.
A=[-2 1 1;0 2 0;-4 1 3]; [V,D]=eig(A)
V=
-0.7071 -0.2425
0
0
-0.7071 -0.9701
0.3015 0.9045 0.3015
1 2
有否解? 2x1 x2 2x3 2x4 3
MATLAB初步
A=[1 -2 3 -1;3 -1 5 -3;2 1 2 -2];%输入系数矩阵A的值
% first,input the coefficient matrix A
matlab线性代数例题[大全5篇]
![matlab线性代数例题[大全5篇]](https://img.taocdn.com/s3/m/76804a34a7c30c22590102020740be1e650ecc80.png)
matlab线性代数例题[大全5篇]第一篇：matlab线性代数例题《数学实验》在线习题3 Matlab程序设计部分一.分析向量组a1=[1T2a23=]-,-T[a31T=2，0],a4=[1-2-1]T，a5=[246]T的线性相关性，找出它们的最大无关组，并将其余向理表示成最大无关组的线性组合。
解，a1=[1 2 3]';a2=[-1-2 0]';a3=[0 0 1]';a4=[1-2-1]';a5=[2 4 6]';A=[a1,a2,a3,a4,a5];[R,S]=rref(A)r=length(S)R =1.0000 0 0.3333 02.0000 0 1.0000 0.3333 0 0 0 0 0 1.0000 0S =4r =线性相关 a1,a2,a3,a4,a5 最大无关组是a1,a2,a4 其余向量的线性组合是a3=1/3a1+1/3a2 a5=2a1二.计算行列式x13D4=x23x33x43x12y1x22y2x32y3x42y4x1y12x2y22x3y32x 4y42y13y23y3323的值。
其中[1解，syms x1 x2 x3 x4 y1 y2 y3 y4 xxxy43x4]=[2357],[y1y2y3y4]=[4567]。
D=[x1^3 x1^2*y1 x1*y1^2 y1^3;x2^3 x2^2*y2 x2*y2^2 y2^3;x3^3 x3^2*y3 x3*y3^2 y3^3;x4^3 x4^2*y4 x4*y4^2 y4^3];d=det(D)x1=2;x2=3;x3=5;x4=7;y1=4;y2=5;y3=6;y4=7;eval (d)d = ans =153664 三.已知向量a={1,-1,0}，b={-1,0,-1}，求向量a与b的夹角的度数。
解，a=[1-1 0];b=[-1 0-1];x=a.*b;x1=sum(x,2);x2=norm(a);x3=norm(b);y=x1/(x2*x3)y1 =acos(y)y =-0.5000y1 =2.0944四.已知线性方程组clear⎧2x1-x2+3x3+2x4=0⎪9x-x+14x+2x=1⎪1234⎨⎪3x1+2x2+5x3-4x4=1⎪⎩4x1+5x2+7x3-10x4=2，求系数矩阵的秩和方程组的通解。
用MATLAB做线性代数实验
2
0
, 2
5 3
, 3
1
3
, 4
1
4
, 5
1
2
。
3
6
0Hale Waihona Puke 73【程序如下】：
% (1)
A=[1 2 1 3;4 -1 5 6;1 -3 -4 7;1 2 1 1]' r=rank(A) [R,IP]=rref(A) % (2) A=[1 2 0 2 1;-2 -5 1 -1 1;0 -3 3 4 2;P3 6 0 -7 3] r=rank(A) [R,IP]=rref(A)
例如：
已知
A
1 3
2 4
,
B
1 1
2 0
，解矩阵方程
(1)
AX
B , (2) XA B 。
MATLAB 程序如下：
A=[1 2;3 4];
B=[1 2;-1 1];
X1=inv(A)*B % AX=B or
X1=A\B
X2=B*inv(A) % XA=B
X2=A/B
将 p(x) 分解为最简分式之和 q( x)
[p,q]=residue(a,b,r) 将简单分式之和合并为有理分式
例如，将有理分式
f
(x)
x2 x3 2x2 3x2
分解为最简分式之和的程序如下：
p=[1 2];
q=[1 2 3 2];
[a,b,r]=residue(p,q)
输出：a =
-0.2500 - 0.4725i
p=[1 -6 11 -6];
基于MATLAB的线性代数应用案例
Advances in Applied Mathematics 应用数学进展, 2019, 8(3), 424-429Published Online March 2019 in Hans. /journal/aamhttps:///10.12677/aam.2019.83048Application Examples of Linear AlgebraBased on MATLABWei Yang1, Shuping Gao2, Bing Han3, Huaichen Chen31School of Physics and Optoelectronic Engineering, Xidian University, Xi’an Shaanxi2School of Mathematics and Statistics, Xidian University, Xi’an Shaanxi3School of Electornoc Engineering, Xidian University, Xi’an ShaanxiReceived: Feb. 14th, 2019; accepted: Feb. 27th, 2019; published: Mar. 6th, 2019AbstractThis paper gives three examples of the application of linear algebra. It includes rigid body plane motion, information retrieval and bacteria conversion. MATLAB software is used for numerical calculation and drawing.KeywordsLinear Algebra, Application of Linear Algebra, MATLAB基于MATLAB的线性代数应用案例杨威1，高淑萍2，韩冰3，陈怀琛31西安电子科技大学物理与光电工程学院，陕西西安2西安电子科技大学数学与统计学院，陕西西安3西安电子科技大学电子工程学院，陕西西安收稿日期：2019年2月14日；录用日期：2019年2月27日；发布日期：2019年3月6日摘要本文给出了线性代数三个精彩应用案例，包括刚体平面运动、情报检索及细菌转换，并利用MATLAB软件进行数值计算和绘图分析。
利用Matlab进行线性代数问题求解的方法与案例
利用Matlab进行线性代数问题求解的方法与案例引言线性代数是数学的一个重要分支，广泛应用于工程、物理、计算机科学等领域。
而Matlab作为一种功能强大的数值计算软件，提供了各种实用的工具和函数，可以方便地解决线性代数问题。
本文将介绍一些常用的线性代数问题求解方法，并通过具体的案例来展示Matlab在实际应用中的效果。
一、线性方程组的求解线性方程组是线性代数中最基础的问题之一。
Matlab提供了多种求解线性方程组的函数，如“backslash”操作符（\）和“linsolve”函数等。
下面通过一个实例来说明Matlab的线性方程组求解功能。
案例：假设有以下线性方程组需要求解：2x + 3y - 4z = 53x - 2y + z = 8x + 5y - 3z = 7在Matlab中输入以下代码：A = [2 3 -4; 3 -2 1; 1 5 -3];b = [5; 8; 7];x = A\b;通过以上代码，我们可以得到线性方程组的解x=[1; -2; 3]。
这表明在满足以上方程组的条件下，x=1，y=-2，z=3。
可以看出，Matlab在求解线性方程组时，使用简单且高效。
二、矩阵的特征值和特征向量求解矩阵的特征值和特征向量也是线性代数中的重要概念。
利用特征值和特征向量可以得到矩阵的许多性质和信息。
在Matlab中，我们可以通过“eig”函数来求解矩阵的特征值和特征向量。
案例：假设有一个2x2矩阵A，需要求解其特征值和特征向量。
在Matlab中输入以下代码：A = [2 3; 1 4];[V, D] = eig(A);通过以上代码，我们可以得到矩阵A的特征向量矩阵V和特征值矩阵D。
具体结果如下：特征向量矩阵V = [0.8507 -0.5257; 0.5257 0.8507]特征值矩阵D = [1.5858 0; 0 4.4142]由结果可知，矩阵A的特征向量矩阵V和特征值矩阵D可以提供有关该矩阵的很多信息，如相关线性变换、对称性等。
使用Matlab求解线性代数问题（一）
使用Matlab求解线性代数问题一、用行列式求解线性方程组x1+ x2–2*x3 = -35*x1–2*x2+ 7*x3 = 222*x1–5*x2+ 5*x3 = 4解法1：使用矩阵求解方程组在Matlab 7.1命令行界面中输入（A表示系数矩阵，b表示向量）A = [1 1 -2; 5 -2 7; 2 -5 5]; % 系数矩阵b = [-3; 22; 4]; % 向量format rat; % 设置结果的显示形式为分数x = A\b; % 计算结果x = x' % 转置以行向量形式显示界面显示结果为x = 47/56 163/56 27/8即方程的解为x1 = 47/56; x2 = 163/56; x3 = 27/8解法2：使用Cramer法则求解方程组在Matlab 7.1命令行界面中输入D（系数行列式）D = [1 1 -2; 5 -2 7; 2 -5 5]; % 系数矩阵DD1 = [b D(:,2) D(:, 3)]; % D1D2 = [D(:,1) b D(:, 3)]; % D2D3 = [D(:, 1) D(:,2) b]; % D3format rat; % 设置结果的显示形式为分数x1 = det(D1)/det(D)x2 = det(D2)/det(D)x3 = det(D3)/det(D)输出结果为x1 = 47/56x2 = 163/56x3 = 27/8即方程的解为：x1 = 47/56; x2 = 163/56; x3 = 27/8二、求排列65872134的逆序数，并确定其奇偶性。
解：排列的逆序数定义为：如在1,2,3,…,n的一个全排列（s1 s2 s3…s n）中，有i<j时，s i >s j，这时s i、s j违反了自然顺序，就说它们构成了一个逆序。
排列（s1 s2 s3…s n）中逆序的总数称为该序列的逆序数。
根据定义计算逆序数的Matlab程序如下：a = [6 5 8 7 2 1 3 4]; % 将排列看作行向量Num = 0; % 逆序数初始值设为0Len = length(a); % 行向量中元素的个数for i = 1:1:(Len-1)n = length(find( a(i+1: Len) < a(i))); % 计算第i个元素% 的逆序数disp(strcat('第', num2str(i), '个元素', num2str(a(i)), '的逆序数为', num2str(n) ));Num = Num + n; %增加序列逆序数的个数end;disp(strcat('排列', num2str(a), '的逆序数为', num2str(Num)))将上述Matlab命令输入Matlab窗口命令行，执行结果为：第1个元素6的逆序数为5第2个元素5的逆序数为4第3个元素8的逆序数为5第4个元素7的逆序数为4第5个元素2的逆序数为1第6个元素1的逆序数为0第7个元素3的逆序数为0排列6 5 8 7 2 1 3 4的逆序数为19即排列65872134的逆序数为19，由于19为基数，所以排列65872134为逆排列。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Advances in Applied Mathematics 应用数学进展, 2019, 8(3), 424-429Published Online March 2019 in Hans. /journal/aamhttps:///10.12677/aam.2019.83048Application Examples of Linear AlgebraBased on MATLABWei Yang1, Shuping Gao2, Bing Han3, Huaichen Chen31School of Physics and Optoelectronic Engineering, Xidian University, Xi’an Shaanxi2School of Mathematics and Statistics, Xidian University, Xi’an Shaanxi3School of Electornoc Engineering, Xidian University, Xi’an ShaanxiReceived: Feb. 14th, 2019; accepted: Feb. 27th, 2019; published: Mar. 6th, 2019AbstractThis paper gives three examples of the application of linear algebra. It includes rigid body plane motion, information retrieval and bacteria conversion. MATLAB software is used for numerical calculation and drawing.KeywordsLinear Algebra, Application of Linear Algebra, MATLAB基于MATLAB的线性代数应用案例杨威1，高淑萍2，韩冰3，陈怀琛31西安电子科技大学物理与光电工程学院，陕西西安2西安电子科技大学数学与统计学院，陕西西安3西安电子科技大学电子工程学院，陕西西安收稿日期：2019年2月14日；录用日期：2019年2月27日；发布日期：2019年3月6日摘要本文给出了线性代数三个精彩应用案例，包括刚体平面运动、情报检索及细菌转换，并利用MATLAB软件进行数值计算和绘图分析。
杨威 等关键词线性代数，线性代数应用，MATLABCopyright © 2019 by author(s) and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License (CC BY)./licenses/by/4.0/1. 引言线性代数是一门实用性很强的基础课程，航空、建筑、计算机、经济、金融等等相关学科都需要大量的线性代数知识。
但是，在目前国内教材里，关于线性代数的应用内容涉及很少，甚至没有涉及。
随着矩阵阶数的增加，线性代数中各种运算的运算量将会剧增，所以利用计算机软件实现线性代数运算是当今线性代数教学改革的必要内容。
国内大多教材缺少了这方面内容，虽然个别新教材加入了MATLAB 软件内容，但其与线性代数知识分离，并没有很好地把线性代数理论与MATLAB 软件深度融合。
本文分别介绍了三个精彩的线性代数应用案例，并用MATLAB 软件进行计算和绘图，并对结果进行了分析讨论。
2. 刚体的平面运动用矩阵X 来描述平面坐标系中的一个闭合图形(即“刚体”)，X 的列向量表示图形n 个顶点的坐标。
为了实现刚体的平移及旋转运算，给矩阵X 添加元素值都为1的一行，使矩阵X 的形状为3 × n 。
若有矩阵：121001001c M c =，cos sin 0sin cos 0001t t R t t − = ，且：1Y MX =，2Y RX = 可以证明，矩阵1Y 是刚体X 沿x 轴正方向平移1c ，沿y 轴正方向平移2c 后的结果；矩阵2Y 是刚体X 以坐标原点为中心逆时针旋转t 弧度的结果[1]。
例1. 表1给出大写字母A 的11个顶点坐标值，其中最后一列与第一列相同，表示一个完整的闭合图形。
对图形A 先逆时针旋转3π4，然后向上移动30，再向右移动20，并绘制移动前后的图形。
Table 1. Vertex coordinates of letter A 表1. 字母A 的顶点坐标x 0 4 6 10 8 5 3.5 6.1 6.5 3.2 2 0 y141411664.54.5解：构造刚体矩阵0461085 3.5 6.1 6.5 3.22001414001166 4.5 4.500111*********X= ， 旋转矩阵3π3πcos sin 0443π3πsin cos 044001R −=及平移矩阵10200130001M = 。
Open Access杨威 等在MATLAB 的M 文件编辑器中编写程序l1.m % 刚体的平面运动 close allX=[0,4,6,10,8,5,3.5,6.1,6.5,3.2,2,0;0,14,14,0,0,11,6,6,4.5,4.5,0,0;ones(1,12)]; % 构造原刚体矩阵 M=[1,0,20;0,1,30;0,0,1]; %构造平移转矩阵 R=[cos(3*pi/4),-sin(3*pi/4),0;sin(3*pi/4),cos(3*pi/4),0;0,0,1]; %构造旋转矩阵Y3=M*R*X; %通过线性变换计算移动后的刚体矩阵 plot(X(1,:),X(2,:)); % 绘制原来刚体 hold on axis equalfill(Y3(1,:),Y3(2,:),'black'); % 绘制旋转及平移后刚体 grid on hold off在MATLAB 命令窗口中输入：l1 绘制图形如图1所示。
Figure 1. Letter A before and after moving 图1. 字母A 移动前后图形3. 情报检索现代情报检索技术是在矩阵理论的基础上发展起来的。
情报中心的数据库中存放着大量文件，我们希望从中搜索到与自己特定关键词相匹配的文件。
假设数据库中包含了n 个文件，而搜索所用的关键词有m 个，把关键词排序，我们就可以把数据库表示为m × n 矩阵A ，称矩阵A 为搜索矩阵。
若第i 个关35302520151050-10 -5 0 5 10 15 20 25 30杨威 等键词出现在第j 个文件中，则矩阵A 的第i 行第j 列元素为1，否则为0。
用于搜索的关键词清单可用m 维列向量x 表示，称为关键词搜索向量。
如果关键词清单中第i 个关键词在搜索列中出现，则x 的第i 个元素就为1，否则就为0。
当确定了搜索矩阵A 和搜索向量x 后，进行线性变换运算：T y A x =，其中向量y 即为检索结果[1]。
例2. 若数据库包含以下四个书名：线性代数，线性代数及其应用，线性代数与解析几何，矩阵代数及其应用。
而有五个搜索关键词：几何，代数，线性，矩阵，应用。
请写出搜索矩阵A 。
如果某读者输入的关键词为“代数，矩阵”，请写出关键词搜索向量x ，并计算搜索结果向量y ，最后说明搜索结果向量的含义。
解：当第i 个关键词出现在第j 本书名上时，矩阵A 的元素i j a 就等于1，否则就等于0。
表2给出了搜索矩阵构造表。
Table 2. Search matrix table of Information retrieval 表2. 情报检索搜索矩阵表关键词 书线性代数线性代数及其应用线性代数与解析几何矩阵代数及其应用几何 0 0 1 0 代数 1 1 1 1 线性 1 1 1 0 矩阵 0 0 0 1 应用11根据表2得搜索矩阵为00101111111000010101=A ，关键词搜索向量为01010=x ，用MATLAB 计算搜索结果向量T y A x =。
在MATLAB 命令窗口输入：A=[0,0,1,0;1,1,1,1;1,1,1,0;0,0,0,1;0,1,0,1] x=[0;1;0;1;0] y=A'*x 计算结果为： y = 1 1 1 2y 的各个分量表示各书与关键词搜索向量x 的匹配程度，即y 的第k 个分量表示第k 本书所含关键词的个数，y 的第四个分量是2，说明第四本书包含了两个关键词，故第四本数匹配程度最高。
4. 细菌转换人口迁移问题、基因突变问题、细菌转换问题等都属于线性代数的同一类应用实例，此类问题的解杨威 等决方法是先构造转换矩阵，然后根据相似对角化来求矩阵的高次幂，最后利用MATLAB 软件计算和分析结果[1]。
例3. 李博士培养了一罐细菌，在这个罐子里存放着A 、B 、C 三类不同种类的细菌，最开始A 、B 、C 三种细菌分别有108、2 × 108、3 × 108个。
它们每天都要发生类型转化，转化情况如下：A 类细菌一天后有5%的变为B 类细菌、15%的变为C 细菌；B 类细菌一天后有30%的变为A 类细菌、10%的变为C 类细菌；C 类细菌一天后有30%的变为A 类细菌、20%的变为B 类细菌。
请问一周后李博士的A 、B 、C 类细菌各有多少个？随着时间的增加，细菌数量会不会保持一个稳定状态？解：表3给出了第n 天与第n + 1天三种细菌相互变换情况。
Table 3. Table of conversion of three bacteria 表3. 三种细菌相互变换情况表A nB nC n A n +1 0.8A n 0.3B n 0.3C n B n +1 0.05A n 0.6B n 0.2C n C n +10.15A n0.1B n0.5C n用向量n x 来表示第n 天三种细菌的个数，即n n n n A x B C=；用向量1n x +来表示第n +1天三种细菌的个数，即1111n n n n A x B C ++++ = ，那么根据表3有矩阵关系：1110.80.30.30.050.60.20.150.10.5n n n n n n A A B B C C +++ = ，设0.80.30.30.050.60.20.150.10.5P = ，则有：277650x Px P x P x ==== ，其中8011023x =。
在MATLAB 命令窗口输入： p=[0.8,0.3,0.3;0.05,0.6,0.2;0.15,0.1,0.5] x0=10^8*[1;2;3] x7=p^7*x0 计算结果为 x7 =357968750 112555355 1294758957x 的三个元素分别为一周后李博士A 、B 、C 三种细菌数量。