Как построить нелинейную систему дифференциальных уравнений в MATLAB?

Question

На самом деле, я задаю тот же вопрос снова

Dx = у

Dy = -k * уи^3 + 9,8 * соз (т)

inits = ('х (0) = 0, y (0) = 0 ')

Это дифференциальные уравнения, которые я хотел построить.

Во-первых, я попытался решить дифференциальное уравнение, а затем построить график.

Dsolve ('Dx = у', 'Dy = -k * у-х^3 + 9,8 * соз (т)', inits)

, как это, однако, не было явное решение этой системы.

теперь я застрял :(

, как вы можете построить эту систему, не решая уравнений?

Крис Тейлор решил мой вопрос, откладывая у против т и х против т графов.

Однако Я хотел бы знать, как построить у против й графы (это право? Я означающий граф с вертикальной осью у, а горизонтальная ось х)

Пожалуйста, помогите мне

Answer 1

Используя код от Криса Тейлора answer, как показано ниже:

k = 1; 
f = @(t,x) [x(2); -k * x(2) - x(1)^3 + 9.8 * cos(t)]; 
tspan = [0, 10]; 
xinit = [0, 0]; 
[t, y] = ode45(f, tspan, xinit); 

дает вам значения х и у как функции т от 0 до 10, хранящихся в первом и втором столбцах у соответственно. Таким образом, вы можете построить у (х), просто вызывая

plot(y(:,2),y(:,1));