diff --git a/psdplot.m b/psdplot.m index b51dab4..564020e 100755 --- a/psdplot.m +++ b/psdplot.m @@ -3,6 +3,7 @@ filename = argv(){1}; col = eval(argv(){2}); mult = eval(argv(){3}); +data.rate = eval(argv(){4}); if length(col) == length(mult) figure @@ -11,7 +12,6 @@ if length(col) == length(mult) cc = 'bkcgmry'; for i = [1:length(col)] data.freq = load(filename)(:,col(i)).*mult(i); - data.rate = 1; [p, f] = pwelch(data.freq, [], 0.95, [], data.rate.*mult(i), 'onesided'); semilogx(f, 10*log10(p), cc(mod(i, length(cc)))) leg{i} = strcat(filename, ' col', num2str(col(i)));