请教一个傅立叶变换的问题 - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>EE - 电子工程

请教一个傅立叶变换的问题

请教一个傅立叶变换的问题# EE - 电子工程

x*b2015-01-31 08:01

1 楼

Amex offer codes for Gold Premier Card bonus Match 75K points
http://goo.gl/RwURG
目前是没有75000 points的链接的，只能先申请再match offer了。可以用上面的链接
，也可以用其他的任何链接，都是可以match的。
http://goo.gl/RwURG
对于已经申请了Amex Gold Premier Card的用户，请注意赶紧发信match 75000points
。对于没有这张卡的用户也可以考虑持有一下了，75K points千万不要换成gift card
，可以在二手版以1.5-1.55比例换成cash, 大概也值个1000刀吧。
申请拿到卡后，请在security mail里面发送信件要求match offer, 以下是offer id号
, 具体match offer的内容在下面两个论坛中有详细记录。
Amex Personal Gold Premier Card Match offer:
New Bonus ID 6939 gold premier 75K spend $1000 within 3 months
6661 gold premier 75K spend $1000 within 3 months, apply between July 1st
and October 31
5894 gold premier 75K spend $500 within 3 months
5985 gold premier 75K spend $1000 within 3 months
以上是75000points
http://www.mitbbs.com/article_t/Money/31581117.html
http://www.mitbbs.com/article_t/Money/31602041.html
http://slickdeals.net/forums/showthread.php?t=3095414
New Bonus ID: 6664 这个目前似乎不能match，还是早下手的好，目前只能拿
75000points了
The offer will give you 10,000 points for 15+ transactions you make. 50,000
point maximum. Make 15 transactions on your card in a calendar month and get
10,000 bonus points
两个加一起就是125K的bonus了，记得两个bonus要分开联系amex拿
http://www.mitbbs.com/article_t/Money/31604839.html
以下是我match offer的回信
Hi **,
Thank you for your email regarding the "Premier Rewards Gold - Earn 75,000
Membership Rewards Bonus Points" promotion.
You will be happy to know that I have linked the offer 6661 to your card
account. As per this offer, you will earn 75,000 Membership Rewards bonus
points after spending $1,000.00 within the first three months of
Cardmembership.
Therefore, I would request you to write back to me after meeting the spend
requirement and I will be happy to issue the remaining 65,000 points on your
card account.
**, thank you for choosing American Express since 200* and we look forward
to serving you for many more years to come!
Have a wonderful day.
Sincerely, **** Email Servicing Team American Express Interactive Services

a*32015-01-31 08:01

2 楼

请教一下大家，我原来在北京上学，所以去年签证是在北京签的，但我家在黑龙江，我
这次回去递签是要去北京使馆辖区的中信银行递签，还是沈阳使馆辖区的中信？

h*n2015-01-31 08:01

3 楼

狮子座与12星座的速配指数
狮子座VS.白羊座
这是一个非常火辣的配对，因为双方都拥有燃烧的热情，所以对别人来说这是热烈得让
人受不了的一对，因为两人对事情有非常雷同的看法，所以算是相当速配的一个组合了!
过招秘诀：狮子怎么说都比羊羊来得冷静一些，加上他们天生的群众魅力，羊羊可是崇
拜得很哪!当然两个人在一起时狮子发号施令的可能性是比较大的，羊羊基于对英雄的尊
敬，通常是照听不误的，加上两个人的想法和个性都很合，所以通常都能够顺顺利利地
这么发展下去。
狮子只要表现出自己的王者风度就可以迷死羊羊了，毕竟羊羊是相当天真的星座，不管
心里想什么都会自然而然地表现出来，所以羊羊的表现是很主动的，反而狮子就显得比
较冷静，而这种有点神秘的气质也是羊羊欣赏的特色，所以狮子只要保持自己尊贵的气
质就行了。
狮子座VS.金牛座
这是个很容易互相斗到死的星座配对，两人的想法差了不止十万八千里远，加上一样无
与伦比的固执，真是让每个人都捏把冷汗，只希望两人到最后还能留个全尸就能偷着笑
了!
过招秘诀：想要牛牛乖乖听话，那狮子得学习一下商人的本领，因为金牛座就是天生的
商人，而且非常懂得在商言商的规矩。如

