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

0

问题是这Tp是一个三维数组。我需要调用Tp(:,:,ii)与当前场景对应的。连同这个并定义Wc为一个单元，我解决了我的问题。

for ii = 1:10
    Tp(:,1,ii) = diff(sign(diff([M(1,2,ii);M(:,2,ii)]))) < 0;
    Wc{:,:,ii} = M(Tp(:,:,ii),1,ii);
    L{:,ii} = diff(cell2mat(Wc(ii)));
end
于 2013-10-22T17:44:37.580 回答