应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 求高手解答Matlab求解微分方程的初始值和边界值问题
51/1返回列表
查看: 2831  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

akang9

铜虫 (初入文坛)

[求助] 求高手解答Matlab求解微分方程的初始值和边界值问题

为了一项工程模拟，有四个微分方程，有三个已知初始值，另一个只有终值，怎么编程求解？我也试过对未知初始值的方程假设一个初始值，用ode45求解验证计算的终值与实际值是否吻合的方法，换了很多假设值，都得不到合理的结果，甚至曲线趋势与实际完全相反。我该怎么办？

比如：
f(1)=-(y(1)-20)*3.8*pi+5*y(4);
f(2)=0.1*y(2)/(0.003*100);
f(3)=y(2)/(0.003*y(1));
f(4)=0.55*(1-y(4))/(y(3)/60)*exp(-4*log(2)*(y(1));

已知初值和终值为
y(1)0=220
y(2)0=14
y(3)0未知，y(3)终值为1000
y(4)0=0
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-02-02 20:54:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

akang9

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by akang9 at 2013-02-04 21:49:13
dingd您好！您给我的code直接能够在命令窗口运行吗？

直接运行的话，跳出这样的错误信息：
??? Undefined function or method 'Variable' for input arguments of type 'char'.

我的版本是R2010b。...

dingd您给我的代码是在1stOpt中的代码吧。
但是下载了1.5Pro运行了，报错信息为

x=[0中变量定义有误！
3]中变量定义有误！
y3=[中变量定义有误！
1000]中变量定义有误！
赞一下 回复此楼
5楼2013-02-04 22:16:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

dingd

铁杆木虫 (职业作家)
ODEFunction y1'=-(y1-20)*3.8*pi+5*y4;
            y2'=0.1*y2/(0.003*100);
            y3'=y2/(0.003*y1);
            y4'=0.55*(1-y4)/(y3/60)*exp(-4*ln(2)*y1);
[/code]

边值估算:
        y3(x=0): 11.2471932817486

t1.jpg
赞一下 回复此楼
2楼2013-02-03 09:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
两点边值问题，Matlab有个BVP4C专门求解这种问题，有点麻烦，不如1stOpt来的直观：
CODE:
Variable  x=[0,3],y1=220,y2=12,y3=[,1000],y4=0;
Plot y1,y3[y2];
ODEFunction y1'=-(y1-20)*3.8*pi+5*y4;
            y2'=0.1*y2/(0.003*100);
            y3'=y2/(0.003*y1);
            y4'=0.55*(1-y4)/(y3/60)*exp(-4*ln(2)*y1);
赞一下 回复此楼
3楼2013-02-03 09:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

akang9

铜虫 (初入文坛)
引用回帖:
3楼: Originally posted by dingd at 2013-02-03 09:20:28
两点边值问题，Matlab有个BVP4C专门求解这种问题，有点麻烦，不如1stOpt来的直观：

Variable  x=,y1=220,y2=12,y3=,y4=0;
Plot y1,y3;
ODEFunction y1'=-(y1-20)*3.8*pi+5*y4;
            y2'=0.1*y2/(0.003 ...

dingd您好！您给我的code直接能够在命令窗口运行吗？

直接运行的话，跳出这样的错误信息：
??? Undefined function or method 'Variable' for input arguments of type 'char'.

我的版本是R2010b。
赞一下 回复此楼
4楼2013-02-04 21:49:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