小编BBD*_*Sys的帖子

对于已知具有季节性或日常模式的数据,我希望使用傅立叶分析来进行预测.在运行时间序列数据fft后,我获得系数.如何使用这些系数进行预测？

我相信FFT假设它接收的所有数据构成一个周期,那么,如果我只是使用ifft重新生成数据,我也在重新生成我的函数的延续,那么我可以将这些值用于未来的值吗？

简单地说:我运行fft为t = 0,1,2,.. 10然后在coef上使用ifft,我可以使用重新生成的时间序列为t = 11,12,... 20？

当我调用rnorm单个值作为均值时,很明显会发生什么:从Normal(10,1)生成一个值.

y <- rnorm(20, mean=10, sd=1)

但是,我看到传递给整个矢量的例子rnorm(或rcauchy等等); 在这种情况下,我不确定R机器到底做了什么.例如:

a = c(10,22,33,44,5,10,30,22,100,45,97)
y <- rnorm(a, mean=a, sd=1)

有任何想法吗？

Emacs Mx编译没有看到.bashrc中设置的任何别名.如果我使用Mx shell然后键入别名,它没关系.我尝试从/ etc/profile,从〜/ .profile,〜/ bash_env采购.bashrc,我能想到的任何东西都无济于事.

我在Emacs 23和Ubuntu 11上.我使用/ usr/bin/emacs%F从桌面按钮启动emacs.

我试图将包含\运算符的代码从Matlab(Octave)转换为Python.示例代码

B = [2;4]
b = [4;4]
B \ b

这有效并产生1.2作为答案.使用此网页

http://mathesaurus.sourceforge.net/matlab-numpy.html

我翻译为:

import numpy as np
import numpy.linalg as lin
B = np.array([[2],[4]])
b = np.array([[4],[4]])
print lin.solve(B,b)

这给了我一个错误:

numpy.linalg.linalg.LinAlgError: Array must be square

为什么Matlab \与B的非方矩阵一起工作？

对此有何解决方案？

我怎样才能genfromtxt返回list自动检索的列名names=True？当我做:

data = np.genfromtxt("test.csv",names=True,delimiter=",",dtype=None)
print data['col1']

它打印整个列值col1.

但是,我需要遍历所有列名称.我怎样才能做到这一点？

我尝试了data.keys()各种其他方法,但无论返回什么genfromtxt似乎都不是字典兼容的对象.我想我可以list自己传递列名,但从长远来看,这对我来说是不可维护的.

有任何想法吗？

我想将我的数据帧转换为矩阵,将单个因子列扩展为多个,并根据因子分配1/ 0.例如

C1 C2 C3
A  3  5
B  3  4
A  1  1

应该变成类似的东西

C1_A C1_B C2 C3
1      0  3  5
0      1  3  4
1      0  1  1

我怎么能在R中这样做？我试过了data.matrix,as.matrix但没有回复我想要的东西.它们为单个因子列分配"整数"值,没有扩展.

有没有人有绘制椭圆体的示例代码？matplotlib现场有一个球体,但椭圆体没有.我想策划

x**2 + 2*y**2 + 2*z**2 = c

其中c是一个定义椭球的常数(如10).我尝试了这meshgrid(x,y)条路线,重新设计了等式,所以z就在一边,但这sqrt是一个问题.该matplotlib球体例子使用角度u,v,但我不知道如何工作,对于椭球.

有没有人知道Matlab/Octave bwdist()函数的Python替换？对于给定矩阵,此函数将每个单元的欧几里德距离返回到最接近的非零单元.我看到了一个Octave C实现,一个纯Matlab实现,我想知道是否有人必须在ANSI C(它不包括任何Matlab/Octave头,因此我可以很容易地从Python集成)或纯Python中实现它.

我提到的两个链接如下:

C++

Matlab M-File

作为测试,Matlab代码/输出看起来像这样:

bw= [0   1   0   0   0;
     1   0   0   0   0;
     0   0   0   0   1;
     0   0   0   0   0;
     0   0   1   0   0]

D = bwdist(bw)

D =

   1.00000   0.00000   1.00000   2.00000   2.00000
   0.00000   1.00000   1.41421   1.41421   1.00000
   1.00000   1.41421   2.00000   1.00000   0.00000
   2.00000   1.41421   1.00000   1.41421   1.00000
   2.00000   1.00000   0.00000   1.00000   2.00000

我在Python中测试了一个推荐的distance_transform_edt调用,它给出了这个结果:

从scipy import ndimage导入numpy as np

a = np.array(([0,1,0,0,0],
              [1,0,0,0,0],
              [0,0,0,0,1],
              [0,0,0,0,0],
              [0,0,1,0,0]))

res = …

我正在研究光流,基于这里的讲义和互联网上的一些示例,我编写了这个Python代码.

所有代码和示例图像也都存在.对于大约4-5个像素的小位移,计算的矢量方向似乎很好,但矢量的幅度太小(这就是为什么我必须在绘制它们之前将u乘以v乘以3).

这是因为算法的限制,还是代码中的错误？上面分享的讲义也说动作需要很小"你,v小于1个像素",也许这就是原因.这种限制的原因是什么？

我想用PIL向用户显示图像,当用户点击此图像的任何位置时,我想要调用def onmousedown(x,y).我会在这个函数中做一些额外的事情.我怎么能在PIL中这样做？

谢谢,