Python中文
首页
教程
问答
标签
搜索
登录
注册
重采样,插值矩阵
回答此问题可获得
20
贡献值,回答如果被采纳可获得
50
分。
<p>我正试图插入一些数据以便绘图。例如,给定N个数据点,我希望能够生成一个“平滑”图,由10*N左右的插值数据点组成。</p> <p>我的方法是生成一个N乘10*N的矩阵,并计算原始向量和我生成的矩阵的内积,得到一个1乘10*N的向量。我已经算出了我想用来插值的数学公式,但是我的代码相当慢。我对Python还不太熟悉,所以我希望这里的一些专家能给我一些想法,让我试着加速我的代码。</p> <p>我认为问题的一部分是生成矩阵需要10*N^2调用以下函数:</p> <pre><code>def sinc(x): import math try: return math.sin(math.pi * x) / (math.pi * x) except ZeroDivisionError: return 1.0 </code></pre> <p>(这个<a href="http://en.wikipedia.org/wiki/Whittaker%E2%80%93Shannon_interpolation_formula" rel="noreferrer">comes from sampling theory</a>。本质上,我试图从它的样本中重新创建一个信号,并将其放大到更高的频率。)</p> <p>矩阵由以下各项生成:</p> <pre><code>def resampleMatrix(Tso, Tsf, o, f): from numpy import array as npar retval = [] for i in range(f): retval.<a href="https://www.cnpython.com/list/append" class="inner-link">append</a>([sinc((Tsf*i - Tso*j)/Tso) for j in range(o)]) return npar(retval) </code></pre> <p>我正在考虑把任务分解成更小的部分,因为我不喜欢记忆中的N^2矩阵。我可能会把‘重采样矩阵’变成一个生成器函数,然后一行一行地执行内部产品,但我认为在开始对内存中的内容进行分页之前,这不会大大加快代码的速度。</p> <p>提前谢谢你的建议!</p>
0 条评论
分类:
Python问答
请先
登录
后评论
默认排序
时间排序
1 个回答
匿名
1天前
擅长:python、mysql、java
<p>如果您想以一种非常通用和快速的方式插值数据,样条曲线或多项式是非常有用的。Scipy有Scipy.interpolate模块,非常有用。你可以在官方网页上找到<a href="http://docs.scipy.org/doc/scipy/reference/tutorial/interpolate.html" rel="nofollow noreferrer">many examples</a>。</p>
请先
登录
后评论
针对此问题:
更多的回答
关注
89
关注
收藏
1
收藏,
216
浏览
网友 提问于 2天前
相关Python问题
为什么在使用strptime时会出现未进行转换的数据错误?
9 回答
为什么在使用strptim时会出现这个datetime日期错误
1 回答
为什么在使用StyleFrame时索引列的标题不显示sf.至excel()?
1 回答
为什么在使用sum()函数时会发生“int”对象不可调用的错误?
2 回答
为什么在使用sympy.dsolve时会得到“'list'对象没有属性'func'”?
5 回答
为什么在使用tabla时会得到一个空的数据帧?
1 回答
为什么在使用tensorboard时需要add_graph()的第二个参数?
9 回答
为什么在使用TensorFlow Lite转换YOLOv4时,推断时间/大小没有改进?有什么可能的改进吗?
8 回答
为什么在使用Tensorflow加载训练批时会出现内存泄漏?
3 回答
为什么在使用tensorflow时会收到警告/错误(使用函数API,但未实现错误)
5 回答
为什么在使用tetpyclient发出POST请求时出现403错误?
9 回答
为什么在使用TextBlob时会出现HTTP错误?
1 回答
为什么在使用TFIDF时出现错误“IndexError:list index out of range”pyspark.ml.feature?
6 回答
为什么在使用timedelta格式化之后,我在python中的日期是错误的?
9 回答
为什么在使用timeit或exec函数时,函数中的变量不会在提供的全局命名空间中搜索?
5 回答
为什么在使用tkinter时不能使用复选框?
2 回答
为什么在使用todoistpythonapi时会返回这个奇怪的ID?
6 回答
为什么在使用TQM时,在调整图像大小时,处理时间会有很大的差异?
10 回答
为什么在使用Tweepy下载用户时间线时收到错误消息
6 回答
为什么在使用twitter帐户登录Django应用程序时重定向127.0.0.1:8000?
5 回答