% PlotLimitCycleWithAllIntialCons % % Script delivers phase space plots with all intial conditions % % Inputs: % PhaseSpace Data % Geneva Mottet UAF 12-22-18 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LoadPhaseSpace load('IntialConditions') load('LimitCycleVector','LCDataHtotal','LCDataPtotal') Place = 400000; UseCondition = 1; %% No Migration Limit Cycle subplot(1,2,1) % For legend plot(LCDataHtotal,LCDataPtotal, 'r', 'LineWidth',3) hold on plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 .75 .25]) hold on UseCondition = 1; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)), 'Color',[0 1 1]) hold on UseCondition = 2; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[1 1 0]) hold on UseCondition = 3; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 1 0]) hold on UseCondition = 4; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 1 .25]) hold on UseCondition = 5; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 .75 .25]) hold on UseCondition = 6; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 .75 .5]) hold on UseCondition = 7; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 .5 0]) hold on UseCondition = 8; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 .25 0]) hold on UseCondition = 9; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[.25 .25 0]) hold on UseCondition = 10; plot(Version4HAMigHp0Pp0(UseCondition,1:Place),abs(Version4PAMigHp0Pp0(UseCondition,1:Place)),'Color',[0 .75 .25]) hold on UseCondition = 1; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 2; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 3; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 4; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 5; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 6; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 7; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 8; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 9; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 10; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 1; plot(LCDataHtotal,LCDataPtotal, 'r', 'LineWidth',3) hold on plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 2; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 3; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 4; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 5; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 6; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 7; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 8; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 9; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 10; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on title('Phase Space for Ten Intial Conditions In Patch A') xlabel('Parasite Population') ylabel('Host Population') legend('Limit Cycle','Intial Conditions','Final Positions', 'Trajectories') set(gca,'TickDir','out'); box off; set(gcf,'color','white'); subplot(1,2,2) UseCondition = 1; BurnIn = 3000; TimeSteps = 5000; MaxPopHost = max(Version4HAMigHp0Pp0(UseCondition, BurnIn:TimeSteps)); MaxPopHostVector = MaxPopHost*ones(TimeSteps); MaxPopPara = max(Version4PAMigHp0Pp0(UseCondition, BurnIn:TimeSteps)); MaxPopParaVector = MaxPopPara*ones(TimeSteps); plot(Version4LCSAHp0Pp0(UseCondition,1:TimeSteps), Version4HAMigHp0Pp0(UseCondition,1:TimeSteps),'Color',[0 0 1]) %plot(Version4LCSAHp0Pp0(UseCondition,1:TimeSteps), Version4HBMigHp0Pp0(UseCondition,1:TimeSteps),'Color',[0 0 .5]) hold on plot(Version4LCSAHp0Pp0(UseCondition,1:TimeSteps), Version4PAMigHp0Pp0(UseCondition,1:TimeSteps),'Color',[1 0 0]) hold on plot(Version4LCSAHp0Pp0(UseCondition,1:TimeSteps), MaxPopParaVector, '--', 'Color', [0 0 0]); hold on plot(Version4LCSAHp0Pp0(UseCondition,1:TimeSteps), MaxPopHostVector, '--', 'Color', [0 0 0]); %plot(Version4LCSAHp0Pp0(UseCondition,1:TimeSteps), Version4PBMigHp0Pp0(UseCondition,1:TimeSteps),'Color',[.5 0 0]) title('Time Series for One Intial Condition for Patch A') xlabel('Time Steps') ylabel('Population') legend('Host','Parasite', 'Population at Limit Cycle') set(gca,'TickDir','out'); box off; set(gcf,'color','white'); %% NO Tragectories figure UseCondition = 1; plot(LCDataHtotal,LCDataPtotal, 'y', 'LineWidth',3) hold on plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 1; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 2; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 3; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 4; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 5; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 6; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 7; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 8; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 9; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 10; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PAMigHp0Pp0(UseCondition,1)), '*r') hold on UseCondition = 1; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 2; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 3; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 4; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 5; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 6; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 7; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 8; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 9; plot(Version4HAMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on UseCondition = 10; plot(Version4HBMigHp0Pp0(UseCondition,1),abs(Version4PBMigHp0Pp0(UseCondition,1)), '*m') hold on plot(LCDataHtotal,LCDataPtotal, 'y', 'LineWidth',3) hold on UseCondition = 1; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 2; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 3; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 4; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 5; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 6; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 7; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 8; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 9; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 10; plot(Version4HAMigHp0Pp0(UseCondition,Place),abs(Version4PAMigHp0Pp0(UseCondition,Place)), 'ob','MarkerSize',6) hold on UseCondition = 1; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 2; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 3; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 4; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 5; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 6; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 7; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 8; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 9; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on UseCondition = 10; plot(Version4HBMigHp0Pp0(UseCondition,Place),abs(Version4PBMigHp0Pp0(UseCondition,Place)), 'og','MarkerSize',6) hold on title('Phase Space for Ten Intial Conditions for Both Ptaches') xlabel('Parasite Population') ylabel('Host Population') legend('Limit Cycle','Patch A Intial Conditions','Patch B Intial Conditions', 'Patch A Final Position', 'Patch A Final Position') set(gca,'TickDir','out'); box off; set(gcf,'color','white');