搜狗做题网matlab中ode45括号的意思,导数?
来源：学生作业帮 编辑：搜狗做题网作业帮 分类：综合作业 时间：2024/08/11 15:46:02
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)=.
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)=.
首先介绍一下函数用法
[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;
希望你能看明白,祝你学习愉快!
[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;
希望你能看明白,祝你学习愉快!