Итак -
рабочий скрипт:
clear();
figure
[X,Y] = meshgrid(0:0.1:1,1:15); % values of x and y at which we want to evaluate
for x = 1:15
for y = 1:11
% Z(x,y) = ppscq(X(x,y), Y(x,y));
Z(x,y) = ppscq(X(x,y), Y(x,y));
end
end
surf(X,Y,Z);
обёртка для предпериода(здесь мы генерируем последовательности ,чтобы затем определить их предпериод):