fitnessfuc.m

有限差分法（crank_nicholson finite difference ) 求解色谱模型

function scores = fitnessfuc(pop)

scores= zeros(size(pop,1),1); %Pre-allocate y
for i = 1:size(pop,1)
    a0=pop(i,1);
    b0=pop(i,2);
    c00=pop(i,3);
    d00=pop(i,4);
    cal=nonlinear_CrankNicolson(a0,b0,c00,d00,8.11);
    expe1=[10,20,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76;0,0.0015,0.003375,0.00725,0.021,0.06475,0.1285,0.231,0.3285,0.41225,0.506,0.596,0.68975,0.7335,0.8135,0.84225,0.86725,0.8885,0.921,0.93225,0.936,0.95725,0.946,0.96475,0.9685,1.001,0.961,0.9785,0.96475,0.97225,0.9685,0.9685,0.9935,0.9935,0.96475,0.9935,0.9785,0.99225,0.9935,0.9685,0.96475,0.9685,0.99975,0.98225,1.001,0.9685,0.9785,1.001,1.001;]';
    calculate1(:,2)=interp1(cal(:,1),cal(:,2),expe1(:,1),'spline');
    scores1=sum((calculate1(:,2)-expe1(:,2)).^2);


    cal=nonlinear_CrankNicolson(a0,b0,c00,d00,10.41);
    expe2=[10,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75;0,0,0.001536983669549,0.001248799231508,0.000480307396734,0,-0.002497598463016,-0.000384245917387,0.000480307396734,0.004226705091258,0.028242074927954,0.10028818443804,0.206916426512968,0.33852065321806,0.424975984630163,0.551777137367915,0.643996157540826,0.728530259365994,0.790009606147935,0.832276657060519,0.878386167146974,0.889913544668588,0.912968299711816,0.939865513928915,0.962920268972142,0.970605187319885,0.966762728146014,0.966762728146014,0.962920268972142,0.970605187319885,0.98597502401537,0.993659942363112,0.98597502401537,0.989817483189241,0.999423631123919,0.993659942363112,0.993659942363112,1.001344860710855,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.98597502401537,0.997502401536984,0.997502401536984,0.997502401536984,0.997502401536984,0.997502401536984,0.997502401536984;]';
    calculate2(:,2)=interp1(cal(:,1),cal(:,2),expe2(:,1),'spline');
    scores2=sum((calculate2(:,2)-expe2(:,2)).^2);


    cal=nonlinear_CrankNicolson(a0,b0,c00,d00,19.62);
    expe3=[10,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75;0,0.000101936799185,0.000101936799185,0.000713557594292,0.000713557594292,0.000713557594292,0.003363914373089,0.039245667686035,0.140163098878695,0.281855249745158,0.399082568807339,0.525484199796126,0.649847094801223,0.749745158002039,0.790519877675841,0.853720693170234,0.888379204892966,0.892456676860347,0.914882772680938,0.931192660550459,0.94750254841998,0.943425076452599,0.95769622833843,0.9638124362895,0.96585117227319,0.95973496432212,0.992354740061162,0.976044852191641,0.982161060142711,0.980122324159021,0.974006116207951,0.95973496432212,0.984199796126402,0.986238532110092,0.987257900101937,0.985728848114169,0.985219164118246,0.967889908256881,0.990316004077472,0.976044852191641,0.986238532110092,0.978083588175331,0.985219164118246,0.986238532110092,0.988277268093782,0.971967380224261,0.974006116207951,0.983690112130479,0.984709480122324,0.998470948012232,0.992354740061162,0.991335372069317,0.988277268093782,0.986238532110092,0.993374108053007,0.990316004077472,0.997451580020387,0.996432212028542,0.996432212028542,0.983180428134556,0.985219164118246,0.996432212028542;]';
    calculate3(:,2)=interp1(cal(:,1),cal(:,2),expe3(:,1),'spline');
    scores3=sum((calculate3(:,2)-expe3(:,2)).^2);


    cal=nonlinear_CrankNicolson(a0,b0,c00,d00,29.6);
    expe4=[10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51;0,-0.000033783783784,0,0.000067567567568,0.000067567567568,0,0.000844594594595,0.001756756756757,0.001993243243243,0.026993243243243,0.121925675675676,0.273277027027027,0.42597972972973,0.562466216216216,0.677331081081081,0.782736486486486,0.82597972972973,0.858412162162162,0.920574324324324,0.934087837837838,0.939493243243243,0.948952702702703,0.973277027027027,0.971925675675676,0.970574324324324,0.988141891891892,0.98679054054054,0.981385135135135,0.990844594594595,0.998952702702703,0.988141891891892,1.007060810810811,1.000304054054054,0.990844594594595,0.990844594594595,0.996925675675676,0.997601351351351,0.994898648648649,0.996587837837838,0.997601351351351,0.999290540540541,0.997939189189189;]';
    calculate4(:,2)=interp1(cal(:,1),cal(:,2),expe4(:,1),'spline');
    scores4=sum((calculate4(:,2)-expe4(:,2)).^2);
    
    cal=nonlinear_CrankNicolson(a0,b0,c00,d00,49.3);
    expe6=[15,0.003488844;16,0.061561866;17,0.20801217;18,0.366227181;19,0.514705882;20,0.636409736;21,0.730527383;22,0.741075051;23,0.794624746;24,0.90010142;25,0.90821501;26,0.935801217;27,0.924442191;28,0.937423935;29,0.973123732;30,0.971501014;31,1.005578093;32,0.997464503;33,0.981237323;34,0.973123732;35,0.989350913;36,0.989350913;37,1.012068966;38,0.981237323;39,0.973123732;40,0.977991886;41,1.008823529;42,0.987728195;43,1.000709939;44,0.976369168;45,0.973123732;46,1.004563895;47,1.005578093;48,1.005578093;49,1.005578093;];
    calculate6(:,2)=interp1(cal(:,1),cal(:,2),expe6(:,1),'spline');
    scores6=sum((calculate6(:,2)-expe6(:,2)).^2);
    
    scores(i)=scores1+scores3+scores4+scores6;
 
end