From 063f74dea989afa8da89696faa00ed06b9271f93 Mon Sep 17 00:00:00 2001
From: marwan <>
Date: Tue, 24 Oct 2006 14:36:06 +0000
Subject: [PATCH] minor bug: different lengths of vectors z0 and z1
---
crqad_big.m | 15 +++++++++++----
1 file changed, 11 insertions(+), 4 deletions(-)
diff --git a/crqad_big.m b/crqad_big.m
index a4b081a..914210a 100644
--- a/crqad_big.m
+++ b/crqad_big.m
@@ -51,9 +51,9 @@ function out=crqad_big(varargin)
% Y.Lp = Lp
% Y.Lm = Lm
%
-% Examples: a = sin(0:.1:80) + randn(1,801);
-% b = sin(0:.1:80) + randn(1,801);
-% crqad_big(a,b,3,15,.1,100,'fan')
+% Examples: a = sin(0:.1:800) + randn(1,8001);
+% b = sin(0:.1:800) + randn(1,8001);
+% crqad_big(a,b,3,15,.1,100,'euc')
%
% See also CRQA, CRQAD, CRP, CRP2, CRP_BIG, DL, TT.
%
@@ -71,6 +71,9 @@ function out=crqad_big(varargin)
% $Revision$
%
% $Log$
+% Revision 5.1 2006/10/09 07:49:03 marwan
+% initial check in
+%
%
%
% This program is part of the new generation XXII series.
@@ -334,12 +337,14 @@ clear z z0 z1
z=diff(X); z0 = [];
if ~isempty(find(~(z-1))),z0(:,1)=find(~(z-1));else,z0(1:length(X))=0;end,
if ~isempty(find(~(z+1))),z1=find(~(z+1));else,z1(1:length(X))=0;end
+ z1 = z1(:); z0 = z0(:);
if z0(1)>z1(1)
z0(2:end+1,1)=z0(1:end);z0(1)=0;
end
if length(z0)>length(z1), z0(end)=[]; end
- l=sort(z1-z0); l1=l(find(l>lmin));
+ l=sort(z1-z0);
+ l1=l(find(l>lmin));
end
DET(w-j+1)=sum(l1)/sum(X);
L(w-j+1)=mean(l1);
@@ -489,6 +494,8 @@ clear z z0 z1
if z0(1)>z1(1)
z0(2:end+1,1)=z0(1:end);z0(1)=0;
end
+ z1 = z1(:); z0 = z0(:);
+
if length(z0)>length(z1), z0(end)=[]; end
l=sort(z1-z0); l1=l(find(l>lmin));
end
--
GitLab