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

first hg version after svn

+function c = mirtimes(a,b)
+
+d = get(a,'Data');
+e = get(b,'Data');
+f = cell(1,length(d));
+p = get(a,'Pos');
+q = get(b,'Pos');
+r = cell(1,length(d));
+l = get(a,'Name');
+m = get(b,'Name');
+n = cell(1,length(l));
+for i = 1:length(d)
+    f{i} = cell(1,length(d{i}));
+    for j = 1:length(d{i})
+        pij = p{i}{j}(:,1,1);
+        qij = q{i}{j}(:,1,1);
+        pq = pij; %union(pij,qij);
+            % The sampling of the product is the sampling of the first
+            % operand.
+        dij = d{i}{j}; %interp1(pij,d{i}{j},pq);
+        eij = interp1(qij,e{i}{j},pq);
+        dij = max(dij,0);
+        eij = max(eij,0);
+        eij = eij./repmat(max(eij)+1e-16,[size(eij,1),1,1,1]);
+            % The second operand of mirtimes is scaled from 0 to 1.
+            % In this way, the range of value of the first operand is kept.
+        f{i}{j} = dij.*eij;
+        [x y] = find(isnan(f{i}{j}));
+        x = unique(x);
+        pq(x) = [];
+        f{i}{j}(x,:) = [];
+        r{i}{j} = repmat(pq,[1 size(p{i}{j},2) size(p{i}{j},3)]);
+    end
+    if strcmpi(l{i},m{i})
+        n{i} = l{i};
+    else
+        n{i} = [l{i} '*' m{i}];
+    end
+end
+c = set(a,'Pos',r,'Data',f,'Name',n,...
+          'Title',[get(a,'Title') ' * ' get(b,'Title')]);