淡水介电常数计算程序

淡水介电常数是科学研究中经常遇到的一个进行计算的问题，现在把计算程序代码列出，需要的话可以自行编译运行。

    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;