| 查看: 2240 | 回复: 4 | |||
[交流]
【求助完毕】用ode逆向求解微分方程问题
|
|
在Matlab中用ode45函数求解微分方程,已知终端时刻的值y(tf),要反向积分求初值y(0),请问高手应该注意什么?方程中含有随时间变化的参数。 例如: tspan=[5 0]; %%%%% (时间需 反向) dy1dt=-f*y(1)+g*y(2); %%%%%%% f 和g都是随时间变化的函数 dy2dt=f^2*y(1); %%% f 和 g 可以 数值 求得,需要实时 代入微分方程中求解 [ Last edited by cenwanglai on 2011-4-3 at 11:31 ] |
» 猜你喜欢
材料复试求调剂
已经有8人回复
一志愿0703化学招61最终排名62化学求调剂
已经有17人回复
考研调剂
已经有5人回复
085600材料与化工301分求调剂院校
已经有33人回复
284求调剂
已经有9人回复
269求调剂
已经有5人回复
285求调剂
已经有6人回复
283求调剂
已经有17人回复
292求调剂
已经有8人回复
070300化学 求调剂
已经有11人回复
» 本主题相关价值贴推荐,对您同样有帮助:
请高手帮忙求解一个常微分方程的题目
已经有6人回复
帮忙求解一道微分方程的题目
已经有7人回复
【求助】matlab 求解微分方程中的未知参数
已经有20人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
【求助】MATLAB 有限差分法(FDM)求解偏微分方程
已经有22人回复
【求助】积分微分方程matlab求解
已经有6人回复
【求助】Python 中 求解微分方程 odeint 的用法
已经有11人回复
【求助】matlab怎么求解偏微分方程组啊,先谢谢了
已经有13人回复
【求助】求教这种微分方程的详细求解过程
已经有14人回复
【求助】如何用matlab求解这个微分方程?【已解决】
已经有3人回复
» 抢金币啦!回帖就可以得到:
★★爆★★0860生物与医药调剂考生看过来,河北大学化学与材料科学学院
+1/213
2026秋季全奖港博招生:硬件安全 / 芯片设计 / AI算力方向博士招生
+1/85
西安电子科技大学通信工程学院招收 “非全日制 调剂生“
+2/64
延安大学国家级高层次人才团队招收材料、化学、化工相关领域调剂研究生数名
+1/41
福建理工大学招收环境、化工、能源动力、新能源科学与工程、材料、机械等学科背景
+1/37
【硕士调剂】【分子医学 10医类 最后机会】2026年青岛大学硕士调剂补招
+1/21
可调声波驱动仪
+1/20
上海应用技术大学姚子建课题组招生,点击就送研究生学历
+1/20
湖州师范大学信息工程学院2026年硕士研究生招生调剂公告
+1/16
福建农林大学资环学院李香真-姚敏杰教授课题组招收硕士调剂生
+1/7
全国重点实验室调剂
+1/5
深圳大学谢和平院士团队招收2026级博士生(AI for Science/计算化学方向,第二轮)
+1/5
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/5
大连工业杰青、长江团队-生物质材料方向招收2026级博士生
+1/4
福建理工大学材料学院招收专业代码08开头学术型硕士研究生
+1/4
西南石油大学招收2026年“申请-考核“博士研究生 方向:不对称催化氢化
+1/3
吉林工程技术师范学院 交叉学科研究院(08/07) 招收调剂
+1/3
【#上海调剂急录#】能接受985联合培养的速来!带你发一区文章!
+1/2
大连工业大学食品交叉科学研究院招收调剂研究生
+1/2
基于质谱从头测序与蛋白语言模型的肿瘤新抗原识别新技术研发 (杭州)
+1/2
在里面不是高手如云吗?就急啊!!
|
用ode(例如 ode45 ; ode113 )求解微分方程,如果已知终端条件,进行反向积分时,除了把积分区间设置之外,还需要注意什么。今天编写了一个程序怎么结果都会出错,请高手指点。 因为中间参数随时间变化,不知道问题出在哪儿: ft = linspace(5,0,25); % Generate t for f f = ft.^2 - ft - 3; % Generate f(t) gt = linspace(6,1,25); % Generate t for g g = 3*sin(gt-0.25); % Generate g(t) Tspan = [5 1]; % Solve from t=5 to t=1 % 计算范围修改 IC = -0.1518; % y(t=5) --终端y值 [T Y] = ode45(@(t,y) myode(t,y,ft,f,gt,g),Tspan,IC); % Solve ODE figure(1) plot(T,Y) --------------------------------------------- function dydt = myode(t,y,ft,f,gt,g) %%%%% 中间参数变化 f =interp1(ft,f,t); % Interpolate the data set (ft,f) at time t -----------这儿需要修改吗??如何变? g =interp1(gt,g,t); % Interpolate the data set (gt,g) at time t-----------这儿需要修改吗??如何变? dydt = -f.*y + g; % Evalute ODE at time t----------原微分方程 应该无需改动吧?? %%%%%%%%%% 运算结果错误,量级都差太远,请高手快点救急啊!! ---------------------------------------------------------------- 又试验了一下,貌似如果是常系数、或者方程简单就容易求解,应该改变积分区间和初值就没问题。 只是现在必选考虑变参数情况,请高手帮忙啊…… 分析上面的主要问题出在哪儿?? |
2楼2010-11-19 16:31:45
3楼2010-11-19 23:18:07
4楼2010-11-20 11:24:25
5楼2010-11-23 17:02:39













回复此楼