MIT_线性代数笔记:第 24 讲 马尔可夫矩阵;傅里叶级数

news2024/11/18 9:34:55

目录

  • 马尔可夫矩阵 Markov matrices
  • 傅里叶级数 Fourier series

本讲学习马尔可夫矩阵和傅里叶级数,两者是关于特征值和投影矩阵的应用。

马尔可夫矩阵 Markov matrices

A = [ 0.1 0.01 0.3 0.2 0.99 0.3 0.7 0 0.4 ] A = \begin{bmatrix} 0.1 & 0.01 & 0.3 \\ 0.2 & 0.99& 0.3\\ 0.7 & 0& 0.4 \end{bmatrix} A= 0.10.20.70.010.9900.30.30.4
形如矩阵 A,任何元素非负,且每列的元素加和为 1 的矩阵被称为马尔可夫矩阵。马尔可夫矩阵主要应用在概率领域。将一个马尔可夫矩阵进行方幂运算得到的仍旧是马尔可夫矩阵。

当处理一个微分方程问题时,特征值 0 意味着得到一个稳态。当进行矩阵的方幂运算时,特征值给出稳态的条件包括:

  1. λ 1 λ_1 λ1=1 是特征值之一。
  2. 其它特征值的绝对值都比 1 小, ∣ λ i ∣ \begin{vmatrix} λ_i \end{vmatrix} λi <1。

如我们所知,如果矩阵具有 n 个线性无关的特征向量,则有:
U k = A k U 0 = c 1 λ 1 k X 1 + c 2 λ 2 k X 2 + . . . . . . + c n λ n k X n U_k = A^kU_0 =c_1λ_1^kX_1 + c_2λ_2^kX_2 + ...... +c_nλ_n^kX_n Uk=AkU0=c1λ1kX1+c2λ2kX2+......+cnλnkXn

如果 λ 1 λ_1 λ1=1 并且其他的特征值都小于 1,则系统在 k 增大过程中趋近于 u 0 u_0 u0的分量 c 1 x 1 c_1x_1 c1x1,即给出了一个稳态状况。这里特征向量 x 1 x_1 x1的每一分量都是正的,因此若初始值为正,则最终的稳态也是正的。

Markov 矩阵每一列的元素加和为 1 这个条件,保证了矩阵具有 1 这个特征值。
A − I = [ − 0.9 0.01 0.3 0.2 − 0.09 0.3 0.7 0 − 0.6 ] A - I = \begin{bmatrix} -0.9 & 0.01 & 0.3 \\ 0.2 & -0.09& 0.3\\ 0.7 & 0& -0.6 \end{bmatrix} AI= 0.90.20.70.010.0900.30.30.6
从每一列减去 1,则每列的加和都从 1 变为 0。这时候行向量相加的结构就是 0向量,因此行向量线性相关,矩阵为奇异矩阵。矩阵 A 有特征向量在 A-I 的零空间中,其对应的特征值为 1。回带计算可得 x 1 x_1 x1 = [ 0.6 33 0.7 ] \begin{bmatrix} 0.6\\33\\0.7 \end{bmatrix} 0.6330.7

我们用马尔可夫矩阵来研究人口流动问题。
[ u c a l u m a s s ] t = k + 1 = [ 0.9 0.2 0.1 0.8 ] [ u c a l u m a s s ] t = k \begin{bmatrix} u_{cal} \\ u_{mass} \end{bmatrix}_{t=k+1} = \begin{bmatrix} 0.9 & 0.2 \\ 0.1 & 0.8 \end{bmatrix} \begin{bmatrix} u_{cal} \\ u_{mass} \end{bmatrix}_{t=k} [ucalumass]t=k+1=[0.90.10.20.8][ucalumass]t=k

方程中 u 的分量分别代表加利福尼亚州和马萨诸塞州的人口,矩阵中的每一列中元素代表着人口去留比例,比如第一列 0.9 表示留在加州的人口占加州人口的90%,而 10%进入麻省,第二列中由麻省进入加州的人口占麻省的 20%,而 80%选择留在麻省。可以看到列向量分量的加和为 1 保证了整体人数不会变化,而在这种问题中矩阵也不会出现负的元素。

如果取初值 [ u c a l u m a s s ] 0 \begin{bmatrix} u_{cal} \\ u_{mass} \end{bmatrix}_{0} [ucalumass]0 = [ 0 1000 ] \begin{bmatrix} 0 \\ 1000 \end{bmatrix} [01000]

