Построение графиков с двумя независимыми осями в Matlab

1. Отсутствие аналога hold on («родной» hold on работает не совсем корректно с plotyy). Для того, чтобы добавить более, чем 2 кривые необходимо использовать вот такую конструкцию:

[hAx,hLine1,hLine2] = plotyy([x1',x2'],[y1' y2'],[x3',x4'],[y3',y4']); %Добавляет 4 кривые

Из этой конструкции вытекает неприятная особенность №2:

2. Размерности массивов, заключенных в квадратные скобки должны совпадать, т.к. из них формируются матрицы элементов. На практике такое бывает очень не часто.

3. Оформление серьезно страдает оттого, что нельзя программными методами изменить цвета и типы всех линий подряд, можно форматировать только набор линий, относящихся к конкретной оси (hLine1 и hLine2) — во всяком случае, я не смог. При этом, я не говорю сейчас об изменении параметров руками, т.е. редактированием в окне «figure» — только непосредственно кодом в .m-файле.

Резюмируя вышесказанное: plotyy() не очень хорошо подходит для отображения нескольких наборов графиков для разных осей. Разве что для простеньких зависимостей типа этих:

x = linspace(0,10);
y1 = 200*exp(-0.05*x).*sin(x);
y2 = 0.8*exp(-0.5*x).*cos(10*x);
y3 = 0.2*exp(-0.5*x).*cos(10*x);
y4 = 150*exp(-0.05*x).*sin(x);

X=[x',x'];
Y1=[y1' y4'];
Y2=[y2',y3'];

[hAx,hLine1,hLine2] = plotyy(X,Y1,X,Y2);

% График во весь экран
h = figure('Units', 'normalized', 'OuterPosition', [0 0 1 1]); 

% Настройки шрифта:
F='Times New Roman';
FN='FontName';
FS='FontSize';
l=30;% Размер шрифта
set(gca, FN, F, FS, l)
box on  % Рамка

hold on

%% ПРАВАЯ ОСЬ
yyaxis right

hPlot_1 = plot(time*10^6, T);             % Прямой численный метод
hPlot_2 =plot(time_exp_T*10^6, T_exp);    % Экспериментальное значение
hPlot_3 =plot(time_pk, T_pk);             % Моментный метод

set( hPlot_1, 'LineWidth', 3, 'LineStyle', ':', 'Color',  'k' );
set( hPlot_2, 'LineWidth', 1, 'LineStyle', 'none', 'Color', 'k', 'Marker', 'o' );
set( hPlot_3, 'LineWidth', 3, 'LineStyle', '--', 'Color', 'k' );

% Пределы по у (справа)
ylim([180 270]);

% Подписи осей по у (справа)
yticks([190 200 210 220 230 240 250 260 267])
yticklabels({'190', '200', '210', '220', '230', '240','250' '260', 'it T, K'})

% Цвет осей
set(gca,'xcolor','k'); 
set(gca,'ycolor','k');

%% ЛЕВАЯ ОСЬ
yyaxis left

hPlot_10 = plot(time_P*10^6, P_p0);            % Прямой численный метод
hPlot_11 = plot(time_exp_P_p0*10^6, P_p0_exp); % Экспериментальное значение
hPlot_13 = plot(time_pk, Pp0_pk);              % Моментный метод

set( hPlot_10, 'LineWidth', 3, 'LineStyle', ':', 'Color',  'k' );
set( hPlot_11, 'LineWidth', 1, 'LineStyle', 'none', 'Color', 'k', 'Marker', 'x' );
set( hPlot_13, 'LineWidth', 3, 'LineStyle', '--', 'Color', 'k' );

% Пределы по у (слева)
ylim([0.1 0.8]);

% Подписи осей по у (слева)
yticks([0.2 0.3 0.4 0.5 0.6 0.7 0.77])
yticklabels({'0.2', '0.3', '0.4', '0.5', '0.6', '0.7', 'it p/p_0 '})

% Цвет осей
set(gca,'xcolor','k'); 
set(gca,'ycolor','k');
%% ОФОРМЛЕНИЕ ГРАФИКОВ

% Маркеры и подписи оси Х
xticks([0 20 40 60 80 100 120 140 160 185])
xticklabels({'0' '20' '40' '60' '80' '100' '120' '140' '160', 'it X/U, мc'})

% СТРЕЛКИ
text(20,0.3,'leftarrow', FS, 60, FN, F);
text(140,0.6,'rightarrow', FS, 60, FN, F);
 
text(160,0.24,'1', FS, l, FN, F);
text(60, 0.25,'2', FS, l, FN, F);
text(120, 0.2,'3', FS, l, FN, F);
 
text(140,0.51,'4', FS, l, FN, F);
text(80, 0.52,'5', FS, l, FN, F);
text(120,0.48,'6', FS, l, FN, F);

text(5,0.75,'it D_2O - N_2', FS, l, FN, F);

hold off