24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

watertxf

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 504.8
散金: 42
红花: 1
帖子: 46
在线: 16.5小时
虫号: 1289373
注册: 2011-05-08
性别: MM
专业: 同步辐射技术及其应用

[求助] matlab求解微分方程

要求解方程为 dy/dt+{1+16*[1-exp(-t)-y]}*y=1，要求画出y随t的变化曲线。为什么用dsolve得不出解呢？求帮忙，非常感谢！

y=dsolve('Dy+y+18*y-18*y^2-1=18*y*exp(-t)','t')
返回
Warning: Explicit solution could not be found.
> In dsolve at 194
In Untitled at 4

y =

[ empty sym ]

回复此楼

» 猜你喜欢

三无产品还有机会吗已经有4人回复
投稿返修后收到这样的回复，还有希望吗已经有7人回复
压汞仪和BET测气凝胶孔隙率已经有4人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复

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

1楼2013-11-27 18:23:30

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68121.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

9楼: Originally posted by watertxf at 2013-11-27 19:30:05
t的范围不是1~5，是0~5
为什么我改成tspan=0:0.01:5;后得出的结果成下边的图1了？很明显不对的，t的范围就不对.这是什么原因呢？谢谢！

D_V9`YP`W~TCPD[5EA(@AZT.jpg
...

CODE:

function solve_ode

clear all;clc

y0=0;

tspan=0:0.01:5;

[ts ys]=ode45(@odefun,tspan,y0);

figure(1),

plot(ts,ys,'-ob')

function dydt=odefun(t,y)

dydt=1-(1+16*(1-exp(-t)-y))*y;