如何在python中找到函数的傅里叶变换量？
在Python中，我们可以使用NumPy库来进行傅里叶变换。傅里叶变换是将一个函数从时域（时间域）转换到频域的过程，它将信号分解为一系列正弦和余弦函数的和，得到信号的频谱表示。
要在Python中找到函数的傅里叶变换量，可以按照以下步骤进行：
- 首先，确保已经安装了NumPy库。可以使用以下命令安装NumPy：
pip install numpy- 导入NumPy库：
import numpy as np- 创建一个表示输入函数的数组。假设我们要对一个函数f(x)进行傅里叶变换，首先需要将函数离散化为一组点（数组），可以使用NumPy的linspace函数生成一组等间隔的点：
x = np.linspace(start, end, num_points) # 根据实际情况替换start、end和num_points的值- 计算函数在这些点上的值：
f_x = f(x) # 根据实际情况替换f(x)为待求函数- 进行傅里叶变换：
fft_values = np.fft.fft(f_x)- 可选步骤：计算变换结果的频率轴：
freqs = np.fft.fftfreq(len(x), d) # 根据实际情况替换d为点之间的时间间隔（采样周期）现在，变量fft_values中存储了函数的傅里叶变换量。这些值表示信号在不同频率上的振幅和相位。你可以根据具体的需求进一步处理这些变量。
这里推荐的腾讯云产品是腾讯云音视频处理（https://cloud.tencent.com/product/vod）和腾讯云人工智能（https://cloud.tencent.com/product/ai）产品。这些产品提供了丰富的音视频处理和人工智能功能，可以与Python的傅里叶变换结合使用，实现更多复杂的应用场景。
相关·内容
1回答
标题似乎有点不清楚，但我要告诉你的整个故事，并希望得到你的意见。我试图在给定的图像中检测到某些特定的边缘，尤其是那些比其他边缘更多的边缘。要做到这一点，最好的思想是傅里叶变换，将离散像素域映射到频域。在应用傅里叶变换之前，我必须测量一些想要的边的距离的平均值，然后用傅里叶变换来找出它们的频率。问题是如何将平均距离设置为FFT算法(在Python或Matlab上)。
谢谢
浏览 0提问于2015-09-11得票数 2
回答已采纳
1回答
我有，它代表一幅图像(模拟的输出)，我想用低通和高通滤波器过滤这幅图像的相位，
My代码(仅限低通滤波器)：
%% Calculate the phase
field_phase = angle(field);
%% Low-pass filter
% Gaussian filter
fieldPhase_Filtered = imgaussfilt(field_phase, 2);
输出：
imagesc(fieldPhase_Filtered); colorbar;
1)我的步骤是否正确？
2)在将相位计算为之前，是否需要进行傅里叶变换？
3)如果我有比矩阵更多的矩阵，
浏览 0提问于2017-09-17得票数 0
1回答
解卷积函数
、、、
我有一个图像(例如(7x7x3)和一个过滤器(3x3x3))。我将图像与过滤器转换成(3x3)输出。如果我想做逆运算，并希望它成为图像输出和过滤器。如何使用Numpy在Python中实现此操作？
我不知道我应该使用哪一个运算与过滤器(逆或转置)？
这是我的反褶积代码：
import numpy as np
def deConv(Z, cashe):
'''
deConv calculate the transpoe Convoultion between the output of the ConvNet and the filter
Arguments:
浏览 0提问于2019-04-06得票数 1
回答已采纳
2回答
非周期信号的傅里叶变换
、、、
我试图计算python中非周期数据的Fourier变换，我的意思是非周期性的。通常我们使用的是定期数据：
x=[0,1,2,3,4,5,6,7]
y=[1,0,2,0,1,0,1,0]
可以很容易地得到傅里叶变换和频域，
from scipy.fft import fft, fftfreq, fftshift
f=fftfreq(len(x), dx)
fy=fft(y)
用它我们可以很容易地用scipy fft计算傅里叶变换，但是在我们的数据不是完全周期性的情况下，例如，
x=[0,1,1.5,1.7,3,4,5...]
y=[...]
任何人都知道如何使用python获得傅里叶变换(
浏览 23提问于2022-04-21得票数 0
1回答
我目前正在研究一个完全卷积的神经网络(图像输入，图像输出)，我正在尝试实现一个丢失函数，在对2幅图像进行一些操作之前对其进行快速傅里叶变换，代码如下所示
def fourierLoss2(y_actual,y_pred):
actual_fft = tf.signal.rfft3d(y_actual)
pred_fft = tf.signal.rfft3d(y_pred)
lossV=tf.math.real(tf.math.reduce_mean(tf.math.square(actual_fft-pred_fft)))
return lossV
with strategy
浏览 0提问于2020-11-26得票数 1
回答已采纳
1回答
为什么经过快速傅里叶变换(FFT)后，我计算的振幅与原来的相距很远？
我有一个1024点的信号，采样频率为1/120000。我用scipy.fftpack在Python中应用了快速傅里叶变换。当我只绘制正数时，我将计算出的震级按桶数归一化，然后乘以2。
当我的初始信号幅度在64 dB左右时，我得到的很低的幅值小于1。
请看我的密码。
Signal = well.ReadWellData(SignalNDB)
y, x = Signal.GetData(numpy=np)
N = y.size # Number of sample points 1024 ...
T = 1/12000
浏览 4提问于2019-12-16得票数 1
2回答
一种检测特定频率的算法？
、、、
使用Goertzel算法，我检测到iphone麦克风中的某个频率(我使用缓冲区样本)。
它的工作原理，但它有许多奇怪的稳定性问题，当值改变。(在同一设备上的频谱应用程序中，它们是恒定的，但在Goertzel算法中则不然)
我想用另一种方法，在C中，检测某个频率，或者某个频率范围内的能量，(我不知道FFT是否良好和准确，如果是，我需要一个好的算法)。如果你有一个函数，它只获取样本和长度，并在一个频率谱或某一已知频率上返回能量，这会有所帮助。我需要一个严肃的过滤器，也许是二级过滤器。
这是我的Goertzel
float goertzel_mag(int16_t* data ,int SAMPLI
浏览 4提问于2013-01-01得票数 0
回答已采纳
2回答
矩形函数的数值傅里叶变换
、、、
这篇文章的目的是正确理解Python或Matlab上的数值傅里叶变换，并举例说明解析傅里叶变换是众所周知的。为此，本文选择了矩形函数，给出了它的解析表达式及其傅里叶变换，并给出了。
这里是Matlab中的代码
x = -3 : 0.01 : 3;
y = zeros(length(x));
y(200:400) = 1;
ffty = fft(y);
ffty = fftshift(ffty);
plot(real(ffty))
这里是Python中的代码
import numpy as np
import matplotlib.pyplot as plt
x = np.arange(-3,
浏览 2提问于2016-11-18得票数 7
回答已采纳
1回答
我把高斯滤波器和我想要过滤的图像转换成傅里叶空间。在傅里叶空间中应用滤波器的规则与在图像空间中的规则相同吗？(例如，仅用np.convolve进行卷积)？
浏览 10提问于2022-05-26得票数 0
回答已采纳
1回答
矩形脉冲的傅里叶逆
、、、、
我很难用Python解释矩形脉冲的逆傅里叶变换的结果。我正在使用来自库numpy.fft的函数irfft。
这是我的代码：
import numpy as np
import matplotlib.pyplot as plt
n_on = 100
n_off = n_on
Y = np.concatenate(( np.ones(n_on), np.zeros(n_off) ))
y = np.fft.irfft(Y)
plt.figure()
plt.subplot(131)
plt.plot(Y, '.')
plt.title("Y")
plt.s
浏览 0提问于2018-04-19得票数 0
回答已采纳
1回答
我目前的知识是，如果我想过滤一个图像，比方说，在一个ppm图像上的蓝色色调，我会循环通过RGB值并增加蓝色值。
然而，像Instagram这样的应用程序有一个滑块，可以用来调整图像过滤器的范围。我假设程序不会一直循环所有RGB值来修改它们。那么，它是如何这么快实现的呢？
浏览 1提问于2021-12-24得票数 2
1回答
当我在我的PC (20个核)上运行Matlab仿真时，它使用了所有的核。有没有办法限制Matlab的CPU可用性，使它只使用10个核心？
浏览 8提问于2022-06-14得票数 1
回答已采纳
2回答
绘制相位和幅值图像傅里叶
、、、、
如何用MATLAB绘制二维图像的傅里叶变换相位和幅值？我使用angle和abs，然后使用imshow，但是我得到了一个黑色的图像。
fftshift在这个绘图中的作用是什么？
浏览 2提问于2014-03-06得票数 2
回答已采纳
fourier(series,k)返回一个矩阵，其中包含傅里叶级数( cos和sin)中的项，直到K(parameter)级。在任何python模块/库中，R中是否有"fourier“的等效函数？如果没有，那么如何在python中对其进行编码？
浏览 12提问于2022-04-11得票数 0
1回答
我有传感器监控泵的实时数据。我需要在频域上分析泵的占空比.我在快速傅立叶变换方面很生疏，所以我希望找到一个C++或Python，它有一个简单的API接收时间序列数据(即一系列{时间戳，值}实数对)，并返回一系列的频率、振幅和相位。
我发现好几个包裹都是复数系列的。此外，他们的示例代码通常使用随机数据作为输入，这使得它很难与生产代码集成。基于Matlab的解决方案对我来说也很难集成到代码中。
我发现的一个包裹是：
https://github.com/flatironinstitute/finufft
我知道问题在于我缺乏知识，但鉴于传感器数据分析是当今的热门话题，我惊讶地发现这样一个图书馆是如
浏览 0提问于2019-09-08得票数 0
1回答
如何在DirectSound中使用C#实现沥青移位？我在谷歌找不到有用的东西。SetFrequency对我不好，因为这也改变了声音的速度。我不想为每个音高创建一个wav文件，因为这将导致一个非常大的软件。
谢谢!
浏览 2提问于2011-09-22得票数 0
2回答
如何获得图像的低、中、高频分量？
、、、、
我目前正在寻找一些图像处理项目，只是想知道如何获得图像的低频，中和高频成分？例如，正如这张图片所显示的(我是从googling获得的，没有详细描述如何获得这张图片，但可能使用了一些过滤)。
另外，我遇到了使用离散余弦变换(DCT)的帖子，它可以帮助我们获得图像的低频和高频分量。只是想知道如何利用DCT来获得中频分量？
我对过滤也有非常基本的知识。我认为也有高斯高通/低通滤波器可用。以及基于小波的滤波。只是想知道高斯滤波、小波滤波和DCT滤波有什么区别？我应该用哪一种？
浏览 3提问于2015-06-17得票数 2
1回答
我一直在阅读关于傅里叶变换及其重要性的文章。为了通过我自己的形象化来获得一种感觉，我尝试了一些东西。
绘制一个信号及其傅里叶变换
绘制了两个不同频率(但振幅相同)的信号及其傅里叶变换。
绘制了3种不同频率(但振幅相同)的信号及其傅里叶变换
这是我得到的。
一个信号及其傅里叶变换：
两个信号及其傅里叶变换：
三个信号及其傅里叶变换：
注:橙色线代表傅里叶变换，蓝线代表正在处理的信号。
我理解傅里叶变换有助于分解信号，给出构成主信号的成分。换句话说，我们得到主信号所包含的不同频率的数目。但我不明白上述的输出。如果我们考虑第一个信号，其中有一个单一的
浏览 2提问于2019-01-25得票数 1
回答已采纳
1回答
FFT实现会产生故障。
、、、
我在FFT图中发现了一个奇怪的白噪声小故障：
我已经检查了参考程序和噪音文件似乎是好的。这是实现中的一个错误吗？
void four1(float data[], int nn, int isign) {
int n, mmax, m, j, istep, i;
float wtemp, wr, wpr, wpi, wi, theta;
float tempr, tempi;
n = nn << 1;
j = 1;
for (int i = 1; i < n; i += 2) {
if (j >
浏览 3提问于2016-01-29得票数 0
回答已采纳
2回答
实现多个嵌套函数
、、、、
我有一个有趣的方程，我想定义它，它有一个傅里叶变换，里面有两个方程。
我不是要解微分方程，我可以用一个正规函数来表示，L。
这个大方程中的其他方程是
现在，你可以，这些是一些重要的功能相关。我对嵌套函数不太熟悉，也不知道如何按顺序定义什么。当我说得有条不紊的时候，我的意思是计算机正确理解的最有效的方法。
我想从定义我的第一个函数L开始，但是我不知道首先实现什么，第一个函数，傅里叶变换，或者傅里叶变换中的函数，然后是傅里叶变换。任何类似的例子和提示都将不胜感激。
谢谢。
编辑:我不是要求任何人为我编写所有代码，我只是想了解什么才是更有效的。是的，有些值是有意义
浏览 6提问于2015-12-04得票数 0
回答已采纳
1回答
是否可以使用一个数字序列作为一个特征？
例如，使用libsvm数据格式：
<label> <index1>:<value1> <index2>:<value2>
+1 1:123.02 2:1.23 3:5.45,2.22,6.76
+1 1:120.12 2:2.23 3:4.98,2.55,4.45
-1 1:199.99 2:2.13 3:4.98,2.22,6.98
...
对于这类数据有什么特殊的机器学习算法吗？
浏览 0提问于2015-08-04得票数 1
回答已采纳
2回答
C的regex库中是否有类似于python中的findall函数：
re.findall(pattern, input_string)
我有一个字符串，其中包含的文件名如下："path/to/file1.ics,path/file2.ics"，其中包含一个未定的文件数。我希望找到所有的文件名(包括它们的路径)，并将它们放在一个字符串数组中。
我正在使用GNU regex.h库
浏览 0提问于2013-11-30得票数 5
回答已采纳
对于这个问题，我使用来自的数学库
我的目的是在对输入值的前向傅里叶变换的结果进行逆傅里叶变换后，得到我的输入。
当我对输入的正向傅里叶变换的结果进行逆傅里叶变换时，得到正确的输出。
我到底做错什么了？
public void fourierTestTemp(){
double[] input = new double[]{1,0,0,0,0,0,0,66,888,0,0,0,0,0,0,0};//Length = 16
double[] result = new double[input.length];//This double array will hold the re
浏览 10提问于2016-04-23得票数 0
回答已采纳
1回答
我在想量子计算机是否真的破坏了施诺尔的签名方案。Schor的算法是通过量子傅里叶变换来实现的，它揭示了周期时间，从而揭示了phi。然而，对于乘法群mod q，一个素数，每个人都知道循环时间，但这不是问题(不像RSA)。问题在于找到一个特定的值。
与RSA不同，在RSA中，必须考虑一个数字，Schnorr的签名依赖于离散对数问题。当你把15插入Schor's时，当你测量状态时得到8。然后你用这个来计算数字。然而，在Schnorr的，你可能得到16当你插入17，但这是没有帮助的。你试图在多项式环生成器2 mod 17中找出哪一项是9，我很确定循环时间没有帮助，我们已经知道了。
我哪里出问题
浏览 0提问于2023-05-02得票数 -6
回答已采纳
上下文
我正在用Python和Numpy实现2D离散傅里叶变换算法。根据图像处理理论，在进行变换之前，图像的每个强度f(x，y)都需要乘以(-1)^(x + y)，其中x和y是二维阵列中表示图像的强度指数。
尝试过什么，
最明显的方法是使用两个for循环迭代每个强度及其指数。
问题
是否有一个使用Python/Numpy矩阵操作的更优雅/更快的解决方案，还是应该使用两个for循环？
浏览 7提问于2019-10-20得票数 0
回答已采纳
1回答
t-SNE的并行版本
、、、、
有没有带有there算法并行版本的Python库？还是存在多核/并行the算法？
我正在尝试使用t降低词汇表中所有word2vecs的维数(300 D Word2vecs 2d)。
问题:词汇量约为130000，要花很长时间才能进行the。
浏览 0提问于2016-11-04得票数 10
回答已采纳
2回答
高斯核函数的FFT错误
、、、、
我正在尝试生成高斯内核的FFT，以便在以后的过滤中使用。我的理解是，高斯核的FFT应该产生一个在视觉上与原始强度图像相似的幅度。这就是我要说的
这是一个可以重现这张图片的代码片段。
import numpy as np
import cv2
from matplotlib import pyplot as plt
ksize = 50
ksize = ksize*2+1
sigma = 15
fil = cv2.getGaussianKernel(ksize,0)
fil = fil * fil.T
fil_ff
浏览 0提问于2018-09-21得票数 1
1回答
我正在学习反向过滤，我试着对它进行编码，我从网上找到了一些参考资料。每个人都考虑了光学传递函数，这在冈萨雷斯书中是看不到的，我指的是。
% Inverse_Filter_Demo-
clc
clear all
close all
original_image=imread('cameraman.jpg'); %loading the original (un-blurred) image
original_image=double(original_image);
%The blur function (PSF- Point Spread Function) tha
浏览 1提问于2012-04-18得票数 2
1回答
我试图在MATLAB中绘制符号变量，我对此使用了与相同的策略。
这是我的代码，它输出一个空白图：
syms t w
x=exp(-t^2)
h=exp(-t)*heaviside(t)+exp(t)*heaviside(-t)
X=fourier(x,w);
H=fourier(h,w);
right=ifourier( rewrite(X*H, 'exp'),t)
fplot(right,[0 8])
我怎样才能使图表出现呢？
浏览 16提问于2022-05-18得票数 1
回答已采纳
1回答
我在用Python做音频分析。我的最终目标是得到一个频率列表和它们各自的体积，如{ frequency : volume (0.0 - 1.0) }。
我把我的音频数据作为在-1.0和+1.0之间有值的帧的列表。我在这个列表上使用了numpy的傅里叶变换- numpy.fftpack.fft()。但由此产生的数据对我来说毫无意义。
我确实理解傅立叶变换从时间到频域的转换，但并不完全是它的数学工作方式。所以我才不太明白结果。
numpy.fftpack.fft()返回的列表中的值意味着什么？我如何处理它/解释它？
如上文所述，在列表上执行的傅里叶变换的最大/分钟值是多少？
如何才能达
浏览 1提问于2014-05-12得票数 2
回答已采纳
我正在构建一个函数来计算给定组件/子系统的可靠性。为此，我用脚本编写了以下内容：
import math as m
import numpy as np
def Reliability (MTBF,time):
failure_param = pow(MTBF,-1)
R = m.exp(-failure_param*time)
return R
该函数对于我在函数中调用的任何时间值都很好。现在我想调用函数来计算给定数组的可靠性，比如np.linspace(0,24,25).但是，我得到了诸如"Type :只有length-1数组可以转换为Python标量“之类的错误。
浏览 4提问于2017-10-05得票数 3
我最近用Python代替了Matlab，因为我是真实的；Y被Sympy所兴奋。
但现在我遇到了以下问题：
我需要一种用Python可靠地计算连续傅里叶变换的方法。Sympy有问题的解决方案，包括Diracs (三角洲-函数)，因为他们发生，例如，三角函数等。
例如，如果我尝试
fourier_transform(cos(x),x,v)
输出是0，它应该基于狄拉克增量函数。
有谁知道，如果Sympy的这一部分要改进，还是有另一种方法可以用Python解析地找到fourier变换？
感谢您的回答或任何建议！
浏览 2提问于2015-10-19得票数 8
回答已采纳
我有两个想要同步的信号(找到时间延迟)。我已经使用了"ccf“函数并找到了最大值，在下面的文章中：
我想我得做点什么
fft1 <- fft(my.vector1)
fft2 <- fft(my.vector2)
ccf(fft1, fft2, lag.max = 6000, plot = FALSE)
然而，出于效率的原因，我想实现与快速傅立叶变换(FFT)的互相关，就像在其他文章中所建议的那样。我对300.000个样本进行了多次测试(在2000赫兹下采样1.5分钟)，最大延迟时间为-3至3秒。
关于如何在R中这样做，有什么提示吗？
我知道fft和ccf函数，但不知道如
浏览 1提问于2014-12-31得票数 2
回答已采纳
5回答
我有一个周期为T的周期函数，想知道如何获得傅立叶系数的列表。我尝试使用numpy中的模块，但它似乎更专注于傅里叶变换而不是级数。也许是因为缺乏数学知识，但我不知道如何从fft计算傅立叶系数。
感谢您的帮助和/或例子。
浏览 29提问于2010-11-24得票数 33
回答已采纳
首先，使用傅里叶变换是否是识别不同说话人的好方法？我不确定如果说的是不同的话，它是否能识别出一个声音。我知道google和amazon在他们的语音助手中都有语音/说话人识别的功能，但是如果傅里叶变换不能成功的话，还有什么好的方法来做到这一点呢？
我想用神经网络识别声音，为了做到这一点，我首先需要为神经网络获得一个很好的输入，但是仅仅把录音作为输入，我不认为它会工作，因为它是基于频率和时间的。所以我找到了傅里叶变换，现在我尝试用傅里叶变换我的音频文件并绘制它。
我的问题是：
如何在python中用音频输入来绘制傅里叶变换？如果这是可行的，我如何在神经网络中输入傅里叶变换(我想也许给每个神经元一个
浏览 0提问于2018-09-19得票数 1
回答已采纳
2回答
在python中，包含一个或多个模块的目录有时会有__init__.py，因此可以将该目录视为python package，对吗？__init__有什么不同？(另一个问题是，python模块只是一个python代码文件，其中包含相关的(可能是独立于其他文件的)一组类、函数和变量？)
浏览 1提问于2011-10-18得票数 0
回答已采纳
1回答
你如何估计有噪声的正弦波的频率？
、、、、
我得到了一个数据点(时间和振幅)的列表，我用matplotlib绘制了一条噪声正弦曲线。
现在我必须估计出噪声波的频率，我认为我应该使用numpy's函数(np.fft)的傅里叶变换来实现这一点。
我不知道我应该做什么来执行傅里叶变换，一旦我做了，我如何确定哪个频率占主导地位。
浏览 0提问于2019-07-31得票数 0
2回答
MATLAB未能计算此函数的fft：
syms t;
x= 1/(1+t^2);
X= fft(x)
并抛出这个错误：
未定义函数'fft‘类型的输入参数'sym’。
为什么不采用符号变量的傅里叶变换？我认为应该这样做，因为我们可能想作为一个象征性的值来获得结果。
我还试图用非符号变量来解决这个问题。
t= -10:0.01:10;
x= zeros(2001);
x= 1/(1+t.^2);
fft(x)
这次我的错误是：
矩阵尺寸必须一致。
然而，它们具有相同的维度。我的错在哪里？
浏览 12提问于2016-02-06得票数 1
2回答
我试图为Python中的数组定义一个傅里叶变换。。
但是，输出始终是空数组。你能帮我纠正一下我的错误吗？
import numpy as np
def TF(T):
N = len(T)
Tr = np.array([])
for k in range(0,N-1):
Tl = 0
for l in range(0,N-1):
Tk += (T[l])*np.exp((-2j*np.pi*k*l)/N)
np.append(Tr,Tk)
print('Tr =',Tr)
T
浏览 5提问于2019-12-13得票数 2
回答已采纳
3回答
我要做的是，从具有周期模式的x-y点的列表中，计算周期。以我有限的数学知识，我知道傅立叶变换可以做这类事情。
我在写Python代码。
我找到了一个相关的答案，但它使用了一个均匀分布的x轴，即dt是固定的，这对我来说不是这样的。因为我不太理解它背后的数学原理，所以我不确定它在我的代码中是否能正常工作。
我的问题是，它有用吗？或者，在numpy中是否有某种方法已经完成了我的工作？或者，我该怎么做？
编辑：所有的值都是Pythonic float (即双精度)
浏览 6提问于2015-12-23得票数 10
回答已采纳
2回答
STFT中窗口函数的工作原理
、、、、
在信号处理和STFT方面有经验的人能向我解释为什么下面张贴的图像中的窗口函数来自(t- t‘)，因为t是总时间，t’是窗口的宽度？
我无法理解，因为最初，窗口位于t=0，如果一个示例的窗口长度为3，则该窗口将从t=0 -> t=3中跨出，如果一个示例的总时间T= 10，则窗口函数将类似于w(T-3)，即7?！我真的不能理解它，我相信有什么隐藏的东西我无法理解。
请澄清并提供指导。谢谢图像：
浏览 2提问于2014-12-30得票数 1
回答已采纳
1回答
我需要读取音频流并确定生成的超声波如何从音频流中找到特定的声音序列？
起初我认为是DTMF的方向，但后来拒绝了它，因为它是人耳所能听到的。
如果你有任何其他的想法，我很乐意倾听。
浏览 0提问于2011-04-16得票数 1
2回答
我正在绘制正弦波(左列)及其各自的频域表示(右列)：
第一波(振幅: 10；频率: 0.5)的fft表示有点混乱
第二波(振幅: 15；频率: 5.0)看上去与预期完全一致。
第三波是第一波，第二波是总结和继承的问题。
第二频率图在x=5 (频率)、y=15 (振幅)处正好有一个峰值。
为什么在只有一个频率的情况下，第一个频率图有多个峰值？
import numpy as np
import matplotlib.pyplot as plt
def sine(freq, time_interval, rate, amp=1):
w = 2. * np.p
浏览 3提问于2016-03-29得票数 4
我想在一个动态系统上进行集成，就像ODE一样。
x_ddot + d*x_dot + k*x = a*sin(omega*t)
通过一项修改:外力a* sin(omega * t)必须由另一个周期信号代替，这是基于一些测量数据的。这意味着:我需要对我的离散测量数据做曲线拟合，得到的函数必须是周期性的。我有两个解决问题的办法：
1)采用傅里叶变换(numpy.fft)。但是，利用离散傅里叶变换从离散数据中生成连续函数是很困难的。
2)使用曲线拟合，具有类似于a1+a2*sin(omega*t)+a3*sin(2*omega*t)+a4*sin(....的函数其中omega为2*pi/(测量数据的
浏览 19提问于2016-06-29得票数 2
回答已采纳
我对matlab和信号处理非常陌生，但我正在做一个项目，而且我被困住了。我的目标是做一个带通快速傅里叶变换(FFT)滤波器，切断频率为0.5赫兹和3赫兹的实际信号。我的输入信号是一个信号采样ad 405 ad和约35000个样本。我的实现如下所示：
%ZFiltered = My signal
filtered = bandpass(fft(zFiltered), [0.5 3], sampling_frequency);
res = ifft(filtered);
subplot(2,1,1)
plot(zFiltered)
subplot(2,1,2)
plot(abs(res))
浏览 1提问于2021-03-05得票数 0
回答已采纳
1回答
我有一个MATLAB程序，我想要移植到Python。问题是，在它中，我使用内置的spectrogram函数，尽管matplotlib specgram函数看起来是相同的，但在运行这两个函数时得到的结果是不同的。
这是我一直在运行的代码。
MATLAB：
data = 1:999; %Dummy data. Just for testing.
Fs = 8000; % All the songs we'll be working on will be sampled at an 8KHz rate
tWindow = 64e-3; % The window must be long
浏览 1提问于2015-06-28得票数 5
回答已采纳
1回答
在我的情节里，我想要不同颜色的涟漪。所以你能告诉我如何添加这个可能的三维图形的涟漪吗？我的代码如下所示。这是我的形象
我想增加波纹表面的轴波纹:样本图像：
import scipy as np
from mayavi import mlab
x,y = np.mgrid[-3:3:100j, -3:3:100j]
a = 0.3
a0 = 0.2
r1=np.sqrt((x-a)**2 + y**2)
values = np.exp(-r1/a0)
mlab.surf(values, warp_scale='auto')
mlab.outline()
浏览 0提问于2018-09-26得票数 2
回答已采纳
1回答
创建输入图像的谱直方图
、、、、
我正在尝试用Python创建单个输入图像的谱直方图。
这意味着一个直方图，它捕获图像中的空间频率，并在图像中显示它们的分布。
例如，对于一个值(例如全黑色)的图像，我们期望在DC (零)频率上有一个较高的值，而没有任何其他值。
更大的图片是，我有一个降噪算法，作为平滑滤波器的图像。我可以从视觉上看到，它创造了一幅“更平滑-不那么详细”的图像，但我想用一个直方图来图形化地捕捉它，它可以证明使用该算法从图像中过滤出更高的频率。
我从图像上的FFT变换开始，很快就不知道如何以正确的方式导出这样的直方图。
我开始使用的一个示例是以下代码：
import scipy.fftpack as fp
im2
浏览 0提问于2019-08-25得票数 0
我对Python和FFT函数比较陌生。我试图使用numpy.fft.fft()函数将正方形脉冲(一维衍射狭缝函数)转换为正弦函数(一维衍射图)，并使输出图与由方程给出的平方脉冲的解析变换相同：
F(u) = sin(πau)/(πu)
但是，我很难格式化FFT，使它看起来与分析功能相同。这是我的最低工作代码：
a = 40.0 # slit width
N = 5000 # sample points
T = 100.0 # sample spacing
x = np.linspace(-T/2.0,T/2.0,N)
y = np.piecewise(x,[abs(x)>
浏览 5提问于2015-08-03得票数 3
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