H*o2015-01-31 08:01

4 楼

Thanks!

i*r2015-01-31 08:01

5 楼

有2列随机数字（不重复），很大比如2个G，排序之后，我想找出所有共同的数字，比
如A开头可能是：
1
2
3
7
10
74
100
349
。。。
B是：
2
3
6
78
89
90
120
。。。
他们相同的部分是
2
3
所以我需要2边同时往下读，然后不停地track谁大。实际上我的data里面，B总是比A增
加的快，所以我写了2个for loop，对每一个A，看B里面有没有相同的，找到了就记下
位置，作为下次loop的起点。
我不是科班出身，觉得自己的算法很烂（目前发现work不错，没有问题），相信有优化
算法，故请教大家
谢谢！
附上我的code结构，注意不要笑岔气了：
ind3=1;
for ind1=1:length(A)
for ind2=ind3:length(B)
if (A==B)
blahblah
end
ind3=ind2;
end
end

h*92015-01-31 08:01

6 楼

刚刚开始学java,这周就要交网上作业，可是有个地方怎么试都不对，请各位好心人帮
帮忙！
我写了一个getFNSN(FN,SN);，来模块化我的代码，可是应用起来下一行出来的总是原
始值，而不是用户的输入值，试来试去不知道怎样改才好。急死人了！
Result of adding 0.00 and 0.00 is 0.00.
部分相关代码如下：
{
public static void main(String[] args) {

float FN=0;
float SN=0;
{out.printf( "Please enter two numbers to add, separated by a space: ");

getFNSN(FN,SN);
out.printf("Result of adding %5.2f and %5.2f is %5.2f.n",FN, SN, FN+ SN);
}
static void getFNSN (float fn, float sn){
do {Scanner In = new Scanner(System.in);
try {
fn = In.nextFloat();
sn = In.nextFloat();
break;
}
catch (final InputMismatchException e) {
out.println("You have entered an invalid choice. Try again."
);
In.nextLine();
continue;
}
} while (true);}
}
结果:
What would you like to do? 1
Please enter two numbers to add, separated by a space: 2 3
Result of adding 0.00 and 0.00 is 0.00.
Press enter key to continue …

b*l2015-01-31 08:01

7 楼

假如我simulate一个w=1的正弦信号：
x=linspace(-10,10,200);
y=sin(x);
我想用Matlab看一下在频率谱图里面，w=1那正好有个1的值，其他的频率是0.这个怎么
搞啊？
stem(fftshift(fft(y)))
可以吗？
目的：我有一段音频，比如是64kHz和130kHz混合的，但这两个频率我事先不知道。我
想用Maltab里面的傅立叶变换函数画个图看一下这两个峰值在哪。请问这个怎么搞？
多谢。。。

s*f2015-01-31 08:01

8 楼

学习了

D*o2015-01-31 08:01

9 楼

作古式...香炉?

【在 H********o 的大作中提到】

: Thanks!

g*92015-01-31 08:01

10 楼

ind1 = 1
ind2 = 1
while ind1 < length(A) and ind2 < length(B)
if (A==B)
blahblah
ind1++
ind2++
else if (A < B)
ind1++
else
ind2++
end
end

【在 i***r 的大作中提到】

: 有2列随机数字（不重复），很大比如2个G，排序之后，我想找出所有共同的数字，比
: 如A开头可能是：
: 1
: 2
: 3
: 7
: 10
: 74
: 100
: 349

k*g2015-01-31 08:01

11 楼

float is a primitive type. It is called-by-value, so the caller's float
variable cannot be updated by the function even if the function modifies it.
(The function is modifying a copy of the value, not the caller's value.)
You can return the two values as a new float[2], a tuple, or a data object
that contains two float fields.
(A data object would be a reference type, which means that if the function
modify the data object's value, it will be seen by the caller.)

e*y2015-01-31 08:01

12 楼

你的音频fft后的spectrum找峰值就好了。需要自动找峰值的话写个peak finder吧

【在 b**********l 的大作中提到】

: 假如我simulate一个w=1的正弦信号：
: x=linspace(-10,10,200);
: y=sin(x);
: 我想用Matlab看一下在频率谱图里面，w=1那正好有个1的值，其他的频率是0.这个怎么
: 搞啊？
: stem(fftshift(fft(y)))
: 可以吗？
: 目的：我有一段音频，比如是64kHz和130kHz混合的，但这两个频率我事先不知道。我
: 想用Maltab里面的傅立叶变换函数画个图看一下这两个峰值在哪。请问这个怎么搞？
: 多谢。。。

H*o2015-01-31 08:01

13 楼

是香炉，不知道啥时代的。

x*z2015-01-31 08:01

14 楼

common = set(A) & set(B)

【在 i***r 的大作中提到】

: 有2列随机数字（不重复），很大比如2个G，排序之后，我想找出所有共同的数字，比
: 如A开头可能是：
: 1
: 2
: 3
: 7
: 10
: 74
: 100
: 349

h*92015-01-31 08:01

15 楼

Thank you for your answer! I got it a little bit. Could you give me an
example?
Thank you again!

it.

【在 k**********g 的大作中提到】

: float is a primitive type. It is called-by-value, so the caller's float
: variable cannot be updated by the function even if the function modifies it.
: (The function is modifying a copy of the value, not the caller's value.)
: You can return the two values as a new float[2], a tuple, or a data object
: that contains two float fields.
: (A data object would be a reference type, which means that if the function
: modify the data object's value, it will be seen by the caller.)

z*n2015-01-31 08:01

16 楼

你这个采样周期可以出现\omega=1的频谱么？
如果要得到你想要的结果，x必须是以1/2pi为周期吧。我觉得是
没法实现的。。。

【在 b**********l 的大作中提到】

H*o2015-01-31 08:01

17 楼

【在 H********o 的大作中提到】

: 是香炉，不知道啥时代的。

i*r2015-01-31 08:01

18 楼

这个好！没想到loop还可以这样写。。。。
多谢！！

【在 g******9 的大作中提到】

: ind1 = 1
: ind2 = 1
: while ind1 < length(A) and ind2 < length(B)
: if (A==B)
: blahblah
: ind1++
: ind2++
: else if (A < B)
: ind1++
: else

h*92015-01-31 08:01

19 楼

我按照你说的又试了几次，终于弄出来了~~谢谢！

it.

【在 k**********g 的大作中提到】

b*l2015-01-31 08:01

20 楼

多谢。。。我就想知道fft画出的图，那些峰值是就是对应的频率值么？还是这些频率
值的函数什么的（然后还得自己估算实际频率值）？
（比如就是y=sin(x)的fft画出来之后，是就是w=1那正好有个峰，别的地方几乎都是0
么？我试了一下fftshift(fft(y))貌似不是，不过当然很有可能我编错了。。。）

【在 e***y 的大作中提到】

: 你的音频fft后的spectrum找峰值就好了。需要自动找峰值的话写个peak finder吧

c*v2015-01-31 08:01

21 楼

仿古式，对炉子没研究，看这名可能不太老

【在 H********o 的大作中提到】

: Thanks!

i*r2015-01-31 08:01

22 楼

谢谢。不过我的数据不只有数字，我还需要把同一行的其他值打出来
另外你这个common是哪个语言的function？

【在 x*z 的大作中提到】

: common = set(A) & set(B)

b*l2015-01-31 08:01

23 楼

这个还和采样周期有关么？（我觉得只是越多越好，所以就取了个n=200）
当然我面临的具体问题信号已经给定了，不用自己simulate了，所以其实就直接关注怎
么把所含的频率找出来就行了。。
【具体问题是：手机上的1-9按键是由(697,770,853)Hz和(1209,1336,1477)Hz这两个向
量里面的频率各取一个组成的。所以每个数字对应两个频率。
现在我有一个按电话号码的音频文件，想看一下按的是哪几个号码（按的时候比如按0.
5秒，等1秒之类的间隔）。
所以我是不是就应该把音频分成10份（美国的电话号码有10位），按照上面的方法fft
一下，看看峰值在哪两个频率上，就对应是按的几？】
多谢！

D*o2015-01-31 08:01

24 楼

据说这种款有明代的, 大虾这个俺不懂
http://tag.cangdian.com/html/221/Tag_221_pm_p13.html

【在 H********o 的大作中提到】

x*z2015-01-31 08:01

25 楼

谢包子，这个是Python里面的set，我理解你只要找出两个set里面
共同的部份，然后爱怎么处理怎么处理，这个可以省掉你排序的部
份。

【在 i***r 的大作中提到】

: 谢谢。不过我的数据不只有数字，我还需要把同一行的其他值打出来
: 另外你这个common是哪个语言的function？

e*y2015-01-31 08:01

26 楼

哦，我知道你的问题所在了。fft出来的spectrum横轴是normalized frequency也就是
没有单位的。你要用采样频率把它转换为analog frequency。fft是dft而dft出来的
spectrum从-fs/2到fs/2，中间隔了n点，n是fft长度。例如你fft的第i点对应的频率是
i*fs/n.
这些你看看matlab fft文档就明白了。

0.
fft

【在 b**********l 的大作中提到】

: 这个还和采样周期有关么？（我觉得只是越多越好，所以就取了个n=200）
: 当然我面临的具体问题信号已经给定了，不用自己simulate了，所以其实就直接关注怎
: 么把所含的频率找出来就行了。。
: 【具体问题是：手机上的1-9按键是由(697,770,853)Hz和(1209,1336,1477)Hz这两个向
: 量里面的频率各取一个组成的。所以每个数字对应两个频率。
: 现在我有一个按电话号码的音频文件，想看一下按的是哪几个号码（按的时候比如按0.
: 5秒，等1秒之类的间隔）。
: 所以我是不是就应该把音频分成10份（美国的电话号码有10位），按照上面的方法fft
: 一下，看看峰值在哪两个频率上，就对应是按的几？】
: 多谢！

H*o2015-01-31 08:01

27 楼

哈哈！不过这几个字写的真漂亮，虽然不咋认识，后仿估计仿不出这种书法来。700刀
，考虑要不要拿下。

【在 c**v 的大作中提到】

: 仿古式，对炉子没研究，看这名可能不太老

w*w2015-01-31 08:01

28 楼

Memory够的话，hash一下就出来了。不够的话得分批。
#!/usr/bin/env python
smallfile = file("small.txt","r")
largefile = file("large.txt","r")
outfile=file("out.txt","w")
HASHSIZE=10*1024*1024
maxvalue=0
hash={}
while 1:
line=largefile.readline()
if not line:
break
v=int(line.strip())

if v>maxvalue:
hash={}
while 1:
r=smallfile.readline()
if not r:
break
maxvalue=int(r.strip())
hash[maxvalue]=1
if len(hash)==HASHSIZE:
break

if len(hash)==0:
break
if hash.has_key(v):
outfile.write("%s\n"%(v,))

b*l2015-01-31 08:01

29 楼

太给力了。。。貌似您说的和那帮助文档上写的就是我想要的例子。。。我仔细学习一
下。。。尤其是为啥它非要搞2的幂次当长度。。。
非常感谢。。。

【在 e***y 的大作中提到】

: 哦，我知道你的问题所在了。fft出来的spectrum横轴是normalized frequency也就是
: 没有单位的。你要用采样频率把它转换为analog frequency。fft是dft而dft出来的
: spectrum从-fs/2到fs/2，中间隔了n点，n是fft长度。例如你fft的第i点对应的频率是
: i*fs/n.
: 这些你看看matlab fft文档就明白了。
:
: 0.
: fft