Two x-axis (Arrhenius plots)

16 views (last 30 days)
ASC
ASC on 13 Jan 2022
Answered: ASC on 21 Jan 2022
I would like to create an Arrhenius plot. Arrhenius plots are made by plotting k versus 1/T. I have poached a photo of a “typical” Arrhenius plot from the web (see image). Because 1/T is an inconvenient unit it is typical to the temperature (T) on the top x axis.
How do I put in a second (top) x-axis? Note the non-linaear scale.
My variables are k, T and T_reciprocal.

Accepted Answer

Chunru
Chunru on 14 Jan 2022
Edited: Chunru on 15 Jan 2022
Not sure your equation. But you can modify the following to produce the plot:
x = linspace(2.5, 4, 20); % x corrdinate
y = 20.2 -1589.6 * x; % your equation
ax1 = axes;
plot(ax1, x, y);
ax1.Box = 'off'; % update: turn off extra ticks
% create a second axis
ax2 = axes('Position', ax1.Position, 'Color', 'none', 'YTick', [], 'XAxisLocation', 'top', 'YTick', [], ...
'XDir', 'reverse', 'XLim', 1000./flip(ax1.XLim)-273);
TC = (-20:20:120); % ticks at T deg C
xTC = 1000./(TC + 273); % Corresponding x
% map xTc onto TC
xTCmap = (xTC - ax1.XLim(1))*(-diff(ax2.XLim)/diff(ax1.XLim)) + ax2.XLim(2)
xTCmap = 1×8
-18.2569 10.6996 35.7031 57.5112 76.6997 93.7139 108.9035 122.5471
ax2.XTick = xTCmap;
ax2.XTickLabel = TC;
ax2.YAxisLocation = 'right'; % Update: move yaxis of ax2 to right
ax2.YTick = [];
xlabel(ax1, "1000/T (K^{-1})")
xlabel(ax2, "Temperature (^\circ C)");
ax1.Position(4) = ax1.Position(4) - 0.1; % Update: make space on top
ax2.Position(4) = ax2.Position(4) - 0.1;
  4 Comments
Chunru
Chunru on 21 Jan 2022
Can you post your code and make it executable so that I can run here and see what is problem?
You need the following to turn off the extra ticks. When box is on, the x-axis ticks will show up on both bottom and top.
ax1.Box = 'off'; % update: turn off extra ticks

Sign in to comment.

More Answers (1)

ASC
ASC on 21 Jan 2022
I missed this addtional line in your code. Adding
ax1.Box = 'off'; % update: turn off extra ticks
solved the problem. Thank you very much for you help.

Products


Release

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!