freqinfo.m 387 Bytes
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
figure;
w = linspace(0.1, 1000, 10000);
[sv,w] = sigma(LQc,w);
svn = sv/max(sv);
larger = find(svn>0.71); % -3dB
smaller = find(svn<=0.71);
bf = [w(larger(end)) w(smaller(1))];
bsv = [svn(larger(end)) svn(smaller(1))];
step=(bsv(1)-0.71)/(bsv(1)-bsv(2)); 
bw = bf(1)+(bf(2)-bf(1))*step
sample_h=[0.15 0.6]/bw
sample_f=1./sample_h
%csvwrite('sigma.csv', [w svn']);
loglog(w, svn);
grid()