①および②を検出するプログラムを作成してみました．以下順を追って説明してみます．
まず添付していただいたファイルをテーブルとして読み込みます．
%% ファイル読み込み
filename = 'sample file.xlsx'; % ファイル名
T = readtable(filename); % ファイル読み込み
T.Properties.VariableNames{1} = 'Frame'; % １列目の変数名を変更
エクセルの第一列目の変数名だけ’Frame’に変更しています．（＃が変数名に使えないので）
次に，質問にある①ですが，０からの立ち上がりが２回あるので，両方を検出して二番目の立ち上がりの行数を求めることを考えました．また立ち上がりはデータ(GRF)が10以上になるところを検出したいということでしたので，まず10をしきい値として二値化してみます．
%% 2回目の立ち上がり位置のインデックス取得
mask = T.GRF>=10; % しきい値10で2値化
立ち上がり位置の検出をどのように行うか迷ったのですが，ここでは二値化した矩形波(mask)の勾配を計算し，関数findpeaksでピークを検出することで求めました．
（find関数を用いれば，idx = find(T.GRF>=10,1,'first') というようにして最初にデータが10以上になるインデックスを取得することができますが，２回目にデータが10以上になるインデックスを取得する方法がわかりませんでした．他によい方法を御存知の方がおられましたら教えて下さい）
dmask = gradient(mask); % 勾配を計算
[~,idx] = findpeaks(dmask); % 立ち上がり位置検出
idx = idx+1; % 立ち上がり位置は+1進んだところ
条件を満たす立ち上がり位置はfindpeaksによって検出されたピーク位置に＋１したところです．また２回目の立ち上がり位置はidx(2)に対応します．
図で表すと以下のようになります．
figure;
subplot(3,1,1);
plot(T.Frame,T.GRF,'-r',T.Frame(idx(2)),T.GRF(idx(2)),'*b');
ylabel(T.Properties.VariableNames{3});
subplot(3,1,2);
plot(T.Frame,mask,'-r',T.Frame(idx(2)),mask(idx(2)),'*b');
ylabel('binarization')
subplot(3,1,3);
plot(T.Frame,dmask,'-r',T.Frame(idx(2)),dmask(idx(2)),'*b');
ylabel('gradient')
xlabel(T.Properties.VariableNames{1});
青いところが二回目の立ち上がり位置（10を超える）です．対応するFrameは
fprintf('2回目の波が現れる始まりのFrame(具体的には0から10以上になった瞬間のFrame)は%dです\n',T.Frame(idx(2)));
で確認してください．
次に②についてですが，こちらはfindpeaks関数で負のピーク位置を求めると良いと思います．
%% ２回目に最小になるインデックス取得
peakTH = -0.8; % ピーク検出用のしきい値
[~,minidx] = findpeaks(-T.VerticalPstn,'MinPeakHeight',peakTH); % 2つの負のピーク検出
figure;
plot(T.Frame,T.VerticalPstn,'-r');
hold on
plot(T.Frame(minidx(2)),T.VerticalPstn(minidx(2)),'*b');
xlabel(T.Properties.VariableNames{1});
ylabel(T.Properties.VariableNames{2});
具体的には，データの正負を反転させて，振幅のしきい値を指定して求めています．
２回目に最小値を取るインデックスはminidx(2)に対応します．
図で表すと以下のようになります．
青の＊が２回目に最小になる位置です．
対応するFrameは
fprintf('VerticalPstnの2回目の波が来たときの最小値が現れるFrameは%dです\n',T.Frame(minidx(2)));
で確認できます．
以上参考になりましたでしょうか？一つの方法として参考にしていただけましたら幸いです．
