From c5afc09bc49c90b7b00211b1b5faf4050f3c7a76 Mon Sep 17 00:00:00 2001
From: Fredrik Bagge Carlson <fredrikb@control.lth.se>
Date: Sun, 6 Sep 2015 21:27:11 +0200
Subject: [PATCH] dev tf

---
 src/transfer_functions.jl | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)
 create mode 100644 src/transfer_functions.jl

diff --git a/src/transfer_functions.jl b/src/transfer_functions.jl
new file mode 100644
index 0000000..d862c4a
--- /dev/null
+++ b/src/transfer_functions.jl
@@ -0,0 +1,17 @@
+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))
-- 
GitLab