學點數學(2)-特征函數

特征函數

  • 1.數列特征方程
  • 2.矩陣特征方程
  • 3.微分方程特征方程
  • 4.積分方程特征方程

特征方程是為研究相應的數學對象而引入的一些等式,這些等式描述了特定對象的特性。依據研究的對象不同,特征方程包括數列特征方程、矩陣特征方程、微分方程特征方程、積分方程特征方程等等。

1.數列特征方程

數列的特征方程 可以用于求導 數列通項公式
數列的二階線性推導:
an+1=pan+qan?1a_{n+1}=pa_n+qa_{n-1}an+1?=pan?+qan?1?
由待定系數法,對比下列等比公式:
an+1?tan=s(an?tan?1)a_{n+1}-ta_n=s(a_n-ta_{n-1})an+1??tan?=s(an??tan?1?)
可得:
{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+1??t1?an?an+1??t2?an??==?(a2??t1?a1?)s1n?1?...(3)(a2??t2?a1?)s2n?1?...(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?=(t2??t1?)s1?a2??t1?a1??s1n?+(t2??t1?)s1?a2??t2?a1??s2n?=c1?s1n?+c2?s2n?
即只需知道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=0s2?ps?q=0

此方程即定義為數列二階線性推到的特征方程

以上推導過程,總結由數列的二次遞推式an+1=pan+qan?1a_{n+1}=pa_n+qa_{n-1}an+1?=pan?+qan?1? 求解 通項公式的步驟:
step1:求解數列特征方程s2?ps?q=0?>(s1,s2)s^2-ps-q=0->(s_1,s_2)s2?ps?q=0?>(s1?,s2?)
step2:由初始條件求解c1,c2c_1,c_2c1?,c2?
step3: 通項公式為an=c1s1n+c2s2na_n=c_1s_1^n+c_2s_2^nan?=c1?s1n?+c2?s2n?

參考資料:
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?,x?2,...,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?=er1?x,y2?=er2?x為方程的兩個解,這兩個解的各種線性組合也是微分方程的解,所以,微分方程的通解為:
y=c1er1x+c2er2x...(3.5)y=c_1e^{r_1x}+c_2e^{r_2x}...(3.5)y=c1?er1?x+c2?er2?x...(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=c1?er1?x+c2?er2?x...(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?λi?1?ψi?(x)ψi?(y)

具體內容詳見下面資料。(本部分還不會,逃!總有一天我會回來的

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

本文來自互聯網用戶投稿,該文觀點僅代表作者本人,不代表本站立場。本站僅提供信息存儲空間服務,不擁有所有權,不承擔相關法律責任。
如若轉載,請注明出處:http://www.pswp.cn/news/445120.shtml
繁體地址,請注明出處:http://hk.pswp.cn/news/445120.shtml
英文地址,請注明出處:http://en.pswp.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 …