24小时热门版块排行榜    

Znn3bq.jpeg
查看: 705  |  回复: 13
当前主题已经存档。

yy3780

木虫 (著名写手)

[交流] 【求助】方程组本征值问题

方程组可以化为182×182阶矩阵。其中w,s,delta是已知的。w,s,delta,g是已知的。
l,m,n可以从1~9。
{c0,d0,c1,c2,...,c9,d1,d2,...,d9,c11,c12,...........,c99,d11,d12,.......,d99}
可以将方程组化为矩阵,编程序求方程组本征值,得到能量本征值。
自己的困难是怎样表示矩阵,将矩阵的元素的值,其行列都要表示出来。。麻烦大家帮我一下。
如果阶数更大的,这样的方程组,有简单的方法吗?
谢谢大家帮忙啊。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yy3780(金币+5,VIP+0): 8-6 09:48
yy3780(金币+5,VIP+0): 8-6 14:28
182×182,这么大,不知道能不能算得来。
怎样表示矩阵?用二维数组A(i,j)表示所有矩阵元,当然要先定义数组为182×182的了。
求矩阵本征值的话,可以直接调用内部程序,当然如果是用CVF,安装时再选择库文件。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
2楼2009-08-06 07:04:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yy3780

木虫 (著名写手)

要把矩阵元素表示出来感觉比较麻烦,还要确定行和列
3楼2009-08-06 09:49:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yy3780(金币+40,VIP+0):嗯,谢谢argo,只有这样子了,呵呵 8-6 14:27
呵呵,不写出矩阵元就想解本征值?个人认为好象不太可能吧。
如果矩阵元有某些规律,可以不必一个一个写,否则,好象没有什么好办法吧。
计算机只是帮用户做些体力劳动而已,个人理解。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
4楼2009-08-06 10:48:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水


小木虫(金币+0.5):给个红包,谢谢回帖交流
不是吧,这样回答就给我50金币,太多了吧。实在不好意思呀。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
5楼2009-08-06 15:26:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yy3780

木虫 (著名写手)

没有其它好办法了吧,就自己做些苦力活,一个一个写出矩阵元了
谢谢指教argo呀
金币我自个再挣,哈哈
6楼2009-08-06 15:39:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
spc08(金币+1,VIP+0):感谢积极参与,欢迎常来 8-6 17:57
建议你写程序之前,认真看一下各个矩阵元之间有没有什么规律,直觉应该有很多是为0的,其他有些元素应该可以写成一个或几个通式,如果确实如此,就可以通过循环语句让计算机半自动完成。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
7楼2009-08-06 17:04:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yy3780

木虫 (著名写手)

嗯,谢谢您的建议。
能否再请教您一下,好吗?
这是一个实对称矩阵。上三角元素和下三角元素是一样的。再写矩阵元素时,可以只写出矩阵的上三角元素,再加上主对角线上的元素,然后解本征值吗?fortran里有现成的程序可以求这样只写出上面元素的矩阵本征值吗?
谢谢指教呀

[ Last edited by yy3780 on 2009-8-6 at 18:20 ]
8楼2009-08-06 18:17:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yy3780(金币+10,VIP+0):谢谢snoopyzhao,我看看下 8-6 18:32
spc08(金币+1,VIP+0):感谢交流,欢迎常来 8-7 09:07
看一下 blas, linpack 等标准库中有没有,感觉应该是有的。
http://www.netlib.org/blas/
http://www.netlib.org/linpack/
9楼2009-08-06 18:23:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水


小木虫(金币+0.5):给个红包,谢谢回帖交流
应该是要全部写出来吧,既然知道上三角元素与下三角元素是一样的,那为什么不直接写出来?这个完全可以由程序自动完成的呀。
如果我没有理解错的话,应该就是A(i,j)=A(j,i)吧。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
10楼2009-08-06 22:31:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yy3780 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 6/300 2026-05-19 17:23 by kopsia
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 3/150 2026-05-19 13:56 by xtlilibin
[教师之家] 上海大学实验技术岗位非升即走 +10 嘻嘻哈哈乐呵呵 2026-05-15 11/550 2026-05-19 10:03 by 嘻嘻哈哈乐呵呵
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 同样的基金本子,换个专家直接从C变A! +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见