ちょっと美しくないですがやってみました。結果としては・・
こんな感じ。
散布データを内挿して等高線図を描くとどうも綺麗に描けなかったので、散布データを三角形の集合として取り扱い、それぞれの三角形に色を付ける方法を取りました。コードのほとんどが散布データを無理やり三角形の集合として（強度も含め）まとめなおすことに費やされています。
もしもともと FEM の結果であれば各要素の情報（座標点・連結・強度）として出力していればもう少し処理は楽かもしれません。
さらに言うと、もう PDE Toolbox を使って強度計算しちゃえばもっと楽です。
https://jp.mathworks.com/products/pde.html
以下コードです。
頂戴したデータはさらにデータ点数削減して、mat ファイルに纏めました（添付しています）。
Excelファイルからの読み込み方はとりあえずおいておいて、まずは描画に集中します。
% 各点で三角形を作り patch オブジェクトとして色を付ける方法を試します。
% データロード
% X : 1001x1 double
% Y : 1001x327 double
% Strength : 1001x327 double
load data_small.mat
% X座標点をYと同じサイズの行列にします。
XX = repmat(X,1,327);
% FEM の結果(?)ですし、各点からドローネ三角形分割します。
DT = delaunayTriangulation(XX(:),Y(:));
% ドローネ三角形分割の情報と Strength （色）情報を
% 結合するために table 型変数を活用します。
dataTable = table([XX(:),Y(:)],Strength(:),...
'VariableNames',{'Points','Strength'});
% 重複している点があるので削除します。
dataTable = unique(dataTable);
% ドローネ三角形分割の座標点情報
triTable = table(DT.Points,'VariableNames',{'Points'});
% ドローネ三角形分割の座標点と一致する位置にある
% 元データの Strength を検索して結合
triColorTable = join(triTable,dataTable,'Keys','Points');
% それぞれの三角形を patch オブジェクトして描画
% 各三角形の頂点に色を指定し、patch の面は各頂点の色
% から内挿するよう設定（'interp')
patch('Faces',DT.ConnectivityList,...
'Vertices',DT.Points,...
'EdgeColor','none','FaceColor','interp',...
'CData',triColorTable.Strength);
