run2.m

模式识别工具箱。非常丰富的底层函数和常见的统计识别工具

data = load('5_1_c.txt');
tc = data2tc(data,'WKModel', 1,100);


lower = [1e-8 1e-8 ]
upper = [1e-6 1e-6]
x0 = [1e5 ;1e5];
a0 = [1e-6 1e-6 1e-6 1e-6]

%diffga(fname,VTR,XVmin,XVmax,NP,itermax,F,CR,strategy,refresh, ode_ind, x0, u, y, x, xmask);

[x,f,nf] = diffga('ODEObj',1e-3, lower, upper,60, 20, 0.5, 0.8, 2, 2, 17, x0, tc.x, tc.y, a0, [1 2]);

xhat = [x 1e-6 1e-6];

TestCase = odesim(17, tc.x ,  x0, xhat,200)

yhat = (TestCase.y(:,1) - TestCase.y(:,2))*xhat(1);

figure
plot(yhat);
hold on;
plot(tc.y,'r')

figure
plot(tc.x,'k');
hold on;
plot(tc.y,'r')

param = [1e5 1e5 1e-6 1e-6 1e-6 1e-6]
parambox(1,:) = [1e2 1e2 1e-8 1e-8 1e-8 1e-8]
parambox(2,:) = [1e6 1e6 1e-1 1e-1 1e-1 1e-1]
paramrange = [3 4 5 6];
MethodId = 1;

[x_hat,xr] = WKPARFIT(1,1,'c',1,100,17,param,parambox,paramrange,MethodId);