PMV-PPD程序，精确度很高。

% PMV Calculation Matlab Code Edited by Shen Wei（Loughborough University）
% Referenced from ‘Human Thermal Environments‘ by Professor Ken Parsons and
% Thermal Comfort Tool Source Code

% Variable definition
Met=58.15; Clo=0.5; Ta=25; Tmrt=25; v=0.5; pa=1014; W=0;   % 1Met=58.15W/m2
RMW=Met-W;
Tolerance=0.00015;        % Accuracy definition

% Calculate FCl Value 
FCL=1.05+0.1*Clo;       % Clothing Area Factor
if（Clo<0.5）
    FCL=1.0+0.2*Clo;
end

% Calculate TCL Value
%% First guess for surface temperature
TAA=Ta+273;
TRA=Tmrt+273;
TCLA=TAA+（35.3-Ta）/（3.5*（Clo+0.1））;  
XN=TCLA/100;       %TCLA=Tcl+273
XF=XN;

%% Compute Surface Temperature of Clothing by Successive Substitution Iterations
FCIC=Clo*0.155*FCL;
P1=FCIC*TAA;
P2=FCIC*3.96;
P3=FCIC*100;
P4=308.7-0.028*RMW+P2*（TRA/100）^4;

nIterations=0;
XF=XN;

while（（nIterations<150））
    XF=（XF+XN）/2;
    
%% HC Calculation
    HCF=12.1*（