Очевидно 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^3Rmin и 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, и мы пытаемся адаптировать ее для более удобной обработки данных метод. Но я геохимик, поэтому в принципе бесполезен. Любая помощь приветствуется! Я могу включить больше подробностей, если какая-то добрая душа готова помочь.