One option is to use the cellfun function (keeping ‘two’ as a cell array). The other option is to use the cell2mat function to convert ‘two’ to a numeric array.
I use cellfun here —
clear;clc
LD = load('3sn_Varying.mat');
u = LD.u;
two = LD.two;
[m,n] = size(two);
MSE = cellfun(@(x)mean((u-x).^2,2), two);
RMSE = sqrt(MSE);
%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Arrange in descending order
%%%%%%%%%%%%%%%%%%%%%%%%%%%
MSE=sort(MSE,'descend');
RMSE=sort(RMSE,'descend');
bestRMSE=min(RMSE)
%%%%%%%%%%%%%%%%%%%%
% Plotting RMSE
%%%%%%%%%%%%%%%%%%%%
figure
x=1:m;
semilogy(x,RMSE,'linewidth',2);
xlabel('\bf Independent Runs')
ylabel('\bf RMSE')
title('\bf RMSE vs Runs');
.
