Is there the function of "DEBYE" in MATLAB ?

18 views (last 30 days)
ADNAN KIRAL
ADNAN KIRAL on 22 Jan 2021
Answered: Jens on 5 Feb 2022
Hi
I need the "DEBYE model". how can I plot that in matlab ? any suggestion ?
thanks in advance
  7 Comments
Walter Roberson
Walter Roberson on 22 Jan 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.

Sign in to comment.

Answers (1)

Jens
Jens on 5 Feb 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));

Categories

Find more on Programming in Help Center and File Exchange

Products


Release

R2020a

Community Treasure Hunt

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

Start Hunting!