Qurak

Nd solve delay differential equations comparison and contrast with odes

ట్యుటోరియల్ 8 దశలు, ఒకే సెషన్‌లో క్రమంగా నడుస్తాయి. ప్రతి దశా ఇంజిన్‌పై అమలు చేయబడింది, కింది అవుట్‌పుట్ అది ఇచ్చినదే.

దశ 1
Manipulate[
	Module[{sol = NDSolve[{x'[t] == x[t - 1](1 - x[t]), x[t /; t ≤ 0] == ϕ}, x, {t, -2, 2}]},
	Plot[Evaluate[x[t] /. First[sol]], {t, -2, 2}]], {ϕ, {Exp[t], Cos[t], 1 - t, 1 - Sin[t]}}]
అవుట్‌పుట్
Manipulate[Module[{sol = NDSolve[{Derivative[1][x][t] == x[t - 1]*(1 - x[t]), x[t /; t <= 0] == ϕ}, x, {t, -2, 2}]}, Plot[Evaluate[x[t] /. First[sol]], {t, -2, 2}]], {ϕ, {E^t, Cos[t], 1 - t, 1 - Sin[t]}}]
దశ 2
Manipulate[
	Module[{T = 50, sol, x, t}, sol = First[x /. NDSolve[{x'[t] == a x[t](1 - x[t - 1]), x[t /; t ≤ 0] == 0.1}, x, {t, 0, T}]];
	If[pp, ParametricPlot[{sol[t], sol[t - 1]}, {t, 1, T}, PlotRange -> {{0, 3}, {0, 3}}],
	Plot[sol[t], {t, 0, T}, PlotRange -> {{0, 50}, {0, 3}}]]], {{pp, False, "Plot in Phase Plane"}, {False, True}}, {{a, 1}, 0, 2}]
అవుట్‌పుట్
Manipulate[Module[{T = 50, sol, x, t}, sol = First[x /. NDSolve[{Derivative[1][x][t] == a*x[t]*(1 - x[t - 1]), x[t /; t <= 0] == 0.1}, x, {t, 0, T}]]; If[pp, ParametricPlot[{sol[t], sol[t - 1]}, {t, 1, T}, PlotRange -> {{0, 3}, {0, 3}}], Plot[sol[t], {t, 0, T}, PlotRange -> {{0, 50}, {0, 3}}]]], {{pp, False, Plot in Phase Plane}, {False, True}}, {{a, 1}, 0, 2}]
దశ 3
sol1 = First[NDSolve[{x'[t] == Sin[x[t - 20]], x[t /; t ≤ 0] == .0001}, x, {t, 0, 500}]];
sol2 = First[NDSolve[{x'[t] == Sin[x[t - 20]], x[t /; t ≤ 0] == .00011}, x, {t, 0, 500}]];

అవుట్‌పుట్ లేదు - ఈ దశ తదుపరి దానికి ఏదో సిద్ధం చేస్తుంది.

దశ 4
Plot[Evaluate[x[t] /. {sol1, sol2}], {t, 0, 500}]
అవుట్‌పుట్
-Graphics-
దశ 5
Grid[Table[sol = First[NDSolve[{x'[t] == Sin[x[t - τ]], x[t /; t ≤ 0] == .1}, x, {t, 100τ, 200τ}, MaxSteps -> Infinity]];
	{ParametricPlot[Evaluate[{x[t - 1], x[t]} /. sol], {t, 101τ, 200τ}].
	Plot[Evaluate[x[t] /. sol], {t, 100τ, 200τ}]}, {τ, 4.9, 5.1, .1}]]
అవుట్‌పుట్
Grid[{{-Graphics- . -Graphics-}, {-Graphics- . -Graphics-}, {-Graphics- . -Graphics-}}]
దశ 6
Block[{λ = 1 / 2, μ = -1, T = 25}, Plot[Evaluate[First[x[t] /. NDSolve[{x'[t] == λ x[t] + μ x[t - 1], x[t /; t ≤ 0] == 1 - t}, x, {t, 0, T}]]], {t, 0, T}, PlotRange -> All]]
అవుట్‌పుట్
-Graphics-
దశ 7
Block[{λ = -7 / 2, μ = 4, T = 25}, Plot[Evaluate[First[x[t] /. NDSolve[{x'[t] == λ x[t] + μ x[t - 1], x[t /; t ≤ 0] == 1 - t}, x, {t, 0, T}]]], {t, 0, T}, PlotRange -> All]]
అవుట్‌పుట్
-Graphics-
దశ 8
Manipulate[Module[{T = 25, x, t}, Plot[Evaluate[First[x[t] /. NDSolve[{x'[t] == λ x[t] + μ x[t - 1], x[t /; t ≤ 0] == 1 - t}, x, {t, 0, T}]]], {t, 0, T}, PlotRange -> All]], {λ, -5, 5}, {μ, -5, 5}]
అవుట్‌పుట్
Manipulate[Module[{T = 25, x, t}, Plot[Evaluate[First[x[t] /. NDSolve[{Derivative[1][x][t] == λ*x[t] + μ*x[t - 1], x[t /; t <= 0] == 1 - t}, x, {t, 0, T}]]], {t, 0, T}, PlotRange -> All]], {λ, -5, 5}, {μ, -5, 5}]

వాడిన ఫంక్షన్లు

సంబంధిత రెసిపీలు

అన్ని రెసిపీలు · ఫంక్షన్ రిఫరెన్స్ · దీన్ని MCP క్లయింట్ నుండి వాడండి