Two x-axis (Arrhenius plots)

16 次查看(过去 30 天)
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.

采纳的回答

Chunru
Chunru 2022-1-14
编辑:Chunru 2022-1-15
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;
  5 个评论
Chunru
Chunru 2022-1-21
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
Walter Roberson
Walter Roberson 2023-1-2
Dustin Mangus comments to @Chunru about the Answer:
The right side of the upper x-axis is on a moving scale based on the TC array that you define. It is adjustable and seem to only be accurate at the y-axis. This means that for all higher temperatures (further to the left) are inaccurate unless TC is adjusted precisely. I give this code a 3/10.

请先登录,再进行评论。

更多回答(1 个)

ASC
ASC 2022-1-21
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.

类别

Help CenterFile Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by