Описание тега differentialequations.jl

Вопросы с тегом
1 ответ

Проблема с diffrentialEquations.jl

Я совсем новичок в Юлии, и в настоящее время я учусь, как решать дифференциальные уравнения с ним. Я пытался запустить простой готовый код Кристофера Ракауцкаса, но получил ошибку. Код можно найти здесь. Я также напишу это здесь: using DifferentialE…
10 май '17 в 08:02
1 ответ

Решить уравнение теплопроводности с ненулевыми БК Дирихле с неявным линейным решателем Эйлера и сопряженного градиента?

Многие пользователи спрашивают, как решить уравнение теплопроводности, u_t = u_xx, с ненулевыми BC Дирихле и с сопряженными градиентами для внутреннего линейного решателя. Это общая упрощенная проблема PDE перед переходом к более сложным версиям пар…
06 фев '19 в 02:13
0 ответов

Юлия Дифференциальные уравнения: не удается решить проблему MonteCarloProblem, когда работает "ручное" решение (ошибка преобразования)

Я работаю над немного большим проектом, в котором я активно использую типы MonteCarloProblem DifrentialEquations.jl. Обычно все работает нормально, но я столкнулся со странным случаем, в котором я не могу разобраться, потому что ручная итерация по к…
04 дек '18 в 10:26
1 ответ

Выполнение проблемы ODEP, пока не будет подан сигнал

