[交流] 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:
不过我还是不明白，我联系你了，能加你吗？

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 );