Question

matlab中ode45括号的意思，导数？

function dy = vdp1000(t,y)
dy = zeros(2,1); % a column vector
dy(1) = y(2);
dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1);

[t,Y] = ode15s(@vdp1000,[0 3000],[2 0])

plot(t,Y(:,1),'-o')

请问ode45中dy(1) ，y（2）等 括号是导数的意思吗?是不是matlab自己就能识别括号加数字是导数的意思？还有就是要是要是我需要使用循环数组 ，下面的格式对吗？满意的话可加分

for i=1:5
dy(1)=zeros(2,1);

dy(i)(1)=y(i)(2);

dy(i)(2)=。。。。。。

Answer 1

首先介绍一下函数用法
[TOUT,YOUT] = ODE45(ODEFUN,TSPAN,Y0)
TOUT和YOUT是最后得到的自变量和对应的函数值
ODEFUN是你需要求解的微分方程或方程组
TSPAN是你的微分方程组的自变量的范围
Y0是方程组的初值

括号不是导数的意思，方程组是一个独立的函数，它的写法是首先定义导数，然后计算每个导数的值
例如微分方程组为
dA/dt=-k1*A+k2*B
dB/dt=k1*A-k2*B-k3*B

则函数的写法为
function dydt = react(t,y)
dydt = zeros(size(y));
k1 = 5; k2 = 2; k3 = 1;
A = y(1);
B = y(2);
C = y(3);
% 求解导数
dydt(1) = -k1*A + k2*B;
dydt(2) = k1*A - (k2+k3)*B;

希望你能看明白，祝你学习愉快！