注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mister.Hu(巷里人家)

Go abroad!

 
 
 

日志

 
 
关于我

A campus photograph palyer,an enthusiastic reader,a solitary writer,a future traffic engineer.

网易考拉推荐

Simulation in escalator[Part 3]  

2016-03-05 16:06:41|  分类: Senior |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
function [RM]=Get_RM(TRP,PM,SDM)
%To get the integrated path to solve the conflict.
global Wid;
RM=cell(length(TRP),9);
for k=1:length(TRP)
    for s=1:length(TRP{k})
        A=(TRP{k}(s)-PM(k));
        if A==0%keep static
            [InPM] = Get_InPM(PM(k));
            RM(k,s)={reshape(InPM,1,numel(InPM))};
        else
            switch SDM(k)
                %SDM=1
                case 1
                    InRM=[PM(k)-4*Wid-4 PM(k)-3*Wid-4 PM(k)-2*Wid-4 PM(k)-Wid-4 PM(k)-4 PM(k)+Wid-4;
                        PM(k)-4*Wid-3 PM(k)-3*Wid-3 PM(k)-2*Wid-3 PM(k)-Wid-3 PM(k)-3 PM(k)+Wid-3;
                        PM(k)-4*Wid-2 PM(k)-3*Wid-2 PM(k)-2*Wid-2 PM(k)-Wid-2 PM(k)-2 PM(k)+Wid-2;
                        PM(k)-4*Wid-1 PM(k)-3*Wid-1 PM(k)-2*Wid-1 PM(k)-Wid-1 PM(k)-1 PM(k)+Wid-1;
                        PM(k)-4*Wid PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k) PM(k)+Wid;
                        PM(k)-4*Wid+1 PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1 PM(k)+Wid+1;
                        PM(k)-4*Wid+2 PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2 PM(k)+Wid+2;
                        PM(k)-4*Wid+3 PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3 PM(k)+Wid+3];
                    switch A
                        case -1
                            RM(k,s)={reshape(InRM(4:8,2:5),1,numel(InRM(4:8,2:5)))};
                        case -2
                            RM(k,s)={reshape(InRM(3:8,2:5),1,numel(InRM(3:8,2:5)))};
                        case -3
                            RM(k,s)={reshape(InRM(2:8,2:5),1,numel(InRM(2:8,2:5)))};
                        case -4
                            RM(k,s)={reshape(InRM(1:8,2:5),1,numel(InRM(1:8,2:5)))};
                        case Wid
                            RM(k,s)={reshape(InRM(5:8,2:6),1,numel(InRM(5:8,2:6)))};
                        case -Wid
                            RM(k,s)={reshape(InRM(5:8,1:5),1,numel(InRM(5:8,1:5)))};
                        case Wid-1
                            RM(k,s)={reshape(InRM(4:8,2:6),1,numel(InRM(4:8,2:6)))};
                        case -Wid-1
                            RM(k,s)={reshape(InRM(4:8,1:5),1,numel(InRM(4:8,1:5)))};
                    end
                    %SDM=5
                case 5
                    InRM=[PM(k)-4*Wid PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k) PM(k)+Wid;
                        PM(k)-4*Wid+1 PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1 PM(k)+Wid+1;
                        PM(k)-4*Wid+2 PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2 PM(k)+Wid+2;
                        PM(k)-4*Wid+3 PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3 PM(k)+Wid+3;
                        PM(k)-4*Wid+4 PM(k)-3*Wid+4 PM(k)-2*Wid+4 PM(k)-Wid+4 PM(k)+4 PM(k)+Wid+4;
                        PM(k)-4*Wid+5 PM(k)-3*Wid+5 PM(k)-2*Wid+5 PM(k)-Wid+5 PM(k)+5 PM(k)+Wid+5;
                        PM(k)-4*Wid+6 PM(k)-3*Wid+6 PM(k)-2*Wid+6 PM(k)-Wid+6 PM(k)+6 PM(k)+Wid+6;
                        PM(k)-4*Wid+7 PM(k)-3*Wid+7 PM(k)-2*Wid+7 PM(k)-Wid+7 PM(k)+7 PM(k)+Wid+7;];
                    switch A
                        case 1
                            RM(k,s)={reshape(InRM(1:5,2:5),1,numel(InRM(1:5,2:5)))};
                        case 2
                            RM(k,s)={reshape(InRM(1:6,2:5),1,numel(InRM(1:6,2:5)))};
                        case 3
                            RM(k,s)={reshape(InRM(1:7,2:5),1,numel(InRM(1:7,2:5)))};
                        case 4
                            RM(k,s)={reshape(InRM(1:8,2:5),1,numel(InRM(1:8,2:5)))};
                        case Wid
                            RM(k,s)={reshape(InRM(1:4,2:6),1,numel(InRM(1:4,2:6)))};
                        case -Wid
                            RM(k,s)={reshape(InRM(1:4,1:5),1,numel(InRM(1:4,1:5)))};
                        case Wid+1
                            RM(k,s)={reshape(InRM(1:5,2:6),1,numel(InRM(1:5,2:6)))};
                        case -Wid+1
                            RM(k,s)={reshape(InRM(1:5,1:5),1,numel(InRM(1:5,1:5)))};
                    end
                    %SDM=3
                case 3
                    InRM=[PM(k)-7*Wid-1 PM(k)-6*Wid-1 PM(k)-5*Wid-1 PM(k)-4*Wid-1 PM(k)-3*Wid-1 PM(k)-2*Wid-1 PM(k)-Wid-1 PM(k)-1;
                        PM(k)-7*Wid PM(k)-6*Wid PM(k)-5*Wid PM(k)-4*Wid PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k);
                        PM(k)-7*Wid+1 PM(k)-6*Wid+1 PM(k)-5*Wid+1 PM(k)-4*Wid+1 PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1;
                        PM(k)-7*Wid+2 PM(k)-6*Wid+2 PM(k)-5*Wid+2 PM(k)-4*Wid+2 PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2;
                        PM(k)-7*Wid+3 PM(k)-6*Wid+3 PM(k)-5*Wid+3 PM(k)-4*Wid+3 PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3;
                        PM(k)-7*Wid+4 PM(k)-6*Wid+4 PM(k)-5*Wid+4 PM(k)-4*Wid+4 PM(k)-3*Wid+4 PM(k)-2*Wid+4 PM(k)-Wid+4 PM(k)+4];
                    switch A
                        case -Wid
                            RM(k,s)={reshape(InRM(2:5,4:8),1,numel(InRM(2:5,4:8)))};
                        case -2*Wid
                            RM(k,s)={reshape(InRM(2:5,3:8),1,numel(InRM(2:5,3:8)))};
                        case -3*Wid
                            RM(k,s)={reshape(InRM(2:5,2:8),1,numel(InRM(2:5,2:8)))};
                        case -4*Wid
                            RM(k,s)={reshape(InRM(2:5,1:8),1,numel(InRM(2:5,1:8)))};
                        case -Wid+1
                            RM(k,s)={reshape(InRM(2:6,4:8),1,numel(InRM(2:6,4:8)))};
                        case -Wid-1
                            RM(k,s)={reshape(InRM(1:5,4:8),1,numel(InRM(1:5,4:8)))};
                        case 1
                            RM(k,s)={reshape(InRM(2:6,5:8),1,numel(InRM(2:6,5:8)))};
                        case -1
                            RM(k,s)={reshape(InRM(1:5,5:8),1,numel(InRM(1:5,5:8)))};
                    end
                    %SDM=7
                case 7
                    InRM=[PM(k)-3*Wid-1 PM(k)-2*Wid-1 PM(k)-Wid-1 PM(k)-1 PM(k)+Wid-1 PM(k)+2*Wid-1 PM(k)+3*Wid-1 PM(k)+4*Wid-1;
                        PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k) PM(k)+Wid PM(k)+2*Wid PM(k)+3*Wid PM(k)+4*Wid;
                        PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1 PM(k)+Wid+1 PM(k)+2*Wid+1 PM(k)+3*Wid+1 PM(k)+4*Wid+1;
                        PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2 PM(k)+Wid+2 PM(k)+2*Wid+2 PM(k)+3*Wid+2 PM(k)+4*Wid+2;
                        PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3 PM(k)+Wid+3 PM(k)+2*Wid+3 PM(k)+3*Wid+3 PM(k)+4*Wid+3;
                        PM(k)-3*Wid+4 PM(k)-2*Wid+4 PM(k)-Wid+4 PM(k)+4 PM(k)+Wid+4 PM(k)+2*Wid+4 PM(k)+3*Wid+4 PM(k)+4*Wid+4];
                    switch A
                        case Wid
                            RM(k,s)={reshape(InRM(2:5,1:5),1,numel(InRM(2:5,1:5)))};
                        case 2*Wid
                            RM(k,s)={reshape(InRM(2:5,1:6),1,numel(InRM(2:5,1:6)))};
                        case 3*Wid
                            RM(k,s)={reshape(InRM(2:5,1:7),1,numel(InRM(2:5,1:7)))};
                        case 4*Wid
                            RM(k,s)={reshape(InRM(2:5,1:8),1,numel(InRM(2:5,1:8)))};
                        case Wid-1
                            RM(k,s)={reshape(InRM(1:5,1:5),1,numel(InRM(1:5,1:5)))};
                        case Wid+1
                            RM(k,s)={reshape(InRM(2:6,1:5),1,numel(InRM(2:6,1:5)))};
                        case -1
                            RM(k,s)={reshape(InRM(1:5,1:4),1,numel(InRM(1:5,1:4)))};
                        case 1
                            RM(k,s)={reshape(InRM(2:6,1:4),1,numel(InRM(2:6,1:4)))};
                    end
                    %SDM=2
                case 2
                    InRM=[PM(k)-5*Wid-2 PM(k)-4*Wid-2 PM(k)-3*Wid-2 PM(k)-2*Wid-2 PM(k)-Wid-2 PM(k)-2;
                        PM(k)-5*Wid-1 PM(k)-4*Wid-1 PM(k)-3*Wid-1 PM(k)-2*Wid-1 PM(k)-Wid-1 PM(k)-1;
                        PM(k)-5*Wid PM(k)-4*Wid PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k);
                        PM(k)-5*Wid+1 PM(k)-4*Wid+1 PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1;
                        PM(k)-5*Wid+2 PM(k)-4*Wid+2 PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2;
                        PM(k)-5*Wid+3 PM(k)-4*Wid+3 PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3];
                    switch A
                        case -Wid-1
                            RM(k,s)={reshape(InRM(2:6,2:6),1,numel(InRM(2:6,2:6)))};
                        case -2*Wid-2
                            RM(k,s)={reshape(InRM(1:6,1:6),1,numel(InRM(1:6,1:6)))};
                        case -1
                            RM(k,s)={reshape(InRM(2:6,3:6),1,numel(InRM(2:6,3:6)))};
                        case -2
                            RM(k,s)={reshape(InRM(1:6,3:6),1,numel(InRM(1:6,3:6)))};
                        case -Wid
                            RM(k,s)={reshape(InRM(3:6,2:6),1,numel(InRM(3:6,2:6)))};
                        case -2*Wid
                            RM(k,s)={reshape(InRM(3:6,1:6),1,numel(InRM(3:6,1:6)))};
                        case -Wid-2
                            RM(k,s)={reshape(InRM(1:6,2:6),1,numel(InRM(1:6,2:6)))};
                        case -2*Wid-1
                            RM(k,s)={reshape(InRM(2:6,1:6),1,numel(InRM(2:6,1:6)))};
                    end
                    %SDM=6
                case 6
                    InRM=[PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k) PM(k)+Wid PM(k)+2*Wid;
                        PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1 PM(k)+Wid+1 PM(k)+2*Wid+1;
                        PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2 PM(k)+Wid+2 PM(k)+2*Wid+2;
                        PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3 PM(k)+Wid+3 PM(k)+2*Wid+3;
                        PM(k)-3*Wid+4 PM(k)-2*Wid+4 PM(k)-Wid+4 PM(k)+4 PM(k)+Wid+4 PM(k)+2*Wid+4;
                        PM(k)-3*Wid+5 PM(k)-2*Wid+5 PM(k)-Wid+5 PM(k)+5 PM(k)+Wid+5 PM(k)+2*Wid+5];
                    switch A
                        case Wid+1
                            RM(k,s)={reshape(InRM(1:5,1:5),1,numel(InRM(1:5,1:5)))};
                        case 2*Wid+2
                            RM(k,s)={reshape(InRM(1:6,1:6),1,numel(InRM(1:6,1:6)))};
                        case 1
                            RM(k,s)={reshape(InRM(1:5,1:4),1,numel(InRM(1:5,1:4)))};
                        case 2
                            RM(k,s)={reshape(InRM(1:6,1:4),1,numel(InRM(1:6,1:4)))};
                        case Wid
                            RM(k,s)={reshape(InRM(1:4,1:5),1,numel(InRM(1:4,1:5)))};
                        case 2*Wid
                            RM(k,s)={reshape(InRM(1:4,1:6),1,numel(InRM(1:4,1:6)))};
                        case Wid+2
                            RM(k,s)={reshape(InRM(1:6,1:5),1,numel(InRM(1:6,1:5)))};
                        case 2*Wid+1
                            RM(k,s)={reshape(InRM(1:5,1:6),1,numel(InRM(1:5,1:6)))};
                    end
                    %SDM=4
                case 4
                    InRM=[PM(k)-5*Wid PM(k)-4*Wid PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k);
                        PM(k)-5*Wid+1 PM(k)-4*Wid+1 PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1;
                        PM(k)-5*Wid+2 PM(k)-4*Wid+2 PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2;
                        PM(k)-5*Wid+3 PM(k)-4*Wid+3 PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3;
                        PM(k)-5*Wid+4 PM(k)-4*Wid+4 PM(k)-3*Wid+4 PM(k)-2*Wid+4 PM(k)-Wid+4 PM(k)+4;
                        PM(k)-5*Wid+5 PM(k)-4*Wid+5 PM(k)-3*Wid+5 PM(k)-2*Wid+5 PM(k)-Wid+5 PM(k)+5];
                    switch A
                        case -Wid+1
                            RM(k,s)={reshape(InRM(1:5,2:6),1,numel(InRM(1:5,2:6)))};
                        case -2*Wid+2
                            RM(k,s)={reshape(InRM(1:6,1:6),1,numel(InRM(1:6,1:6)))};
                        case 1
                            RM(k,s)={reshape(InRM(1:5,3:6),1,numel(InRM(1:5,3:6)))};
                        case 2
                            RM(k,s)={reshape(InRM(1:6,3:6),1,numel(InRM(1:6,3:6)))};
                        case -Wid
                            RM(k,s)={reshape(InRM(1:4,2:6),1,numel(InRM(1:4,2:6)))};
                        case -2*Wid
                            RM(k,s)={reshape(InRM(1:4,1:6),1,numel(InRM(1:4,1:6)))};
                        case -Wid+2
                            RM(k,s)={reshape(InRM(1:6,2:6),1,numel(InRM(1:6,2:6)))};
                        case -2*Wid+1
                            RM(k,s)={reshape(InRM(1:5,1:6),1,numel(InRM(1:5,1:6)))};
                    end
                    %SDM=8
                case 8
                    InRM=[PM(k)-3*Wid-2 PM(k)-2*Wid-2 PM(k)-Wid-2 PM(k)-2 PM(k)+Wid-2 PM(k)+2*Wid-2;
                        PM(k)-3*Wid-1 PM(k)-2*Wid-1 PM(k)-Wid-1 PM(k)-1 PM(k)+Wid-1 PM(k)+2*Wid-1;
                        PM(k)-3*Wid PM(k)-2*Wid PM(k)-Wid PM(k) PM(k)+Wid PM(k)+2*Wid;
                        PM(k)-3*Wid+1 PM(k)-2*Wid+1 PM(k)-Wid+1 PM(k)+1 PM(k)+Wid+1 PM(k)+2*Wid+1;
                        PM(k)-3*Wid+2 PM(k)-2*Wid+2 PM(k)-Wid+2 PM(k)+2 PM(k)+Wid+2 PM(k)+2*Wid+2;
                        PM(k)-3*Wid+3 PM(k)-2*Wid+3 PM(k)-Wid+3 PM(k)+3 PM(k)+Wid+3 PM(k)+2*Wid+3];
                    switch A
                        case Wid-1
                            RM(k,s)={reshape(InRM(2:6,1:5),1,numel(InRM(2:6,1:5)))};
                        case 2*Wid-2
                            RM(k,s)={reshape(InRM(1:6,1:6),1,numel(InRM(1:6,1:6)))};
                        case -1
                            RM(k,s)={reshape(InRM(2:6,1:4),1,numel(InRM(2:6,1:4)))};
                        case -2
                            RM(k,s)={reshape(InRM(1:6,1:4),1,numel(InRM(1:6,1:4)))};
                        case Wid
                            RM(k,s)={reshape(InRM(3:6,1:5),1,numel(InRM(3:6,1:5)))};
                        case 2*Wid
                            RM(k,s)={reshape(InRM(3:6,1:6),1,numel(InRM(3:6,1:6)))};
                        case Wid-2
                            RM(k,s)={reshape(InRM(1:6,1:5),1,numel(InRM(1:6,1:5)))};
                        case 2*Wid-1
                            RM(k,s)={reshape(InRM(2:6,1:6),1,numel(InRM(2:6,1:6)))};
                    end
            end
        end
    end
    %     RM(k,(length(TRP{k})+1):9)={[]};
