f = figure;
m = 2;
n = 3;
tiles = tiledlayout(f, m, n);
for i = 1:m  %rows
    for j = 1:n  %cols
        ax(i,j) = nexttile;
        plot(linspace(0,1,10)*i, linspace(2,1,10)*j);
        linkaxes(ax(:,j), 'y');
        setappdata(ax(i,j), 'XLim_listeners', linkprop(ax(i,:),'XLim')); %using linkaxes(ax(i, :), 'x'); 
    end % for all cols
end % for all rows

@Mark, the zoom function calls resetplotview (undocumented). It, in turn, resets all modes (X/YLimMode) to auto. That is why the listeners do not trigger.

To solve the issue, link these properties too:

linkprop([a1 a2],{'XLimMode','YLimMode'}); % trigger linking on double click

x = 1:5;
y = randn(1,5);

figure;
a1 = subplot(1,2,1);
p1 = plot(x,y);

a2 = subplot(1,2,2);
p2 = plot(y,x);

addlistener(a1,'XLim','PostSet',@(~,~)set(a2,'YLim',get(a1,'XLim'))); % link a1.XLim to a2.YLim
addlistener(a1,'YLim','PostSet',@(~,~)set(a2,'XLim',get(a1,'YLim'))); % link a1.YLim to a2.XLim
% and vice versa
addlistener(a2,'XLim','PostSet',@(~,~)set(a1,'YLim',get(a2,'XLim'))); % link a1.XLim to a2.YLim
addlistener(a2,'YLim','PostSet',@(~,~)set(a1,'XLim',get(a2,'YLim'))); % link a1.YLim to a2.XLim

%Use case, Zoom about 4 Times in and zoom out per doubleclick an you will see, that it then doesn't update the axes

@Yaroslav – you are correct, I was misled by the wording of linkprop‘s documentation:

“…maintain the same values for the corresponding properties of different graphics objects.”

Serves me right to trust the builtin documentation! – I should have known better… 🙂

@Yair — I am afraid that your solution doesn’t answer @Amin’s question. It, actually, behaves almost exactly as

linkaxes([hAxes1 hAxes2], 'xy');

@Amin — If you want to link two different properties, you must specify the listeners explicitly. For example,

addlistener(hAxes1,'XLim','PostSet',@(~,~)set(hAxes2,'YLim',get(hAxes1,'XLim'))); % link hAxes2.YLim to hAxes1.XLim
addlistener(hAxes1,'YLim','PostSet',@(~,~)set(hAxes2,'XLim',get(hAxes1,'YLim'))); % link hAxes2.XLim to hAxes1.YLim
addlistener(hAxes2,'XLim','PostSet',@(~,~)set(hAxes1,'YLim',get(hAxes2,'XLim'))); % link hAxes1.YLim to hAxes2.XLim
addlistener(hAxes2,'YLim','PostSet',@(~,~)set(hAxes1,'XLim',get(hAxes2,'YLim'))); % link hAxes1.XLim to hAxes2.YLim

The linkprop function won’t do, since it links the same property between two different handle objects.

@Amin – I believe that you can link multiple properties/handles for example:

linkprop([hAxes1,hAxes2], {'XLim','YLim'})

aught to link hAxes1.XLim to hAxes2.YLim

One question though. Is it possible to link two different properties? For example ‘x-axis’ of an ax to ‘x-axis’ and ‘y-axis’ of another ax?

The solution you proposed with linkprop has a small issue: it does not restore the linking after saving and reloading. On the other hand, linkaxes recovers its linking correctly.

A simple workaround is to use Matlab’s own solution: “MCOS graphics cannot rely on custom machinery in hgload to restore linkaxes. Instead, create a matlab.graphics.internal.LinkAxes to wrap the linkprop which will restore the linkaxes when it is de-serialized.” In simple words,

% Link the relevant axes
restorable_linkprop = @(varargin)matlab.graphics.internal.LinkAxes(linkprop(varargin{:}));
%
setappdata(ax(1,1), 'YLim_listeners', restorable_linkprop(ax(1,:),'YLim'));
setappdata(ax(2,1), 'YLim_listeners', restorable_linkprop(ax(2,:),'YLim'));
setappdata(ax(1,1), 'XLim_listeners', restorable_linkprop(ax(:,1),'XLim'));
setappdata(ax(1,2), 'XLim_listeners', restorable_linkprop(ax(:,2),'XLim'));