用matlab编写pmv方程

nibam11

我最近在学matlab，想用matlab编写PMV热舒适性方程，公式如下：
PMV=(0.303*exp((-0.036)*M)+0.028)*((M-W)-3.05*10^(-3)*(5733-6.99*(M-W)-Pa)-0.42*((M-W)-58.15)-1.7*10^(-5)*M *(5867-Pa)-0.0014*M*(34-ta)-3.96*10^(-8)*fcl*((tcl+273)^4-(tr+273)^4)-fcl*hc*(tcl-ta))；
其中tcl,Pa,hc,fcl,与以下公式有关：
tcl=35.7-0.028*(M-W)-I_clo.*(3.96*10^(-8).*fcl.*((tcl+273).^4-(tr+273).^4)+fcl.*hc.*(tcl-ta))；
Pa=phia*exp(16.6536-4030.183/(ta+235));
hc和fcl是分段函数：
hc:    if (2.38*(tcl-ta).^(1/4))>(12.1*sqrt(va))
              hc=2.38*(tcl-ta).^(1/4);
         elseif 2.38*(tcl-ta).^(1/4)<12.1*sqrt(va)
              hc=12.1*sqrt(va);
fcl:    if I_clo<=0.078
           fcl=1+1.290*I_clo;
        elseif I_clo>0.078
           fcl=1.05+0.645*I_clo;


式中：ta、va、tr、phia是输入变量，I_clo、M、W为恒定值。
请问如何实现PMV值随着各输入量的变化而变化的功能（例如在其他三个输入变量一定的情况下，PMV随着ta变化的情况）？