how to solve nonlinear coupled dgl second order
14 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Christian Dieterich
am 29 Mär. 2018
Kommentiert: Christian Dieterich
am 30 Mär. 2018
Hello everyone,
Does anybody know how to solve the following two differential equotations using ODE45?
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/171652/image.png)
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
Akzeptierte Antwort
Abraham Boayue
am 30 Mär. 2018
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;
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/180696/image.png)
0 Kommentare
Weitere Antworten (2)
Siehe auch
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!