24小时热门版块排行榜    

查看: 1478  |  回复: 10

xzz1007

金虫 (小有名气)

[求助] 求助:解矩阵的有关问题(对上一贴中的问题进一步描述了一下)

大侠你好,
我是学化学的,对于数学的学习,仅限于大学学的那点,至今已经近十年了,所以,很多问题都不明白,莫见笑!!
已经整理成word,于附件中![ 来自小组 减肥吧 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 三.doc
  • 2013-07-04 16:52:46, 37.5 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

阳光总在风雨后,再坚持坚持吧.....
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★ ★ ★ ★
xzz1007: 回帖置顶 2013-07-05 10:26:16
xzz1007: 金币+5, ★★★★★最佳答案 2013-07-05 10:26:22
见图片
求助:解矩阵的有关问题(对上一贴中的问题进一步描述了一下)
快照2.jpg

3楼2013-07-04 20:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
问题 (1) q1, q2 ,… q3N-6 是什么,
         (2)  R1, R2 ,… R3N 是什么
         (3)   G1, G2 ,… G3N 是什么 Gibs 自由能?
         (4)  V 是什么的函数 V =V (q,R)? 还是, q 是一组独立变量, R是另一组独立变量; 用了 q 做自变量,R也就是q的函数,或用了R 做自变量,q也就是R的函数?
我的化学还不如你的数学,可能。
2楼2013-07-04 18:39:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzz1007

金虫 (小有名气)

引用回帖:
2楼: Originally posted by pippi6 at 2013-07-04 18:39:17
问题 (1) q1, q2 ,… q3N-6 是什么,
         (2)  R1, R2 ,… R3N 是什么
         (3)   G1, G2 ,… G3N 是什么 Gibs 自由能?
         (4)  V 是什么的函数 V =V (q,R)? 还是, q 是一组独立变量, R是另 ...

N代表的是分子中的原子个数

假设分子为非线性,那么它的自由度也就是3N-6个

R是分子的笛卡尔坐标,也就是直角坐标,一个原子有一组(x,y,z),也就有三个R,那么N个原子也就有N组(x,y,z),因此也就3N个R

q是内坐标,也就是使用分子中的键长(N-1个)、键角(N-2个)、以及二面角(N-3个),来描述一个分子的,一共有3N-6个q

G是笛卡尔坐标系下的梯度

g是内坐标下的梯度,即所求

V是势能,V可以是V =V (q),也可以是V =V (R),现在就是想找到R与q的关系,好做微分,这个是我个人的想法!
阳光总在风雨后,再坚持坚持吧.....
4楼2013-07-05 09:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzz1007

金虫 (小有名气)

引用回帖:
3楼: Originally posted by pippi6 at 2013-07-04 20:00:16
见图片

快照2.jpg

非常感谢您的帮助!!
基本上能看懂!
u是个什么?具体应该是什么样的呢?它在里面的作用是什么呢?
式9中,这么多矩阵相乘,我怎么有点晕呢?
阳光总在风雨后,再坚持坚持吧.....
5楼2013-07-05 09:43:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
5楼: Originally posted by xzz1007 at 2013-07-05 09:43:45
非常感谢您的帮助!!
基本上能看懂!
u是个什么?具体应该是什么样的呢?它在里面的作用是什么呢?
式9中,这么多矩阵相乘,我怎么有点晕呢?...

你word文档里的

u is a diagonal matrix with the reciprocals of the atomic masses on the diagonal.

在数学里,就是权重方阵。在每一个变量上选择权重,一种重要性的量化表示。在你这里的问题,原文献建议用原子质量倒数。为什么这么取,我就不清楚了。你也许更有感觉。
6楼2013-07-05 10:11:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzz1007

金虫 (小有名气)

引用回帖:
6楼: Originally posted by pippi6 at 2013-07-05 10:11:47
你word文档里的

u is a diagonal matrix with the reciprocals of the atomic masses on the diagonal.

在数学里,就是权重方阵。在每一个变量上选择权重,一种重要性的量化表示。在你这里的问题,原文献建 ...

那wilson B矩阵中的微分我该怎么求呢?
u是对角非零的3N*3N矩阵,意思是非对角全是0,对角为原子质量倒数?
阳光总在风雨后,再坚持坚持吧.....
7楼2013-07-05 10:25:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzz1007

金虫 (小有名气)

引用回帖:
6楼: Originally posted by pippi6 at 2013-07-05 10:11:47
你word文档里的

u is a diagonal matrix with the reciprocals of the atomic masses on the diagonal.

在数学里,就是权重方阵。在每一个变量上选择权重,一种重要性的量化表示。在你这里的问题,原文献建 ...

我还有个问题的,我为什么不能直接用公式6,直接求g呢?
阳光总在风雨后,再坚持坚持吧.....
8楼2013-07-05 14:37:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
8楼: Originally posted by xzz1007 at 2013-07-05 14:37:52
我还有个问题的,我为什么不能直接用公式6,直接求g呢?...

因为超定。 如果你有 10个方程,5个变量,怎么解啊?方程数多于未知量数,叫超定(overdetermined)。方程数多少于未知量数,叫不定(underdetermined)。因为超定,所以用最小二乘法。
9楼2013-07-05 14:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzz1007

金虫 (小有名气)

引用回帖:
9楼: Originally posted by pippi6 at 2013-07-05 14:59:53
因为超定。 如果你有 10个方程,5个变量,怎么解啊?方程数多于未知量数,叫超定(overdetermined)。方程数多少于未知量数,叫不定(underdetermined)。因为超定,所以用最小二乘法。...

我怎么觉得我是3N个方程,3N-6个未知量呢??
大侠能否,推荐下,我目前应该看点什么书,或者学点什么呢?
现在的问题,还是B矩阵,和u矩阵,怎么办呀?
阳光总在风雨后,再坚持坚持吧.....
10楼2013-07-05 15:18:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xzz1007 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 18/900 2026-08-25 14:50 by Weed114114
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 +4 ziyangfang 2026-08-25 4/200 2026-08-25 14:41 by ZJTJZ
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +6 yuleib84 2026-08-25 7/350 2026-08-25 14:40 by 六两废铜
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 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
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见