Math 板


LINE

想请问下列递回方程式是否有解的分类与规律: ------------------------------------------------------------- 令 x_n为一实数列, n>=0 y_(-1), y_(-2),..., y_(-N) 为N个初始值 a_1, a_2,..., a_N 为N个实数, a_N != 0 定义递回关系式 y_n = (Σ_{k=1~N} a_k*y_(n-k) ) + x_n, n>=0 想请问有没有关於 怎样的a_k与x_n的条件 会导致 y_n有怎样的行为, 特别是n→∞时的行为 或是条件更强一点, 加入" x_n = 0 for n>=M, some M " 会有更好的结果吗? (我是做实验时发现n很大时感觉y_n会趋近周期函数, 不知道是特例还是怎样) ------------------------------------------------------------- 因为是工作上遇到的, 对於推导不太要求, 想知道有没有确定的分类 目前尝试写成矩阵的方式, 可是得不到什麽定性的结论: (A^t代表A的transpose, "*"只是乘法而已) 令 Y_n:=[y_n, y_(n-1),..., y_(n-N+1)]^t X_n:=[x_n, 0, ..., 0]^t [a_1, a_2, ..., a_(N-1), a_N] [ 1 , 0 , ..., 0 , 0 ] A :=[ 0 , 1 , ..., 0 , 0 ] [ ..........................] [ 0 , 0 , ..., 1, , 0 ] 则 Y_n = A*Y_(n-1) + X_n, n>=0 一直拆下去可以得到, Y_n = Σ_{k=0~n} A^(n-k)*X_k 之後把A对角化(只要a_N != 0都可以对角化)得到 A = P*D*P^(-1), P可逆 代回去得到 Y_n = P*(Σ_{k=0~n} D^(n-k)*P^(-1)*X_k), n>=0 我就卡在这了, 即便讨论一些eigenvalue的值也得不出什麽结果 (把P过来Y_n这边後又导致了y_n被污染了, 拿不到y_n) 或是加入更强条件" x_n = 0 for n>=M, some M ", 我们得到了: Y_n = P*(Σ_{k=0~M} D^(n-k)*P^(-1)*X_k), n>=M 此时summation只有到M, 所以只要所有eigenvalue的绝对值都<1, 那Y_n就收敛到0 但是这个case是trivial的 可是如果只有某些eigenvalue的绝对值<1, 我又得不到什麽结果 最後回想到硕班ODE的一阶线性微分方程也是写成类似的形式 然後就能藉由eigenvalue去分类解会: (1) 发散到无穷 / 收敛到一点 / 在某个有界区域回绕 (2) 轨迹是椭圆 / 轨迹是值线 所以觉得原题目(ODE的离散化?)应该也会有相应的结果吧@@? 只是好几年没碰这些了, 只剩些印象... 再请板友们帮忙, 谢谢~ ----------------------------------------------------------------- <Update> 线代启示录网站有写关於x_n = q_n * s^n, q_n是多项式的情形 因此general的x_n应该不可能归类了 只是我工作遇到的x_n有compact support, 应该还是有机会(? --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.164.101.117 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Math/M.1651132986.A.52A.html
1F:推 LPH66 : Um, 数列的 compact support?04/28 18:31
2F:→ LPH66 : 你是想讲只有有限个 n 有非零 x_n?04/28 18:31
3F:→ LPH66 : (就是你中间的 " x_n = 0 for n>=M, some M ")04/28 18:32
4F:推 LPH66 : 是说啦, 如果只有有限项 x_n 非零, 要探讨 n→∞ 时04/28 18:36
5F:→ LPH66 : 的行为的话就只看 x_n 是 0 的「最後区域」就好04/28 18:37
6F:→ LPH66 : 那这其实就成了齐次线性递回, 一般的理论就能用了04/28 18:37
嗨L大, 我刚刚边骑机车边想时突然发现改变一下初始值就变成原本的case了XDDD 耍笨了 如果x_n=0 for all n>=M, 那考虑y_n, n>=M就像你说的是其次线性了@@ 另外请问一下齐次线性的解有什麽分类讨论吗? 就像我文中说的ODE会讨论解的有界/无界, 收敛/发散...这些行为 比如 y_n = a_1*y_(n-1) + a_2*y_(n-2) + a_3*y_(n-3) 我是想知道有没有可能存在a_1,a_2,a_3跟初始值y_0, y_1,y_2 使得当n→∞时, y_n趋近於周期数列(即来回震荡), 但是又不要y_n = (-1)^n这种极端的 比如y_n = -1, -1/2, 0, 1/2, 1, 1/2, 0, -1/2, -1.... 目前证出: 当只有a_1,a_2两个系数时, 不可能 (2x2矩阵的对角矩阵很容易得证) 而且一直造例子发现如果y_n震荡都会是(-1)^n这种一上一下, 不会有中间态的情形... 谢谢回答~ ※ 编辑: znmkhxrw (59.102.225.191 台湾), 04/28/2022 21:11:17
7F:推 LPH66 : 齐次线性递回去看它的特徵方程04/28 23:09
8F:→ LPH66 : 它的解会跟特徵方程的根有关04/28 23:10
9F:→ LPH66 : 在没有重根的状况会是以其根为底的指数之线性组合04/28 23:10
10F:→ LPH66 : 那如果根出现复根, 经过一些复数的推导之後04/28 23:15
11F:→ LPH66 : (尤拉公式与棣美弗) 会得到 cos 与 sin 表示的结果04/28 23:16
12F:→ LPH66 : 这就能解释你观察到某些状况的震荡性了 04/28 23:16
喔喔喔!!! 想起来了, 难怪我一直用实根做实验都不会有中间态的震荡 谢谢~~~ ※ 编辑: znmkhxrw (59.102.225.191 台湾), 04/28/2022 23:33:24 ※ 编辑: znmkhxrw (59.102.225.191 台湾), 04/28/2022 23:33:44
13F:推 Vulpix : 看到的第一印象是用生成函数,也就是z转换。因为那05/01 13:15
14F:→ Vulpix : 个convolution实在太抢眼了。实际做起来也挺快的。05/01 13:15
大学没修离散数学, 所以对生成函数这些名词以及递回方程做法仅有个大概方向 谢谢回覆~ ※ 编辑: znmkhxrw (114.137.102.213 台湾), 05/01/2022 16:18:03







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灯, 水草

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

TOP