Оптимизация значений параметров с использованием нелинейных наименьших квадратов в R (с интегралами) - PullRequest
Новая
       61

Оптимизация значений параметров с использованием нелинейных наименьших квадратов в R (с интегралами)

0 голосов
/ 06 февраля 2020

Очевидно R (и математика) любитель. Я работал более 10 часов, пытаясь заставить это работать, поэтому я решил попытаться опубликовать здесь как снимок.

У меня есть данные, собранные из эксперимента с двумя переменными: Iq и q. Эти данные являются линейными при построении графика в пространстве журнала. Я пытаюсь найти для двух других переменных, por и r, в следующем уравнении: 

Iq=SLD^2*(por/Vra)*integral{Rmin to Rmax}((Vr)^2*f(r)*F dr)

Где:

  • SLD = известная константа
  • por = неизвестно
  • Vra=integral{0 to Inf}(Vr*f(r)dr)
  • Vr=(4/3)*pi*r^3
  • Rmin и Rmax = известные константы
  • f(r)=((r^-(1+fd))/(Rmin^(-fd) - Rmax^(-fd))/fd)
  • r = неизвестно
  • fd = известная константа
  • F=(3*(sin(q*r)-q*rcos(q*r))/(q*r)^3)^2

Я пытался сделать много попыток, но не могу обернуть мой мозг вокруг переменных внутри переменных в код. Раньше эта проблема решалась в программе решателя Excel, которая оптимизировала значения параметров с использованием нелинейных наименьших квадратов, которые работают только в (imo) Windows 95 Excel, и мы пытаемся адаптировать ее для более удобной обработки данных метод. Но я геохимик, поэтому в принципе бесполезен. Любая помощь приветствуется! Я могу включить больше подробностей, если какая-то добрая душа готова помочь.

...