Здраствуйте, в рамках своей работы я столкнулся с моделирование процесса переноса сигнала между нейронами человеческого мозга. Моделирую это процесс следующим уравнением
![$i\hbar \frac{\partial \psi(x,t)}{\partial t} = -\frac{\hbar^2}{2m} \frac{\partial^2 \psi(x,t)}{\partial x^2} + V(x,t) |\psi(x,t)|^2 \psi(x,t)$ $i\hbar \frac{\partial \psi(x,t)}{\partial t} = -\frac{\hbar^2}{2m} \frac{\partial^2 \psi(x,t)}{\partial x^2} + V(x,t) |\psi(x,t)|^2 \psi(x,t)$](https://dxdy-01.korotkov.co.uk/f/c/e/5/ce57e346503ab7ed66747ac3e58e7b0b82.png)
![$\psi(x,t)$ $\psi(x,t)$](https://dxdy-03.korotkov.co.uk/f/6/5/3/65388eac83fde2cd79458efe154770d582.png)
- представляет собой комплексную величину, которая содержит информацию об амплитуде и фазе сигнала. Амплитуда дает представление о величине сигнала в точке x в момент времени t.
![$|\psi(x,t)|^2$ $|\psi(x,t)|^2$](https://dxdy-04.korotkov.co.uk/f/b/b/a/bbae63ed883a2c2fcb8ab995e22325bb82.png)
- абсолютное значение есть плотность сигнала или проще говоря, вероятность наличия сигнала в точке x в момент времени t
![$V(x,t)$ $V(x,t)$](https://dxdy-03.korotkov.co.uk/f/a/6/4/a64748703bdd906bdf08d292f21db5ec82.png)
- потенциальная энергия, зависящая от координаты x и времени t. В контексте моей модели это может быть потенциал мембраны аксона, который меняется в зависимости от локального окружения и времени.
Это уравнение я преобразовываю методом разделения переменных,
![$\psi(x,t)=\varphi(x)\chi(t), \chi(t)=Ae^\(i\omega t\)$ $\psi(x,t)=\varphi(x)\chi(t), \chi(t)=Ae^\(i\omega t\)$](https://dxdy-02.korotkov.co.uk/f/d/0/e/d0e1f5c8c1c78de20c8bb6e7ac5e00d482.png)
, подставив в первое уравнение получаю следующее уравнение
![$E\varphi(x) = -\frac{\hbar^2}{2m} \frac{\partial^2 \varphi(x)}{\partial x^2} + V(x,t) A^2|\varphi(x)|^2 \varphi(x)$ $E\varphi(x) = -\frac{\hbar^2}{2m} \frac{\partial^2 \varphi(x)}{\partial x^2} + V(x,t) A^2|\varphi(x)|^2 \varphi(x)$](https://dxdy-04.korotkov.co.uk/f/7/0/3/70350c8287723d64a36689ae451eed4382.png)
A - нормированное значение амплитуды и ровно 1.
Данное уравнение я планировал решать численно методом Рунге-Кутты 4-го порядка, разложив вторую производную, методом замены переменной тем самым перейдя к производной первого порядка.
Однако при решение я столкнулся с проблемой отсутствия моделей для задание потенциала
![$V(x,t)$ $V(x,t)$](https://dxdy-03.korotkov.co.uk/f/a/6/4/a64748703bdd906bdf08d292f21db5ec82.png)
, где имеется зависимость от координаты, есть популярная модель Ходжкина-Хаксли, а также модель Ижикевича, но они в своем стандартном виде не подходят так как зависят только от времени. Можете подсказать каким образов я могу задать уравнение для потенциала с зависимостью и от времени и от координаты.
Так же для решение численным методом, мне нужно задать начальные значения для функции
![$\varphi(x)$ $\varphi(x)$](https://dxdy-02.korotkov.co.uk/f/5/1/7/51733a6ece5add8ff3a19f7275dd196d82.png)
, но не могу подобрать функцию для её задания. Можете и с этим вопросом подсказать?