应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 应用数学与力学 » 求教一个证明数值积分递推公式无条件稳定的问题
61/1返回列表
查看: 1282  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ggamis

至尊木虫 (职业作家)

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

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

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

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

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

ps：我可以证明B-（h/2）A的范数大于B+（h/2）A的范数，不知道有没有用啊！！
回复此楼

» 猜你喜欢
1楼 2016-03-10 21:59:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

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

» 本帖已获得的红花（最新10朵）

ggamis
善恶到头终有报,人间正道是沧桑.
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的回帖

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的回帖

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的回帖

ggamis

至尊木虫 (职业作家)
引用回帖:
5楼: Originally posted by wurongjun at 2016-03-11 13:03:49
不是说的数值积分!
放大矩阵由24式决定,你好好看文章吧!...

谢谢。

文章我和实验室的同学看了好多遍了。

审稿人让我分析我提出的数值算法的稳定性分析。 我想通过证明放大矩阵的范数小于1来证明。

但是参考的文献中，没有关于算法稳定性的证明，所以很尴尬啊。

麻烦您能帮忙看看，方程（24）的放大矩阵的特征值的模怎么就是1了。系统矩阵A,B也给出来了，h是积分的步长，也就是要确定的值。
赞一下 回复此楼
6楼2016-03-11 13:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ggamis 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