则经过一次迁徙 [ u c a l u m a s s ] 1 \begin{bmatrix} u_{cal} \\ u_{mass} \end{bmatrix}_{1} [ucalumass]1 = [ 0.9 0.2 0.1 0.8 ] \begin{bmatrix} 0.9 & 0.2 \\ 0.1& 0.8 \end{bmatrix} [0.90.10.20.8] [ 0 1000 ] \begin{bmatrix} 0 \\ 1000 \end{bmatrix} [01000] = [ 200 800 ] \begin{bmatrix} 200 \\ 800 \end{bmatrix} [200800]

为了获取长时间后的人口分布,我们需要了解矩阵的特征值和特征向量。因为这是马尔可夫矩阵,所以有一个特征值 1,则另一个特征值为

0.9+0.8-1=0.7。可以求得 x1= [ 2 1 ] \begin{bmatrix} 2\\1 \end{bmatrix} [21],x2= [ 1 − 1 ] \begin{bmatrix} 1\\-1 \end{bmatrix} [11]。从 [ 2 1 ] \begin{bmatrix} 2\\1 \end{bmatrix} [21],可知最后的稳态为加州人口 2/3,麻省人口 1/3。

通解为: u k = c 1 [ 2 1 ] + c 2 ( 0.7 ) k [ − 1 1 ] u_k = c_1 \begin{bmatrix} 2 \\ 1 \end{bmatrix} + c_2(0.7)^k \begin{bmatrix} -1 \\ 1 \end{bmatrix} uk=c1[21]+c2(0.7)k[11] ,可以从 u0解得 c1=1000/3,c2=2000/3。

傅里叶级数 Fourier series

如果有一组标准正交基 q 1 , q 2 … … q n q_1,q_2……q_n q1q2……qn为则任意向量 v 可以写成:

v = x 1 q 1 + x 2 q 2 … … x n q n v = x_1q_1+x_2q_2……x_nq_n v=x1q1+x2q2……xnqn
因为当 i,j 不相等时有 q i T q j q_i^Tq_j qiTqj=0。因此有 q i T v = x 1 q i T q 1 + x 2 q i T q 2 … … x n q i T q n = x i q_i^Tv=x_1q_i^Tq_1+x_2q_i^Tq_2……x_nq_i^Tq_n=x_i qiTv=x1qiTq1+x2qiTq2……xnqiTqn=xi
我们得到了分量 xi的公式: x i = q i T v x_i=q_i^Tv xi=qiTv
因为 v = [ q 1 . . . q n ] [ x 1 . . . x n ] v=\begin{bmatrix} q_1&...& q_n\end{bmatrix}\begin{bmatrix} x_1\\...\\ x_n\end{bmatrix} v=[q1...qn] x1...xn ,即 v = Q x v=Qx v=Qx,所以 x = Q − 1 v x =Q^{-1}v x=Q1v

因为 Q 为正交矩阵,所以有 Q − 1 = Q T Q^{-1}=Q^T Q1=QT x = Q T v x=Q^Tv x=QTv,这与我们之前得到的 x i = q i T v x_i=q_i^Tv xi=qiTv 完全相同。这里给出了求分量的思路就是用空间的一组标准正交基去点乘目标向量,利用其标准正交的性质得到所求。

标准正交是此处的核心概念。而傅里叶级数也是在这个概念上构建的。我们可以对任意函数做傅里叶展开,得到表达式:
f ( x ) = a 0 + a 1 c o s x + b 1 s i n x + a 2 c o s 2 x + b 2 s i n 2 x + . . . f(x) = a_0 + a_1cosx +b_1sinx +a_2cos2x +b_2sin2x +... f(x)=a0+a1cosx+b1sinx+a2cos2x+b2sin2x+...

与之前的有限个标准正交向量组成的正交矩阵不同,这个空间是无限维,它的一组基是 1,cosx,sinx,cos2x,sin2x……
此处的正交概念与 Rn空间不同,点积的概念也不同。
向量 V T W = v 1 w 1 + v 2 w 2 + . . . + v n w n 向量 V^TW =v_1w_1+v_2w_2+...+v_nw_n 向量VTW=v1w1+v2w2+...+vnwn
函数 f T g = ∫ 0 2 π f ( x ) g ( x ) d x 函数 f^Tg =\int^{2π}_0{f(x)g(x)}{\rm d}x 函数fTg=02πf(x)g(x)dx
在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1359109.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Python爬虫实战技巧:如何在爬取过程中动态切换代理IP

