24小时热门版块排行榜    

查看: 3068  |  回复: 28

Yefee

新虫 (初入文坛)

引用回帖:
19楼: Originally posted by pippi6 at 2015-01-15 05:53:51
你这个矩阵第一行和最后一行写的不对,应该是Neumann边界条件,不能是内点方程...

请问我的思路对吗?
21楼2015-01-15 10:35:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
14楼: Originally posted by Yefee at 2015-01-14 22:38:59
这个x是什么呢?...

简单的 说 不论你有什么数值格式 原方程(忽略边值)等价于 离散该方程得到的就是的形式;这里T空间离散得到的向量; 从而对应的就是微分方程(无限维)的特征值的有限维近似
22楼2015-01-15 12:29:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yefee

新虫 (初入文坛)

引用回帖:
22楼: Originally posted by mathstudy at 2015-01-15 12:29:49
简单的 说 不论你有什么数值格式 原方程(忽略边值)等价于 \frac{d^2T}{dx^2}-b(x)T=\lambda T 离散该方程得到的就是A T=\lambda T的形式;这里T空间离散得到的向量; 从而\lambda对应的就是微分方程(无限维)的特征值 ...

恩,尝试做了一下真是这样的!
感谢!
23楼2015-01-15 13:52:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yefee

新虫 (初入文坛)

引用回帖:
22楼: Originally posted by mathstudy at 2015-01-15 12:29:49
简单的 说 不论你有什么数值格式 原方程(忽略边值)等价于 \frac{d^2T}{dx^2}-b(x)T=\lambda T 离散该方程得到的就是A T=\lambda T的形式;这里T空间离散得到的向量; 从而\lambda对应的就是微分方程(无限维)的特征值 ...

还有个问题,当n=1000与n=500的时候,模拟的特征函数振幅差别很大。
请问产生这种结果的原因是什么?自己这方面的知识实在捉急啊
左图n=1000,右n=500;b(x)=1;
偏微分方程的特征根问题
24楼2015-01-15 15:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by Yefee at 2015-01-14 22:07:13
可以讲的详细点吗?
我的思路是这样的
...

把A矩阵第一列和最后一列去掉,第一行第三列的1-〉2,倒数一行的倒数第三列的1改为2试一试
长期招收博士生,参见http://fsl-unsw.com
25楼2015-01-15 18:35:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
24楼: Originally posted by Yefee at 2015-01-15 15:38:29
还有个问题,当n=1000与n=500的时候,模拟的特征函数振幅差别很大。
请问产生这种结果的原因是什么?自己这方面的知识实在捉急啊
左图n=1000,右n=500;b(x)=1;
...

不太懂  楼主的意思,特征函数  模拟的不好还是?有可能的是差分格式的问题  矩阵太大
26楼2015-01-15 19:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yefee

新虫 (初入文坛)

引用回帖:
26楼: Originally posted by mathstudy at 2015-01-15 19:38:22
不太懂  楼主的意思,特征函数  模拟的不好还是?有可能的是差分格式的问题  矩阵太大...

对啊,前几个特征函数差别不应该太大啊
可是500与1000直接振幅差的很多,形状倒是一致的
27楼2015-01-15 21:27:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaolongsb

新虫 (初入文坛)

数学是很神奇的
28楼2015-01-15 22:46:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
27楼: Originally posted by Yefee at 2015-01-15 21:27:20
对啊,前几个特征函数差别不应该太大啊
可是500与1000直接振幅差的很多,形状倒是一致的...

特征函数可以相差一个常数
29楼2015-01-15 23:40:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Yefee 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 7/350 2026-08-25 03:24 by BZKMTicpDhFj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见