end
% RM((length(TRP)+1):TNUM,9)={[]};
% RM(cellfun(@(x) ~isempty(x)&&any(isnan(SM(x))),RM,'UniformOutput',true))={[]};
end

function [RM] = RM_Refine(fsub,RM,SM,PM)
%to refine the RM.
if fsub~=0
    for k=1:9 %to refine the RM.
        for s=fsub:length(PM)
            if ~isempty(RM{s,k})
                if any(isnan(SM(RM{s,k})))%Delete the NaN
                    RM{s,k}=[];
                else
                    PMA=PM(s+1:end);
                    InPM=reshape(Get_InPM(PMA),1,numel(Get_InPM(PMA)));
                    if ~isempty(intersect(RM{s,k},InPM))
                        RM{s,k}=[];
                    end
                end
            end
        end
    end
end
end

function [FRM]=Solve(PP,fsub,RM,PARM)
FRM=cell(length(PP),1);
if fsub~=0
    for k=1:9
        if ~isempty(RM{fsub,k})
            FRM{fsub,1}=PARM{fsub,k};
            break
        end
    end
    if length(PP)>=fsub+1
        for k=fsub+1:length(PP)%k means the main man now.
            for s=1:9%s can indicate the different destination of the main man.
                K=1;
                for v=fsub:(k-1) %the main man should compare with the man before him.
                    if isempty(RM{k,s})||~isempty(intersect(RM{k,s},reshape(Get_InPM(FRM{v,1}'),1,numel(Get_InPM(FRM{v,1}')))))
                        K=0;
                        break
                    end
                end
                if K
                    FRM{k,1}=PARM{k,s};
                    break
                end
            end
        end
        
    end
