pmvf.m

由Fanger教授提出PMV－PPD评价指标广泛用于室内热环境评价

%PMVF Fanger教授PMV－PPD指标计算程序
%
% 大纲:	[PMVOUT]=pmvf(PMVPAR)
%			
% 描述:
%  PMVPAR 是一个包含7个变量用以计算PMV和PPD的计算函数 
% 	PMVPAR(:,1): 人体能量代谢率 [W/m^2]
% 	PMVPAR(:,2): 人体对外所做机械功 [W/m^2]
% 	PMVPAR(:,3): 平均辐射温度 [℃]     
% 	PMVPAR(:,4): 空气温度 [℃]
% 	PMVPAR(:,5): 相对湿度    [0 < Rh < 1]
% 	PMVPAR(:,6): 服装热租 [clo]    
% 	PMVPAR(:,7): 空气流速 [m/s]
%
%  PMVOUT 输出一个二维数组，第一列为PMV计算结果
%  第二列为PPD的计算结果
% 例如:
%
%  PMVVAR=[58.2 0 20 20 0.5 1 0.2;58.2 0 20 20 0.5 1 0.3];
%  [PMVOUT]=pmvf(PMVVAR)
%
%  PMVOUT =
%
%   -1.1337   32.0718
%   -1.2965   40.0881
%   





function q=pmvf(in)

insize=size(in); 

q=zeros(insize(1),2);

if insize(2)==7
 q=zeros(insize(1),1);
  for i=1:insize(1)
  pm=in(i,:);
  uu11=pm;
  pm(6)=pm(6)*0.155;
   if pm(6) <= 0.078,
     fclpmv=1+1.29*pm(6);
   else
     fclpmv=1.05+0.645*pm(6);
   end
  pm(8)=fclpmv;
  x=fminsearch('pmveqf',[30 30]',[],pm);
  
  tclpmv=x(1);
  hcpmv=x(2);
  pm=uu11;
  q1=(0.303*exp(-0.036*pm(1))+0.028);
  q(i,1)=q1*(pm(1)-pm(2)-3.05e-3*(5733-6.99*(pm(1)-pm(2))-pm(5)*psatf(pm(4)))-0.42*(pm(1)-pm(2)-58.15)-1.7e-5*pm(1)*(5867-pm(5)*psatf(pm(4))) -0.0014*pm(1)*(34-pm(4))-3.96e-8*fclpmv*( (tclpmv+273)^4 -(pm(3)+273)^4 ) - fclpmv*hcpmv*(tclpmv-pm(4)));
  q(i,2)=ppdf(q(i,1)); 
end

else
  disp('Wrong inputmatrix. Inputmatrix must be: n-by-7 ')
  disp('CALCULATION INTERRUPTED !')
  q=[]

end