| 查看: 2590 | 回复: 0 | ||
[求助]
求助各位MATLAB微分方程互相迭代怎么做？
|
|
我要做两个微分方程的互相迭代，程序如下 h0=[0.5 0;0 -0.5]; h1=[0 0.5;0.5 0]; odefun1=@(t,psi) -i*(h0+cos(0.5*t)*h1)*psi;%求初始psi tspan1=[0:0.1:10]; psi0=[1;0]; [t,psi]=ode45(odefun1,tspan1,psi0); psi.' %因为我的版本是R2010，这个psi求出来理论上是列向量，但是2010求出来是行向量，之前已在论坛发帖，以后热心 %网友在R14版本测试为列向量。 for i1=1:1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% temp_1=[]; temp_1=h1*psi.'; for i0=1:size(temp_1,2) val_1=temp_1(:,i0); odefun2=@(t,chi) -i*(h0+h1*imag(dot(chi,val_1)))*chi;%这部分为第一个微分方程，chi为所求变量，因为要迭代所以我在开头使用了for循环， tspan2=[10:-0.1:0]; % psi为上面所求出来的psi chi0=[0,0;0,1]*psi(end,  .';[t,chi]=ode45(odefun2,tspan2,chi0); end chi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% temp_2=[]; temp_2=chi.'; for i0=1:size(chi,2) val_1=temp_2(:,i0); odefun3=@(t,psi) -i*(h0+h1*imag(dot(val_1,h1*psi)))*psi;%这部分是迭代中的第二个微分方程，所求变量为psi，把上一步的chi带入这一步 tspan3=[0:0.1:10]; %求出psi，再把这一步求得的psi带入上一步求出chi，在带入这一步求psi，如此迭代 psi0=[1;0]; [t,psi]=ode45(odefun3,tspan3,psi0); end psi %%%%%%%%%%%%%%%%%%%%%%% for i0=1:size(chi,1) u(i0)=imag(chi(i0, *h1*psi(i0,.');end end %%%%%%%%%%%%%%%%%% 问题：在最上面的部分，求得的psi是61个列向量，但是带入迭代的第一个方程中得到的解和理论上出入很大。研究后我发现在matlab求得的chi结果中并不是拿来一个psi值带入得到一个chi解，得到的只是psi中的最后一个解带入的结果。求各位指点，在这部分卡了好久了。另，迭代中的第一个微分方程由于给出的是终值而不是初值，所以我把时间设为[10:-0.1:0]，不知道对这种知道终值的微分方程这样是否欠妥，但是我用最上面的微分方程反演试验了一下，结果是正确的。 |
回复此楼» 猜你喜欢
0703求调剂383分
已经有9人回复
-
求调剂
已经有8人回复
-
化学调剂求助
已经有7人回复
-
一志愿哈尔滨工业大学085600英一数二337分求调剂
已经有12人回复
-
085600，320分求调剂
已经有16人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有12人回复
-
0703化学调剂325分
已经有10人回复
-
085600，321分求调剂
已经有12人回复
-
化学357分，考研调剂
已经有12人回复
-
一志愿南昌大学，085600，344分求调剂
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
请教matlab如何增加迭代次数
已经有11人回复
- matlab如何求解多元多次方程，包括绘图
已经有7人回复
- 偏微分方程Matlab求解
已经有7人回复
- 用matlab 求解参数方程组
已经有8人回复
- 求助：matlab求解一个方程组，其中一个方程有三角函数
已经有10人回复
- 请教大家微分方程与动力系统大家都在做什么
已经有7人回复
- Matlab求解多元多次方程组问题
已经有7人回复
- matlab中利用LM算法解非线性方程组的问题
已经有5人回复
- 关于matlab求解非线性指数方程组出现问题
已经有3人回复
- 分数阶微分方程的求解（用差分法，matlab编程求解）
已经有9人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- Matlab求解偏微分方程组
已经有22人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
- 【求助】解常微分方程组
已经有7人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
