Commit c5afc09b authored by Fredrik Bagge Carlson's avatar Fredrik Bagge Carlson
Browse files

dev tf

parent e5371d18
using DSP
function tfest(y,u)
Cyu = xcorr(y,u)
Cuu = xcorr(u,u)
Syu = welch_pgram(Cyu)
Suu = welch_pgram(Cuu)
return H = TF(Syu.power./Suu.power,Syu.freq)
end
N = 200000;
u = randn(N);
y = filt(ones(30),30,u);
H = tfest(y,u);
loglog(H.f,abs(H.P))
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment