摘要: 淡水介电常数是科学研究中经常遇到的一个进行计算的问题,现在把计算程序代码列出,需要的话可以自行编译运行。 function result=debye(f,temp) %compute dielectric constant of wate...
淡水介电常数是科学研究中经常遇到的一个进行计算的问题,现在把计算程序代码列出,需要的话可以自行编译运行。
function result=debye(f,temp) %compute dielectric constant of water by debye model %edited by zhengxingming 2010.12.06 dcw_inf=4.9; %dielectric constant of water under infinite high frequency dcw_zero = 88.045 - 0.4147 * temp + 6.295 * 0.0001 * temp * temp + 1.075 * 0.00001 * temp * temp * temp; tau = 0.11109 - 0.003824 * temp + 0.00006938 * temp * temp - 5.096 * 0.0000001 * temp * temp * temp; dcw_real = dcw_inf + (dcw_zero - dcw_inf) ./ (1.0 + (tau.*f.*tau.*f)) ; dcw_imag = tau .* f .* (dcw_zero - dcw_inf) ./ (1.0 + (tau.*f).^2); hold on; loglog(f,dcw_real,'r-',f,dcw_imag,'r:') result=dcw_real+dcw_imag*j;