Question

查找波高和波长的最大值

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循环中完成吗？请随时给我发电子邮件或其他明智的人询问更多信息。

Answer 1

问题是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

尺寸问题

1 个答案: