Описание тега lorenz-system

Система Лоренца - это система обыкновенных дифференциальных уравнений (уравнения Лоренца, обратите внимание, что это не Лоренц), впервые изученная Эдвардом Лоренцем.
2 ответа

3D-печать аттрактора Лоренца с использованием MATLAB

Я использую MATLAB для построения аттрактора Лоренца, и мне было интересно, как я могу экспортировать координаты XYZ в файл для 3D-печати! У меня проблемы с переходом от координат XYZ к поверхности (я должен даже сделать это?). Кроме того, даже дума…
21 июл '16 в 01:56
1 ответ

Gnuplot - график с непрерывно меняющимся цветом линии - график с переходом цвета / переходом

Я использую GNU plot (gnuplot) для построения некоторых выходных данных от интеграции странного аттрактора ( см. Википедию) - и, поскольку это странный аттрактор, линия "вращается вокруг себя в трехмерном пространстве". Это затрудняет понимание того…
0 ответов

Python Nolds: как получить правильное значение для системы Лоренца

В настоящее время я использую следующий код для генерации серии Лоренца def generate(x, stop=10000, s=10, b=8/3, r=28): def lor(v): return np.array([s * (v[1] - v[0]), v[0] * (r - v[2]) - v[1], v[0] * v[1] - b * v[2]]) ret = [] step = 0.1 xtemp = x.…
05 окт '18 в 14:45
1 ответ

Каков наиболее описательный способ построения системы Лоренца?

Я (численно) решаю систему Лоренца, используя разные методы. Я строю это с использованием matplotlib, но я хотел бы найти способ лучше различать точки. Например: Давайте предположим, что точки, которые будут нанесены на график, хранятся в массиве a …
01 ноя '16 в 10:51
1 ответ

Система Лоренца в MATLAB- создание симуляторов и фильмов

Я пытаюсь смоделировать траектории в системе Лоренца в MATLAB, в настоящее время использую следующий код - clear all clf; clc; % Solution [t1,x1] = ode45('g',[0 30],[0;2;0]); [t2,x2] = ode45('g2',[0 30],[0;2.001;0]); [C,h] = size(x2); ang = 0; for j…
19 мар '16 в 18:50
1 ответ

Интеграция хаотических систем

У меня есть 3 уравнения следующим образом u'= -a*(u-v) v'= c*u-v-u*w w'= -b*w+u*v a=5.0 b=0.9 and c=8.2 Я пытаюсь использовать scipy.integrate.odeint решить от t = 0 до t=10 Мои начальные условия: u(0)=0, v(0)=1.0 и w(0)=2.0 Я не могу много полезных…
04 дек '15 в 15:12
3 ответа

Константы Рунге-Кутты, расходящиеся для системы Лоренца?

Я пытаюсь решить систему Лоренца, используя метод Рунге Кутты 4-го порядка, где dx/dt=a*(y-x) dy/dt=x(b-z)-y dx/dt=x*y-c*z Поскольку эта система не зависит от времени, возможно, эту часть итерации игнорировать, поэтому я просто имею dX=F(x,y,z) def …
1 ответ

Использование разных наборов моментов времени (параметр t) для встроенной функции python scipy ODEINT

Я работаю по системе Лоренца. Я использовал встроенную функцию scipy.integrate.odeint для интеграции, как это было предложено в Википедии 1. Система Лоренца имеет три переменные: x, y, z. Когда я сравнил эволюцию x для двух систем Лоренца с одинаков…
12 дек '19 в 13:59
3 ответа

Аттрактор Лоренца с питоном Рунге-Кутта

Здравствуйте, мне нужно запрограммировать функцию Python для решения дифференциальных уравнений Лоренца с использованием 2-го класса Рунге-Кутта сигма =10, r=28 и b=8/3 с начальными условиями (x,y,z)=(0,1,0) это код, который я написал, но он выдает …
1 ответ

В Python: как составить бифуркационную диаграмму системы Лоренца при изменении значения параметра?

