我想通过反向快速傅立叶变换找到时间置换曲线

Question

import numpy as np
import matplotlib.pyplot as plt

N=1000; T=10; dt=T/N
t=np.arange(0,N*dt,dt)

#system properties
wn = 4 * np.pi
w=0.2*wn
Tn = 2*np.pi/wn
fn = 1/Tn

m = 1
Z = 0.05
k = wn*wn * m
c = 2*m*wn*Z
######################################################
P = np.sin(wn * t)

#Frequency response
L = len(P)
a = np.arange(L)
freq = a * (1/dt) / L
freq = freq[range(int(L/2))]

Pf = np.fft.fft(P)/L
Pf = Pf[range(int(L/2))]
FFT = 2*abs(Pf)
#########################################################
# transfer function
H=(1/k)*(1/(1+(w/wn)**2)+1j*(2*Z*(w/wn)))

Uf=2*Pf*H

# for i in range(0,int(N/2)):
#     Uf[N-i-1] = np.conjugate(Uf[i+1])

ut= np.real(np.fft.ifft(Uf)*(1/dt))

plt.plot(t, ut)
plt.show()

我不知道这个 ValueError ValueError：x 和 y 必须具有相同的第一维度，但具有形状 (1000,) 和 (500,) 我想通过快速傅里叶逆变换找到时间-位移曲线

Answer 1

plt.plot（）功能只是抱怨t和ut没有相同的长度。

此行pf = pf [range（int（l/2））]将pf从长度n到n /2。这会影响uf，因此ut，这就是为什么您会收到错误消息的原因。我不确定您在此命令中要做什么...您是否试图过滤某些频率？无论哪种方式，再看一下此命令都可能给您想要的东西。