Is there the function of "DEBYE" in MATLAB ?
16 次查看(过去 30 天)
显示 更早的评论
Hi
I need the "DEBYE model". how can I plot that in matlab ? any suggestion ?
thanks in advance
7 个评论
Walter Roberson
2021-1-22
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.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/497273/image.png)
回答(1 个)
Jens
2022-2-5
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));
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!