| 查看: 2489 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请高手帮忙求解一个常微分方程的题目
已经有6人回复
- 帮忙求解一道微分方程的题目
已经有7人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】Python 中 求解微分方程 odeint 的用法
已经有11人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求教这种微分方程的详细求解过程
已经有14人回复
- 【求助】如何用matlab求解这个微分方程？【已解决】
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
美国肯塔基大学Dr. Sheng Tong博士后招聘（空间可控的基因组编辑与肿瘤免疫治疗）
+1/94
-
招收 哈工大深圳﹣大湾区大学联合培养博士生1名，2027年入学
+1/94
-
高校教师 科研教学岗
+1/79
-
【推免】】中科院国家级人才团队招收27级推免硕士生 （材料+机械）
+1/71
-
求助邢其毅第四版基础有机化学PDF版
+1/44
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
+1/32
-
Biofunctional Materials青年编委全球招募，邀您共建学术阵地！
+1/15
-
聚智生物材料，共赴学术新程｜Biofunctional Materials青年编委全球招募
+1/9
-
哈尔滨工业大学（深圳）吴晓军教授课题组招收2027级“申请-考核”博士研究生
+2/8
-
怀柔实验室与中国矿业大学（北京）氢能与燃料电池课题组招收2027年联培博士生
+1/6
-
中国石油大学（华东）生物质绿色化学转化和生物基功能高分子团队诚招27级推免生
+1/5
-
关于征集高能质子与氦离子辐照实验用户需求的公告
+1/4
-
【推免招生】招收材料、化工、环境、冶金类推免生
+1/4
-
浙江大学Kemal Celebi教授课题组招收2027年秋季博士生
+1/3
-
无人机方向 博士后与青年教师招聘公告
+1/3
-
北理工集成电路杰青团队 | 诚招科助理
+1/3
-
宁波诺丁汉大学招收2027年春季入学CO2电催化方向博士生
+1/3
-
浙江大学Said Mikki教授课题组招生2027秋博士生-电子科学与技术
+1/2
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
-
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题，Scopus期刊
+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
