Զ庯

:

Щ£ܸ,ҪʹԶ庯ʵӦ,ʱɿ·ʽ.


:

1׼

ʹVB.NETдԶ庯Ŀ¼µ"Calculator\DefineFunction\VBNET"ļ,ɲοԶ庯

2

2.1

ֺ֧:IntegralCortezIntegralDiffIntegralGaussIntegralGaussAdaptIntegralInfiniteIntegralInfinite2IntegralRombergIntegralSimpsonIntegralxyIntegralxyz

:

sin(x)*cos(x+5)/(4+x*x)*sin(sin(x)*cos(x+5)/(4+x*x)+x)x[0,1]ϵĶ

(1):"Calculator\DefineFunction\VBNET"½Ϊ"Ftest.vb"ļ,ļ/*

Public Function Ftest(ByVal x as Double) as Double
dim y,z as Double
y=sin(x)*cos(x+5)
z=y/(4+x*x)
Return z*sin(z+x)
End Function

*/

 

(2):ִ

Integral(Ftest,0,1)//سõ½
ans =
[ 0.05562043407933 ]

2.2

ֺ֧:{nlinFit2<\nlinFit2>}{nlinFitDE<\nlinFitDE>}{nlinFitGA<\nlinFitGA>}{nlinFitPSO<\nlinFitPSO>}{nlinFitPSO2<\nlinFitPSO2>}{nlinFitPSO3<\nlinFitPSO3>}

ע:ڱдʱ,Ĳһһ˳ıƶӦ,


:

//֪һ,ݱʽsin(z(1))*a+b*z(2)+c*z(3)-d+z(4)=0,z(1)z(2)z(3)z(4)ֱӦzݵĵ1234,Ҫϳϵabcd

z =
[ 1.00000000000000    0.84147098480789    0.54030230586814   -0.11354116913921
  2.00000000000000    1.81859485365136   -0.83229367309428   -2.98857287563957
  3.00000000000000    0.42336002417960   -2.96997748980134    8.57093735278360
  4.00000000000000   -3.02720998123171   -2.61457448345445    26.5074268849209
  5.00000000000000   -4.79462137331569    1.41831092731613    30.5547959392623
  6.00000000000000   -1.67649298919356    5.76102171990220    10.3420277278667
  7.00000000000000    4.59890619103152    5.27731578040313   -18.9578735381232
  8.00000000000000    7.91486597298705   -1.16400027046891   -27.4422548545962
  9.00000000000000    3.70906636717581   -8.20017235696209   -0.69668553794993
  10.0000000000000   -5.44021110889370   -8.39071529076452    40.8716652807862 ]


(1):"Calculator\DefineFunction\VBNET"½Ϊ"Ftest.vb"ļ,ļ/*

Public Function Ftest(ByVal a as Double,ByRef b as Double,ByRef c as Double ,ByRef d as Double,ByRef z1 as Double,ByRef z2 as Double,ByRef z3 as Double,ByRef z4 as Double) as Double
Return sin(z1)*a+b*z2+c*z3-d+z4
End Function

*/

b=Zeros(10,1);zn="z1,z2,z3,z4";xn="a,b,c,d";nlinFit2(Ftest,xn,zn,b,z)//سõ½
ans =
[ 4.99999999548462    4.00000000015479    0.99999999989246    7.99999999914101    7.4121897549E-17 ]

//,һҪעԶ庯Ftestղ˳ϱ(xn),Ա(zn)


2.3Է

ֺ֧:{nlinSolve<\nlinSolve>}{nlinSolveGA<\nlinSolveGA>}{nlinSolvePSO<\nlinSolvePSO>}{nlinSolvePSO2<\nlinSolvePSO2>}{nlinSolvePSO3<\nlinSolvePSO3>}

ע:ж,һ.

:

·/*

sin(x)*8+4*y^2-z=10
x^2-cos(z)*7=1.235
log(x^2+4)*7-y*z+7=12.78

*/


(1):"Calculator\DefineFunction\VBNET"½Ϊ"Ftest.vb"ļ,ļ/*

Public Function Ftest(ByVal x as Double,ByVal y as Double,ByVal z as Double) as Double
dim f1,f2,f3 as Double
f1=sin(x)*8+4*y^2-z-10
f2=x^2-cos(z)*7-1.235
f3=log(x^2+4)*7-y*z+7-12.78
Return f1^2+f2^2+f3^2
End Function

*/

 

(2):ִ

x = "x,y,z";//ע,xyz뺯Ftest(x,y,z).Ϊ x="x,z,y"뺯ͱFtest(x,z,y),ʱһҪע

nlinSolveGA(Ftest,x)//سõ½
ans =
[ 0.28575236344791   -0.89414919069016   -4.54687684265185    1.5899325674E-16 ]

2.4ֵ֧

{DiffFirstOrder<\DiffFirstOrder>}{DynamicWaveU<\DynamicWaveU>}{DynamicWaveV<\DynamicWaveV>}