1

我正在尝试模拟一个非常用于分析磁性材料的模型。最初，我用python写了这段代码，但没有得到结果；然后，我决定使用 Mathematica 进行相同的模拟 [ 1 ] 并且它似乎工作正常，至少它显示了与预期相似的曲线，但是当它被转换为 python 时（使用 lmfit 包和 scipy.integrate.quad）它只显示一条线。我还应该考虑什么才能让它在 python 中正常工作？

数学代码

这是我所拥有的：

from lmfit import minimize, Minimizer, Parameters, Parameter, report_fit
from numpy import loadtxt, vectorize, sqrt, log, inf, exp, pi, tanh,     linspace
from scipy.integrate import quad
from scipy.constants import codata
import matplotlib.pyplot as plt
from mpmath import coth


T1 = linspace(1, 300, 100)

def PDF(Dmag, v):
#    v = params.valuesdict()
    pdf=(1/((Dmag*1e-9)*v['w']*sqrt(2*pi)))*exp(-((log((Dmag*1e-9)/Dmag))**2)/(2*v['w']**2))*1e-9
    return pdf

def mZFC(Dmag, args):
    T1, params = args
    v  = params.valuesdict()
    V  = ((pi/6)*(Dmag*1e-9)**3)
    nu = v['v0']*exp(-(V*v['Kv'])/(v['kB']*T1))
    dt = v['kB']*T1**2/((V*v['Kv'])*v['vT']);
    return (v['B']*(v['Ms']*V)**2/(3*(V*v['Kv'])))*(exp(-nu*dt) +  ((V*v['Kv'])/(v['kB']*T1))*(1 - exp(-nu*dt)))*PDF(Dmag, v)

def curveZFC(T1, params):
    return params['NT']*quad(mZFC, 0.0, inf, [T1, params])[0]

vcurveZFC = vectorize(curveZFC, excluded=set([1]))


kB = codata.value('Boltzmann constant')

params = Parameters()
params.add('Dm' , value= 3.2e-9 , vary= True)
params.add('w'  , value= 0.26    , vary= True)
params.add('Ms' , value= 1350000   , vary= True)
params.add('T'  , value= 300     , vary= False)
params.add('kB' , value= kB      , vary= False)
params.add('NT' , value= 1.7e12  , vary= True)
#m_ZFC params
params.add('B'   , value= 0.0050, vary= False)
params.add('v0'  , value= 1e10 , vary= False)
params.add('Kv'  , value= 178000  , vary= False)
params.add('vT'  , value= 0.05 , vary= False)


y2init = vcurveZFC(T1, params)

plt.plot(T1, y2init, 'r*')
plt.show()

任何帮助，将不胜感激。谢谢！

4

0 回答 0