how to solve nonlinear coupled dgl second order

17 次查看（过去 30 天）
Christian Dieterich
评论： Christian Dieterich 2018-3-30
Hello everyone,
Does anybody know how to solve the following two differential equotations using ODE45?
My Problem is, that i don't know how to rewrite the phi'' and x'' during the transformation in a system of dgl's first order.
Thanks

请先登录，再回答此问题。

采纳的回答

Abraham Boayue
Abraham Boayue 2018-3-30
Here is an example that I have obtained using ode45; the solution seems quite logical since we a dealing with sine waves.
clear variables
close all
N = 500;
L = 5; g = 9.81; mw = 2; mk = 5; Fan = 4; FR = 2;
q1 = -1/L;
q2 = -g/L;
F1 = -(mk*L)/(mw + mk);
F2 = -F1;
S = (Fan - FR)/(mw + mk);
F = @(t,y) [ y(2) ;
(1./(1-q1*F1*(cos(y(1)).^2))).*(0.5*q1*F2*sin(2*y(1)).*y(2).^2+...
q1*S*cos(y(1)) + q2*sin(y(1)))];
t0 = -2*pi;
tf = 2*pi;
tspan = t0:(tf-t0)/(N-1):tf;
ic = [0 0];
[t,y] = ode45(F, tspan, ic);
figure
plot(t,y(:,1),'-o')
hold on
plot(t,y(:,2),'-o')
a = title('\theta vs \theta_{prime}');
legend('\theta','\theta_{prime}');
set(a,'fontsize',14);
a = ylabel('y');
set(a,'Fontsize',14);
a = xlabel('t [-2\pi 2\pi]');
set(a,'Fontsize',14);
xlim([t0 tf])
grid
grid minor;

更多回答（2 个）

Abraham Boayue
Abraham Boayue 2018-3-30
Can you provide some values for the constants?
  1 个评论
Christian Dieterich
Christian Dieterich 2018-3-30
Hello Abraham Boayue,
First of all, Thanks for your help. Let's say that mw is 1000 kg and mw is 10 kg. L is 1 m and the initial condition of x' und x are 0, the initial condition of phi' is 0 and the initial condition of phi is 0.4 rad.

请先登录，再进行评论。

Abraham Boayue
Abraham Boayue 2018-3-30

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by