0

我想编写一个快速的 matlab 代码来解决参数与空间和时间相关的 ODE 系统。有没有办法以矢量方式解决它?以下是代码的一部分:

 function dM=testfun(T,M,B1)

   Mx=M(1);
   My=M(2);
  Mz=M(3);


dM=[My*B1(3)-Mz*B1(2);Mz*B1(1)-Mx*B1(3);Mx*B1(2)-My*B1(1)]+...
    [-Mx/T2;-My/T2;(1-Mz)/T1];

here B1 is space-time dependent variable 

主程序:

   clc;clear all;
   tmin=0;
   dt=rand(5,1);
   tsteps=10;
   B1x=rand(5,1);% time-dependent
   jx=length(x);
   z=-10:10;
   jz=length(z);
   M0=repmat([0; 0; 1],1,jz);
  for p=1:jz
       B1=[B1x;zeros(1,jx);10*z(jz)*ones(1,jx)]; % 3rd term is space dependent
     for pp=1:jx
       tspan=linspace(tmin,tmin+dt,tsteps);
  [t,M(:,:,pp,p)]=ode45 (@(t,M) testfun(t,M,B1(:,jj)),tspan,M0(:,p));
      tmin=tmin+dt;
    end
end

有没有办法对代码进行矢量化?时间依赖性存在一个问题——matlab 有一个内置函数 interp1 用于处理时间依赖性 ODE,但它需要太多时间。在这里查看 如何求解具有时间相关参数的常微分方程系统 我的代码将基于这个基本方程进行开发,重要的是要节省一些时间。有什么出路吗?先感谢您。

代码的向量化可以在下面 的 Matlab 代码向量化中找到,涉及每次迭代时的 ODE 求解器

下面的帖子也很有帮助。 如何求解具有时间相关参数的常微分方程组

4

0 回答 0