[mathjiang]

楼主试试下述语句——求数值解。
function results = qiufangcheng
RH = 20:1:100; %% 湿度
Tdb = 18.0:0.1:85.0; %% 干球温度
r = length(RH); c = length(Tdb);
results = zeros(r,c);
for l = 1:r
for k = 1:c
results(l,k) = fzero(@(x)myfun(x,RH(l),Tdb(k)),1);
end
end

function y = myfun(x,RH,Tdb)
y = exp( (16.78.*x-116.9)./(x+237.3) ) - 0.0669*(Tdb-x) - (RH/100)*exp((16.78.*Tdb-116.9)/(Tdb+237.3));
end
end