带有积分项的微分方程数值解，高手进

乱世野鬼 · 发表于 2013-10-16 16:21:19

方程的形式和边界条件如图所示，其中A=8.04*10^6, B=-3.6*10^9, D=1.2*10^-8 m=4*10^-8 k=8.85*10^-12 L=10^-6 y是x的函数, 求x=[0,10^-6] ，时y的值，(计算步长可选10^-8)。求相应或者相似的Matlab程序，或者加我QQ492777972一起讨论，谢谢

TBE_Legend · 发表于 2013-10-16 16:21:20

TBE_Legend 发表于 2013-10-22 13:04
我顺便那你的问题测试下手上软件解微积分方程的能力。

放心，如有结果立马发给你。 ...

代码已经短消息发给你。

这是结果：

乱世野鬼 · 发表于 2013-10-17 13:28:12

方程是Dy''(x)，两个点点错位置了。。。

TBE_Legend · 发表于 2013-10-18 13:46:04

x=[0,10^-6]

y[0]=?

y[-L/2]=?

TBE_Legend · 发表于 2013-10-18 13:48:38

到底是 0 到 L 还是 -L/2 到 L/2

TBE_Legend · 发表于 2013-10-18 13:57:46

试了下含个积分，不是问题，问题是：你的边界应该不合适，求出来结果全部为0

乱世野鬼 · 发表于 2013-10-20 10:52:35

TBE_Legend 发表于 2013-10-18 13:57
试了下含个积分，不是问题，问题是：你的边界应该不合适，求出来结果全部为0 ...

能不能留个QQ，我请具体请教下，这个公式是有点问题

TBE_Legend · 发表于 2013-10-20 13:30:06

你把你正确的式子贴出来，写清楚，我给你求。

TBE_Legend · 发表于 2013-10-20 13:32:11

代码也会给你。

乱世野鬼 · 发表于 2013-10-20 15:02:09

TBE_Legend 发表于 2013-10-20 13:32
代码也会给你。

我做了个WORD把原始的方程和系数都写清楚，您看看能不能求

乱世野鬼 · 发表于 2013-10-20 15:02:50

TBE_Legend 发表于 2013-10-20 13:32
代码也会给你。

http://pan.baidu.com/s/1kTdpQ 网盘的地址

TBE_Legend · 发表于 2013-10-20 23:32:15

好的，明天抽空给你做做。

TBE_Legend · 发表于 2013-10-20 23:38:00

好的，明天抽空给你做做。

TBE_Legend · 发表于 2013-10-21 20:59:15

做了下，方程还是奇异，你试试无量钢化，系数太大或太小了。

再有就是你文档中的单位统一吗？

TBE_Legend · 发表于 2013-10-21 21:02:53

做了下，方程还是奇异，你试试无量钢化，系数太大或太小了。

再有就是你文档中的单位统一吗？

TBE_Legend · 发表于 2013-10-21 21:03:52

你改改我再给你试试

乱世野鬼 · 发表于 2013-10-22 11:22:09

本帖最后由乱世野鬼于 2013-10-22 11:38 编辑

TBE_Legend 发表于 2013-10-21 21:03
你改改我再给你试试

不好意思，是我自己没说清楚。我再确定一下

乱世野鬼 · 发表于 2013-10-22 11:55:43

TBE_Legend 发表于 2013-10-21 21:03
你改改我再给你试试

我重新又把系数确定了一下，都换算到了国际单位，word已传到网盘。麻烦你帮我再算下，如果取定是系数的问题，能不能把程序发我，我自己再调试。谢谢。http://pan.baidu.com/s/1rELjV

TBE_Legend · 发表于 2013-10-22 13:04:02

乱世野鬼发表于 2013-10-22 11:55
我重新又把系数确定了一下，都换算到了国际单位，word已传到网盘。麻烦你帮我再算下，如果取定是系数的问 ...

我顺便那你的问题测试下手上软件解微积分方程的能力。

放心，如有结果立马发给你。

乱世野鬼 · 发表于 2013-10-23 08:35:28

TBE_Legend 发表于 2013-10-16 16:21
代码已经短消息发给你。

这是结果：

十分感谢

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

带有积分项的微分方程数值解，高手进

最佳答案

本帖子中包含更多资源