Ä¿µÄ£º½âÓÉE1¡¢E2¡¢E3×é³ÉµÄ·ÇÏßÐÔ·½³Ì×é¡£Çó½â±äÁ¿Îªq,Tf1,Tf2
function S=myfun(q,Tf1,Tf2)
%ÊäÈë³£Êý
Ts=33
Ta=20
a=4
a=0.001*a
%µÚÒ»¸ö·½³ÌE1µÄ±í´ïʽ£º
b = (2.438714 + 0.7784798 * 0.01 * (0.5 * (Ts + Tf1)) - 0.17553068 * 10 ^ (-5) * (0.5 * (Ts + Tf1)) ^ 2) * 10 ^ (-2)
hcon=b/a
hr=( 1 / (1 / 0.97 + 1 / 0.97 - 1))* 5.67 * ((Ts + 273) ^ 4 - (Tf1 + 273) ^ 4) / (Ts - Tf1) * 10 ^ (-8)
hm=hcon+hr
E1=q-hm*(Ts-Tf1)
%µÚ¶þ¸ö·½³ÌE2µÄ±í´ïʽ
E2=q-(Tf1-Tf2)/(0.155*0.2325)
%µÚÈý¸ö·½³ÌE3µÄ±í´ïʽ
%b1,Cp,p,d,v¾ùÊÇÓëTf2Óйصĺ¯Êý
b1 = (2.438714 + 0.7784798 * 0.01 * (0.5*(Tf2+Ta)) - 0.17553068 * 10 ^ (-5) * (0.5*(Tf2+Ta)) ^ 2) * 10 ^ (-2)
Cp = (1005.28 - 0.260338 * 0.1 *(0.5*(Tf2+Ta)) + 0.6370071 * 10 ^ (-3) * (0.5*(Tf2+Ta)) ^ 2) * 10 ^ (-3)
p =1.2926 - 0.00463 * (0.5*(Tf2+Ta)) + 1.2619 * 10 ^ (-5) * (0.5*(Tf2+Ta)) ^ 2
d=1 / ((0.5*(Tf2+Ta)) + 273)
v = (13.213 + 0.09115128 * (0.5*(Tf2+Ta)) + 0.8758729 * 10 ^ (-4) * (0.5*(Tf2+Ta)) ^ 2) * 10 ^ (-6)
%¸ù¾Ýb1,Cp,p,d,vµÃ³öhcov
hcov=0.15*((9.81*d*b1^2*p*Cp*(Tf2-Ta))/v)^(1/3)
Hr1 =0.97* 5.67 * ((Tf2 + 273) ^ 4 - (Ta + 273) ^ 4) / (Tf2 - Ta) * 10 ^ (-8)
ht=hcov+Hr1
E3=q-ht*(Tf2-Ta)
%¹¹½¨myfunº¯Êý
S=[E1,E2,E3]
½¨Á¢myfunº¯Êýºó£¬½«Æä±£´æΪmyfun.m¡£È»ºóÔÚÃüÁî´°¿ÚÖÐÔËÐÐ
>> [q,Tf1,Tf2]=fsolve('myfun',[40,20,20])
½á¹û³öÏÖÒÔÏ´íÎó£º
??? Input argument 'Tf1' is undefined.
Error in ==> G:\matlab\work\myfun.m
On line 9 ==> b = (2.438714 + 0.7784798 * 0.01 * (0.5 * (Ts + Tf1)) - 0.17553068 * 10 ^ (-5) * (0.5 * (Ts + Tf1)) ^ 2) * 10 ^ (-2)
Error in ==> G:\matlab\toolbox\optim\fsolve.m
On line 161 ==> fuser = feval(funfcn{3},x,varargin{:});
Çë¸÷λ¸ßÊÖÖ¸µã£¬ÎÒÊÔÁ˺ܶà´ÎҲûÓгɹ¦¡£Çë°ÑÄãÃÇÇó½âµÄ½á¹ûÒ²¸½ÉÏÀ´°É
[ Last edited by nono2009 on 2009-11-1 at 15:21 ] |