²é¿´: 737  |  »Ø¸´: 4

f2012

гæ (СÓÐÃûÆø)

[ÇóÖú] matlab ÓïÑÔ´íÎó

x0=[0.0350 0.16565];
x=fsolve('mnsi',x0)



ÒýÓÃmÎļþʱ±¨´íError: File: mnsi.m Line: 1 Column: 14
The expression to the left of the equals sign is not a valid target for an
assignment.


mÎļþfunction fx=mnsi(x)
T=1700;
a=[-45584.1,-39280.9,-10069.8,34484.37];
b=[28.50974,6.881164,2.113303,17.83834];
c=[0.5;0.5];
k=zeros(1,4);

ÓÐʲôÎÊÌâ Âé·³°ïæ¿´¿´
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
2Â¥2013-10-14 20:52:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

f2012

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
2Â¥: Originally posted by ÔÂÖ»À¶ at 2013-10-14 20:52:26
xÔÚÄÄÀïÄØ?

function fx=mnsi(x)
T=1700;
a=[-45584.1,-39280.9,-10069.8,34484.37];
b=[28.50974,6.881164,2.113303,17.83834];
c=[0.5;0.5];
k=zeros(1,4);
for i=1:1:4;
k(i)=exp(a(i)/T+b(i));
end
fx(1)=c(1)*x(1)- c(2)*x(2)+(3* c(1)- c(2))*k(1)*((x(1)^3)* x(2))+ (5* c(1)- 3*c(2))*k(2)*(x(1)^5)* (x(2))^3)+( (c(1)- c(2))*k(3)*(x(1)* x(2)+ (11* c(1)- 19*c(2))*k(4)*(x(1)^11)* (x(2))^19);
fx(2)=x(1)+x(2)+ k(1)*((x(1)^3)* x(2)+ k(2)*(x(1)^5)* (x(2))^3)+ k(3)*(x(1)* x(2)+ k(4)*(x(1)^11)* (x(2))^19);
Õâ¸öÊÇÈ«²¿µÄmÎļþ ÎÒ²»¶® ÐÂÊÖ¡£¡£¡£¡£¡£
3Â¥2013-10-14 21:01:02
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
f2012: ½ð±Ò+10, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸ 2013-10-15 20:17:22
ûʲô´óÎÊÌ⣬¾ÍÊÇ´ýÇó½â·½³Ìûд¶Ô¡£ÒÔÏÂËùÓгÌÐò¸´ÖƽøÒ»¸ömÎļþÔËÐм´¿É¡£Äã¿ÉÒÔ¼ì²éÒ»ÏÂfx(1)дµÃ¶Ô²»¶Ô¡£
CODE:
function solve_nonlEqs
format long
x0=[0.0350 0.16565];
x=fsolve(@mnsi,x0)
fval=mnsi(x)

function fx=mnsi(x)
T=1700;
a=[-45584.1,-39280.9,-10069.8,34484.37];
b=[28.50974,6.881164,2.113303,17.83834];
c=[0.5;0.5];
k=zeros(1,4);
for i=1:1:4;
k(i)=exp(a(i)/T+b(i));
end
fx(1)=c(1)*x(1)- c(2)*x(2)+(3* c(1)- c(2))*k(1)*((x(1)^3)* x(2))+ (5* c(1)- 3*c(2))*k(2)*(x(1)^5)* x(2)^3+(c(1)- c(2))*k(3)*(x(1)* x(2)+ (11* c(1)- 19*c(2))*k(4)*(x(1)^11)* (x(2))^19);
fx(2)=x(1)+x(2)+ k(1)*((x(1)^3)* x(2)+ k(2)*(x(1)^5)* (x(2))^3)+ k(3)*(x(1)* x(2)+ k(4)*(x(1)^11)* (x(2))^19);

½á¹û£º
x =

  1.0e-010 *

   0.159854118332060   0.159828336749788


·½³Ì²Ð²îfval =

  1.0e-010 *

   0.000012890791136   0.319682455081905
MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
4Â¥2013-10-15 09:04:27
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

f2012

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
4Â¥: Originally posted by ÔÂÖ»À¶ at 2013-10-15 09:04:27
ûʲô´óÎÊÌ⣬¾ÍÊÇ´ýÇó½â·½³Ìûд¶Ô¡£ÒÔÏÂËùÓгÌÐò¸´ÖƽøÒ»¸ömÎļþÔËÐм´¿É¡£Äã¿ÉÒÔ¼ì²éÒ»ÏÂfx(1)дµÃ¶Ô²»¶Ô¡£
function solve_nonlEqs
format long
x0=;
x=fsolve(@mnsi,x0)
fval=mnsi(x)

function fx=mn ...

Èç¹û·½±ãµÄ»°   Äܲ»ÄÜ  °ïÎÒ¿´Ï Õ⼸¸ö·½³Ì ÓÃʲô·½·¨½â±È½ÏºÃ  N1+N2+N3+N4+N5+N6 =1
aN1-bN2+(3a-b)N3+(5a-3b)N4+(a-b)N5+(11a-19b)N6=0
N3=K1*( N2^3)*N1
N4=K2*( N1^5)*(N2^3)
N5=K3*N1*N2
N6=K4* (N1^11)*(N2^19)

lnK1=45584.0871/T-28.50974
lnK2=39280.9406/T-6.8812
lnK3=10069.7385/T-2.1133
lnK4=34484.3661/T-17.838

T ¿ÉÒÔÈÎÒ⸳ֵÇó N1  N2      a b
ÖµÊÇÒÑÖªµÄa+b=1   ¿ÉÒÔËæÒ⸳ֵ
5Â¥2013-10-15 20:22:42
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ f2012 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸Öйúº£Ñó´óѧ£¬ÉúÎïѧ£¬301·Ö£¬Çóµ÷¼Á +5 1ËïÎò¿Õ 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[¿¼ÑÐ] 0817 »¯Ñ§¹¤³Ì 299·ÖÇóµ÷¼Á ÓпÆÑо­Àú ÓжþÇøÎÄÕ +18 rare12345 2026-03-18 18/900 2026-03-19 22:36 by ѧԱ8dgXkO
[¿¼ÑÐ] 085601²ÄÁϹ¤³Ìר˶Çóµ÷¼Á +10 Ľº®mio 2026-03-16 10/500 2026-03-19 15:26 by ¶¡¶¡*
[¿¼ÑÐ] 266Çóµ÷¼Á +5 ÑôÑôÍÛÈû 2026-03-14 10/500 2026-03-19 15:08 by ÑôÑôÍÛÈû
[¿¼ÑÐ] 317Çóµ÷¼Á +3 Éê×ÓÉêÉê 2026-03-19 6/300 2026-03-19 14:16 by Éê×ÓÉêÉê
[¿¼ÑÐ] Ò»Ö¾Ô¸ Î÷±±´óѧ £¬070300»¯Ñ§Ñ§Ë¶£¬×Ü·Ö287£¬Ë«·ÇÒ»±¾£¬Çóµ÷¼Á¡£ +3 ³¿»èÏßÓëÐǺ£ 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á£¬Çó¸÷λÀÏʦÊÕÁô +10 ÇïÓÐľ±± 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜר˶317Çóµ÷¼Á£¨0854¶¼¿ÉÒÔ£© +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[¿¼ÑÐ] 328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑо­Àú +3 ÉúÎ﹤³Ìµ÷¼Á 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[¿¼ÑÐ] »¯Ñ§¹¤³Ì321·ÖÇóµ÷¼Á +15 ´óÃ×·¹£¡ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[¿¼ÑÐ] 0703»¯Ñ§Çóµ÷¼Á ×Ü·Ö331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[¿¼ÑÐ] 274Çóµ÷¼Á +5 ʱ¼äµã 2026-03-13 5/250 2026-03-17 07:34 by ÈÈÇéɳĮ
[¿¼ÑÐ] 11408 Ò»Ö¾Ô¸Î÷µç£¬277·ÖÇóµ÷¼Á +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[¿¼ÑÐ] [µ¼Ê¦ÍƼö]Î÷ÄϿƼ¼´óѧ¹ú·À/²ÄÁϵ¼Ê¦ÍƼö +3 ¼â½ÇСºÉ 2026-03-16 6/300 2026-03-16 23:21 by ¼â½ÇСºÉ
[¿¼ÑÐ] 304Çóµ÷¼Á +3 ÂüÊâ2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[¿¼ÑÐ] 321Çóµ÷¼Á +5 ´óÃ×·¹£¡ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[¿¼ÑÐ] 0703 ÎïÀí»¯Ñ§µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-13 5/250 2026-03-16 10:50 by ÎÒ¿ÉÒÔÉϰ¶µÄ¶ÔÂ
[¿¼ÑÐ] 0856ר˶279Çóµ÷¼Á +5 ¼ÓÓͼÓÓÍ£¡? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[¿¼ÑÐ] 289Çóµ÷¼Á +4 ÕâôÃû×ÖÕ¦Ñù 2026-03-14 6/300 2026-03-14 18:58 by userper
[¿¼ÑÐ] 266Çóµ÷¼Á +4 ѧԱ97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û