clc

clear all

close all

%ixd 5 sensore esterno
%ixd 8 sensore termi

start_val=1000;
end_val=1100;
str=1;
Dati_ch = getappdata(0,'Dati_ch'); 
Dati_id = getappdata(0,'Dati_id'); 

Dati=Dati_id;
%mi manca y 1 la recupero fondendo udati con y1
Udati_norm=[Dati(1,1).t.y_norm  Dati(1,1).t.Udati_norm];
Udati_org=Udati_norm;

[r c]=size(Udati_norm);


for i=1:length(Dati)
    
    tmp=[ Udati_norm(1:start_val,i); zeros(r-start_val,1) ];
    
%     if ( i==1 || i==2 || i==3 || i==4 || i==5 )
%     else
%         Udati_norm(:,i)=tmp;
%     end
    
    if ( i==11 )
        Udati_norm(:,i)=tmp;
    else
    end
    
    
end

%

for j=1:length(Dati)
       Model(j).mod=Dati(j).t.sseh_yy_norm(2,4).Mnpf;
end
   
U=[];


for i=start_val:end_val
    
    i
    
    for j=1:length(Dati)
        U(j).U=cancel_y_from_U(Udati_norm,j);
        my_y=Udati_norm(str:i,j);
        my_dati=U(j).U(str:i,:);
        Dati(j).iddata=iddata(my_y,my_dati);
        %Dati(j).iddata=iddata(Udati_norm(str:i,j),U(j).U(str:i,:));
        y(j).iddata=predict(Model(j).mod,Dati(j).iddata,1);
        if Udati_norm(end,j)==0 && Udati_norm(end-1,j)==0 && Udati_norm(end-2,j)==0 && Udati_norm(end-3,j)==0 
            Udati_norm(i+1,j)=y(j).iddata.y(i);
        end
    end
 
end

disp('fine')

%end_val=end_val-1;

%%

close all

for i=1:length(Dati)
    v=var(Dati(i).t.y_org)^(1/2);
    m=mean(Dati(i).t.y_org);

    figure(i)
    hold on
    
    plot(Udati_norm(str:end_val,i),'b')
    plot(Udati_org(str:end_val,i),'r')
    plot((Udati_norm(str:end_val,i)-Udati_org(str:end_val,i)),'g')
    
    %plot(Udati_norm(str:end_val,i)*v+m,'b')
    %plot(Udati_org(str:end_val,i)*v+m,'r')
    %plot((Udati_norm(str:end_val,i)-Udati_org(str:end_val,i))*v,'g')
    %axis([(start_val-100) (end_val-100) 0 50])
    hold off
    somma=sum(Udati_norm(1:end_val,i)-Udati_org(1:end_val,i));
    if somma==0
        close(i)
    end
end



%print(gcf, '-depsc', 'figure\bode_filter.eps')
