bmarechal / octave_scripts

Download as
Browse Code ยป

Commit 2f13b73569577e6baf7085cc4b9b9ac366408254

Authored by bmarechal
1 parent 45d0ff1fae
Exists in master

add multiple column option

Showing 1 changed file with 16 additions and 5 deletions Inline Diff

#!/usr/bin/octave-cli --persist 1 1 #!/usr/bin/octave-cli --persist
2 2
filename = argv(){1}; 3 3 filename = argv(){1};
col = eval(argv(){2}); 4 4 col = eval(argv(){2});
mult = eval(argv(){3}); 5 5 mult = eval(argv(){3});
6 6
data.freq = load(filename)(:,col).*mult; 7 7 if length(col) == length(mult)
data.rate = 1; 8 8 figure
9 9 hold all
ad = allan(data, 2.^(0:nextpow2(length(data.freq))), strsplit(filename, '/'){end}, 1); 10 10 grid on
input("Press to continue..."); 11 11 cc = 'bkcgmry';
12 for i = [1:length(col)]
13 data.freq = load(filename)(:,col(i)).*mult(i);
14 data.rate = 1;
15 [ad, S, err, tau] = allan(data, 2.^[0:nextpow2(length(data.freq))-3]./data.rate, strcat(strsplit(filename, '/'){end}, num2str(i)), 0);
16 loglogerr(tau, ad, err, cc(mod(i, length(cc))))
17 leg{i} = strcat(filename, ' col', num2str(col(i)));
18 hold on
19 end
20 legend(leg)
21 input("Press to continue...");
22 end
exit 12 23 exit
13 24