0

波高と波長の最大値を求める

dwcL01ただし、数値波動モデルからの出力dwcL10の配列です。<3001x2 double>

私のスクリプトの一部:

%% Plotting results from SWASH
% Examination of phase velocity on deep water with different number of layers
% Wave height 3 meters, wave peroid 8 sec on a depth of 30 meters
clear all; close all; clc;
T=8;
L0=1.56*T^2;

%% Loading results tabels.
load dwcL01.tbl; load dwcL02.tbl; load dwcL03.tbl; load dwcL04.tbl;
load dwcL05.tbl; load dwcL06.tbl; load dwcL07.tbl; load dwcL08.tbl;
load dwcL09.tbl; load dwcL10.tbl;
M(:,:,1) = dwcL01; M(:,:,2) = dwcL02; M(:,:,3) = dwcL03; M(:,:,4) = dwcL04;
M(:,:,5) = dwcL05; M(:,:,6) = dwcL06; M(:,:,7) = dwcL07; M(:,:,8) = dwcL08;
M(:,:,9) = dwcL09; M(:,:,10) = dwcL10;

%% Finding position of wave crest using diff and sign.
for ii=1:10
    Tp(:,1,ii) = diff(sign(diff([M(1,2,ii);M(:,2,ii)]))) < 0;
    Wc(:,:,ii) = M(Tp,1,ii);
    L(:,ii) = diff(Wc(:,1,ii))
end

ループ

for ii=1:10
    Tp(:,1,ii) = diff(sign(diff([M(1,2,ii);M(:,2,ii)]))) < 0;
    Wc(:,:,ii) = M(Tp,1,ii);
    L(:,ii) = diff(Wc(:,1,ii))
end

ii = 1次のエラーを取得するために正常に動作しますii = 2

Index exceeds matrix dimensions.

Error in mkPlot (line 19)
   Wc(:,:,i) = M(Tp,:,i);

異なるセットアップで同じ数の波頭を持たないでください。当然、M(Tp,1,ii)異なる寸法になります。この問題を回避するにはどうすればよいですか? for ループで実行できますか? お気軽に私に電子メールを送信するか、他の方法で詳細を尋ねてください。

4

1 に答える 1