fit_various

d = read_arc('27-jan-2010:00:12:00','27-jan-2010:00:20:00'); %cygA
d = read_arc('27-jan-2010:00:39:30','27-jan-2010:00:40:30'); %cygA

0001 %d = read_arc('27-jan-2010:00:12:00','27-jan-2010:00:20:00'); %cygA
0002 %d = read_arc('27-jan-2010:00:39:30','27-jan-2010:00:40:30'); %cygA
0003 d = read_arc('27-jan-2010:02:46:21','27-jan-2010:02:46:26'); %casA
0004 
0005 
0006 d = reshapeRx(d);
0007 d = reshapeAz(d);
0008 d = interpAz(d);
0009 
0010 az = d.antenna0.servo.interpAz;
0011 rx = d.antenna0.receiver.dataNew(:,1);
0012 
0013 plot(az,rx,'.');
0014 
0015 beta0(1) = min(rx)-max(rx);
0016 minrx=find(rx==min(rx));
0017 beta0(2) = az(minrx(1));
0018 beta0(3) = 1;
0019 beta0(4) = 1;
0020 beta0(5) = mean(rx);
0021 
0022 [betanew]=nlinfit(az,rx,@gaussfit,beta0)
0023 
0024 azcorrected = az-betanew(2);
0025 rxcorrected = rx + betanew(4).*az + betanew(5);
0026 figure
0027 plot(azcorrected,rxcorrected,'.');
0028 
0029 betanew2(1)=betanew(1);
0030 betanew2(2) = 0;
0031 betanew2(3) = 1;
0032 betanew2(4) = 1;
0033 betanew2(5) = 0;
0034 
0035 [betanew2]=nlinfit(azcorrected,rxcorrected,@gaussfit,betanew2)
0036 
0037 hold
0038 plot(azcorrected,gaussfit(betanew2,azcorrected),'r');
0039 
0040 grid on
0041 
0042 betanew2(3)*1.386