目录 前言 第一步&#xff1a;获取代理IP列表 第二步&#xff1a;测试代理IP的可用性 第三步&#xff1a;动态切换代理IP 总结 前言 在进行爬虫开发的过程中&#xff0c;有时候需要使用代理IP来访问目标网站&#xff0c;以避免被封IP或者降低访问频率的限制。本文将介绍如…

菜单显示博客页

结果图 难点及其实现 滚动条自定义&#xff1a; .card::-webkit-scrollbar { width: 4px;} .card::-webkit-scrollbar-thumb { border-radius: 8px; box-shadow: inset 0 0 5px rgba(0, 0, 0, 0.2); background: #666666; } .card::-webkit-scrollbar-thumb:hover { backgr…

c# windows10大小端试

测试代码&#xff1a; unsafe public void ceshi() {byte[] by BitConverter.GetBytes(0x12345678);Debug.WriteLine(" byte[0] 0x" by[0].ToString("x2"));Debug.WriteLine(" byte[1] 0x" by[1].ToString("x2"));Debug.WriteLi…

MinIO 批处理框架添加了对过期时间的支持

您现在可以使用 MinIO 批处理框架执行 S3 删除操作&#xff0c;以通过单个 API 请求删除大量对象。借助 MinIO 批处理框架&#xff0c;可以快速轻松地在 MinIO 部署中执行重复或批量操作&#xff0c;例如批量复制和批量密钥轮换。MinIO 批处理框架处理所有手动工作&#xff0c;…

【自动化测试总结】优点、场景、流程、项目人员构成

一、自动化测试的概念 以程序测试程序&#xff0c;以代码代替思维&#xff0c;以脚本的运行代替手工测试&#xff0c;可以大大提高工作测试的效率。 二、自动化测试的优点 1.回归测试更为方便&#xff0c;可靠。自动化测试最主要的任务和特点&#xff0c;特别是在程序修改比较…

Minitab 20安装包下载及安装教程

Minitab 20下载链接&#xff1a;https://docs.qq.com/doc/DUmNYVGxtUnZkWUpk 1.选中下载好的安装包&#xff0c;鼠标右键解压到”Minitab 20“文件夹 2.选中Setup&#xff0c;鼠标右击选择“以管理员身份运行” 3.点击“下一步” 4.点击“下一步” 5.勾选我接受许可协议中的条…

oracle语法学习

oracle语法学习 1.备份表 create table bd_psndoc_temp as select * from bd_psndoc2.还原表 drop table bd_psndoc; create table bd_psndoc as select * from bd_psndoc_temp3.查询表的前5条记录 select * from bd_psndoc_temp where rownum<54.从一个表中复制所有的列…

stm32学习笔记:TIM-输出比较

四部分讲解内容&#xff0c;本文是第二部分 输出比较主要用于PWM波形的控制电机&#xff08;驱动电机的必要条件&#xff09; 1、定时器基本定时&#xff0c;定一个时间&#xff0c;然后让定时器每隔一段时间产生一个中断&#xff0c;来实现每隔一个固定时间执行一段程序的目…

HarmoryOS Ability页面的生命周期

接入穿山甲SDK app示例&#xff1a; android 数独小游戏 经典数独休闲益智 广告接入示例: Android 个人开发者如何接入广告SDK&#xff0c;实现app流量变现 Ability页面的生命周期 学习前端&#xff0c;第一步最重要的是要理解&#xff0c;页面启动和不同场景下的生命周期的…

老师们居然这样把期末考试成绩发给家长——群发成绩系统的妙用

在现代社会&#xff0c;随着科技的发展和教育的进步&#xff0c;家校之间的沟通方式也在不断更新和完善。其中&#xff0c;群发成绩系统作为一种新型的教育信息化工具&#xff0c;在期末考试成绩的传递方面发挥了重要作用。老师们通过群发成绩系统&#xff0c;将期末考试成绩快…

合合TextIn团队发布 - 文档图像多模态大模型技术发展、探索与应用

