Math 板


LINE

请问一下, 我想证明: 线性实数常系数齐次递回方程如果其每个特徵根的绝对值都小於1 则解都趋近於0 我知道有个证明是找到解空间的基底的长相, 然後就得证 但是这个证法要完整地写出来很复杂, 因为特徵多项式的解有很多case要讨论: (1) 重根/无重根, 重根的话又要看general eigenspace大小 (2) 实根/虚根 粗糙的来说, 如果根是实根λ, 解会是λ^n*(c_0+c_1*n+c_2*n^2+...)这种形式 根是虚根λ, 解会是|λ|^n(cos+sin...)这种形式 最复杂的大概是同时是虚根而且又有重根 而只要|λ|<1, 上述的解都会趋近於0 as n→+∞ --------------------------------------------------- 有监於以上, 我才好奇有没有比较简短的证法 因为有些其他性质我确实用方程式本身就能证出来, 不用透过解的长相 我才觉得这篇的问题应该有机会有其他证法 以下补充严格的数学描述: ================================================================== 令 a_1~a_N€R, a_N !=0 y_n 满足 y_n + a_1*y_(n-1) + ... + a_N*y(n-N) = 0, for all n >= N p(x):= x^N + a_1*x^(N-1) + ... + a_N 为特徵多项式 若 over C来看, p(x)所有的根的绝对值都小於1 则 lim_{n→+∞} y_n = 0 ================================================================== 谢谢帮忙~ --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 59.102.225.191 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Math/M.1670079252.A.C4B.html
1F:推 Vulpix : 可是特徵根都比1小这件事,你也很难用原方程写啊。 12/03 23:00
嗨V大~是的, 我就是自己试了一些方向都没办法用到特徵根 甚至我直接假设y_n极限存在为L, 也只是得到(1+a_1+...+a_N)*L = 0 完全没有特徵根资讯 因此我才想问是不是有其他方式
2F:→ Vulpix : 直接用通解的公式可行,而且没有你说的那麽复杂。 12/03 23:00
有不复杂的讨论方式可以分享吗!? 我记得之前线性ODE时, 2x2是最简单的, 3x3就要讨论几种case了 而老师跟课本举例也只到3x3, 因为general的NxN有太多根的情形了... 如果有简单写出通解的证明我也想知道, 谢谢!
3F:推 Vulpix : 每一项也就λ^n乘上n的多项式而已。都会收敛到0。 12/03 23:23
4F:→ Vulpix : 其实重数、虚实都不用顾虑,全部一起做就好。 12/03 23:24
"每一项也就λ^n乘上n的多项式" 这结论对我来说就是要分析每个λ的重数才能得到 还是V大你意思是, 只要λ是特徵多项式的m重根, 那就存在m个数列 y_k_n = (λ^n)*(n^(k-1)), k=1~m 是线性独立的 如果λ是虚数再换成cos, sin即可? 如果是这样的话确实很快耶! 我文中说要讨论general eigenspace说很复杂的意思是 我把线性方程转成矩阵形式来看 Y_n = A*Y_(n-1), 其中A是矩阵 而A的特徵多项式就跟差分方程的特徵多项式一样 因此接下来要对A做Jordan form时就会需要讨论每个λ的情形 看来直接用原差分方程做会避开很多矩阵会遇到的麻烦耶!
5F:推 wohtp : y_n = A y_(n-1) + x_n 这是你的递回式? 12/04 21:44
6F:→ wohtp : 有 x 在怎麽会收敛? 12/04 21:44
w大是说我上面回V大的那串吧 那我笔误了 我这篇都是讨论齐次的 我改一下 谢谢
7F:推 Vulpix : 并不须要真的知道重数,反正这只是some function of 12/04 22:53
8F:→ Vulpix : A and λ。 12/04 22:53
9F:推 Vulpix : 这种矩阵已经是自己的rational canonical form了。 12/04 22:56
10F:→ Vulpix : 所以每个eigenvalue只有一个Jordan block。 12/04 22:57
11F:→ Vulpix : 然後回你回应第一段:是。理由就是上面这段。 12/04 22:59
嗨V大, 非常谢谢你~了解了! 听你的解释确实给了【原本我在乎A的Jordan form还要讨论很多】一个不用讨论的理由 A因为长相的关系, 导致其Jordan form很简单(每个eigenvalue只有一个Jordan block) 所以不管直接证线性独立, 或是真的找A的Jordan form, 都是简单的事 再次感谢~
12F:推 alan23273850: 我想借串问为何一元高次递回式可以用 ODE 通解,二 12/05 01:53
13F:→ alan23273850: 元一次递回式可以用 power series 去解 ~ 12/05 01:53
嗨a大, 你的一元高次跟二元一次分别是? 线性递回的 y_n = a_1*y_(n-1) + ... + a_N*y_(n-N) + x_n 就是你的一元高次吗? 如果是的话, 我是不知道如何解释"为什麽" 单纯就是令 Y_n:=[y_n y_(n-1)...y_(n-N)]^t, X_n = [x_n 0...0] 的话 原方程很容易证明等价於Y_n = A*Y_(n-1) + X_n
14F:推 alan23273850: 所以他跟 ODE 是有关系的吗??? 12/05 16:31
是很常听到一阶线性常系数ODE的离散型就是差分方程啦 不过细节上还是有差, 但是也是有相似的地方 比如这篇的差分方程写成矩阵: Y_n = A*Y_(n-1) + X_n 而一阶线性ODE: Y(t)' = A*Y(t) + X(t) 通解表达式都会跟求A的eigenvalue有关 不过ODE的A是任意矩阵, 差分方程的A是很多0的矩阵 所以差分方程的A解起来不像ODE那样要讨论一堆A的eigenvalue的重数问题 ※ 编辑: znmkhxrw (59.102.225.191 台湾), 12/05/2022 18:03:12
15F:推 wohtp : 微分也是线性算子,很多结果可以直接套在ODE上面用 12/06 09:49
16F:→ wohtp : 其实就是向量空间维度变成无限这样 12/06 09:51
17F:推 alan23273850: 感谢楼上的开导! 12/06 15:39







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:Soft_Job站内搜寻

TOP