0

根据要求的完整测试代码：

import numpy as np
import matplotlib.pyplot as plt

# Numtopad determines number of padding points
numtopad  = 512

# Define axis
x = np.arange(numtopad)
y = x[:,np.newaxis]

# Offsets which are zero
x0 = 256*0
y0 = 256*0

# Exponentially decaying function in 2D
f = np.exp( -((y-y0) + (x-x0))/(10))

# Fourier transform above function and move zero frqeuencies to center of graph
f2 = np.fft.fft(f,axis=0)
f2 = np.fft.fft(f2,axis=1)
f2 = np.fft.fftshift(f2,axes=0)
f2 = np.fft.fftshift(f2,axes=1)

Delta_t = x[1]-x[0]

# Define a frequency
freq_t = np.fft.fftfreq(numtopad,d = Delta_t)
freq_offset = 200

E1 = freq_t + freq_offset
E2 = freq_t + freq_offset

# plt.contourf(abs(f2))
plt.contourf(E1,E2,abs(f2))
4

1 回答 1

0

您能否提供完整的代码，因为图片不可用，以确保我明白目的？

如果我正确理解了您的问题，则您的数组 E1 和 E2 以 0 为中心：[-0.5,...,0.5] 而函数 f 是以 256 为中心的高斯。您应该更改函数 f 以正确放置尊重到您的数组 E1 和 E2 或标准化您的数组 X，Y：

import numpy as np
import matplotlib.pyplot as plt

numtopad = 512
x = np.arange(numtopad)
y = x[:,np.newaxis]

x0 = 256
y0 = 256
f = exp( -((y-y0)**2 + (x-x0)**2)/9000)

X,Y = numpy.meshgrid(x,y)
X = ((X)/512.)-0.5
Y = ((Y)/512.)-0.5

fig = plt.figure()
axe = fig.add_subplot(111)
axe.contourf(X, Y, abs(f))
fig.show()

如果您只想重新调整数据，您可以使用此代码，甚至可以使用此代码生成 x 和 y（但您必须更改 f）：

x = numpy.linspace(-0.5,0.5,512) #512p linearly spaced [-0.5,0.5]
X,Y = numpy.meshgrid(x,x) #2-D meshgrid on the box [-0.5,0.5]
f = exp( -((X)**2 + (Y)**2))
axe.contourf(X, Y, f)
于 2015-06-11T07:52:47.747 回答