comparison toolboxes/MIRtoolbox1.3.2/MIRToolbox/@mirpattern/mirpattern.m @ 0:e9a9cd732c1e tip

first hg version after svn

+function varargout = mirpattern(orig,varargin)
+%   p = mirpattern(a)
+
+        period.key = 'Period';
+        period.type = 'Boolean';
+        period.when = 'After';
+        period.default = 0;
+    option.period = period;
+        
+specif.option = option;
+     
+varargout = mirfunction(@mirpattern,orig,varargin,nargout,specif,@init,@main);
+
+
+function [x type] = init(x,option)
+if not(isamir(x,'mirpattern'))
+    x = mirsimatrix(x);
+end
+type = 'mirpattern';
+
+
+function p = main(orig,option,postoption)
+if not(isamir(orig,'mirpattern'))
+    b = get(orig,'Branch');
+    fp = get(orig,'FramePos');
+    pp = get(orig,'PeakPos');
+    for i = 1:length(b)
+        bi = b{i}{1}{1};
+        pi1 = sort(pp{i}{1}{bi(1,1)});
+        pi2 = sort(pp{i}{1}{bi(end,1)}); 
+        p.pattern{1}.occurrence{1}.start = ...
+            fp{i}{1}(1,bi(1,1)) - mean(fp{i}{1}(1:2,pi1(bi(1,2))));
+        p.pattern{end}.occurrence{2}.start = ...
+            fp{i}{1}(1,bi(1,1));
+        p.pattern{end}.occurrence{1}.end = ...
+            fp{i}{1}(1,bi(end,1)) - mean(fp{i}{1}(1:2,pi2(bi(end,2))));
+        p.pattern{end}.occurrence{2}.end = ...
+            fp{i}{1}(1,bi(end,1));
+    end
+    p = class(p,'mirpattern',mirdata(orig));
+end
+if postoption.period
+    for i = 1:length(p.pattern)
+        poi = p.pattern{i}.occurrence;
+        if poi{1}.end > poi{2}.start
+            poi{1}.end = poi{2}.start;
+            cycle = poi{1}.end - poi{1}.start;
+            ncycles = floor((poi{2}.end-poi{2}.start)/cycle)+2;
+            poi{ncycles}.end = poi{2}.end;
+            poi{2}.end = poi{2}.start + cycle;
+            for j = 2:ncycles-1
+                poi{j}.end = poi{j}.start + cycle;
+                poi{j+1}.start = poi{j}.end;
+            end
+        end
+        p.pattern{i}.occurrence = poi;
+    end
+end