Итак, я видел закодированное решение своего вопроса в Mathematica, но, имея очень слабое понимание математики, я еще не смог его воспроизвести. Это то, что я пытаюсь сделать с Python: https://mathematica.stackexchange.com/questions/159211/how-to-mak…
11 июл '20 в 00:04
1 ответ

Как складывать заполненные области вместо того, чтобы перекрывать их с помощью ggplot?

У меня есть график кривой Лоренца, который я заполнил факторными переменными (мужчина и женщина). Это было сделано достаточно просто, и перекрытие не было проблемой, потому что было всего два фактора. Wage %>% ggplot(aes(x = salary, fill = gender…
14 ноя '20 в 16:23
0 ответов

Внедрение системы Пекоры-Кэрролла Лоренца в Python?

У меня есть следующая система, которая представляет собой две системы Лоренца, соединенные вместе: def coupled_dynamics_Lorenz(x,t,sigma=10.0, beta=8.0/3.0, rho=28.0, k=0.0): xdot = np.zeros(6) xdot[0] = sigma*(x[1] - x[0]) + k*(x[3] - x[0]) xdot[1]…
22 окт '20 в 08:44
0 ответов

Как встроить трехмерный анимированный сюжет с помощью tkinter?

Я новичок в использовании tkinter и хотел создать небольшой пользовательский интерфейс для практики. Я выбрал построение системы Лоренца, позволяя пользователю настраивать входные параметры. Мне удалось подключить диаграмму и поля ввода, и я могу по…
0 ответов

Проблема в запуске кода ОДУ дробного порядка

Не удалось запустить пакет fodeint | Численное интегрирование обыкновенных дифференциальных уравнений дробного порядка для системы Лоренца. В этом пакете мне нужно запустить caputoEuler(a, f, y0, tspan); где а = 0,98 https://github.com/mattja/fodein…
18 янв '22 в 23:22
1 ответ

Как решить модель Лоренца 96 методом Рунге – Кутты?

Я написал приведенный ниже код для решения модели Лоренца 96 с использованием метода Рунге-Кутты, но результаты ненадежны, как видно на рисунке ниже: Правильное соотношение между тремя переменными выглядит следующим образом: Как я могу изменить код,…
1 ответ

RuntimeWarning с уравнениями Лоренца

Я пытаюсь выполнить уравнения Лоренца в питоне (я следую упражнению 8.3 от Марка Ньюмана - Вычислительная физика (2012, независимая издательская платформа CreateSpace)) Я уже получил графику, и все выглядит "правильно". Это, вероятно, математическая…
1 ответ

Лоренц Аттрактор

У меня есть вопрос, связанный с функцией numpy empty(). Я пытаюсь построить график хаоса, тем временем я использую функцию np.empty() для создания пустых массивов, состоящих из 20000 элементов. Код отлично работает, когда я использую функцию empty(…
06 июл '22 в 08:40
0 ответов

Решите и постройте уравнения Лоренца для двух различных начальных условий и двух значений rho в julia

Моя цель — решить уравнения Лоренца и построить их так, как показано на рисунке. У меня есть два разных начальных условия [x0, 1, 0] и x0 = 0, тогда x0 =1* 10^-5, два значения rho равны ρ = 14 и ρ = 28. Я написал этот код в julia, но получаю сообщен…
01 мар '23 в 05:49
0 ответов

Какие «реальные» задачи можно решить с помощью Lorenz Attractor?

Я создал блокнот с анализом аттрактора Лоренца и теперь ищу практические задачи, которые можно с его помощью решить. Причина, по которой мне это нужно, заключается в том, что без каких-либо практических задач мой ноутбук не приносит ничего инновацио…
10 фев '23 в 12:12
0 ответов

фазовый анализ уравнений Лоренца

Я пытаюсь построить анализ фазовой плоскости уравнений Лоренца. Например седло,спираль,узлы. это мой код x = np.linspace(-40,40,100) y = np.linspace(-60,50,100) z = np.linspace(-100,100,100) rho=28 beta=8/3 sigma=10 xx,yy = np.meshgrid(x,y) yy,zz = …
20 дек '22 в 20:58