1，解方程组

rank 函数判定解的情况

NULL

rref-------简化的行阶梯形矩阵（Gauss-Jordan 消去法）

(1)唯一解

（2）唯一解

（3）无解

（4）求通解

使用rref求解：

(4) 特征值与特征向量

设a=(1,2,3),b=(2,4,3), 分别计算a./b, a.\b, a/b, a\b, 分析结果的意义。

1,对应元素分别运算a除以b

2,对应元素分别运算b除以a

3,A/B是A乘B的逆矩阵的意思（/x=乘以x的逆）

4,B\A是B的逆矩阵乘A的意思。（x\=x的逆乘以）

1，解方程组

2. 用矩阵除法解下列线性方程组，并判断解的意义，用矩阵乘法验算。

rank 函数判定解的情况

1）当方程组的系数矩阵的秩与方程组增广矩阵的秩相等且均等于方程组中未知数个数n的时候，方程组有唯一解

2）当方程组的系数矩阵的秩与方程组增广矩阵的秩相等且均小于方程组中未知数个数n的时候，方程组有无穷多解

3）当方程组的系数矩阵的秩小于方程组增广矩阵的秩的时候，方程组无解

NULL

z=null(a) {返回 A 的零空间的标准正交基}

>>a矩阵*x=0

>>解x

rref-------简化的行阶梯形矩阵（Gauss-Jordan 消去法）

用途为解方程组

(1)唯一解

>> A=[4 1 -1;3 2 -6;1 -5 3];b=[9;-2;1];>> rank(A), rank([A,b])
//[A,b]为增广矩阵, 求rank的目的是为了判断解的意义，注意不要写成rank(A,b)。ans =3ans =3
//可见方程组唯一解>> x=A\bx =2.38301.48942.0213

（2）唯一解

>> A=[4 -3 3;3 2 -6;1 -5 3];b=[-1;-2;1];>> rank(A), rank([A,b])ans =3ans =3   %可见方程组唯一解>> x=A\bx =-0.4706-0.29410

（3）无解

>> A=[4 1;3 2;1 -5];b=[1;1;1];>> rank(A), rank([A,b])ans =2ans =3  %可见方程组无解>> x=A\bx =0.3311-0.1219 %最小二乘近似解

（4）求通解

clear;clca=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1 2 3]';%注意b的写法ra=rank(a),rab=rank([a,b])%结果：ra = 3rab = 3%ra==rab<4说明有无穷多解x0=a\b%结果：x0 =1010%x0是一个特解v=null(a)%结果：v =-0.62550.6255-0.20850.4170%{v}为齐次方程ax=0的基础解系，从而x=cv+x0(c为任意常数)为方程的通解

使用rref求解：

>> rref([a,b])ans =1.0000         0         0    1.5000    1.00000    1.0000         0   -1.5000         00         0    1.0000    0.5000    1.0000%前4列分别对应x1-x4, 通解x=k*[-1.5,1.5,-0.5,1]'+[1,0,1,0]'

4. (人口流动趋势) 对城乡人口流动作年度调查，发现有一个稳定的朝向城镇流动的趋势，每年农村居民的5%移居城镇而城镇居民的1%迁出，现在总人口的20%位于城镇。假如城乡总人口保持不变，并且人口流动的这种趋势继续下去，那么

(1) 一年以后住在城镇人口所占比例是多少？两年以后呢？十年以后呢？

(2) 很多年以后呢？

(3) 如果现在总人口70%位于城镇，很多年以后城镇人口所占比例是多少？

(4) 计算转移矩阵的最大特征值及对应的特征向量，与问题(2)(3)有何关系？

clc,clear
city=0.2;country=0.8;
y=10;
for i=1:yif i==1 | i==2 |i==10icity=city*0.99+country*0.05country=country*0.95+city*0.01end
end

很多年之后

city=0.2;country=0.8;
t=0;
y=100000
for i=1:yif country==country*0.95+city*0.01breakelsei,city=city*0.99+country*0.05,country=country*0.95+city*0.01end
end

i =

584

city =

0.8387

country =

