Is there the function of "DEBYE" in MATLAB ?

27 visualizaciones (últimos 30 días)
ADNAN KIRAL
ADNAN KIRAL el 22 de En. de 2021
Respondida: Jens el 5 de Feb. de 2022
Hi
I need the "DEBYE model". how can I plot that in matlab ? any suggestion ?
thanks in advance
  7 comentarios
Walter Roberson
Walter Roberson el 22 de En. de 2021
X = logspace(-3,2);
Y = arrayfun(@debye1, X);
semilogx(X, Y)
However, this function has no parameters such as or W so it is not clear that it is the same model that you need.
ADNAN KIRAL
ADNAN KIRAL el 23 de En. de 2021
many thanks. Appriciate that.

Iniciar sesión para comentar.

Respuestas (1)

Jens
Jens el 5 de Feb. de 2022
Hello,
you mean the Debye-function, specific heat in 3D?
you can easily program it yourself. Either by numerical integration of the definition equatin, or a rational approximant:
R=8.314462618; % molar gas constant in J/mol/K CODATA 2018, exact.
xn=Theta./T;
I=@(x) x.^4.*exp(-x)./(exp(-x)-1).^2;
DEBYE=R*9*(1./xn).^3.*integral(I,0,xn,'reltol',1e-9,'abstol',1e-11);
or
N=[0.0279254827 -0.0358437761 1.05224663 17.2857105 64.7520511 226.68446];
D=[0.000119451046 -0.000210411972 0.00744002583 0.0346090463 0.337538084 ...
1.4272431 9.53997783 21.5840170 75.5614867]; % R. J. Goetsch+ (2011)
res1=R*polyval(N,Tn)./polyval(D,Tn);
res=reshape(res1,size(T));

Categorías

Más información sobre Programming en Help Center y File Exchange.

Productos


Versión

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by