Movatterモバイル変換

[0]ホーム
URL:

Skip to content

Navigation Menu

Sign in
Appearance settings

Search code, repositories, users, issues, pull requests...

Provide feedback

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly

 Sign up
Appearance settings

Frequency response of Downsampled transfer function in Z domain#1126

roccafrancesco started this conversation inGeneral
Discussion options

Hello Everybody
I wrote a simple function that convert a z transfer function(num,den) in frequency response object when the outputl is decimated M fold.

#num: transfer function numerator
#den: transfer function denominator
#M: M fold decimator
#Tc: sample time
#points: points of frequency response

`def` downSampling(num,den,M,Tc,points):    omega = np.logspace(np.log10(0.0001),np.log10(0.99*np.pi/(Tc)),num=points,endpoint=True)    denFreq=np.zeros(omega.size,dtype=complex)    numFreq=np.zeros(omega.size,dtype=complex)     bode=np.zeros(omega.size,dtype=complex)     M1=float(M)        for k in range(M):                zeta=np.exp(1.0j*omega*Tc/M)*np.exp(-2.0j*np.pi*k*Tc/M)        denFreq=denFreq+np.polyval(den, zeta)        numFreq=numFreq+np.polyval(num, zeta)          bode=bode+(numFreq/denFreq)    bode=bode/M           sys=ct.FrequencyResponseData(bode,omega,smooth=True)    return sys

Im not really sure if i did everithig well.
Could someone check and correct eventually the function?
Thank you

You must be logged in to vote

Replies: 0 comments

Sign up for freeto join this conversation on GitHub. Already have an account?Sign in to comment
Category
General
Labels
None yet
1 participant
@roccafrancesco
[8]ページ先頭
©2009-2025 Movatter.jp