sgod.m

来自「用matlab实现利用统计混沌方法解决非线性系统时间序列预测的问题」· M 代码 · 共 59 行

function sgod(filename,contrast);

if nargin<2,
        contrast = 2.0;
end;

% figure
figure;

% coord max
xmax = 0;
ymax = 0;
emax = 0;

% SST lecture fichier filename.sst
[sst,rsst,csst] = ld(strcat(filename,'.sst'));
if rsst*csst==0, oksst=0; else, oksst=1; end;

if oksst==1,
        % SGD lecture fichier filename.sgd
        [sgd,rsgd,csgd] = ld(strcat(filename,'.sgd'));
        if rsgd*csgd==0, oksgd=0; else, oksgd=1; end;
        
        % coord max
        if oksgd==1,
                xmax = max(xmax,sst(max(sgd(:,1))));
                ymax = max(ymax,sst(max(sgd(:,2))));
        end;
        emax = max(emax,max(sst(:,2)));
end;

if oksst==1 & oksgd==1,
        % SGD sous figure
        colormap('gray');
        
        % SGD pas
        step = sst(2,1) - sst(1,1);
        xmax = xmax + 3*step;
        ymax = ymax + 3*step;
        
        % SGD polygone d'arriere plan
        patch([0 0 xmax xmax],[0 ymax ymax 0],+contrast);
        axis([0 xmax 0 ymax]);
        
        % SGD lecture du god
        for cnt = 1:rsgd,
                % polygone
                patch(sst(1+sgd(cnt,1))+[0 0 step step],sst(1+sgd(cnt,2))+[0 step step 0],-sgd(cnt,3));
        end;
        
        shading flat;
        
        % ligne diagonale
        line([0 min(xmax,ymax)],[0 min(xmax,ymax)]);
        
        % SGD ornements
        axis equal;
        axis off;
end;