| ²é¿´: 694 | »Ø¸´: 4 | |||
| µ±Ç°Ö÷ÌâÒѾ´æµµ¡£ | |||
cooler8395½ð³æ (³õÈëÎÄ̳)
|
[½»Á÷]
¡¾ÇóÖú¡¿MATLAB»Ø¹éwilson·½³Ì²ÎÊý
|
||
|
Âé·³¸ßÊÖ£¬ÇóÖúһϻعéwilson·½³Ì²ÎÊýµÄmatlab³ÌÐò¡£×Ô¼º±àµÄ²»ÖªµÀÄÄÀïµÄÔÒò£¬Îó²îºÜ´ó¡£ % wilson »Ø¹é²ÎÊý µ÷Óøñʽ-----[x,fval]=fminsearch(@wilson,[100 100])µ¥´¿Ðη¨Çó×îСֵ function f = wilson(x) A12=x(1); A21=x(2); YY=0; YYY=0; FFF=0; X1=[0.0593 0.1306 0.2660 0.3842 0.4708 0.6256 0.7178 0.8214 0.9037 0.9397]; Y1=[0.3285 0.4902 0.6234 0.6831 0.7125 0.7524 0.7832 0.8349 0.8864 0.9215]; for i=1:length(X1) i x1=X1(i);%ÌáÈ¡±¾´ÎÊÔÑéµãµÄʵÑéÖµ x2=1-x1; % x1,x2 ÊÔÑéÖµ y1=Y1(i); y2=1-y1; % y1,y2 ÊÔÑéÖµ A12; A21; Lnr1=-log(x1+A12*x2)+x2*(A12/(x1+A12*x2)-A21/(x2+A21*x1)); Lnr2=-log(x2+A21*x1)+x1*(A21/(x2+A21*x1)-A12/(x1+A12*x2)); r1=exp(Lnr1); r2=exp(Lnr2); A1=4.1973;A2=3.4337; B1=1575.0;B2=1413.0; C1=-34.29;C2=-44.25; t=70; T=t+273.15; P10=exp(A1-B1/(C1+t)); P20=exp(A2-B2/(C2+t)); P=x1*r1*P10+x2*r2*P20 % P-×Üѹ£»P10-ÎïÖÊ1µÄ±¥ºÍÕôÆûѹ£»P20-ÎïÖÊ2µÄ±¥ºÍÕôÆûѹ. µ¥Î»£ºKpa y1cal=x1*r1*P10/P y2cal=x2*r2*P20/P %------------------------------ Y=abs(y1-y1cal); %Y-±¾´Î¼ÆËãµãµÄ¾ø¶ÔÎó²î YY=YY+Y; %YY-ÒѼÆËãµãµÄ¾ø¶ÔÎó²îÖ®ºÍ %------------------------------- MM=abs((y1-y1cal)./y1); YYY=YYY+MM; %------------------ m=abs(y1cal-y1); M(1,i)=m; %----------------------- FF=(y1cal-y1).^2+(y2cal-y2).^2; FFF=FFF+FF;%Ä¿±êº¯Êý %a=log(r1/r2); %b=log(r1)/x2^2; % c=-log(r2)/x1^2; end f=FFF; fata=YY/i%ƽ¾ùÆ«²î sigma=YYY/i%ƽ¾ùÏà¶ÔÆ«²î sigmaM=max(M)%×î´óÆ«²î [ Last edited by cooler8395 on 2009-8-3 at 20:45 ] |
» ²ÂÄãϲ»¶
282Çóµ÷¼Á
ÒѾÓÐ22È˻ظ´
283Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
265Çóµ÷¼Á
ÒѾÓÐ21È˻ظ´
288Çóµ÷¼Á Ò»Ö¾Ô¸¹þ¹¤´ó ²ÄÁÏÓ뻯¹¤
ÒѾÓÐ10È˻ظ´
µ÷¼Á0855-288
ÒѾÓÐ5È˻ظ´
µç×ÓÐÅÏ¢µ÷¼Á½»²æÑ§¿ÆÓÐÍÆ¼öÂð
ÒѾÓÐ6È˻ظ´
266Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
²ÄÁÏ¿ÆÑ§Ó빤³Ì¿¼ÑÐ
ÒѾÓÐ10È˻ظ´
283·Ö²ÄÁÏÓ뻯¹¤Çóµ÷¼Á
ÒѾÓÐ29È˻ظ´
Ò»Ö¾Ô¸¹þ¶û±õ¹¤Òµ´óѧ085600Ó¢Ò»Êý¶þ337·ÖÇóµ÷¼Á
ÒѾÓÐ4È˻ظ´
![]() ![]() |
2Â¥2009-08-03 16:46:07
cooler8395
½ð³æ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 1455.9
- Ìû×Ó: 41
- ÔÚÏß: 6.4Сʱ
- ³æºÅ: 679287
- ×¢²á: 2008-12-22
3Â¥2009-08-03 18:44:12
¶¬ÌìÀïµÄ½¾Ñô
Ìú¸Ëľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 12 (СѧÉú)
- ½ð±Ò: 7179.7
- É¢½ð: 10
- Ìû×Ó: 387
- ÔÚÏß: 318Сʱ
- ³æºÅ: 439463
- ×¢²á: 2007-10-13
- ÐÔ±ð: GG
- רҵ: ÀíÂۺͼÆË㻯ѧ
4Â¥2009-11-10 16:29:19
¶¬ÌìÀïµÄ½¾Ñô
Ìú¸Ëľ³æ (ÕýʽдÊÖ)
- Ó¦Öú: 12 (СѧÉú)
- ½ð±Ò: 7179.7
- É¢½ð: 10
- Ìû×Ó: 387
- ÔÚÏß: 318Сʱ
- ³æºÅ: 439463
- ×¢²á: 2007-10-13
- ÐÔ±ð: GG
- רҵ: ÀíÂۺͼÆË㻯ѧ
¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
cooler8395(½ð±Ò+10,VIP+0): 11-10 16:30
nono2009(½ð±Ò-10,VIP+0):Ç벻Ҫͨ¹ýÂí¼××ªÒÆ½ð±Ò£¡¿Û³ý×ªÒÆ½ð±Ò¡£ 11-10 16:50
nono2009(½ð±Ò+10,VIP+0):¾½âÊÍ˵Ã÷£¬½ð±ÒÍË»¹¡£µ«²»¹ÄÀøÍ¨¹ýÕâÖÖ·½Ê½½±Àø½ð±Ò£¬Ï²»ÎªÀý¡£ 11-11 11:24
cooler8395(½ð±Ò+10,VIP+0): 11-10 16:30
nono2009(½ð±Ò-10,VIP+0):Ç벻Ҫͨ¹ýÂí¼××ªÒÆ½ð±Ò£¡¿Û³ý×ªÒÆ½ð±Ò¡£ 11-10 16:50
nono2009(½ð±Ò+10,VIP+0):¾½âÊÍ˵Ã÷£¬½ð±ÒÍË»¹¡£µ«²»¹ÄÀøÍ¨¹ýÕâÖÖ·½Ê½½±Àø½ð±Ò£¬Ï²»ÎªÀý¡£ 11-11 11:24
| ˽ÏÂÁÄÁİɣ¬ÄãµÄÎÊÌâÎÒÄܽâ¾ö |
5Â¥2009-11-10 16:29:54














»Ø¸´´ËÂ¥