3Â¥: Originally posted by cooooldog at 2014-09-18 09:34:19
Â¥Ö÷ÄãʲôÒâ˼?
ÎÒ´ð·ÇËùÎÊÂð?

clear all
close all
clc
data=[0.981,0.013,0.291;1.043,1.161,0.091;0.765,-0.123,0.221;0.469,-3.268,0.215;0.381,-4.265,0.208;0.798,0.618,0.224;1.006,-0.645,0.172;0.534,-4.199,0.203;0.95,0.017,0.254;0.525,-3.26,0.196;0.965,0.499,0.079;0.995,0.166,0.114;1.011,-0.146,0.215;0.977,2.121,0.211;0.645,-0.699,0.178;0.539,0.192,0.256;0.719,0.053,0.13;1.125,-0.177,0.200;0.49,-2.974,0.21;1.152,-0.375,0.122;1.007,-0.644,0.227;0.391,-1.051,0.192;0.714,2.019,0.05;0.612,2.222,0.069;0.573,-0.32,0.178;1.053,-0.965,0.226;1.04,0.069,0.364;0.697,0.719,0.401;1.008,0.264,0.487;0.444,-0.233,0.144;0.926,-0.737,0.11;1.08,4.225,0.009;0.583,-2.278,0.241;0.598,3.001,0.351;0.73,1.093,0.158;0.956,0.337,0.074;0.742,0.742,0.1;1.577,-0.577,0.224;1.053,1.704,0.109;0.542,1.215,0.204;0.594,-0.87,0.116;1.16,-0.319,0.233;0.556,2.2,0.21;0.861,-1.008,0.139;0.555,-1.49,0.186;0.378,-0.128,0.186;1.183,-0.164,0.17;0.357,0.042,0.194;0.756,-0.838,0.24;1.304,-0.836,0.083;0.97,-0.337,0.192;1.068,-0.052,0.199;0.485,-1.242,0.219;0.492,-1.479,0.167;0.405,-2.649,0.246;0.63,-0.894,0.104;0.925,0.701,0.5;0.709,1.412,0.064;0.818,0.422,0.138;0.604,-3.035,0.177;0.627,-2.62,0.173;0.526,-1.677,0.264;0.404,-1.063,0.191;0.885,-1.58,0.112;1.224,0.43,0.329;1.135,3.134,0.338;0.701,0.658,0.135;0.317,-1.63,0.275;0.502,-2.067,0.234;0.732,-1.317,0.274;0.494,-1.433,0.156;0.368,0.101,0.165;0.357,-2.017,0.161];
a=data(:,1);
b=data(:,2);
c=data(:,3);
s1=0.6;
s2=0.8;
f1=@(x)mean(c+(1-c)./(1+exp(1.7*a.*(b-x))))-s1
f2=@(x)mean(c+(1-c)./(1+exp(1.7*a.*(b-x))))-s2
[x1,fval,exitflag,output]=fzero(f1,0);
[x2,fval0,exitflag0,output0]=fzero(f2,0);
disp([' s = 0.6, theta =',num2str(x1,'%10.8f')])
disp([' s = 0.8, theta =',num2str(x2,'%10.8f')])

3Â¥: Originally posted by cooooldog at 2014-09-18 09:34:19
Â¥Ö÷ÄãʲôÒâ˼?
ÎÒ´ð·ÇËùÎÊÂð?