aparat.m

different kind of measurements

% prelucrarea datelor obtinute la aparat
clc	
clear
j=1:1:1601;
incarc1=load (get(findobj(gcbf,'tag','fisierul1'),'string')); %incarc fisierul dorit
incarc2=load (get(findobj(gcbf,'tag','fisierul2'),'string'));
incarc3=load (get(findobj(gcbf,'tag','fisierul3'),'string'));

n1=str2double(get(findobj(gcbf,'tag','numar1'),'string')); % numarul de puncte semnificative (201, 401, sau 801)
n2=str2double(get(findobj(gcbf,'tag','numar2'),'string'));
n3=str2double(get(findobj(gcbf,'tag','numar3'),'string'));

jos1=str2double(get(findobj(gcbf,'tag','lim1'),'string')); % domeniul de frecvente in Mhz
sus1=str2double(get(findobj(gcbf,'tag','lim2'),'string'));
jos2=str2double(get(findobj(gcbf,'tag','lim3'),'string'));
sus2=str2double(get(findobj(gcbf,'tag','lim4'),'string'));
jos3=str2double(get(findobj(gcbf,'tag','lim5'),'string'));
sus3=str2double(get(findobj(gcbf,'tag','lim6'),'string'));

for j=1:1601
c1(j)=incarc1(j)+i*incarc1(j+1601); 
c2(j)=incarc2(j)+i*incarc2(j+1601);
c3(j)=incarc3(j)+i*incarc3(j+1601);
end

for j=1:n1
    c11(j)=c1(j);
    c71(j)=abs(c11(j))^2;     
end    

for j=1:n2
    c12(j)=c2(j);
    c72(j)=abs(c12(j))^2;     
end    

for j=1:n3
    c13(j)=c3(j);
    c73(j)=abs(c13(j))^2;     
end    
 

c21=ifft(c11,n1); %raspunsul in timp
c22=ifft(c12,n2);
c23=ifft(c13,n3);

c31=abs(c21); % modulul
c32=abs(c22);
c33=abs(c23);

% pentru camera reala
R=str2double(get(findobj(gcbf,'tag','raport'),'string'));
LOS1=str2double(get(findobj(gcbf,'tag','vad1'),'string'));
LOS2=str2double(get(findobj(gcbf,'tag','vad2'),'string'));
LOS3=str2double(get(findobj(gcbf,'tag','vad3'),'string'));

c31R=c31/(R^0.5);% pentru camera reala (mai mare)
c32R=c32/(R^0.5);
c33R=c33/(R^0.5);

maxim1=max(c31)
for i= 1:length(c31)
    if c31(1,i)==maxim1
       dep1=i
   end
end
if LOS1 ==1
    c31R(1,dep1)=c31R(1,dep1)*(R^0.5)
end  
zerouri1=zeros(1,R*dep1-dep1)
cR31R = [zerouri1 c31R(1:length(c31R)-length(zerouri1))];

maxim2=max(c32)
for i= 1:length(c32)
    if c32(1,i)==maxim2
       dep2=i
   end
end
if LOS2 ==1
    c32R(1,dep2)=c32R(1,dep2)*(R^0.5)
end  
zerouri2=zeros(1,R*dep2-dep2)
cR32R = [zerouri2 c32R(1:length(c32R)-length(zerouri2))];

maxim3=max(c33)
for i= 1:length(c33)
    if c33(1,i)==maxim3
       dep3=i
   end
end
if LOS3 ==1
    c33R(1,dep3)=c33R(1,dep3)*(R^0.5)
end  
zerouri3=zeros(1,R*dep3-dep3)
cR33R = [zerouri3 c33R(1:length(c33R)-length(zerouri3))];

c61=20*(log10(abs(c11)));
c62=20*(log10(abs(c12)));
c63=20*(log10(abs(c13)));

x1=1:1:n1;

f1=jos1*10^6+(x1-1).*((sus1-jos1)*10^6)/(n1-1);
t1=(x1-1).*1/((sus1-jos1)*10^6);

x2=1:1:n2;

f2=jos2*10^6+(x2-1).*((sus2-jos2)*10^6)/(n2-1);
t2=(x2-1).*1/((sus2-jos2)*10^6);

x3=1:1:n3;

f3=jos3*10^6+(x3-1).*((sus3-jos3)*10^6)/(n3-1);
t3=(x3-1).*1/((sus3-jos3)*10^6);

figure(1)

    subplot(3,3,1); plot(t1,c31);
    ylabel((get(findobj(gcbf,'tag','masura1'),'string'))); %axes
    xlabel(' Timp '); %axes
    axis([min(t1) max(t1) min(c31) max(c31)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie1'),'string'))));
    
    subplot(3,3,4); plot(f1,c61);
    ylabel((get(findobj(gcbf,'tag','masura4'),'string'))); %axes
    xlabel(' Frecv '); %axes
    axis([jos1*10^6 sus1*10^6 min(c61) max(c61)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie1'),'string'))));

    subplot(3,3,2); plot(t2,c32);
    ylabel((get(findobj(gcbf,'tag','masura2'),'string'))); %axes
    xlabel(' Timp '); %axes
    axis([min(t2) max(t2) min(c32) max(c32)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie2'),'string'))));
    
    subplot(3,3,5); plot(f2,c62);
    ylabel((get(findobj(gcbf,'tag','masura5'),'string'))); %axes
    xlabel(' Frecv '); %axes
    axis([jos2*10^6 sus2*10^6 min(c62) max(c62)]); 
    title(sprintf((get(findobj(gcbf,'tag','scrie2'),'string'))));

    subplot(3,3,3); plot(t3,c33);
    ylabel((get(findobj(gcbf,'tag','masura3'),'string'))); %axes
    xlabel(' Timp '); %axes
    axis([min(t3) max(t3) min(c33) max(c33)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie3'),'string'))));
    
    subplot(3,3,6); plot(f3,c63);
    ylabel((get(findobj(gcbf,'tag','masura6'),'string'))); %axes
    xlabel(' Frecv '); %axes
    axis([jos3*10^6 sus3*10^6 min(c63) max(c63)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie3'),'string'))));
    
    subplot(3,3,7); plot(t1,cR31R);
    ylabel((get(findobj(gcbf,'tag','masura1'),'string'))); %axes
    xlabel(' Timp '); %axes
    axis([min(t1) max(t1) min(c31) max(c31)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie1'),'string'))));
    
    subplot(3,3,8); plot(t2,cR32R);
    ylabel((get(findobj(gcbf,'tag','masura2'),'string'))); %axes
    xlabel(' Timp '); %axes
    axis([min(t2) max(t2) min(c32) max(c32)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie2'),'string'))));
    
    subplot(3,3,9); plot(t3,cR33R);
    ylabel((get(findobj(gcbf,'tag','masura3'),'string'))); %axes
    xlabel(' Timp '); %axes
    axis([min(t3) max(t3) min(c33) max(c33)]);
    title(sprintf((get(findobj(gcbf,'tag','scrie3'),'string'))));