Matlab实现 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*(v)^0.5; HCN=2.38*abs(100*XF-TAA)^0.25; if(HCF>HCN) HC=HCF; % Convective Heat Transfer Coefficient else HC=HCN; end XN=(P4+P1*HC-P2*XF^4)/(100+P3*HC); nIterations=nIterations+1; if ((nIterations>1) & (abs(XN-XF)<Tolerance)) break; end end if(nIterations<150) TCL=100*XN-273; % Surface Temperature of Clothing % Compute the Predicted Mean Vote (PMV) PM1=3.96*FCL*(XN^4-(TRA/100)^4); PM2=FCL*HC*(TCL-Ta); PM3=0.303*exp(-0.036*Met)+0.028; PM4=0.0; if(RMW>58.15) PM4=0.42*(RMW-58.15); end BMV=RMW-3.05*0.001*(5733-6.99*RMW-pa); CMV=-PM4-1.7*0.00001*Met*(5867-pa)-0.0014*Met*(34-Ta)-PM1-PM2; PMV1=PM3*(BMV+CMV); % Calculate the Predicted Percentage Dissatisfied (PPD) PPD=100-95*exp(-0.03353*PMV1^4-0.2179*PMV1^2); else PMV1=999; PPD=100; end disp(PMV1); disp(PPD);