合合信息TextIn&#xff08;Text Intelligence&#xff09;团队在2023年12月31日参与了中国图象图形学学会青年科学家会议 - 垂直领域大模型论坛。在会议上&#xff0c;丁凯博士分享了文档图像大模型的思考与探索&#xff0c;完整阐述了多模态大模型在文档图像领域的发展与探索…

在pycharm中执行 os.makedirs 提示用户名或密码不正确

问题&#xff1a;在pycharm中运行脚本&#xff0c;在 \10.0.21.249\share 共享目录下创建目录提示错误 发现&#xff1a;手动在该目录下创建目录没有问题。 解决方法&#xff1a; 切换到cmd 命令行运行该脚本成功创建 猜测&#xff1a;感觉应该是pycharm中使用的用户名和密码存…

医疗设备产生谐波对配电系统危害有哪些

1.对电网的影响 电网中的谐波会增加系统损耗&#xff0c;使设备发热&#xff0c;影响设备使用寿命。此外当并联的无功补偿电容器容抗与系统感抗匹配时&#xff0c;将发生n次并联谐振&#xff0c;造成电容器组的过电压和过电流。当基波频率为时&#xff0c;谐振频率&#xff0c…

【重磅新品】小眼睛科技推出紫光同创盘古系列FPGA开发板套件,盘古200K开发板,紫光同创PG2L200H,Logos2系列

FPGA&#xff0c;即现场可编程门阵列&#xff0c;作为可重构电路芯片&#xff0c;已经成为行业“万能芯片”&#xff0c;在通信系统、数字信息处理、视频图像处理、高速接口设计等方面都有不俗的表现。近几年&#xff0c;随着国家战略支持和产业发展&#xff0c;国产FPGA迎来迅…

LLM(九)| 使用LlamaIndex本地运行Mixtral 8x7大模型

欧洲人工智能巨头Mistral AI最近开源Mixtral 8x7b大模型&#xff0c;是一个“专家混合”模型&#xff0c;由八个70亿参数的模型组成。Mistral AI在一篇博客文章&#xff08;https://mistral.ai/news/mixtral-of-experts/&#xff09;介绍了Mixtral 8x7b&#xff0c;在许多基准上…

代码随想录算法训练营第十五天| 二叉树 513. 找树左下角的值 112. 路径总和 106.从中序与后序遍历序列构造二叉树

513. 找树左下角的值 层序遍历 本题用层序遍历可以直接秒了&#xff0c;直接提取每一层中最左边的元素&#xff08;i0&#xff09;&#xff0c;然后保存到最后一层即可。 class Solution { public:int findBottomLeftValue(TreeNode* root) {queue<TreeNode*> que;int…

【金猿人物展】东方金信董事长兼总经理王伟哲:价值化、智能化驱动下的大数据,治理仍需深化...

‍ 王伟哲 本文由东方金信董事长兼总经理王伟哲撰写并投递参与“数据猿年度金猿策划活动——2023大数据产业年度趋势人物榜单及奖项”评选。 大数据产业创新服务媒体 ——聚焦数据 改变商业 2023年&#xff0c;数据要素、公共数据授权运营成为热点话题&#xff0c;也是激活数据…

Flume基础知识(二):Flume安装部署

1. Flume 安装部署 1.1 安装地址 &#xff08;1&#xff09;Flume 官网地址&#xff1a;Welcome to Apache Flume — Apache Flume &#xff08;2&#xff09;文档查看地址&#xff1a;Flume 1.11.0 User Guide — Apache Flume &#xff08;3&#xff09;下载地址&#xf…

Halcon计算封闭区域(孔洞)的面积area_holes

Halcon计算封闭区域&#xff08;孔洞&#xff09;的面积 除了可以用area_center 算子计算区域的面积以外&#xff0c;在Halcon中还可以使用area_holes算子计算图像中封闭区域&#xff08;孔洞&#xff09;的面积。该面积指的是区域中孔洞部分包含的像素数。一个区域中可能不只…

深入了解Swagger注解:@ApiModel和@ApiModelProperty实用指南

在现代软件开发中&#xff0c;提供清晰全面的 API 文档 至关重要。ApiModel 和 ApiModelProperty 这样的代码注解在此方面表现出色&#xff0c;通过增强模型及其属性的元数据来丰富文档内容。它们的主要功能是为这些元素命名和描述&#xff0c;使生成的 API 文档更加明确。 Api…