cdfenlei.m

来自「signal and system answer」· M 代码 · 共 55 行

clear
load que1
cdfen1nihe=ones(140,60)*1000;
cdfen2nihe=ones(111,60)*1000;
cdfen3nihe=ones(106,60)*1000;
cdnihe=ones(356,60)*1000;
for i=1:1665
    j=fix(A(i,1));k=fix(A(i,2)+1);
    if isnan(j)==0&isnan(k)==0&k>0
       cdnihe(j,k)=A(i,3);
    end
end
pcd=zeros(356,3);xiabiao=zeros(356,60);
for i=1:356
    for j=1:60
        if cdnihe(i,j)~=1000
            xiabiao(i,j)=j;
        end
    end
end
pcd1=zeros(139,3);xiabiao1=zeros(139,60);
pcd2=zeros(111,3);xiabiao2=zeros(111,60);
pcd3=zeros(106,3);xiabiao3=zeros(106,60);
j=1;
k=1;
l=1;
for i=1:356
    p=1;
    while cdnihe(i,p)==1000
        p=p+1;
    end
    if cdnihe(i,p)>=0&cdnihe(i,p)<=50
       cdfen1nihe(j,:)=cdnihe(i,:);
       xiabiao1(j,:)=xiabiao(i,:);
       j=j+1;
    end
    p=1;
    while cdnihe(i,p)==1000
        p=p+1;
    end
    if cdnihe(i,p)>50&cdnihe(i,p)<=120
        cdfen2nihe(k,:)=cdnihe(i,:);
        xiabiao2(k,:)=xiabiao(i,:);
        k=k+1;
    end
    p=1;
    while cdnihe(i,p)==1000
        p=p+1;
    end
    if cdnihe(i,p)>120
        cdfen3nihe(l,:)=cdnihe(i,:);
        xiabiao3(l,:)=xiabiao(i,:);
        l=l+1;
    end
end