how to plot a gragh with 2 x axes (top and bottom) with different lables and single left y axis

Step_size = [5 10 15 20 25 30 35 40]; % bottom x axis

FMC_models = [7650000 3830000 2550000 1920000 1530000 1280000 1100000 960000];
RIA_models = [3689 1853 1232 935 747 624 545 463];
PMA_models = [26 26 26 26 26 26 26 26];
PSA_models = [10000 10000 10000 10000 10000 10000 10000 10000];

%% need to use semilog y axis

Initial_trial_force = [0 500 1000 1500 2000 2500 3000 3500]; % top x axis

FMC_models_1 = [960000 840000 710000 600000 460000 340000 210000 90000];
RIA_models_1 = [463 442 398 350 297 244 171 66];

Answer 1

[Ameer Hamza]

Try this

Step_size = [5 10 15 20 25 30 35 40]; % bottom x axis
FMC_models = [7650000 3830000 2550000 1920000 1530000 1280000 1100000 960000];
RIA_models = [3689 1853 1232 935 747 624 545 463];
PMA_models = [26 26 26 26 26 26 26 26];
PSA_models = [10000 10000 10000 10000 10000 10000 10000 10000];

%% need to use semilog y axis
Initial_trial_force = [0 500 1000 1500 2000 2500 3000 3500]; % top x axis
FMC_models_1 = [960000 840000 710000 600000 460000 340000 210000 90000];
RIA_models_1 = [463 442 398 350 297 244 171 66];

labels1 = {'FMC models', 'RIA models', 'PMA models', 'PSA models'};
labels2 = {'FMC models 1', 'RIA models 1'};

figure();
ax_bottom = axes();
ax_bottom.YScale = 'log';
hold(ax_bottom);
plot(ax_bottom, Step_size, FMC_models, ...
    Step_size, RIA_models, ...
    Step_size, PMA_models, ...
    Step_size, PSA_models);
legend(ax_bottom, labels1, 'Location', 'northeast');

ax_top = axes();
ax_top.XAxisLocation = 'top';
ax_top.YAxis.Visible = 'off';
ax_top.Color = 'none';
ax_top.YScale = 'log';
ax_top.YLim = ax_bottom.YLim;
hold(ax_top);
plot(ax_top, Initial_trial_force, FMC_models_1, ...
    Initial_trial_force, RIA_models_1, 'LineStyle', '--');

legend(ax_top, labels2, 'Location', 'northwest', 'Color', 'none');

[Akshay Kumar]

Thank you Ameer Hamza .....

here is a problem while exporting the figure with axes label, two top axes appeaers. Futher i need to add box around. help if possible.

thankyou

[Ameer Hamza]

Can you show me the code which you used to add the xlabels? You can also add this line after plotting on both axes

linkprop([ax_bottom, ax_top],{'Units','Position','ActivePositionProperty'});

[Akshay Kumar]

Dear Ameer Hamza, code os as below

Step_size = [5 10 15 20 25 30 35 40]; % bottom x axis
FMC_models = [7650000 3830000 2550000 1920000 1530000 1280000 1100000 960000];
RIA_models = [3689 1853 1232 935 747 624 545 463];
PMA_models = [26 26 26 26 26 26 26 26];
PSA_models = [10000 10000 10000 10000 10000 10000 10000 10000];

%% need to use semilog y axis
Initial_trial_force = [0 500 1000 1500 2000 2500 3000 3500]; % top x axis
FMC_models_1 = [960000 840000 710000 600000 460000 340000 210000 90000];
RIA_models_1 = [463 442 398 350 297 244 171 66];

labels1 = {'FMC', 'RIA', 'PMA', 'PSA'};
labels2 = {'FMC', 'RIA'};

figure();
ax_bottom = axes();
ax_bottom.YScale = 'log';
hold(ax_bottom);
plot(ax_bottom, Step_size, FMC_models,'-k>','linewidth',1.3,'markersize',4,'MarkerEdgeColor','k','MarkerFaceColor','k')
set(gca, 'XDir','reverse')
hold on
plot(ax_bottom, Step_size, RIA_models,'-bo','linewidth',1.3,'markersize',4,'MarkerEdgeColor','b','MarkerFaceColor','b')
plot(ax_bottom, Step_size, PMA_models,'-rs','linewidth',1.3,'markersize',4,'MarkerEdgeColor','r','MarkerFaceColor','r')
plot(ax_bottom, Step_size, PSA_models,'-gd','linewidth',1.3,'markersize',4,'MarkerEdgeColor','g','MarkerFaceColor','g')

ylim([0 10^8])

xlabel(ax_bottom,'Step size (kN)','FontWeight','Bold')
ylabel(ax_bottom,'No of models','FontWeight','Bold')
legend(ax_bottom, labels1, 'Location', 'northeast');
box on

ax_top = axes();
ax_top.XAxisLocation = 'top';
ax_top.YAxis.Visible = 'off';
ax_top.Color = 'none';
ax_top.YScale = 'log';
ax_top.YLim = ax_bottom.YLim;
hold(ax_top);
plot(ax_top, Initial_trial_force, FMC_models_1, '--k>','linewidth',1.3,'markersize',4,'MarkerEdgeColor','k','MarkerFaceColor','k')
plot(ax_top, Initial_trial_force, RIA_models_1, '-bo','linewidth',1.3,'markersize',4,'MarkerEdgeColor','b','MarkerFaceColor','b')


xlabel(ax_top,'Initial Trial T (kN)','FontWeight','Bold')
legend(ax_top, labels2, 'Location', 'northwest', 'Color', 'none');

[Ameer Hamza]

Are you using an older version of MATLAB? I don't see this issue in R2020a; however, I remember that this used to happen in R2019b. Adding the following line should fix the problem in R2019b.

linkprop([ax_bottom, ax_top],{'Units','Position','ActivePositionProperty'});

[Akshay Kumar]

Yes I am using R2019b and it works.....

Thank you Genius.....

Ameer Hamza

[Ameer Hamza]

I am glad to be of help.