qiuzhongheliaoxiaozhibiao.m

signal and system answer

clear
load data1
clear A
a=zeros(1239,5);b=zeros(1251,5);c=zeros(1254,5);d=zeros(1292,5);
j=1;k=1;l=1;m=1;
for i=1:5036
    if B(i,2)==1
        a(j,:)=B(i,:);
        j=j+1;
    end
     if B(i,2)==2
        b(k,:)=B(i,:);
        k=k+1;
     end
     if B(i,2)==3
        c(l,:)=B(i,:);
        l=l+1;
     end
     if B(i,2)==4
        d(m,:)=B(i,:);
        m=m+1;
    end
end
k=1;j=1;l=0;i=1;r=0;
while i<=1239
    if a(i,1)==k&r==0
        l=l+1;
    end
    if a(i,1)==k
        r=1;
        chafen1(l,j)=a(i,5);
        j=j+1;
        i=i+1;
    else
        k=k+1;
        j=1;
        r=0;
    end
end
k=1;j=1;l=0;i=1;r=0;
while i<=1251
    if b(i,1)==k&r==0
        l=l+1;
    end
    if b(i,1)==k
        r=1;
        chafen2(l,j)=b(i,5);
        j=j+1;
        i=i+1;
    else
        k=k+1;
        j=1;
        r=0;
    end
end
k=1;j=1;l=0;i=1;r=0;
while i<=1254
    if c(i,1)==k&r==0
        l=l+1;
    end
    if c(i,1)==k
        r=1;
        chafen3(l,j)=c(i,5);
        j=j+1;
        i=i+1;
    else
        k=k+1;
        j=1;
        r=0;
    end
end
k=1;j=1;l=0;i=1;r=0;
while i<=1292
    if d(i,1)==k&r==0
        l=l+1;
    end
    if d(i,1)==k
        r=1;
        chafen4(l,j)=d(i,5);
        j=j+1;
        i=i+1;
    else
        k=k+1;
        j=1;
        r=0;
    end
end
for i=1:325
    for j=1:6
        qq(i,j)=(exp(chafen1(i,j+1))-exp(chafen1(i,1)))/(exp(chafen1(i,1))-1);
    end
end
for i=1:324
    for j=1:6
        rr(i,j)=(exp(chafen2(i,j+1))-exp(chafen2(i,1)))/(exp(chafen2(i,1))-1);
    end
end
for i=1:330
    for j=1:8
        ss(i,j)=(exp(chafen3(i,j+1))-exp(chafen3(i,1)))/(exp(chafen3(i,1))-1);
    end
end
for i=1:330
    for j=1:6
        tt(i,j)=(exp(chafen4(i,j+1))-exp(chafen4(i,1)))/(exp(chafen4(i,1))-1);
    end
end
for i=1:325
    for j=1:6
        if isnan(qq(i,j))==1
            qq(i,j)=0;
        end
        if isinf(qq(i,j))==1
            qq(i,j)=0;
        end
         if qq(i,j)==-1
            qq(i,j)=0;
        end
    end
end
for i=1:324
    for j=1:6
          if isnan(rr(i,j))==1
            rr(i,j)=0;
          end
          if isinf(rr(i,j))==1
            rr(i,j)=0;
        end
         if rr(i,j)==-1
            rr(i,j)=0;
        end
    end
end
for i=1:330
    for j=1:8
          if ss(i,j)==-1
            ss(i,j)=0;
          end
          if isinf(ss(i,j))==1
            ss(i,j)=0;
        end
         if isnan(ss(i,j))==1
            ss(i,j)=0;
        end
    end
end
for i=1:330
    for j=1:6
          if isinf(tt(i,j))==1
            tt(i,j)=0;
          end
          if isnan(tt(i,j))==1
            tt(i,j)=0;
        end
         if tt(i,j)==-1
            tt(i,j)=0;
        end
    end
end
zz=zeros(1,6);xx=zeros(1,6);cc=zeros(1,8);dd=zeros(1,6);
for i=1:325
    for j=1:6
       
         if qq(i,j)~=0
            zz(j)=zz(j)+1;
        end
    end
end
for i=1:324
    for j=1:6
        
         if rr(i,j)~=0
           xx(j)=xx(j)+1;
        end
    end
end
for i=1:330
    for j=1:8
          if ss(i,j)~=0
            cc(j)=cc(j)+1;
          end
         
    end
end
for i=1:330
    for j=1:6
          
         if tt(i,j)~=0
           dd(j)=dd(j)+1;
        end
    end
end