Возможно ли это за время ODEProblem быть открытым, т. е. чтобы проблема выполнялась до получения внешнего сигнала (возможно, с помощью дискретного обратного вызова с соответствующим условием)? Интуитивно это означало бы дать что-то вроде tspan = (0.…
07 фев '18 в 23:41
1 ответ

Юлия 0.0.6 LoadError на ODE решить

Я пытаюсь решить простую модель Lotka Volterra using DifferentialEquations f2 = @ode_def_nohes LVtest begin dx = x*(1.0 - A*x - B*y) dy = y*(rho - C*x - D*y) end A B C D rho u0 = [1;1] tspan = (0.0,300.0) p = [0.2,0.3,0.35,0.2,0.25] prob = ODEProble…
22 авг '18 в 20:14
1 ответ

Юлия Плотс; Как я могу увеличить количество образцов / точек данных?

При решении дифференциальных уравнений и построении графиков результатов, как я могу увеличить количество отображаемых точек данных? я имею using DifferentialEquations using Plots function lorenz(du,u,p,t) du[1] = 10.0*(u[2]-u[1]) du[2] = u[1]*(28.0…
27 окт '18 в 23:44
1 ответ

Арифметика произвольной точности в JIT-скомпилированных функциях

Когда я использую numba внутри python, я знаю, что если я попытаюсь jit-скомпилировать функции, которые имеют произвольные значения с плавающей точкой (mpmath) внутри своих циклов, он не сможет скомпилироваться в режиме nopython, и его скорость буде…
29 май '18 в 15:19
1 ответ

Кусочно-дифференциальное уравнение с использованием Юлии

Я новичок в Юлии, я хотел бы решить эту систему \frac{dx}{dt} = k1y \ \frac{dy}{dt} = k2y+I где k1 и k2 - постоянные параметры. однако, I=0, когда y,0 или Ky в противном случае, где k является постоянным значением. Я следовал учебному пособию по ODE…
05 июл '18 в 16:55
1 ответ

Не удалось выполнить первый пример дифференциально-алгебраических уравнений

Я использую JuliaPro v0.6.0.1 и JunoIDE, которые я пытался применить diffrentialEquations.jl. Чтобы запустить первый пример, я добавил пакет DifferentialEquations с последующим оператором using. На следующем шаге я скопировал первый пример: f(t,u) =…
05 окт '17 в 12:33
1 ответ

Генерация функции в цикле для дифференциальных уравнений в Юлии

Принимая пример Лоренца ( JuliaDiffEq GitHub) function lorenz(t,u,du) du[1] = 10.0*(u[2]-u[1]) du[2] = u[1]*(28.0-u[3]) - u[2] du[3] = u[1]*u[2] - (8/3)*u[3] end Если бы мне пришлось генерировать уравнения в цикле, я бы попытался объединить строки в…
22 дек '17 в 07:26
1 ответ

Использование пакета diffrentialEquations в Julia для решения матрицы ODE

Я хотел бы решить: [\ mathbf {M} \ ddot {\ mathbf {U}} + \ mathbf {C} \ dot {\ mathbf {U}} + \ mathbf {K} \ mathbf {U} = \ mathbf {P} (t)] Или в форме пространства состояний: [\ Точка {\mathbf{Y}}= F (\mathbf{Y}, т)] где: [\ mathbf {Y} = \ left [\ b…
19 июл '17 в 14:22
2 ответа

Имитировать прыгающий мяч?

Можно ли создать простую модель прыгающего мяча, используя решатели уравнения Юлии? Я начал с этого: using ODE function bb(t, f) (y, v) = f dy_dt = v dv_dt = -9.81 [dy_dt, dv_dt] end const y0 = 50.0 # height const v0 = 0.0 # velocity const startpos …
13 окт '16 в 22:14
1 ответ

Ошибка метода при использовании пакета diffrentialEquations.jl Julia

Я пытаюсь решить дифференциальное уравнение ode45 с помощью пакета DifferentialEquation.jl, но получаю ошибку метода. using DifferentialEquations M = 400; m = 35; C = 3e3; c = 300; K = 50e3; k = 200e3; A = 0.05; L = 0.5; vh = 13.9 MM = [M 0; 0 m] # …
13 сен '17 в 20:10
1 ответ

Параболические PDE в Юлии

Я пытаюсь решить параболическое дифференциальное уравнение в частных производных численно, используя Джулию, но я не могу найти доступную документацию, которая может помочь. Вот пример: t, x являются 1-мерным вещественным веществом. Я хочу решить дл…
17 янв '18 в 23:17
1 ответ

Решить систему ОДУ с чтением во внешнем форсировании

В Юлии я хочу решить систему ОДУ с внешним принуждением g1(t), g2(t) лайк dx1(t) / dt = f1(x1, t) + g1(t) dx2(t) / dt = f2(x1, x2, t) + g2(t) с принуждениями, считанными из файла. Я использую это исследование для изучения Юлии и пакета diffrentialEq…
22 мар '18 в 12:39
1 ответ

Прекращение интеграции оды после некоторого теста сходимости в DiffEqPhysics Юлия

Недавно я начал знакомиться с пакетом DiffEqPhysics и diffrentialEquations Джулии. Мне интересно, есть ли способ (например, функции обратного вызова) завершить решатель с определенным условием сходимости на выходе. Например, если результат, полученн…
27 май '18 в 04:15
1 ответ

Используя Дифференциальные уравнения: вы не обновляете

Я считаю, что в этом коде есть ошибка. Для краткости я просто напишу функцию, которая определяет ODE function clones(du,u,p,t) (Nmut,f) = p # average fitness phi = sum(f.*u) # constructing mutation kernel eps = 0.01 Q = Qmatrix(Nmut,eps) # defining …
29 сен '18 в 23:05
1 ответ

Юлия Дифференциал Equations.jl скорость

Новичок в Юлии, пытающейся проверить скорость решения ODE. Я использовал уравнение Лоренца в учебнике using DifferentialEquations using Plots function lorenz(t,u,du) du[1] = 10.0*(u[2]-u[1]) du[2] = u[1]*(28.0-u[3]) - u[2] du[3] = u[1]*u[2] - (8/3)*…
26 ноя '17 в 22:50
1 ответ

Как исправить "LoadError: DimensionMismatch (" не может широковещательный массив иметь меньше измерений ")"

Я хотел бы решить следующие два связанных дифференциальных уравнения численно: d/dt Phi_i = 1 - 1/N * \sum_{j=1}^N( k_{ij} sin(Phi_i - Phi_j + a) d/dt k_{ij} = - epsilon * (sin(Phi_i - Phi_j + b) + k_{ij} с определенными начальными условиями phi_0 (…
11 фев '19 в 11:13
1 ответ

Анимация решения для ODE в Юлии

У меня есть код Джулии: using DifferentialEquations using Plots using ParameterizedFunctions plotly() lorenz = @ode_def Lorenz begin dx = σ*(y-x) dy = ρ*x-y-x*z dz = x*y-β*z end σ = 10. β = 8./3. ρ => 28. u0 = [1., 5., 10.] tspan = (0., 2.) prob …
14 фев '17 в 00:00