%file per stampare e salvare le pzmap in maniera sequenziale

clc
%clear all
close all

Dati_id = getappdata(0,'Dati_id'); 

Dati=Dati_id;

sys=[];

j=1;
p=4;
n=6;

%%


dir_pzmap=['pzmap' ];
%mkdir(dir_pzmap);
%
j=0;

for n=1:6
    for p=1:4

        dir_np=[ 'par___n_' num2str(Dati(1).t.sseh_yy_norm(p,n).n_) '___p_' num2str(Dati(1).t.sseh_yy_norm(p,n).p_) ];
        mkdir(dir_pzmap,dir_np);
            
            for j=1:length(Dati)
    
                dir_sens=['sens_' num2str(j) ];
                mkdir([ dir_pzmap '/' dir_np ],dir_sens);

                for i=1:(3*length(Dati)-1)

                    a.a=Dati(j).t.sseh_yy_norm(p,n).Mnpf.a;
                    b(i).b=Dati(j).t.sseh_yy_norm(p,n).Mnpf.b(i,:);
                    sys(i).sys=tf(Dati(j).t.sseh_yy_norm(p,n).Mnpf.b(i,:),Dati(j).t.sseh_yy_norm(p,n).Mnpf.a);

                    

                    if sum(b(i).b)==0
                       %close(i) 
                    else
                       figure(i)
                       pzmap(sys(i).sys)
                       file_name_eps=[ dir_pzmap '\'  dir_np '\'  dir_sens '\eps_zpmap_ft_'  num2str(i)   '.eps'];
                       file_name_jpg=[ dir_pzmap '\'  dir_np '\'  dir_sens '\jpg_zpmap_ft_'  num2str(i)   '.jpg'];
                       print(gcf, '-depsc', file_name_eps );
                       saveas(gcf,file_name_jpg, 'jpg')
                       close(i)
                    end

                end
        end
    end
end

%%

clc
close all

k=0;

disp(['ricerca poli instabili di sseh_yy_norm' ])
disp([' ' ])
for j=1:length(Dati)
    
    for n=1:6

        for p=1:4

        
            poli=abs(roots(Dati(j).t.sseh_yy_norm(p,n).Mnpf.a));

            indx_pol=find(poli>1);

            for ii=1:length(indx_pol)
                
                k=k+1;
                
                
                disp(['sensore usato = ' num2str(j) ])
                disp(['n = ' num2str(Dati(j).t.sseh_yy_norm(p,n).n_) ])
                disp(['p = ' num2str(Dati(j).t.sseh_yy_norm(p,n).p_) ])
                disp(['modulo =' num2str(poli(indx_pol(ii))) ])
                disp(' ')
                disp(' ')

            end

        end

    end

end

disp(['numero di poli non stabili trovati è = ' num2str(k) ])

%%
clc
close all

k=0;

disp(['ricerca poli instabili di sseh_ny_norm' ])
disp([' ' ])
for j=1:length(Dati)
    
    for n=1:6

        for p=1:4

        
            poli=abs(roots(Dati(j).t.sseh_ny_norm(p,n).Mnpf.a));

            indx_pol=find(poli>1);

            for ii=1:length(indx_pol)
                
                k=k+1;
                
                
                disp(['sensore usato = ' num2str(j) ])
                disp(['n = ' num2str(Dati(j).t.sseh_ny_norm(p,n).n_) ])
                disp(['p = ' num2str(Dati(j).t.sseh_ny_norm(p,n).p_) ])
                disp(['modulo =' num2str(poli(indx_pol(ii))) ])
                disp(' ')
                disp(' ')

            end

        end

    end

end

disp(['numero di poli non stabili trovati è = ' num2str(k) ])

%%

clc
close all
% clear all

Dati_id = getappdata(0,'Dati_id'); 

Dati=Dati_id;

k=0;

disp(['ricerca poli instabili di ssglar_norm' ])
disp([' ' ])


for j=1:length(Dati)
            

            poli=abs(roots(Dati(j).t.ssglar_norm.mod.Mnpg.a));

            indx_pol=find(poli>1);

            for ii=1:length(indx_pol)
                
                k=k+1;
               
                disp(['sensore usato = ' num2str(j) ])
                
                %poli(indx_pol(ii))
                disp(['modulo =' num2str(poli(indx_pol(ii))) ])
                disp(' ')
                disp(' ')

            end

end

disp(['numero di poli non stabili trovati è = ' num2str(k) ])



%%

clc
close all

Dati_id = getappdata(0,'Dati_id'); 

Dati=Dati_id;

disp(['ricerca poli instabili di glar_norm Mgl' ])
disp([' ' ])

k=0;

for j=1:length(Dati)
            

            poli=abs(roots(Dati(j).t.glar_norm.mod(1,1).Mgl.a));

            indx_pol=find(poli>1);

            for ii=1:length(indx_pol)
                
                k=k+1;
               
                disp(['sensore usato = ' num2str(j) ])
                
                %poli(indx_pol(ii))
                disp(['modulo =' num2str(poli(indx_pol(ii))) ])
                disp(' ')
                disp(' ')

            end

end

disp(['numero di poli non stabili trovati è = ' num2str(k) ])


%%

clc
close all

Dati_id = getappdata(0,'Dati_id'); 

Dati=Dati_id;

disp(['ricerca poli instabili di glar_norm Mgls' ])
disp([' ' ])

k=0;
for jj=1:2
    for j=1:length(Dati)
                
                

                poli=abs(roots(Dati(j).t.glar_norm.mod(jj).Mgls.a));

                indx_pol=find(poli>1);

                for ii=1:length(indx_pol)

                    k=k+1;
                    
                    disp(['sensore usato = ' num2str(j) ])
                    disp(['n = ' num2str(Dati(j).t.glar_norm.mod(jj).n) ])

                    disp(['modulo =' num2str(poli(indx_pol(ii))) ])
                    disp(' ')
                    disp(' ')

                end

    end
end

disp(['numero di poli non stabili trovati è = ' num2str(k) ])
