【答案】应助回帖

» 猜你喜欢

基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有81人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有13人回复
散金币祈福已经有94人回复
青基已中已经有1人回复
散金币祈福已经有185人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复

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

Matlab求解偏微分方程组已经有22人回复
请问如何求解二元一阶微分方程组已经有24人回复
matlab数值求解边界条件微分方程组已经有7人回复
LABVIEW四阶龙格库塔法求解一阶微分方程组已经有12人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】这个微分方程组simulink求解到底错在哪里？已经有3人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】急请微分方程高手看下这个方程组！已经有14人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【分享】PDECOL-经典数值求解非线性偏微分方程组源程序已经有142人回复

hehe

1楼 2011-11-10 21:09:57

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★
ben_ladeng(金币+1): 谢谢参与 2011-11-11 00:27:25
_aqing(金币+1): 试过好多解法，未能凑效。包括数值解。求详细解法。 2011-11-14 16:13:10

matlab和mathematic求不出解析解，尝试数值解法。

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2楼2011-11-10 22:03:10

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

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

【答案】应助回帖

★
_aqing(金币+1): 愿闻其详。 2011-11-14 16:14:24
dbb627(金币+1): 感谢参与 2011-11-16 16:05:33

用编程迭代法求解

生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`

3楼2011-11-14 09:56:17

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

zhaozhh02

铜虫 (小有名气)

应助: 4 (幼儿园)
金币: 48.1
帖子: 156
在线: 32.7小时
虫号: 546980
注册: 2008-04-16
性别: GG
专业: 粒子物理学和场论

★
dbb627(金币+1): 感谢参与 2011-11-16 16:05:03

这个可以解，Matlab的bvp4c或者bvp5c都行的，有个专门的bvp4c的教程你搜下。

不是很麻烦的，但是要给边界条件和猜测解

4楼2011-11-14 16:39:26

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

zhaozhh02

铜虫 (小有名气)

应助: 4 (幼儿园)
金币: 48.1
帖子: 156
在线: 32.7小时
虫号: 546980
注册: 2008-04-16
性别: GG
专业: 粒子物理学和场论

忘了说了，只能数值求解。这个方程组是非线性的，极难有解析解。

5楼2011-11-14 16:41:31

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

引用回帖:

1楼: Originally posted by _aqing at 2011-11-10 21:09:57:
这个微分方程组，实在解不出。求大家给个意见或参考。能给出解的定重谢。

可以给出代码，但是你需要先给出些参数的值。M c和 b的值和初值条件

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6楼2011-11-14 17:03:25

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

★ ★ ★
cenwanglai(金币+3): 幸苦了~ 2011-11-17 20:35:57

我随便取了些数据，你参考下吧

CODE:

 M=200;c=2;b=2;

 odefun=@(t,x)[M*(-3/(8*pi*x(2)))*(17*c*x(1).^3/x(2).^2+3*b);M*(3/(16*pi*x(1)))*(11*c*x(1).^3/x(2).^2-b)];

[t,x]=ode45(odefun,[0 10],[2 3]);

plot(t,x)

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

7楼2011-11-15 09:49:17

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

zhaozhh02

铜虫 (小有名气)

应助: 4 (幼儿园)
金币: 48.1
帖子: 156
在线: 32.7小时
虫号: 546980
注册: 2008-04-16
性别: GG
专业: 粒子物理学和场论

对，如果是初值问题的ode45也可以解

8楼2011-11-16 14:45:12

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

_aqing

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 561.3
帖子: 377
在线: 115.7小时
虫号: 1151611
注册: 2010-11-20
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

5楼: Originally posted by zhaozhh02 at 2011-11-14 16:41:31:
忘了说了，只能数值求解。这个方程组是非线性的，极难有解析解。

给出参数大家帮忙算算
具体参数
b=34.5086
M=0.0048
c=1.6*10^8*(ln(8.16*L/r)-1)
初值
l=1.35*10^(-4)
r=110*10^(-7)

hehe

9楼2011-11-17 17:14:38

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

引用回帖:

9楼: Originally posted by _aqing at 2011-11-17 17:14:38:
给出参数大家帮忙算算
具体参数
b=34.5086
M=0.0048
c=1.6*10^8*(ln(8.16*L/r)-1)
初值
l=1.35*10^(-4)
r=110*10^(-7)

c=1.6*10^8*(ln(8.16*L/r)-1) c是个变量吗，随着L和r实时变吗？