Refprop9.0工质热物理性质计算程序

NIST的著名热物理性质计算软件，可用Matlab,Fortran,Excel,C语言进行调用。压缩包中Readme给出了64位MATLAB调用该软件的办法，并成功验证。

function varargout = refproparray（ varargin ）
%refproparray Interface for refpropm which allows arrays to be entered
%instead of individual data points
%
%Function calls are very similar to refpropm （identical if only a single
%value is to be evaluated） except that vectors can be entered
%
%To use custom mixtures the substances and concentrations must be arranged
%in rows vectors
%Lengths of the values substances concentrations must be the same length
%or of length 1.
%
%Programmed by Scott Wujek Creative Thermal Solutions Fall 2013

%% convert known numeric input variables to numbers
propVal1 = cell2mat（varargin（3））;
propVal2 = cell2mat（varargin（5））;

%% determine vector lengths
L1=length（propVal1）;
L2=length（propVal2）;
L3=size（char（varargin（6））1）;
if nargin>6
    L4=size（cell2mat（varargin（nargin））1）;
    concentrations=cell2mat（varargin（nargin））;
else
    L4=NaN;
    concentrations=[];
end
maxsize=nanmax（[L1 L2 L3 L4]）;

%%  Make vectors to all be same length （expand singular entries）
if maxsize>1  && （L1==maxsize || L1==1） && （L2==maxsize || L2==1） && （L3==maxsize || L3==1） && （L4==maxsize || L4==1 || isnan（L4））
    if L1==1
        propVal1=repmat（propVal1 [maxsize1]）; %#ok
    end
    if L2==1
        propVal2=repmat（propVal2 [maxsize 1]）; %#ok
    end
    if L3==1
        if nargin==6
            varargin{6}=repmat（varargin（6） [maxsize1]）;
            
        else
            for inarg=6:nargin-1
                varargin{inarg}=repmat（varargin（inarg） [maxsize1]）;
            end
        end
    end
    if L4==1
        concentrations=repmat（concentrations [maxsize1]）; %#ok
    end
elseif maxsize~=1
    error（‘Input value vectors are of different lengths or not of length 1‘）
end

%% determine number of properties requested make output variable and left side of equation that replicates format of refpropm
propReq = lower（char（varargin（1）））;


if propReq==‘x‘
    varargout = cell（[12]）;
    variables=genvarname（repmat（{‘temp‘} [12]））;
    for i=1:length（variables）
        eval（[char（variables（i）） ‘=NaN（maxsize2）;‘]）;
    end
    string1=‘[temp（arraynum:） temp1（arraynum:）]‘;
else
    varargout = cell（size（propReq））;
    variables=genvarname（repmat（{‘temp‘} [1length（propReq）]））;
    string1=‘[‘;
    for i=1:length（variables）
        eval（[char（variables（i）） ‘=NaN（size（propVal1））;‘]）;
        string1=[string1 char（variables（i））]; %#ok
        if i~=length（variables）
            string1=[string1 ‘（arraynum）‘]; %#ok
        else
            string1=[string1 ‘（arraynum）]‘]; %#ok
        end
    end
end

%% create right side of refprop function call and evaluate function
for arraynum=1:maxsize
    if nargin==6 %this is case for pure substances or .mix files
        string2=‘=refpropm（char（varargin（1））char（varargin（2））propVal1（arraynum）char（varargin（4））propVal2（arraynum）char（varargin{6}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-07-26 09:57  REFPROP9\
     目录           0  2016-07-26 09:57  REFPROP9\Latest Update\
     文件     2272768  2013-12-18 23:22  REFPROP9\Latest Update\REFPRP64.DLL
     文件       60416  2013-12-18 23:21  REFPROP9\Latest Update\REFPRP64_thunk_pcwin64.dll
     文件        3581  2013-12-18 23:22  REFPROP9\Latest Update\refproparray.m
     文件       30356  2013-12-18 23:20  REFPROP9\Latest Update\refpropm.m
     文件       51741  2013-12-18 23:22  REFPROP9\Latest Update\rp_proto64.m
     目录           0  2016-07-26 09:57  REFPROP9\ML2010A\
     文件        4010  2011-08-17 14:58  REFPROP9\ML2010A\refprop_example.m
     文件       20225  2011-10-03 09:24  REFPROP9\ML2010A\refpropm.f90
     文件        3521  2011-10-03 09:19  REFPROP9\ML2010A\refpropm.m
     文件      856064  2011-10-03 09:23  REFPROP9\ML2010A\refpropm.mexw32
     文件     1254400  2011-04-07 10:59  REFPROP9\ML2010A\refpropm.mexw64
     目录           0  2016-07-27 18:08  REFPROP9\REFPROP9\
     文件     4737745  2011-07-26 00:07  REFPROP9\REFPROP9\NIST2309.EXE
     文件     2706944  2016-03-05 10:30  REFPROP9\REFPROP9\REFPRP64.DLL
     文件       60416  2016-03-05 10:30  REFPROP9\REFPROP9\REFPRP64_thunk_pcwin64.dll
     目录           0  2016-07-26 09:57  REFPROP9\REFPROP9\Refprop9 Excel Add-In\
     文件        1136  2012-02-24 02:44  REFPROP9\REFPROP9\Refprop9 Excel Add-In\EffyFixedVolRatio.txt
     目录           0  2016-07-26 09:57  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\
     目录           0  2016-07-26 09:57  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\Refprop Excel Add-In\
     文件      171535  2012-03-22 02:08  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\Refprop Excel Add-In\RefpropExcelAddIn.xlam
     文件       21373  2012-03-22 01:32  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\Refprop Excel Add-In\RefpropExcelAddInCreationNotes20120316.docx
     文件      108473  2012-03-16 23:48  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\Refprop Excel Add-In\RefpropExcelAddInHelp.xlsm
     文件       13393  2012-03-16 21:27  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\Refprop Excel Add-In\RefpropExcelAddInInstallationNotes20120316.docx
     文件        9219  2012-01-27 04:10  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\Refprop Excel Add-In\TryThis20120126.xlsm
     文件       13393  2012-03-16 21:27  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files\RefpropExcelAddInInstallationNotes20120316.docx
     目录           0  2016-07-26 09:57  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files （x86）\
     目录           0  2016-07-26 09:57  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files （x86）\Refprop Excel Add-In\
     文件      156404  2012-03-22 02:14  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files （x86）\Refprop Excel Add-In\RefpropExcelAddIn.xlam
     文件       21373  2012-03-22 01:32  REFPROP9\REFPROP9\Refprop9 Excel Add-In\IfRefprop9InProgram Files （x86）\Refprop Excel Add-In\RefpropExcelAddInCreationNotes20120316.docx
............此处省略11个文件信息