Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#4 |
高级会员
注册日期: 2008-01-02
帖子: 247
声望力: 22 ![]() |
![]()
楼主试试下述语句——求数值解。
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 |
![]() |
![]() |