[½»Á÷] MATLAB±à³Ì£¬·ÇÏßÐÔÄâºÏ¡£

514169Â¥: Originally posted by xiegangmai at 2011-04-09 21:24:40
¼òµ¥ÊÔÁËÊÔ£¬×Ô¼º¸Ä¸Ä³õÖµ£¬ÉèÖÃÒ»ÏÂoptionÊÔÊÔ°É¡£

function qiancai
clear
clc
close
x0 = ;
% x0 = ;
= lsqnonlin( @qiancaifun, x0 )
L = ;

mu = ;

plot( L, mu, 'r-' )
muu =  ( x( 3 ) +  ...

Originally posted by qiancai at 2011-04-06 20:17:51:
¡¾ÇóÖú¡¿MATLAB±à³Ì

´ó¼ÒºÃ£¬ÕâÊÇÎÒÔÚһƪÎÄÏ×ÉÏ¿´µ½µÄÒ»¸ö¹«Ê½£¬È總¼þDOC2ÖÐËùʾ¡£

      ¸ù¾Ý¾ßÌåµÄʵÑéÖµ£¬ÎÒÃÇ¿ÉÖªµÀµ¡¢µp(ÔÚÕâÀï¾ÍÊǵ±[L] =0ʱµÄµÖµ)ÒÔ¼°[L]¡£

ÎÒ ...

Originally posted by xiegangmai at 2011-04-07 23:21:44:
»¹È±ÉÙÈý¸ö²ÎÊý°É¡£

¿ÉÒÔÓÃmatlab×Ô´øµÄlsqnonlinº¯Êý£¬ÓÃ×îС¶þ³Ë·¨½øÐÐÄâºÏ£¬Ó¦¸ÃÄܼÆËã³ök1£¬k2µÄ¡£

Originally posted by qiancai at 2011-04-09 20:39:17:
²»¹ýÎÒ»¹ÊDz»Ã÷°×£¬ÎÒÁªÏµÄãÁË£¬ÄܼÓÄãÂð£¿

function qiancai
clear
clc
close
x0 = [ 6.5 , 4.44,   14.9,  15.2,  16.1 ];
% x0 = [ 0.5681    0.5681   15.0217   15.9145   16.9757 ];
[ x, resnorm ] = lsqnonlin( @qiancaifun, x0 )
L = [0;
0.5;
1;
2;
4;
8;
20;
30;
50;
100];

mu = [14.9;
15.67;
15.98;
16.16;
16.35;
16.52;
16.9;
16.97;
17.07;
17.06];

plot( L, mu, 'r-' )
muu =  ( x( 3 ) + ( x(1) + x(2) ) * L * x( 4 )  + 2 * x(1) * x(2) * L .^ 2  * x( 5 ) ) ./ ( ( x(1) + x(2) ) * L + 2 * x(1) * x(2) * L .^ 2 + 1 );
hold on
plot( L, muu, 'b-*')

function F = qiancaifun( x )
L = [0;
0.5;
1;
2;
4;
8;
20;
30;
50;
100];

mu = [14.9;
15.67;
15.98;
16.16;
16.35;
16.52;
16.9;
16.97;
17.07;
17.06];

F = mu - ( x( 3 ) + ( x(1) + x(2) ) * L * x( 4 )  + 2 * x(1) * x(2) * L .^ 2  * x( 5 ) ) ./ ( ( x(1) + x(2) ) * L + 2 * x(1) * x(2) * L .^ 2 + 1 );