|
用ode45来求解常微分方程数值解,函数定义中如果有参数,那么在调用的时候如何传递呢?我的程序如下,请问要如何修改才能实现我的目标呢?
定义函数文件vdp.m
function dx=vdp(t,x,k10,k12,k21)
dx=zeros(2,1);
dx(1)=-k10*x(1)-k12*x(1)+k21*x(2);
dx(2)=k12*x(1)-k21*x(2);
求解所定义的常微分方程
[t,x]=ode45(@vdp,[t1,tf],[x1,x2]);
其中k10,k12,k21在每次求解的时候,都会有新的参数值需要带入计算,而且每一次的求解,方程的初始值x1,x2也都会改变。
处于这样的目的,我的命令应该怎么写呢
|
|