%?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*（