duffing equation ODE45

48 次查看（过去 30 天）
Ramya Raman
Ramya Raman 2019-3-7
评论： ARPIT KUMAR 2020-11-16
I solved the duffing equation using ODE45. Below is the code.
It works fine and plots as expected. However, I noticed that the difference of 't' is not same. How can I have the 't' value equally spaced?
-----------------------------------------------------
tspan = [0 100];
x0=[1,0];
[t,y]=ode45(@duff1,tspan,x0)
figure(1)
%first value
% subplot(221)
plot(t,y(:,1)); %plot(t,y);
xlabel('Time');
ylabel('State');
title('k3=0,k=1,x0=1');
function result = duff1(t,y)
c = 0.1;
k3=0;
k=1;
result= [y(2); -c*y(2)-k*y(1)-k3*y(1).^3]
end
  1 个评论
ARPIT KUMAR
ARPIT KUMAR 2020-11-16
I did the coding for duffing eq but getting error , I am attaching the code below, kindly help me rectify the error.
clc % Time Response and Phase portrait of Duffing Oscillator global muglobal mu om1 om2 alpha mu=0.8; alpha=0.1 om1=2; om2=2*sqrt(2); [T,Y] = ode45(@ff22,[0 5000],[0.01 -0.3]); ns=length(Y); nm=floor(ns*0.9); figure(1) plot(Y(nm:ns,1),Y(nm:ns,2),'r','linewidth',2) set(gca,'FontSize',15) xlabel('\bfx','Fontsize',15) ylabel('\bf u','Fontsize',15) grid on function dy = ff22(t,y) global mu om1 om2 alpha dy = zeros(2,1); % a column vector dy(1) = y(2); dy(2) =cos(om1)+cos(om2) -y(1)-alpha*y(1)^3-2*mu*y(2); figure(2) subplot(2,1,1) plot(T(nm:ns),Y(nm:ns,1),'linewidth',2) grid on set(gca,'FontSize',15) xlabel('\bf Time','Fontsize',15) ylabel('\bfx','Fontsize',15) subplot(2,1,2) plot(T(nm:ns),Y(nm:ns,2),'linewidth',2) grid on set(gca,'FontSize',15) xlabel('\bf Time','Fontsize',15) ylabel('\bfx','Fontsize',15) end

请先登录，再进行评论。

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

采纳的回答

Star Strider
Star Strider 2019-3-7
Use:
tspan = linspace(0, 100, 50);
to create a ‘tspan’ vector with 50 elements between 0 and 100.

更多回答（1 个）

Daniel Dolan
Daniel Dolan 2019-3-7
Pass an array of tpan points. From the ode45 documentation:
tspan — Interval of integration
vector
Interval of integration, specified as a vector. At minimum, tspan must be a two element vector [t0 tf] specifying the initial and final times. To obtain solutions at specific times between t0 and tf, use a longer vector of the form [t0,t1,t2,...,tf]. The elements intspan must be all increasing or all decreasing.
The solver imposes the initial conditions given by y0 at the initial time tspan(1), then integrates from tspan(1) to tspan(end):
  • If tspan has two elements, [t0 tf], then the solver returns the solution evaluated at each internal integration step within the interval.
  • If tspan has more than two elements [t0,t1,t2,...,tf], then the solver returns the solution evaluated at the given points. However, the solver does not step precisely to each point specified in tspan. Instead, the solver uses its own internal steps to compute the solution, then evaluates the solution at the requested points in tspan. The solutions produced at the specified points are of the same order of accuracy as the solutions computed at each internal step.Specifying several intermediate points has little effect on the efficiency of computation, but for large systems it can affect memory management.
The values of tspan are used by the solver to calculate suitable values for InitialStep and MaxStep:
  • If tspan contains several intermediate points [t0,t1,t2,...,tf], then the specified points give an indication of the scale for the problem, which can affect the value of InitialStep used by the solver. Therefore, the solution obtained by the solver might be different depending on whether you specify tspan as a two-element vector or as a vector with intermediate points.
  • The initial and final values in tspan are used to calculate the maximum step size MaxStep. Therefore, changing the initial or final values in tspan could lead to the solver using a different step sequence, which might change the solution.
Example: [1 10]
Example: [1 3 5 7 9 10]

类别

在 Help Center 和 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