Сетевые пути рассеяния
Создайте две сети вейвлет-рассеяния, обе для сигнала длиной 500. Во второй сети установите OptimizePath значение для true.
sf = waveletScattering('SignalLength',500); sfOpt = waveletScattering('SignalLength',500,'OptimizePath',true);
Получение информации о путях обеих сетей. Определите общее количество путей рассеяния в обеих сетях.
[spaths,npaths] = paths(sf); [spathsOpt,npathsOpt] = paths(sfOpt); str = sprintf('Paths in default network: %d\nPaths in path-optimized network: %d\n',... sum(npaths),sum(npathsOpt)); fprintf(str)
Paths in default network: 65 Paths in path-optimized network: 52
Обе сети имеют два банка фильтров. Визуализируйте пути рассеяния, которые включают вейвлеты во втором наборе фильтров. Создайте направленный график. Для каждого вейвлет-фильтра, который находится, по меньшей мере, на одном пути, маркируйте соответствующий узел как waveletNumber.filterbank. Для каждого пути подключите соответствующие узлы. Использовать функцию помощника helperPlotScatteringGraph для построения графиков. Постройте графики обеих сетей.
scatGraph = helperPlotScatteringGraph(spaths);
plot(scatGraph)
title({'Scattering Paths',['OptimizePath: ',num2str(sf.OptimizePath)]})
figure
scatGraphOpt = helperPlotScatteringGraph(spathsOpt);
plot(scatGraphOpt)
title({'Scattering Paths',['OptimizePath: ',num2str(sfOpt.OptimizePath)]})
Вспомогательные функции
plotScatteringGraph
function dirGraph = helperPlotScatteringGraph(networkPaths) % This function is intended for use only in this example. It may change or % be removed in a future release. path = networkPaths{3}.path; % set to 0 if want to show the multiple paths between 0 and each % first level node mkunique = 1; if mkunique == 1 f1 = path(:,1:2); c = unique(f1,'rows'); else c = path(:,1:2); end p1 = string(c(:,1)); p2 = string(c(:,2)+.1); p3 = string(path(:,2)+.1); p4 = string(path(:,3)+.2); dirGraph = digraph([p1;p3],[p2;p4]); end