24小时热门版块排行榜    

查看: 2983  |  回复: 17

musejianglin

银虫 (小有名气)

引用回帖:
6楼: Originally posted by 月只蓝 at 2013-05-15 09:24:36
你的方程少个边界条件吧,是不是少个"dy/dx(x=某数值)= ”的边界条件? 不然按你目前的边界条件,都是x=0处的初始值,这个微分方程已经不是边值问题了,而是初值问题了。在现有条件下,解该问题用ode可解,程序 ...

谢谢,高手。能再帮我看看么?能不能弄出来?

附件.png

11楼2013-05-16 21:51:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

musejianglin

银虫 (小有名气)

引用回帖:
8楼: Originally posted by somomo91 at 2013-05-15 18:28:25
不必用1stopt, MATLAB 完全可以做到
这里是程序 :
===================================
% main program
close all;  clear all;  clc
global a
a = 0.1;
t0    = 0;   tend=5000;
tspan = t0:0.5:tend;
...

谢谢你了,那能不能再帮我弄弄这个。谢谢了

附件.png

12楼2013-05-16 21:52:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fegg7502: 金币+2, 应助指数+1, 3ks 2013-05-18 14:32:18
musejianglin: 金币+9, 有帮助 2013-05-20 09:04:21
引用回帖:
11楼: Originally posted by musejianglin at 2013-05-16 21:51:04
谢谢,高手。能再帮我看看么?能不能弄出来?

附件.png
...

边界条件还是有问题吧,边界条件顾名思义是将解确定在一定的区间之内,你的边界条件只有右端,没有左端,怎么确定解呢,这样子说不定有无穷多个解。
能不能给出比如"y(r=b)=" 或者"dy/dr(r=b)="(b<a)呢?
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
13楼2013-05-17 08:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

musejianglin

银虫 (小有名气)

引用回帖:
13楼: Originally posted by 月只蓝 at 2013-05-17 08:58:08
边界条件还是有问题吧,边界条件顾名思义是将解确定在一定的区间之内,你的边界条件只有右端,没有左端,怎么确定解呢,这样子说不定有无穷多个解。
能不能给出比如"y(r=b)=" 或者"dy/dr(r=b)=&qu ...

哦,这是个球,前面的条件是在r=a处即球的表面处,而在0<r<a范围内都是0.这样可以了么?
14楼2013-05-17 10:34:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)


fegg7502: 金币+1, 3ks 2013-05-18 14:32:29
引用回帖:
14楼: Originally posted by musejianglin at 2013-05-17 10:34:04
哦,这是个球,前面的条件是在r=a处即球的表面处,而在0<r<a范围内都是0.这样可以了么?...

是y(r=0)=0,还是dy/dr(r=0)=0呢?
不可能说0<r<a范围内都是0,如果这样的话,[0 a]这个区间内原函数的数值全部确定了,都不需要解方程了。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
15楼2013-05-17 10:50:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
14楼: Originally posted by musejianglin at 2013-05-17 10:34:04
哦,这是个球,前面的条件是在r=a处即球的表面处,而在0<r<a范围内都是0.这样可以了么?...

你的原始模型的方程一定要建对,不然解出来的结果都是没有意义的。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
16楼2013-05-17 10:51:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

musejianglin

银虫 (小有名气)

引用回帖:
15楼: Originally posted by 月只蓝 at 2013-05-17 10:50:06
是y(r=0)=0,还是dy/dr(r=0)=0呢?
不可能说0<r<a范围内都是0,如果这样的话,这个区间内原函数的数值全部确定了,都不需要解方程了。...

是 y(r=0)=0.
17楼2013-05-20 08:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

musejianglin

银虫 (小有名气)

引用回帖:
16楼: Originally posted by 月只蓝 at 2013-05-17 10:51:16
你的原始模型的方程一定要建对,不然解出来的结果都是没有意义的。...

是 y(r=0)=0
18楼2013-05-20 08:58:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 musejianglin 的主题更新
信息提示
请填处理意见