Matlab - Algoritmo para reemplazar la funcion sigadd

%Signal Addition 

function [y,n] = sigadd(x1,n1,x2,n2)

%implements y(n)=x1(n)+x2(n)

%--------------------------

%[y,n]= sigadd(x1,n1,x2,n2)

% y = sum sequence over n , which includes n1 and n2

% x1 = first sequence over n1 

% x2 = second sequence over n2 ( n2 can be different from n1 )

%

n=min(min(n1),min(n2)):max(max(n1),max(n2)); % duration of y(n)

y1=zeros(1,length(n)); y2=y1;   %initialization

y1(find((n>=min(n1)) & (n<=max(n1))==1))=x1   %x1 with duration of y

y2(find((n>=min(n2)) & (n<=max(n2))==1))=x2   %x2 with duration of y

y=y1+y2;                                      %sequence addition

function m=max_2(varargin)

if nargin==1

    v=varargin{1};

    m=v(1);

    for i=v

        if i>m

            m=i;

        end

    end

elseif nargin==2

    a=varargin{1};

    b=varargin{2};

    if a>b

        m=a;

    else

        m=b;

    end

end

end

function m=min_2(varargin)

if nargin==1

    v=varargin{1};

    m=v(1);

    for i=v

        if i<m

            m=i;

        end

    end

elseif nargin==2

    a=varargin{1};

    b=varargin{2};

    if a<b

        m=a;

    else

        m=b;

    end

end

end

function ind=find_2(cond)

MC=cond;

k=1;

for i=1:length_2(MC)

    if MC(i)>0

        ind(k)=i;

        k=k+1;

    end

end

end

function L=length_2(v)

L=0;

for i=v

    L=L+1;

end

end

function X=zeros_2(varargin)

if nargin==1

    dim=varargin{1};

    X(dim,dim)=0;

elseif nargin==2

    nf=varargin{1};

    nc=varargin{2};

    X(nf,nc)=0;

end

end