allanplot.m: add detrend option

bmarechal
1 parent 6ce4bbcc4e
Showing 1 changed file with 4 additions and 0 deletions Side-by-side Diff
allanplot.m
@@ -11,6 +11,10 @@
     cc = 'bkcgmry';
     for i = [1:length(col)]
         data.freq = load(filename)(:,col(i)).*mult(i);
+        if eval(argv(){4})(i) == 1
+            printf(strcat(filename, ' col', num2str(col(i)), ' drift removed\n\n'))
+            data.freq = detrend(data.freq);
+        end
         data.rate = 1;
         [ad, S, err, tau] = allan(data, 2.^[0:nextpow2(length(data.freq))-3]./data.rate, strcat(strsplit(filename, '/'){end}, num2str(i)), 0);
         loglogerr(tau, ad, err, strcat(cc(mod(i, length(cc))), '-s'))
...	...	@@ -11,6 +11,10 @@
11	11	cc = 'bkcgmry';
12	12	for i = [1:length(col)]
13	13	data.freq = load(filename)(:,col(i)).*mult(i);
	14	+ if eval(argv(){4})(i) == 1
	15	+ printf(strcat(filename, ' col', num2str(col(i)), ' drift removed\n\n'))
	16	+ data.freq = detrend(data.freq);
	17	+ end
14	18	data.rate = 1;
15	19	[ad, S, err, tau] = allan(data, 2.^[0:nextpow2(length(data.freq))-3]./data.rate, strcat(strsplit(filename, '/'){end}, num2str(i)), 0);
16	20	loglogerr(tau, ad, err, strcat(cc(mod(i, length(cc))), '-s'))