clc
clear all
close all


Original_CV_12 = [35.1143271609051,35.7181802545624,36.0728860426700,36.2891011390115,36.2661002418215,36.2351314682920,36.2965888246979,37.0790134891817,41.3255159565718];
s_12 = [2,3,4,5,6,7,8,9,10,11] ;

Original_CV_24 = [35.1143271609051,35.7182401285409,36.0728249738140,36.2883949133339,36.2614868776259,36.2091677712078,36.1544201857041,36.1091027816732,36.3383604575390,36.3387567251609,36.2630612156123,36.1938100976312,36.1362552937180,36.1817442829706,36.3808381844006,36.3041658486465,36.2343235751926,36.1905213522727,36.2533704695256,37.3324278169106,41.1083382170952];
s_24 = [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23] ;

Original_CV_36 = [35.1143271609051,35.7182401285409,36.0728249738140,36.2883949133339,36.2614868776259,36.2091677712078,36.1544201857041,36.1091029999194,36.3383729010492,36.3387440603568,36.2631328650920,36.1937387218283,36.1363164263031,36.1816829965864,36.3807022294202,36.3033582403434,36.2295957965586,36.1650506820031,36.1146701188477,36.3204593587068,36.3441889894022,36.2678362151368,36.1977708071379,36.1392511311136,36.1608231797204,36.3852385323530,36.3080728077903,36.2346769916751,36.1730245958783,36.1417077470761,36.4821232971275,37.2616770255603,41.2233406435180];
s_36 =  [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];

Original_CV_48 = [35.1143271609051,35.7182401285409,36.0728249738140,36.2883949133339,36.2614868776259,36.2091677712078,36.1544201857041,36.1091029999194,36.3383729010492,36.3387440603568,36.2631328650920,36.1937387218283,36.1363164263031,36.1816829965864,36.3807022294202,36.3033582403434,36.2295957965586,36.1650506820031,36.1146701188477,36.3204593587068,36.3441889894022,36.2678739945695,36.1977810631574,36.1392686603807,36.1607578377804,36.3852979745486,36.3079396618658,36.2338408989068,36.1685566070175,36.1171218627424,36.2997838659511,36.3489980339817,36.2725660258644,36.2018574402932,36.1424234997600,36.1384445762855,36.3901432414005,36.3129715652927,36.2385954137639,36.1731978218362,36.1241806564982,36.3101295197254,36.5196850638789,37.1429829031230,41.3383938981190];
s_48 =  [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];



figure
plot(s_12(1,1:length(s_12)-1) , Original_CV_12 , 'o-m')
hold on
plot(s_24(1,1:length(s_24)-1) , Original_CV_24 , 'o-b')
hold on
plot(s_36(1,1:length(s_36)-1) ,Original_CV_36 , 'o-g')
hold on
plot(s_48(1,1:length(s_48)-1) , Original_CV_48 , 'o-k')


legend('1 cm', '2 cm' , '3 cm' ,'4 cm')

xlim([0 50])
title('CV')
xlabel('node #')
ylabel('CV')
title('original parameters')






%%

fiberDx4_CV_12 = [44.7926401401212,79.0744639211390,84.5565766414119,89.2995507782837,96.3203822577206,108.878717520148,130.141082103786,172.600471572635,263.884085329898];
fiberDx4_CV_24 = [44.4938890897030,77.5185924758998,80.9131042049164,82.0834717990866,82.6654935261160,83.0152605913550,83.2558509938770,83.4459608802068,83.6227191558728,83.8126880687249,84.1205260260764,84.5883587522381,85.3458108060674,86.6289352536422,88.8439195422020,92.7112577098163,99.6672557664323,112.194688492415,134.862725585555,178.355895196515,265.694783601356];
fiberDx4_CV_36 = [44.4938890897030,77.5184754745039,80.9118132005635,82.0803273509163,82.6570124001748,82.9950210199902,83.2142874892818,83.3643418094425,83.4672252935008,83.5229632155042,83.6030208975118,83.6626075732074,83.7085946689461,83.7412124132241,83.7670153860682,83.7863860696661,83.8034808118424,83.8217439250537,83.8464981038695,83.8875594801787,83.9509236509612,84.0778559065781,84.3350394627309,84.7650270729514,85.5010634796052,86.7719924024193,88.9856776699040,92.8601521352243,99.8388281222630,112.389277622706,135.140221376278,178.708577732027,265.897407977757];
fiberDx4_CV_48 = [44.4938890897030,77.5184754745039,80.9118132005635,82.0803273509163,82.6570124001748,82.9950210199902,83.2148054347908,83.3638220007785,83.4672252935008,83.5223690757443,83.6023089963398,83.6624342526624,83.7065761974192,83.7382525461894,83.7609807100302,83.7757955572948,83.7842240451074,83.7876512458049,83.7868311425873,83.7813105385098,83.7629026948665,83.7553397339468,83.7738773846527,83.7883062382875,83.7997244796675,83.8085547688779,83.8160134827152,83.8225312962416,83.8294655229175,83.8401802038977,83.8606887698802,83.8966914254348,83.9554764507596,84.0874701786918,84.3429052059684,84.7726614943375,85.5080522452368,86.7797881564392,88.9935647067439,92.8700292881482,99.8413542692659,112.422685361916,135.133225098301,178.729504824911,266.153636094380];



figure
plot(s_12(1,1:length(s_12)-1) , fiberDx4_CV_12 , 'o-m')
hold on
plot(s_24(1,1:length(s_24)-1) , fiberDx4_CV_24 , 'o-b')
hold on
plot(s_36(1,1:length(s_36)-1) , fiberDx4_CV_36 , 'o-g')
hold on
plot(s_48(1,1:length(s_48)-1) , fiberDx4_CV_48 , 'o-k')


legend('1 cm', '2 cm' , '3 cm' ,'4 cm')

xlim([0 50])
title('CV')
xlabel('node #')
ylabel('CV')
title('fiberDx4')


%%

mean_Original_CV_36 = mean(Original_CV_36(1,10:20))
mean_fiberDx4_CV_36 = mean(fiberDx4_CV_36(1,10:20))