matlab ode45解线性微分方程时的参数函数文件:function xdot=FreeOcillation(t,x,dummy,zeta)xdot=[x(2);-2.0*zeta*x(2)-x(1)];命令窗口:zeta=[0.1 1.0 5.0];tspan=linspace(0,40,400);%生成0-40的四百个线性点for i=1:3[t,x]=ode45('F
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/19 20:35:01
matlab ode45解线性微分方程时的参数函数文件:function xdot=FreeOcillation(t,x,dummy,zeta)xdot=[x(2);-2.0*zeta*x(2)-x(1)];命令窗口:zeta=[0.1 1.0 5.0];tspan=linspace(0,40,400);%生成0-40的四百个线性点for i=1:3[t,x]=ode45('F
matlab ode45解线性微分方程时的参数
函数文件:
function xdot=FreeOcillation(t,x,dummy,zeta)
xdot=[x(2);-2.0*zeta*x(2)-x(1)];
命令窗口:
zeta=[0.1 1.0 5.0];
tspan=linspace(0,40,400);%生成0-40的四百个线性点
for i=1:3
[t,x]=ode45('FreeOcillation',tspan,[1 1],[],zeta(i));
plot(t,x(:,1));
hold on
end
函数文件里的变量dummy和命令窗口的空矩阵是什么意思?
matlab ode45解线性微分方程时的参数函数文件:function xdot=FreeOcillation(t,x,dummy,zeta)xdot=[x(2);-2.0*zeta*x(2)-x(1)];命令窗口:zeta=[0.1 1.0 5.0];tspan=linspace(0,40,400);%生成0-40的四百个线性点for i=1:3[t,x]=ode45('F
函数FreeOcillation(t,x,dummy,zeta)写错了,应该把那个dummy删掉.该参数在函数中并未使用,而且调用的时候也没有提供相应的数据.
命令窗口的空矩阵代表ode求解器的选项.各种求解器的调用格式都是
[T,Y] = solver(odefun,tspan,y0,options,p1,p2...)
因为你的函数需要传递附加的参数zeta,前面四个参数都必须提供.前三个参数分别是函数名、时间范围和初值,第四个是options,用以控制求解过程的一些选项(如误差控制等).把options设为空矩阵表示使用默认选项.