学点数学(2)-特征函数

特征函数

  • 1.数列特征方程
  • 2.矩阵特征方程
  • 3.微分方程特征方程
  • 4.积分方程特征方程

特征方程是为研究相应的数学对象而引入的一些等式,这些等式描述了特定对象的特性。依据研究的对象不同,特征方程包括数列特征方程、矩阵特征方程、微分方程特征方程、积分方程特征方程等等。

1.数列特征方程

数列的特征方程 可以用于求导 数列通项公式
数列的二阶线性推导:
an+1=pan+qan−1a_{n+1}=pa_n+qa_{n-1}an+1=pan+qan1
由待定系数法,对比下列等比公式:
an+1−tan=s(an−tan−1)a_{n+1}-ta_n=s(a_n-ta_{n-1})an+1tan=s(antan1)
可得:
{p=s+t...(1)y=st...(2)\left\{ \begin{aligned} p & = & s+t...(1) \\ y & = & st...(2)\\ \end{aligned} \right. {py==s+t...(1)st...(2)
一般情况下上面方程组有两组不同的解(t1,s1)(t2,s2)(t_1,s_1)(t_2,s_2)(t1,s1)(t2,s2)(特殊情况详见参考资料),依据等比数列通项公式可得吗:
{an+1−t1an=(a2−t1a1)s1n−1...(3)an+1−t2an=(a2−t2a1)s2n−1...(4)\left\{ \begin{aligned} a_{n+1}-t_1a_n& = & (a_2-t_1a_1)s_1^{n-1}...(3) \\ a_{n+1}-t_2a_n& = & (a_2-t_2a_1)s_2^{n-1}...(4)\\ \end{aligned} \right. {an+1t1anan+1t2an==(a2t1a1)s1n1...(3)(a2t2a1)s2n1...(4)
(3)-(4)化简得:
an=a2−t1a1(t2−t1)s1s1n+a2−t2a1(t2−t1)s1s2n=c1s1n+c2s2na_n=\frac{a_2-t_1a_1}{(t_2-t_1)s_1}s_1^n+\frac{a_2-t_2a_1}{(t_2-t_1)s_1}s_2^n=c_1s_1^n+c_2s_2^nan=(t2t1)s1a2t1a1s1n+(t2t1)s1a2t2a1s2n=c1s1n+c2s2n
即只需知道s1,s2s_1,s_2s1,s2,通过初始条件(a1,a2a_1,a_2a1,a2的已知值)求解系数c1,c2c_1,c_2c1,c2即可。

好像求到这里没有数列特征方程什么事啊!那数列的特征方程是怎么回事呢?其实,s1,s2s_1,s_2s1,s2是求数列特征方程的根。

通过求解第一个方程组可以求得s,将(2)式子带入(1)式,消去t,方程组变为求解下面二元一次方程:
s2−ps−q=0s^2-ps-q=0s2psq=0

此方程即定义为数列二阶线性推到的特征方程

以上推导过程,总结由数列的二次递推式an+1=pan+qan−1a_{n+1}=pa_n+qa_{n-1}an+1=pan+qan1 求解 通项公式的步骤:
step1:求解数列特征方程s2−ps−q=0−>(s1,s2)s^2-ps-q=0->(s_1,s_2)s2psq=0>(s1,s2)
step2:由初始条件求解c1,c2c_1,c_2c1,c2
step3: 通项公式为an=c1s1n+c2s2na_n=c_1s_1^n+c_2s_2^nan=c1s1n+c2s2n

参考资料:
https://wenku.baidu.com/view/8824638caef8941ea76e0592.html
https://blog.csdn.net/qq_20340417/article/details/78433961

2.矩阵特征方程

(矩阵的特征方程是我最经常看到的特征方程,在写这篇文章之前。我一直以为特征方程指的就是矩阵的特征方程。咦,闲话不多说)

矩阵的特征方程 可以用于求解矩阵的特征值。

每个事物都具有许多的特征,这些特征用于表示该事物的某些属性,矩阵也不例外地具有许多特征。在一个方阵A的众多特征中,其 特征值与特征向量 是两个比较典型的特征。一个方阵A的特征值λ\lambdaλ与特征向量xxx为满足下式:
Ax=λx...(2.1)Ax=\lambda x...(2.1)Ax=λx...(2.1)

上式的几何含义:对向量xxx做变换AAA,其效果等价于对xxx做伸缩变换,伸缩系数为λ\lambdaλ

特征值和特征向量的定义和几何意义很明确,那么如何求解x,λx,\lambdax,λ
呢?将(2.1)式移项得:
Ax−λx=0−>(A−λI)x=0...(2.2)Ax-\lambda x=0->(A-\lambda I)x=0...(2.2)Axλx=0>(AλI)x=0...(2.2)
如果特征向量x=[x1,x−2,...,xn]x=[x_1,x-2,...,x_n]x=[x1,x2,...,xn]为n维向量,那么对式(2.2)求解xxx的过程 等价于 求解一个n元一次方程组。n元一次方程组有非零解的 充要条件 为:系数矩阵的行列式为0,即:
∣A−λI∣=0...(2.3)|A-\lambda I|=0...(2.3)AλI=0...(2.3)

式子(2.3)即为矩阵A的特征方程。

以上推导过程,总结求解方阵A 特征值与特征多项式的步骤:
step1:求解特征方程:∣A−λI∣=0|A-\lambda I|=0AλI=0
step2:带λ\lambdaλ入(2.2)式求解特征向量xxx

参考资料:https://baike.baidu.com/item/%E7%9F%A9%E9%98%B5%E7%89%B9%E5%BE%81%E5%80%BC/8309765?fr=aladdin

3.微分方程特征方程

微分方程的特征方程可以用于求解 二阶 常系数 齐次 常微分方程。

微分方程指含有未知函数及其导数的关系式,如:
dydx=2x..(3.1)\frac{dy}{dx}=2x..(3.1)dxdy=2x..(3.1)

解微分方程就是找出未知函数,如:
y=x2+Cy=x^2+Cy=x2+C

CCC由微分方程的约束条件确定。

微分方程的阶数由方程中函数的几次导数决定,(3.1)式子为阶微分方程;只含有一个未知数的微分方程为常微分方程。

二阶 常系数 齐次 常微分方程:
y′′+py′+qy=0...(3.2)y''+py'+qy=0...(3.2)y+py+qy=0...(3.2)
只要满足式(3.2)的函数就是其解,可以假设y=erxy=e^{rx}y=erx,求解合适的rrr满足(3.2)式。(为啥设置y=erxy=e^{rx}y=erx的具体原因我忘记了)

y=erxy=e^{rx}y=erx带入(3.2)式子,有:
(r2+pr+q)erx=0...(3.2)(r^2+pr+q)e^{rx}=0...(3.2)(r2+pr+q)erx=0...(3.2)

因为:erxe^{rx}erx恒大于0,所以只有:r2+pr+q=0...(3.4)r^2+pr+q=0...(3.4)r2+pr+q=0...(3.4)

rrr才能使(3.3)式为0。式子(3.4)为二阶 常系数 齐次 常微分方程的特征方程,解特征方程得r1,r2r_1,r_2r1,r2(假定有两个不相等的实数解,其他情况,详见参考资料),那么y1=er1x,y2=er2xy_1=e^{r_1x},y_2=e^{r_2x}y1=er1x,y2=er2x为方程的两个解,这两个解的各种线性组合也是微分方程的解,所以,微分方程的通解为:
y=c1er1x+c2er2x...(3.5)y=c_1e^{r_1x}+c_2e^{r_2x}...(3.5)y=c1er1x+c2er2x...(3.5)

以上推导过程,总结求解 二阶常系数齐次常微分方程 的步骤:
step1:求解特征方程:r2+pr+q=0r^2+pr+q=0r2+pr+q=0
step2:依据初始条件求解系数c1,c2c_1,c_2c1,c2
step3:带入3.5式,二阶常系数齐次常微分方程的通解为:y=c1er1x+c2er2x...(3.5)y=c_1e^{r_1x}+c_2e^{r_2x}...(3.5)y=c1er1x+c2er2x...(3.5)

参考资料:
https://baike.baidu.com/item/%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B/4763?fr=aladdin
https://blog.csdn.net/low5252/article/details/90758604

4.积分方程特征方程

积分方程是含有 对未知函数的积分运算 的方程,与微分方程相对。如果积分方程中只含有位置。
积分方程中有核函数k(x,y)k(x,y)k(x,y),核函数可以展开成特征值与特征函数加权和
k(x,y)=∑i=1∞1λiψi(x)ψi(y)k(x,y)=\sum_{i=1}^{\infty}\frac{1}{\lambda_i}\psi_i(x)\psi_i(y)k(x,y)=i=1λi1ψi(x)ψi(y)

具体内容详见下面资料。(本部分还不会,逃!总有一天我会回来的

参考材料:https://wenku.baidu.com/view/344127efaf45b307e87197f2

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

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

相关文章

GCC如何产生core dump

先决条件1.安装apport(automatically generate crash reports for debugging)2.修改/etc/security/limits.conf文件,使允许core dump,或者用ulimit -c unlimited设置core dump文件的大小为unlimited3.C/C的编译开关-g(…

经典的进程同步问题

经典的进程同步问题 普通版:一类进程作为生产者,生产产品,生产的产品放入一个缓冲区,消费者从缓冲区中取出产品,需要保证生产者不可以向满的缓冲区中添加产品,消费者不可以从空的缓冲区中取出产品。同一时刻…

面试题汇总---深度学习(图像识别,NLP内容)

文章目录1.基本概念1.1 为什么神经网络中深度网络的表现比广度网络表现好?1.2 推导BP算法1.3 什么是梯度消失和梯度爆炸?1.4 常用的激活函数有哪些?1.5 常用的参数更新方法有哪些?1.6 解决过拟合的方法?数据层面模型层…

Linux-2.6.25 TCPIP函数调用大致流程

Linux-2.6.25 TCPIP函数调用大致流程学习目的,随手笔记。函数和文字说明会不断补充更新。Changelog2008.10.08 最近找工作忙。暂时缓缓插口层系统调用sendsys_sendsys_sendtosendtosys_sendtosock_sendmsgsendmsgsys_sendmsgsock_sendmsgwritesys_writevfs_write…

Python(28)-文件,os模块

文件1. 文件2. 文件的基本操作3. 读取文件open()3.1 文件指针: 标记从哪一个位置开始读取数据.3.2 文件的打开方式mode3.3 文件按行读取3.3.1 readline()3.3.2 readlines()4.文件输出f.write(),print()5.文件复制5.1 小文件复制(搬家)5.2 大文件复制&…

IOCP的程序

C代码 #include <winsock2.h> #include <mswsock.h> #include <windows.h> #include <stdio.h> #include <stdlib.h> #include <assert.h> #include "vld.h" #pragma message("automatic link to ws2_32.lib and…

PaperNotes(3)-图像分割-RCNN-FCN-Boxsup

图像分割算法对比小结1.{基本概念}2.{R-CNN}2.1R-CNN 网络结构选择性搜索算法为什么选择SVM作分类器边框回归2.2{R-CNN 训练}2.3{R-CNN实验结果}2.4{R-CNN语义分割}2.5{补充材料}2.5.1{R-CNN建议区域放缩}2.5.2{IOU阈值设置不一样的原因}2.5.3{Bounding-box回归修正}2.6{R-CNN存…

Python模块(3)--PIL 简易使用教程

PIL模块-用与记1.图片导入Image.open()2.图像显示.show()4.查看图片属性.format,.size,.mode3.图像格式转换.convert()4.图像模式“L”&#xff0c;“RGB”,"CYMK"5. 图片旋转.rotate()旋转方式1&#xff1a;旋转不扩展旋转方式2&#xff1a;旋转扩展旋转方式3&#…

日志级别 debug info warn eirror fatal

日志级别 debug info warn eirror fatal 软件中总免不了要使用诸如 Log4net, Log4j, Tracer 等东东来写日志&#xff0c;不管用什么&#xff0c;这些东东大多是大同小异的&#xff0c;一般都提供了这样5个日志级别&#xff1a; Debug Info Warn Error Fatal一个等级比一个高&…

输入输出系统

I/O设备&#xff1a;输入输出和存储功能的设备 I/O设备的分类 按传输的速度&#xff1a; 低速设备&#xff08;如键盘、鼠标、语音输入输出设备&#xff09; 中速设备&#xff08;如行式打印机、激光打印机等&#xff09; 高速设备&#xff08;如磁带机、磁盘机、光盘机等&…

vue2源码解析---v-model双向数据绑定

什么是v-model v-model 是 Vue 中的一个指令&#xff0c;用于实现表单元素与 Vue 实例中数据的双向绑定。这意味着当表单元素的值发生变化时&#xff0c;Vue 实例中的数据也会随之更新 工作原理 生成ast树 本质上是语法糖 结合了v-bind和v-on两个指令 示例代码 new Vue({e…

php收集的精典代码

1. οncοntextmenu"window.event.return&#xff06;#118aluefalse" 将彻底屏蔽鼠标右键 <table border οncοntextmenureturn(false)><td>no</table> 可用于Table 2. <body onselectstart"return false"> 取消选取、防止复制…

python外卷(7)--glob

glob模块1.glob.glob()2.对比os.listdir()glob是python自带的一个操作文件的模块&#xff0c;可用于查找 指定路径 中 匹配的 文件。1.glob.glob() 下面是一个测试文件路径&#xff1a; (base) pppp-System-Product-Name:~/Desktop/test_glob$ tree . ├── a │ ├── 1…

Sublime Text 2配置强大的IDE开发环境,运行java

Sublime Text 2是我无意中发现的、据说十分强大的、便捷的编辑器&#xff0c;许多程序员都投入到Sublime Text 2的怀抱中。 1 配置java开发环境的方法如下&#xff1a; 在jdk安装目录下的bin文件夹下新建一个bat格式的文件&#xff0c;文件命为javacexec.bat。 如果是在Wind…

thinkphp的快捷方法实例化对象

D、F、S、C、L、A、I 他们都在functions.php这个文件家 下面我分别说明一下他们的功能 D&#xff08;&#xff09; 加载Model类 M&#xff08;&#xff09; 加载Model类 A&#xff08;&#xff09; 加载Action类 L&#xff08;&#xff09; 获取语言定义 C&#xff08;&#xf…

Python外卷(8)--pdist, squareform

pdist, squareform1.pdist, squareform使用例子2.通过矩阵的四则运算实现上述pdist, squareformscipy.spatial.distance 距离计算库中有两个函数&#xff1a;pdist, squareform&#xff0c;用于计算样本对之间的欧式距离&#xff0c;并且将样本间距离用方阵表示出来。&#xff…

模拟进程调度

功能 data.h #ifndef _Data_h_ #define _Data_h_#include <stdio.h> #include <stdlib.h> #include <string.h>#define ElemType PCB #define Status int #define OK 1 #define ERROR 0 #define TimeSlice 1 #define Infinity 10 //INT_MAX#define NAME_M…

gdb调试多进程和多线程命令

1. 默认设置下&#xff0c;在调试多进程程序时GDB只会调试主进程。但是GDB&#xff08;>V7.0&#xff09;支持多进程的 分别以及同时 调试&#xff0c;换句话说&#xff0c;GDB可以同时调试多个程序。只需要设置follow-fork-mode(默认值&#xff1a;parent)和detach-on-fork…

python外卷(10)--取整

"取整"那些事1.python 内置函数1.1int()--向下取整1.2round()--四舍五入2.math模块取整函数2.1 math.floor()--向下取整2.2 math.ceil()--向上取整2.3 math.modf()--分别取小数部分和整数部分3.numpy模块取整函数3.1 numpy.floor()--向下取整3.2 numpy.ceil()--向上取…

模拟银行家算法

介绍 data.h #ifndef _Data_h_ #define _Data_h_#include <stdio.h> #include <stdlib.h> #include <string.h>#define ElemType PCB #define Status int #define true 1 #define false 0 #define OK 1 #define ERROR 0 #define RESOURCE_NUM …