fitnessfuc.asv

遗传算法数值法确定色谱吸附等温线（基于线上求解法）

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

    cal=nonlinear_CrankNicolson(a0,b0,c00,0.9571);
    expe1=[0,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,280,290,300,310,320,330,340,350,360,370,380,390,400,410,420,430,440,450,460,470,480,490,500,510,520,530,540,550,560,570,580,590,600,610,620,630,640,650,660,670,680,690,700,710,720,730,740,750,760,770,780,790,800,810,820,830,840,850,860,870,880,890,900,910,920,930,940,950,960,970,980,990,1000;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.000534176840598,0.001292052229342,0.009558695582553,0.037555435531227,0.073268043557524,0.11334686228032,0.156037152087083,0.200175248378846,0.244915363408717,0.289615812922569,0.333780202782033,0.377022394209147,0.419043377149054,0.459614804530551,0.498566545007118,0.535776872891062,0.571164476120965,0.604681797453532,0.636309399605236,0.66605113910532,0.693930008850698,0.719984544878086,0.744265705406849,0.766834184538088,0.787758079111913,0.807110886781127,0.824969788283221,0.841414193013738,0.856524497748372,0.870381059558642,0.88306333903758,0.894649202342125,0.905214358020739,0.914831921312326,0.923572075615896,0.931501829041275,0.938684844099095,0.945181342619739,0.951048055600699,0.956338226341118,0.961101652235656,0.965384756868917,0.96923068405165,0.972679417978123,0.975767912787121,0.97853023466109,0.980997710194382,0.983199081030562,0.985160659544571,0.986906479300657,0.988458450734531,0.989836508476758,0.991058754496758,0.992141598112262,0.993099886595158,0.99394703159794,0.994695127221388,0.995355060768322,0.995936623497351,0.996448598390002,0.996898860456001,0.997294456141625,0.997641682738055,0.997946160475795,0.998212898350019,0.998446353676834,0.998650489559848,0.998828825042813,0.998984477948339,0.999120213985695,0.999238473916817,0.999341419529568,0.999430955579519,0.999508762180195,0.999576317789704;];
    
    calculate1(:,2)=interp1(cal(:,1),cal(:,2),expe1(:,1),'spline');
    scores1=sum((calculate1(:,2)-expe1(:,2)).^2);


    cal=nonlinear_CrankNicolson(a0,b0,c00,10.3487);
    expe2=[0,1,2,3,4,5,6,7,8,9,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,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,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.00012146308232,0.000545791741958,0.002402906258757,0.009014199657928,0.019890761448298,0.034708550832472,0.05316753862804,0.074964529071284,0.099785987611971,0.127307288644951,0.157194667736044,0.189108455748065,0.22270691729396,0.257650315112043,0.293604959366877,0.330247078376994,0.367266392107221,0.40436930252109,0.441281641655474,0.477750938378734,0.513548182670287,0.548469083556389,0.582334829495492,0.614992373341579,0.646314273097104,0.676198127687536,0.704565655009808,0.731361460376666,0.756551550629548,0.780121645424063,0.802075339511243,0.822432167132104,0.841225615487936,0.858501131929614,0.874314163904645,0.888728264129794,0.901813291911061,0.913643731000029,0.924297143602578,0.933852771459217,0.942390291727464,0.94998873211128,0.956725541952129,0.962675817832192,0.967911674896364,0.972501754809783,0.976510862233904,0.979999715906346,0.983024794418623,0.985638285001981,0.987888095123059,0.989817937518722,0.991467457748316,0.992872413926387,0.994064879646719,0.995073471063999,0.99592358653744,0.996637659802681,0.997235404446935,0.997734064181975,0.998148650555142,0.998492170997324,0.998775850106777,0.999009325809039,0.999200839718999,0.999357411075787,0.999484990385266,0.999588602433156,0.999672471904684,0.999740138374868,0.999794553905322,0.999838169045388,0.999873014968064,0.999900764347213,0.999922790302164,0.999940217611874,0.999953961367128,0.999964765622735,0.999973232386677,0.999979845777731,0.99998499521679,0.999988991854049,0.999992084029878,0.999994468870486,0.99999630195097,0.999997706958362,0.999998780523158,0.999999598017142,1.000000218384918,1.000000688009122,1.000001042643038,1.000001309343202,1.000001509368326,1.000001659145593,1.000001770270662,1.000001852406583,1.000001914250099,1.000001958700126,1.0000019925208;];
    calculate2(:,2)=interp1(cal(:,1),cal(:,2),expe2(:,1),'spline');
    scores2=sum((calculate2(:,2)-expe2(:,2)).^2);


    cal=nonlinear_CrankNicolson(a0,b0,c00,78.69487242);
    expe3=[0,1,2,3,4,5,6,7,8,9,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,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.000001599150644,0.000312442813956,0.007241162514979,0.015753992609432,0.091922966952115,0.191378185625752,0.296974907903816,0.40035079541368,0.49677653570943,0.58370793367804,0.660044757665363,0.725650285215433,0.781009643191617,0.82698062555515,0.864613381172096,0.895023490340543,0.91930653905145,0.93848442630971,0.953475289377075,0.965080381066626,0.973982565070926,0.980752282168222,0.985857872873592,0.989678006452769,0.992514665245143,0.994605677496254,0.996136201329438,0.997248856406197,0.998052404539557,0.998629012680618,0.999040206798662,0.99933167028613,0.999537051702207,0.999680945397986,0.999781195731871,0.999850658174799,0.999898530527391,0.999931350951586,0.999953736646212,0.999968928482024,0.999979187469582,0.999986081817246,0.999990693043641,0.999993762874087;];
    calculate3(:,2)=interp1(cal(:,1),cal(:,2),expe3(:,1),'spline');
    scores3=sum((calculate3(:,2)-expe3(:,2)).^2);


    cal=nonlinear_CrankNicolson(a0,b0,c00,103.12);
    expe4=[0,1,2,3,4,5,6,7,8,9,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,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.010333472633825,0.152198517940264,0.383539356380915,0.64463565079519,0.900782738460046,1.131280361714507,1.32591218871218,1.482026185027153,1.601873243793638,1.690381964701319,1.753494574282389,1.797068573506594,1.826261010473235,1.845272395267649,1.857325376260667,1.864773559930178,1.869264559736229,1.87190924359969,1.873431530256012,1.874288617145074,1.874760941621412,1.875015858223429,1.875150671062839,1.875220565360745,1.875256107447634,1.875273841155935,1.875282526183087,1.875286702870442,1.875288676299457,1.875289591737781,1.87529000969744,1.87529019685803,1.875290279286268,1.875290315166796,1.8752903306827,1.875290337470908,1.87529034038014,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884,1.875290341349884;];
    calculate4(:,2)=interp1(cal(:,1),cal(:,2),expe4(:,1),'spline');
    scores4=sum((calculate4(:,2)-expe4(:,2)).^2);
    
    
    scores(i)=scores1+scores2+scores4+scores6;
end