24小时热门版块排行榜    

查看: 1301  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ggamis

至尊木虫 (职业作家)

[求助] 求教一个证明数值积分递推公式无条件稳定的问题 已有1人参与

求教一个证明数值积分递推公式无条件稳定的问题
上图是递推公式。其中h是时间步长,也就是要确定的值。

求教一个证明数值积分递推公式无条件稳定的问题-1
这里是矩阵A,B的表达式,其中M,C,K分别为质量刚度阻尼矩阵,质量阵是正定的,刚度和阻尼矩阵是半正定的

求教一个证明数值积分递推公式无条件稳定的问题-2
作者有一句话,说这个数值积分方法是无条件稳定的,还说相关的放大矩阵的特征值为1。

为什么我看不懂啊!!作者写的这么理所当然,应该是很明显的道理吧!!求各位高人指点呐!!

ps:我可以证明B-(h/2)A的范数大于B+(h/2)A的范数,不知道有没有用啊!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggamis

至尊木虫 (职业作家)

引用回帖:
2楼: Originally posted by wurongjun at 2016-03-11 09:56:29
作者是说这里提出的线性插值方法是无条件稳定的!
放大矩阵特征值的模为1,说明误差是不增长的!所以稳定!

或者说,作者说线性差值方法是无条件稳定的。不能代表说这个数值积分的算法是无条件稳定的吗?

那么所谓的放大矩阵是指的啥啊?
4楼2016-03-11 12:48:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
作者是说这里提出的线性插值方法是无条件稳定的!
放大矩阵特征值的模为1,说明误差是不增长的!所以稳定!

» 本帖已获得的红花(最新10朵)

善恶到头终有报,人间正道是沧桑.
2楼2016-03-11 09:56:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggamis

至尊木虫 (职业作家)

送红花一朵
引用回帖:
2楼: Originally posted by wurongjun at 2016-03-11 09:56:29
作者是说这里提出的线性插值方法是无条件稳定的!
放大矩阵特征值的模为1,说明误差是不增长的!所以稳定!

对对。我们对z(t)做了线性假设(也就是线性差值方法)

请问,放大矩阵是指的什么啊?就是这里不懂啊!

大神不吝赐教啊!!工科的矩阵论和数值分析学得不太好。
3楼2016-03-11 12:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

引用回帖:
4楼: Originally posted by ggamis at 2016-03-11 12:48:42
或者说,作者说线性差值方法是无条件稳定的。不能代表说这个数值积分的算法是无条件稳定的吗?

那么所谓的放大矩阵是指的啥啊?...

不是说的数值积分!
放大矩阵由24式决定,你好好看文章吧!
善恶到头终有报,人间正道是沧桑.
5楼2016-03-11 13:03:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见