How do we plot differential Equation
显示 更早的评论
With theta1(0)=0.1, theta2(0)=0.3, theta3(0)=0.2 , I have to plot the solution in the range of [0.1] all in one graph.
I have no idea how to do this:
diff(theta1)=1+sin(theta2-theta1);
diff(theta2)=1+sin(theta1-theta2)+sin(theta3-theta2);
diff(theta3)=1+sin(theta2-theta3);
ezplot(
Please help
采纳的回答
Ameer Hamza
2020-11-4
编辑：Ameer Hamza
2020-11-4
Use ode45() for estimating a numerical solution
IC = [0.1, 0.3, 0.2];
tspan = [0 1];
[t, thetas] = ode45(@odeFun, tspan, IC);
plot(t, thetas);
legend({'\theta_1', '\theta_2', '\theta_3'}, 'FontSize', 16, 'Location', 'best');
function dthetas = odeFun(t, thetas)
theta1 = thetas(1);
theta2 = thetas(2);
theta3 = thetas(3);
dtheta1 = 1+sin(theta2-theta1);
dtheta2 = 1+sin(theta1-theta2)+sin(theta3-theta2);
dtheta3 = 1+sin(theta2-theta3);
dthetas = [dtheta1; dtheta2; dtheta3];
end
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