0.1677

（3）如果现在总人口70%位于城镇，很多年以后城镇人口所占比例是多少？

（稍作改良）

city=0.7;country=0.3;
t=0;
y=100000
for i=1:yif country-country*0.95+city*0.01<=1e-6breakelsei,city=city*0.99+country*0.05,country=country*0.95+city*0.01end
end

(4) 特征值与特征向量

[V,D]= eig(A) 返回特征值的对角矩阵 D 和矩阵 V，其列是对应的右特征向量，使得 A*V = V*D。

//(4) 计算转移矩阵的最大特征值及对应的特征向量，与问题(2)(3)有何关系？x0=[0.2 0.8]';a=[0.99 0.05;0.01 0.95];x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667

>> [v,e]=eig(a)v =0.9806   -0.70710.1961    0.7071e =1.0000         00    0.9400>> v(:,1)./xans =1.17671.1767    %v(:,1)与x成比例，说明x也是最大特征值1对应的特征向量

MATLAB experiment＜3＞相关推荐

MATLAB experiment＜6＞
1,解非刚性微分方程 [t,y] = ode45(odefun,tspan,y0)(其中 tspan = [t0 tf])求微分方程组 y′=f(t,y) 从 t0 到 tf 的积分,初始条件为 y0 ...
＜binding＞和＜operation＞元素
<binding>和<operation>元素 Binding栏是完整描述协议.序列化和编码的地方,Types, Messages和PortType栏处理抽象的数据内容,而Bi ...
水鱼学习回顾＜1＞
学习回顾 <1> 1.注释 2.选择语句 3.循环语句 4.函数 5.数组 6.操作符 7.关键字 8.其他穿插知识点 1.注释鼠标选中代码块 Ctrl+K+C 快速注释 int num ...
封装自定义的redis切库工具类ByteArrayRedisTemplate，读取byte数组反序列化成List＜Object＞
封装自定义的redis切库工具类ByteArrayRedisTemplate,读取byte数组反序列化成List<Object>(使用lettuce连接池) 代码环境框架:springb ...
mybatis中＜mappers＞，mapperLocations，和MapperScannerConfigurer 用法
mybatis中<mappers> ,mapperLocations,和MapperScannerConfigurer 用法 1. mappers标签在mybatis单独使用时,myba ...
训练时发生的错误：Couldn‘t open shared file mapping: ＜000001910A228862＞, error code: ＜1455＞
今天在运行TTSR时第14个epoch时发生了以下错误: 看到一篇博客有写: RuntimeError: Couldn't open shared file mapping: <torch_37 ...
C++头文件＜functional＞和bind、placeholders占位符使用简单例子
相关博文:C++新特性探究(十五):bind 相关博文:C++头文件<functional>和bind.placeholders占位符使用简单例子相关博文:<Essential C ...
opencv中 .at＜uchar＞（）和.ptr＜uchar＞（）使用方法的区别
opencv中 .at<uchar>()和.ptr<uchar>()使用方法的区别在opencv中,.at<uchar>()和.ptr<uchar>() ...
vector＜int＞ v, vector＜int＞ v(n) 与vector＜int＞ v[n]的区别
一.vector<int> v表示声明一个容器v,并且不给他预定存储空间. 运行: vector<int> v;cout << v.size(); 输出为: 0 二 ...

MATLAB experiment＜3＞

1，解方程组

rank 函数判定解的情况

NULL

rref-------简化的行阶梯形矩阵（Gauss-Jordan 消去法）

(1)唯一解

（2）唯一解

（3）无解

（4）求通解

使用rref求解：

(4) 特征值与特征向量

MATLAB experiment＜3＞相关推荐

最新文章

热门文章

MATLAB experiment＜3＞

1，解方程组

rank 函数判定解的情况

NULL

rref-------简化的行阶梯形矩阵（Gauss-Jordan 消去法）

(1)唯一解

（2）唯一解

（3）无解

（4） 求通解

使用rref求解：

(4) 特征值与特征向量

MATLAB experiment＜3＞相关推荐

最新文章

热门文章

（4）求通解