end
end

function [FSVM,SDM]=Update_Direction(PP,FRM,PM,SDM)
%update the direction of pedes.
global Wid;
FSVM=zeros(length(PP),1);
for k=1:length(PP)
    if ~isempty(FRM{k})
        A=(FRM{k}(end)-PM(k));
        switch SDM(k)
            %SDM=1
            case 1
                switch A
                    case -1
                        FSVM(k)=1;SDM(k)=1;
                    case -2
                        FSVM(k)=2;SDM(k)=1;
                    case -3
                        FSVM(k)=3;SDM(k)=1;
                    case -4
                        FSVM(k)=4;SDM(k)=1;
                    case Wid
                        FSVM(k)=2;SDM(k)=7;
                    case -Wid
                        FSVM(k)=2;SDM(k)=3;
                    case -Wid-1
                        FSVM(k)=2;SDM(k)=2;
                    case Wid-1
                        FSVM(k)=2;SDM(k)=8;
                end
                %SDM=5
            case 5
                switch A
                    case 1
                        FSVM(k)=1;SDM(k)=5;
                    case 2
                        FSVM(k)=2;SDM(k)=5;
                    case 3
                        FSVM(k)=3;SDM(k)=5;
                    case 4
                        FSVM(k)=4;SDM(k)=5;
                    case Wid
                        FSVM(k)=2;SDM(k)=7;
                    case -Wid
                        FSVM(k)=2;SDM(k)=3;
                    case Wid+1
                        FSVM(k)=2;SDM(k)=6;
                    case -Wid+1
                        FSVM(k)=2;SDM(k)=4;
                end
                %SDM=3
            case 3
                switch A
                    case -Wid
                        FSVM(k)=1;SDM(k)=3;
                    case -2*Wid
                        FSVM(k)=2;SDM(k)=3;
                    case -3*Wid
                        FSVM(k)=3;SDM(k)=3;
                    case -4*Wid
                        FSVM(k)=4;SDM(k)=3;
                    case -Wid+1
                        FSVM(k)=2;SDM(k)=4;
                    case -Wid-1
                        FSVM(k)=2;SDM(k)=2;
                    case 1
                        FSVM(k)=1;SDM(k)=5;
                    case -1
                        FSVM(k)=1;SDM(k)=1;
                end
                %SDM=7
            case 7
                switch A
                    case Wid
                        FSVM(k)=1;SDM(k)=7;
                    case 2*Wid
                        FSVM(k)=2;SDM(k)=7;
                    case 3*Wid
                        FSVM(k)=3;SDM(k)=7;
                    case 4*Wid
                        FSVM(k)=4;SDM(k)=7;
                    case Wid-1
                        FSVM(k)=2;SDM(k)=8;
                    case Wid+1
                        FSVM(k)=2;SDM(k)=6;
                    case -1
                        FSVM(k)=1;SDM(k)=1;
                    case 1
                        FSVM(k)=1;SDM(k)=5;
                end
                %SDM=2
            case 2
                switch A
                    case -Wid-1
                        FSVM(k)=2;SDM(k)=2;
                    case -2*Wid-2
                        FSVM(k)=4;SDM(k)=2;
                    case -1
                        FSVM(k)=2;SDM(k)=1;
                    case -2
                        FSVM(k)=3;SDM(k)=1;
                    case -Wid
                        FSVM(k)=2;SDM(k)=3;
                    case -2*Wid
                        FSVM(k)=3;SDM(k)=3;
                    case -Wid-2
                        FSVM(k)=3;SDM(k)=1;
                    case -2*Wid-1
                        FSVM(k)=3;SDM(k)=3;
                end
                %SDM=6
            case 6
                switch A
                    case Wid+1
                        FSVM(k)=2;SDM(k)=6;
                    case 2*Wid+2
                        FSVM(k)=4;SDM(k)=6;
                    case 1
                        FSVM(k)=2;SDM(k)=5;
                    case 2
                        FSVM(k)=3;SDM(k)=5;
                    case Wid
                        FSVM(k)=2;SDM(k)=7;
                    case 2*Wid
                        FSVM(k)=3;SDM(k)=7;
                    case Wid+2
                        FSVM(k)=3;SDM(k)=5;
                    case 2*Wid+1
                        FSVM(k)=3;SDM(k)=7;
                end
                %SDM=4
            case 4
                switch A
                    case -Wid+1
                        FSVM(k)=2;SDM(k)=4;
                    case -2*Wid+2
                        FSVM(k)=4;SDM(k)=4;
                    case 1
                        FSVM(k)=2;SDM(k)=5;
                    case 2
                        FSVM(k)=2;SDM(k)=5;
                    case -Wid
                        FSVM(k)=2;SDM(k)=3;
                    case -2*Wid
                        FSVM(k)=3;SDM(k)=3;
                    case -Wid+2
                        FSVM(k)=3;SDM(k)=5;
                    case -2*Wid+1
                        FSVM(k)=3;SDM(k)=3;
                end
                %SDM=8
            case 8
                switch A
                    case Wid-1
                        FSVM(k)=2;SDM(k)=8;
                    case 2*Wid-2
                        FSVM(k)=4;SDM(k)=8;
                    case -1
                        FSVM(k)=2;SDM(k)=1;
                    case -2
                        FSVM(k)=3;SDM(k)=1;
                    case Wid
                        FSVM(k)=2;SDM(k)=7;
                    case 2*Wid
                        FSVM(k)=3;SDM(k)=7;
                    case Wid-2
                        FSVM(k)=3;SDM(k)=1;
                    case 2*Wid-1
                        FSVM(k)=3;SDM(k)=7;
                end
        end
    end
end
end
  评论这张
 
阅读(16)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016