From 473f7e4ca011d36e4bb1291a8dd8e517299f5962 Mon Sep 17 00:00:00 2001 From: giulio Date: Sat, 30 Oct 2021 16:11:17 +0000 Subject: [PATCH] Added warning for time --- dist/tfest-0.1.11.tar.gz | Bin 0 -> 3075 bytes tfest.egg-info/PKG-INFO | 4 ++-- tfest/tfest.py | 6 +++++- 3 files changed, 7 insertions(+), 3 deletions(-) create mode 100644 dist/tfest-0.1.11.tar.gz diff --git a/dist/tfest-0.1.11.tar.gz b/dist/tfest-0.1.11.tar.gz new file mode 100644 index 0000000000000000000000000000000000000000..d262d66ecbdc3b800c402e7224194e43e4accd10 GIT binary patch literal 3075 zcmV+e4E*ySiwFq*ZGB+^|72-%bT4#fWpi{bFfK7JF)=Q5VR8WN8*6XdI5M9F=m)%xU6Nww-QQc_v8q`!@km9)I;=lli<&LyvX7-S`jt!+{z9Z$17o`7@z0 z;kHwXS{{Mmb-*=qBn*Sdi4|`_*?*rZtfn6+c`!^DE%EB?~*^c!V2V{(tIaql? zY9!Wq5(GG2jnPGnL++6nc}X}EK!Sjrf+dc zc@ytepatLyf%+EvBov5w=*`K+>yxK?_%h1;<>=@Ly$r$XeP8utoR|Fm`|o@~d_PA- zG(t2bUE>C%n>oK&;vn)#!SsTmIS~jW6v1BDehVEp%OxPoG?QVP$QICK$@$c%a`RSs$-2<&nW#txP#QJe-zYI&$ZL*O)Hp$o&{YO$c7PIHHfD&y1A zy_l#jv{aiDr>eVGWCc^o<=UyWa&cO~#>dyzJLnS3u3(HlynXfd>|bxKXFlefdX#|g zc}{ML&!Rw1uycV$!qIq)`lyXgq&KM0z6=EkU2uLUiA$jgB1ts^aHSR63l`CtEI((S zxWh52v;_F43H4ndV}7-w_HB~n?EED{LKqHa4m4@No3EwxqgFHV44a|r!=HuIxR4% zon)W)7Yx>RLs~C#NhG}`c;cx~!yCSdro0caUM5&H^`+1JY%4YMzw?N(`6~NAU1x3n z_m2C0GynGiCT5%Yzis}1Dc`z#{twsYzvDP&{_k)6jT12QU(WwL1^lvY?3(|Mvv&UD z4En>Nng9EM%Rp}cUs-NxmG4ns?w=K_1M3pP!c>NZD@)+{^jd99Opy9;vD$6xTE>z5 zt>^!jI`H}I^Pl74v5EiY{KuUCnDZZV{$tL6%=wQw|1q$y^A<0z*U8fTE6@PuUd=Y0q%^s;6Ns|M0o9=YfXe9DOpHwozg^Lig%7zS5d1dl|I%rMzNmd!7bxj-w~&wb2|BZwqOuF2`B`Nq?Qv%J3oex zr?4ujAsDsV0KOn{x`{HQ4lup;Nllu)%w=Xt9Lu&qMo7MxN}%G*qSG$zNC&i#^-swJ zvO14T0y64H#D7JqbepDJ#i-R(A^@h|G$HdlKmYmi{{s5-|9>10M}71Ek8d>o+3k)! z<$q3hP#XV%#5=>j$^Z76{}q8423~DiK!l5YVd-)ekR>jREX#{oP_7X`w9ocxBJQ@( zjJqk33feYs!Ukc``>O=&v|9Fa!bX(6ZCLU>8}#A(_*2^vjw^hC9M%^ORNDn>Y0=QLBN9`9_6 zYvB|-T$=|xr5x)GIhn&k9n{d7E$g9^|I~3 z$t3N~xRlVtnr79lNCN>|m=RdcOw$khI=oXb;B4q(WXr8MBPd9}u>AYSbv zF3)f6Y9y7-!C|$KJCTa&_^>+Goo+tN?o_p@+wbJv&S{8WRTGR=RZBh8-msB zhAd|bNZ^TD$w6niT$=FC$oUy!cmj)!7|24oPH)2^rM-YvTA1_ngA>=PZ)v~TLcj3)@>f`?EoH@O3M`_G#bc^UO=UglWPG;kskf%59YQ2ph}B;g zV)aNNvg)0zOC=#VsrEX?6z0;4<&3CpgGTK`coV|%U5i?4ZVx{9Ctn_cPS+h(T2fL? z;c@b<^wiH}K_%3Fun0vHaAr*Pl286qUxFF=^)Yfty9Y~t7d&5G9tkGozlF2XM+>_= zR@zoM{Crw?D*8B9?-ctejh@n>>Fh)8{p3}&Ami*Mb@r8P#m>`|?VdFg@ROYBmS*fT zjd@QvIO+#9g`UW)psJlJu9|g(RdrTVkFcuC`}s%lQHZka%I%+`aWm<2P z`9fON&8@c{5+7V?J0RGt7pAX&n;WMX@GfT7=A%MMw`Cq=26FF=bX#~WO*r{liQy7LYlvER-6sk|{52s@`~gh4lC(mTKuU7Fg@`r5UKPmJ2L4 zSL-TRM)~-f_>)?3lvCSYCCzw>N-*ruki6-%R>>!G61oka%_g-Zioh^pI$wy|BklA% z-2W~8-)lejf4hBs|F1jf$&J6RGaL=g{lC3{+Gs0Xld9agGRGr7*7;9<-oMxTA5K4; z|4w(<9m)AW=ypeD{_g|yoz+9@>s?kow07Q7J+uz3d-qZgt?WYTp>^*jsw&K?WtUDt zie0x&RlW3d`_L-;I^dLFES2y1jnYGF^Howww0?(F(ym<|Rj4gFzZPQaJ<&>Mt2aY+ z)gR(6=(-vnbNf@(G49=f0RsjM7%*VKfB^#r3>YwAz<>b*1`HT5V8DO@0|pEjFkrxd R0R!JL{11dYNWcJi001VG7v literal 0 HcmV?d00001 diff --git a/tfest.egg-info/PKG-INFO b/tfest.egg-info/PKG-INFO index 436eb88..7823c98 100644 --- a/tfest.egg-info/PKG-INFO +++ b/tfest.egg-info/PKG-INFO @@ -1,12 +1,12 @@ Metadata-Version: 2.1 Name: tfest -Version: 0.1.10 +Version: 0.1.11 Summary: Transfer function estimation based on frequency response. Home-page: https://github.com/giuliovv/tfest Author: Giulio Vaccari Author-email: io@giuliovaccari.it License: MIT -Download-URL: https://github.com/giuliovv/tfest/archive/refs/tags/v0.1.10-alpha.tar.gz +Download-URL: https://github.com/giuliovv/tfest/archive/refs/tags/v0.1.11-alpha.tar.gz Description: # Tfest Transfer function estimation with frequency response. Only works for SISO systems for now. Python equivalent of MATLAB tfest diff --git a/tfest/tfest.py b/tfest/tfest.py index 234d619..9a3a5ac 100644 --- a/tfest/tfest.py +++ b/tfest/tfest.py @@ -1,5 +1,6 @@ import matplotlib.pyplot as plt import numpy as np +import warnings from matplotlib.mlab import psd, csd from scipy.optimize import minimize @@ -33,13 +34,16 @@ def loss(self, x, nzeros, freq, H): risp /= np.array([np.polyval(poles, s) for s in 1j*freq]) return np.linalg.norm((risp-H).reshape(-1, 1), axis=1).sum() #+ np.abs(x).sum() - def frequency_response(self, method="h1", time=1): + def frequency_response(self, method="h1", time=None): """ method: "fft" or "h1" (default) or "h2" time: time for fft return: frequency response and frequency """ + if time == None: + warnings.warn("Setting default time=1") + time = 1 if method == "fft": u_f = fft(self.u) y_f = fft(self.y)