24小时热门版块排行榜    

查看: 3130  |  回复: 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一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:12 by vZfe6xYu34yj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:01 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:33 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:24 by vZfe6xYu34yj
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 QUjhNVAcOSff 2026-09-13 7/350 2026-09-14 21:21 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见