版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

akang9

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 92.5
帖子: 23
在线: 15.9小时
虫号: 1259352
注册: 2011-04-08
专业: 纤维

[求助] 求高手解答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

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有49人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
三无产品还有机会吗已经有6人回复

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

Matlab求解二阶偏微分方程组，希望能给出相应的求解代码已经有10人回复
分数阶微分方程的求解（用差分法，matlab编程求解）已经有9人回复
如何用matlab求解矩阵系数的二阶微分方程已经有13人回复
求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？已经有5人回复
matlab求解偏微分方程组已经有3人回复
Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点已经有5人回复
高手进来用matlab求解常微分方程，十分感谢！已经有17人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
matlab解微分方程已经有10人回复
Matlab求解偏微分方程组已经有22人回复
偏微分方程数值解有基本理论和一些例子（化工类）MATLAB源程序已经有296人回复
matlab数值求解边界条件微分方程组已经有7人回复
高金求助matlab解微分方程组已经有12人回复
matlab解微分方程组已经有15人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
偏微分方程(PDEs)的MATLAB数值解法已经有262人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复

1楼 2013-02-02 20:54:56

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

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

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +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

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 92.5
帖子: 23
在线: 15.9小时
虫号: 1259352
注册: 2011-04-08
专业: 纤维

引用回帖:

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的回帖

akang9

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 92.5
帖子: 23
在线: 15.9小时
虫号: 1259352
注册: 2011-04-08
专业: 纤维

引用回帖:

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的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

内容已删除

赞一下

回复此楼

6楼2013-02-04 23:34:11

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

相关版块跳转我要订阅楼主 akang9 的主题更新

返回列表