a****l
发帖数: 8211 | 1 I wrote this Matlab script to examine the effect mentioned in previous posts:
x=[1:pi/100:pi*100];
y=sin(2*pi*200*x);
hold on;
fy=log(abs(fft(y)));
plot(fy,'b')
fy=log(abs(fft(y.*y)));
plot(fy,'.k')
fy=log(abs(fft(abs(y))));
plot(fy,'r')
I believe this plot explains the difference pretty well. Very interesting.
Nshaper |
|
b**********l
发帖数: 116 | 2 多谢。。。我就想知道fft画出的图，那些峰值是就是对应的频率值么？还是这些频率
值的函数什么的（然后还得自己估算实际频率值）？
（比如就是y=sin(x)的fft画出来之后，是就是w=1那正好有个峰，别的地方几乎都是0
么？我试了一下fftshift(fft(y))貌似不是，不过当然很有可能我编错了。。。） |
|
s********k
发帖数: 6180 | 3 两个flow分别做FFT的结果见附件，蓝色表示输入的FFT（对应首帖figure1），红色表
示输出的FFT对应首帖figure2） |
|
m****s
发帖数: 18160 | 4 ☆─────────────────────────────────────☆
fft (coool) 于 (Tue Jan 23 01:13:42 2007) 提到:
有这么一说吗？
如果是这样，Overseas版我推荐版二qed接任版主，此人刚从国内回来，而且手段狠一
点。
☆─────────────────────────────────────☆
vox (voice) 于 (Tue Jan 23 01:17:20 2007) 提到:
同走。再重新申请
☆─────────────────────────────────────☆
fft (coool) 于 (Tue Jan 23 01:18:13 2007) 提到:
ft，这个规矩
☆─────────────────────────────────────☆
fft (coool) 于 (Tue Jan 23 01:22:34 2007) 提到:
那我这样行不行，我把我的id转让给我信过的人，这样板斧们也不用下台。当然，如果
这个建议过于ridiculous，大家可以反对
☆── |
|
J********a
发帖数: 5208 | 5 狗屁，做出FFT的就是FFT里头最牛，做出Quick sort就是quick sort 里面最牛。图灵
会quick sort吗？不会。over |
|
N*******d
发帖数: 5641 | 6 他说的就是第一个作出fft的
按你的说法，第一个做fft不过是把计算速度提高了，根本算不得大牛 |
|
N*******d
发帖数: 5641 | 7 是不是我用fft做任何事情都不算原创？因为fft不是我第一个推出来的？ |
|
J********a
发帖数: 5208 | 8 回头，就说土共的车，380公里是不是技术创新？重点就是380公里运营。怎么了？
我看你是读书读傻了，fft推导个屁理论，人家是个算法知道不？还会推fft拉倒吧你 |
|
N*******d
发帖数: 5641 | 9 有人说过高铁是创新了么？不是原创的就不能做了？
说第一是在速度方面，设计制造当然是模仿别人的，你不要故意混淆视听。
用了fft去做了一个最牛的软件，难道因为没有发明fft没有新理论就否定他最牛？ |
|
N*******d
发帖数: 5641 | 10 别人的platform？别人怎么跑不出380？
我给你fft，你能跑出比fft更快的也算你牛 |
|
c********v
发帖数: 1278 | 11 fft还能show off？
况且你都不懂FFT |
|
N*******d
发帖数: 5641 | 12 你学了啥了？皮毛的历史掌故？FFT的背景谁都知道，就是不知道google一下就行了，
你如果还把知道这个历史当作会fft的证据，这说明你根本不懂学算法的目的。
没有humanity这句话暴露了文科生的本质，法西斯德国没有humanity，不妨碍他们有发
达的科技。 |
|
t***h
发帖数: 5601 | 13 发信人: petriv (上完QBQ论坛的目标-移民美国), 信区: Reader
标 题: 关于朱学勤涉嫌抄袭事件的公开信（修订稿）
发信站: 水木社区 (Sat Jan 22 00:40:51 2011), 站内
水木社区读书心得版部分网友就朱学勤博士涉嫌抄袭问题致复旦大学学术规范委员会的
公开信(修订稿)
修订说明：
第二稿公开信根据部分热心网友，尤其是部分参与联署的网友的意见，修改了部分措辞
，修订内容主要有三方面，首先是单纯文字错误或不当表述的修订，其次是对部分称谓
的统一，最后是删除了部分过于激烈的措辞。并没有涉及对公开信主旨思想的修订，亦
请各位已经联署或将要联署的网友审阅。
正文：
水木社区（newsmth.net）“读书心得（Reader）”版是揭露复旦大学博士朱学勤先生
涉嫌抄袭事件的发源地，在这里长期活动的大部分是从事与学术有关职业或对学术问题
感兴趣的网友，我们目睹和参与了朱学勤博士涉嫌抄袭事件从提出、争论、媒体跟踪和
复旦大学学术规范委员会着手调查的全部过程。根据网友“Isaiah”、“留美博士”和
“邢玉思”等人的文章以及大量的分析、辩论，我们认为，从公认... 阅读全帖 |
|
k*****r
发帖数: 21039 | 14 高铁，歼20，就是西南交大和西北工大最好的名片。
米国的MIT二战冷战时代牛逼过，现在绑架来的毛子国，纳粹国科学家死光了，剩下蜜
油自己的人，就只能搞伪科学了。
那个dina ka~~~啥美女教授，麻痹的为了吸引眼球发一片文章，比FFT还快的FFT，草，
整个业界都震精了。
仔细一看，根本不是这么回事儿， 麻痹的糊弄paper跟个car dealer一样的胡吹，也太
鸡巴low了。 |
|
t******l
发帖数: 10908 | 15 所以我支持正反两方用 FFT spectrogram 撕逼。。。骗得过人耳的 post-processing
大概率无法同时也骗过 FFT spectrogram / spectrum，因为素数理论的不对称性。。
。卡拉OK版无数次证明了这点。 |
|
t******l
发帖数: 10908 | 16 但是弃婴这次栽了的是，FFT 切高频的 smoothing 做不到高质量的 smoothing，因为
会发生类似 "伪影" 的现象，因为高频切掉以后，低频信号在时域或空间域发生震荡波
。。。
: 取决于具体算法
: FFT 切高频的 smoothing 无法还原，因为高频信息丢失。
: 但是 2D convolution 算子的 smoothing 可以还原，因为高频信号并未丢失，
所谓模
: 糊只是猴子眼里的 illusion。
|
|
t******l
发帖数: 10908 | 17 数字卷积器本身并不低通滤波。。。FFT 本身并不损失信息，频域大幅度低通滤波器才
大幅度损失信息。
如果卷积器加频域大幅低通滤波器，那确实理论上有可能信息量损失太大而还原不出来。
但频域大幅度低通滤波器的问题是，模糊后的艺术效果不好，比如有震荡波效果啥的。
。。而 P 图软件本质是美丽软件而不是加密软件，其 LPF 不一定是真正的频域大幅低
通滤波器。。。所以有可能出现邓尼兹潜艇医你个妈密码机的漏洞。。。当然这些都是
赌。。。
: 我觉得你说的有问题，虽然弃婴扯到泛函我也不懂，但是卷积在频率域上
就是两
者fft
: 变换的乘积，如果是低通滤波器的话那就是把高频直接截掉了，想完全恢
复过去
是不可
: 能的
: ：高端数学象牙塔把数字卷积叫泛函分析。 |
|
d***u
发帖数: 943 | 18 你把均值改成FFT不就能还原了?
要不是FFT会怎么样？ 有时能但一定会增加基础噪声，图形的基础噪声quantization是
很重要的一方面 |
|
S******8
发帖数: 24594 | 19 哈工大最新一期学报透露，两院院士、中国新体制雷达
奠基人刘永坦领导研制的引力波无源雷达进入测试阶段。
极大地振奋了人心、士气。
中国的新体制雷达在刘院士的带领下，从80年代初就
开始进行了积极的研究，可谓30年磨一剑。
传统雷达依靠发射雷达波，接受被探测物体反馈的信
号实现探测。
但在隐形技术迅猛发展的今天，采用吸波纳米材料、特殊外形设计、波形相干干涉
等技术实现了对飞机、舰艇等武器的隐形，让传统的旧体制雷达陷入英雄无用武之地的
尴尬局面。
严峻的现实，促使新体制雷达的研制进程大大加快。
中国的新体制雷达经过短短20余年，呕心沥血，从无到有，从理论、雷达数据处理软件
、样机小型化方面均取得成熟的科学成果，目前仅有美国在这方面有类似的研究，但距
离实用应该还有不少的距离。
外媒猜测，哈工大雷达的这次泄密是让美军F22面临停产的真正原因。
引力波无源雷达是一种全新概念的雷达。
我们知道，任何有质量的物体都会向外发出引力波，但对引力波的探测是当今世界的一
大科学难题。
只要实现了对引力波的有效探测，世界上就不会存在可以隐形的东西。目前美国最先进
的F22只能实现对传统雷达的隐形，但在“... 阅读全帖 |
|
p*****u
发帖数: 225 | 20 新版本的matlab（好像是2010以后）fft是自动多核的，如果你的瓶颈是fft用C写不会
快多少，可以用profiler看不同函数运行了多长时间和次数。用GPU能块，不过你所有
的东西都要在GPU上实现，不然在GPU内存和系统内存间互相倒数据也会很费时间，特别
是数据很大的时候。 |
|
n*******l
发帖数: 2911 | 21 很有道理。多谢。
我刚刚用profile看了看，2/3的时间是fft和ifft。估计那些部分是没办法了，看看
其他部分能不能省点儿是点儿。
我的程序总共就一个400来行，一层循环，就是每个时间步长计算若干次FFT 和 IFFT，
让我看看能不能把IFFT省掉。 |
|
c****u
发帖数: 584 | 22 换个算法。
我解过积分方程，时域的算法是local的，可以用非等分grid. Matlab 1分钟就能算到1
毫秒，精度到皮妙量级。蘋域FFT的算法是global 的，要用等分的grid. 我们用了640
个CPU, 两天才能算出同样结果。而且FFT 对round off errors非常敏感，还要用
infinity precision package 来算。 |
|
n*******l
发帖数: 2911 | 23 我就是算着玩儿。专业做计算的肯定是cluster，个人用我建议
Mac Pro。我新买的Mac Pro配置如下：
CPU: 2.7GHz 12-core Xeon E5,
Memory: 64GB 1867 MHz DDR3 ECC，
GPU: AMD FirePro D700 6144 MB,
Storage: 1TB flash storage.
我还有去年买的一台iMac:
CPU: 3.5GHz 4-core i7，
Memory: 32GB 1600MHz DDR3,
GPU: Nvidia GTX 775M,
Storage: 1.12TB fusion drive.
我就是用谱方法解方程，主要是Matlab里调用FFT，因为Matlab的矩
阵计算和FFT都是并行的，这两台机器性能差别非常明显。同样的程序，
同样的参数，iMac需要20个小时，Mac Pro只要9个小时。
另外，就我的应用来说，上述iMac的性能和15寸MacBook Pro没什么
区别。 |
|
r*****e
发帖数: 620 | 24 离散傅立叶变换的复杂度是就N^2吧, A是对的
FFT是NlogN
估计你没分清楚DFT和FFT, 错上加错给了个正确答案，没有丢脸:) |
|
d***a
发帖数: 13752 | 25 哈哈，我以前做过一个项目，在胖老师说的对面的位置，不过对方不是搞生物的。我们
最头大的是，已经和对方开会讨论N次，终于把需要说明书搞好，对方签字了，等我们
做出来，对方说这不是他们想要的...
胖老师说的对，我也感觉楼主这些年，从职业角度来说，精力花的不是地方。如果是计
算机科班出身的，第一不会花这么长时间，第二不会把编程看得这么重。编程是工具和
方法。前面有位版友提到FFT。从专业角度来说，FFT的意义是算法上的突破，和编程并
没有什么关系。楼主如果想在computational science上发展，算法上的突破最重要，
或者建立新的计算模型，或者考虑把计算问题map到新的计算机硬件上（GPU, FPGA，特
殊构造的CPU等）。单纯追求编程上的提高，很fun，但不能体现自己的价值。
这样说可能有些不nice，但从职业的长远来说，楼主也需要有人客观冷静地指出问题。
楼主这些年的收获还是很大的，才七年时间，老婆孩子都有了，也做了自己想尝试的事
情。接下来好好规划下职业发展，前景应该会不错。
至于你能做什么？我建议你去做学术界和编程界的接口工作。我很久以前的
时候也能编程，也对... 阅读全帖 |
|
g**e
发帖数: 6127 | 26 面的资深码工职位，遇到某bt俄罗斯老头，于是有了以下题目：
1. what is a Ring in algebraic structures
2. tell me about your understanding about Riemannian geometry
3. write down FFT definition. name a few implementations of FFT algorithms
4. how do you implement taylor expansion
5. what are eigenvalues and eigenvectors, how do you calculate them
6. how do you solve a system linear equations
其它常见IT题目就不说了，说几个不常见的
1. how do you extract information from an udp packet
2. given the following DB table:
symbol, timestamp, bid, ask
A... 阅读全帖 |
|
j****x
发帖数: 149 | 27 fft到频率域，算完后再inverse fft到时间域 |
|
f******n
发帖数: 198 | 28 FFT是比n^2慢，我是对上面那个说CG处理的说的。Box filter的特殊性质决定了你可以
用较小的运算量从相邻的点算出当前的点的数值，所以可以做到n^2。如果是一个不规
则的filter，就要FFT了。 |
|
s****r
发帖数: 278 | 29 来自主题: JobHunting版 - CS真难学 FFT易于实现而已，无它。面试谁都不会去问蝶形图吧，who cares 怎么实现啊。信号
处理关键是把傅立叶级数到傅立叶变换，周期离散到非周期离散，离散到连续，采样定
理，上下采样，前因后果弄明白了，才是信号系统的精髓。至于FFT怎么实现，总有苦
逼去干的。。。I dont care, lol
本科DSP的教材和老师专注于各种变态蝶形图，实在是丢了西瓜捡芝麻。 |
|
t*******r
发帖数: 22634 | 30 俺上次在音乐版就说过，合唱的和声通常不会比唱单音更容易准，只是和声不容易
听出偏差（当然，大三小六要往纯律方向偏，偏的方向没搞错）。。。一堆人跟俺
争。。。
俺被搞急了，拿了给的一个合唱录音里面的一长音做了 FFT spectrogram 分析，
指出有多少个声部，每个声部的音，跟平律差多少，跟纯律差多少。每个声部里面
自己人跟自己人有偏差的话，差多少。。。哼哼。。。俺耳朵不好使，但眼睛看
FFT spectrogram/spectrum 难道还看不出来么。。。 |
|
t*******r
发帖数: 22634 | 31 确定 note 弹对，实在不行你看小孩的键盘，要是再不行。录下来看
FFT 频谱。。。不过这个 FFT 老实说看起来也太累了。。。 |
|
t*******r
发帖数: 22634 | 32 三截棍扒谱，就是用个 keyboard 或 piano 边听磁带边试。
不过这个已经不先进了，俺目前典藏 "捷克式轻机枪扒谱"，就是用免费软件 Audacity
把 mp3 投射成 FFT spectrogram，然后选择战略高地 plot 该地区的 FFT spectrum
，最后用尺子量量唱歌小妞唱偏了多少 cent，唱了些啥效果，伴奏乐器坐在哪里，电
子效果器大概是咋回事，等等。。。当然，核电还是水电的问题，目前捷克式轻机枪无
法分辨。。。 |
|
t*******r
发帖数: 22634 | 33 FFT spectrogram 上，可以看到 singer 的 pitch bend 的细节，用尺子量泛音分布的
变化，也可以推测音色变化的方向和技术。
不过主要是太麻烦，而且开车时没法同时干活，所以我虽然有典藏的捷克式轻机枪，其
实平时并不用。。。当然曾有人跟俺争论合唱准不准的问题，俺一急，端出捷克式轻机
枪，FFT spectrogram / spectrum 上清清楚楚看看他们的某一声部唱成三个偏差不同
的频率峰值，幅值的相对差别一清二楚，一个扫射就把山头削成山谷了。。。 |
|
a****l
发帖数: 8211 | 34 FFT是不适合测量人的说话的声音的,这也就是你说的量的准不容易的原因,因为你根本
就是在用错误的工具测量.
还有前面的谁说的FT,更是没边了,大概以为FT就是更慢更准确的FFT? |
|
t*******r
发帖数: 22634 | 35 唱歌界文盲很多，把旋律记住就行了。
俺走音公鸭嗓自己能看一点儿谱子，如果是器乐需要谱子。。。但唱歌俺从来懒得看谱
或 FFT，直接记忆旋律更方便。
当然老帕多半不会真看不懂谱子，我觉得是多半是嫌麻烦不看，也不去做啥狗屁听音写
音名。
而且很多听音写音名的测试，还不一定在调性环境里。。。就算是 perfect pitch，也
都是在调性环境里更准。。。那种测试很多纯粹就是回字九种写法而已。。。当然要扒
谱另说，不过活人扒不过 FFT 的，除非脑外伤导致 island intelligence 那种。。。
而且有科学家怀疑，没有真正意义的 perfect pitch，其实都是 tonal memory 的一种
。（童子功的 tonal memory 更坚固，其实也算脑畸的一种）。。。从这个角度看，不
识谱的老帕文盲倒可能是正确方向，因为伊可能把小时数都投资在最最实用的 tonal
memory 上了。。。而那些执着于回字九种写法的，可能最后大半是孔乙己了。。。 |
|
t*******m
发帖数: 999 | 36 我也学过FFT，还写过FFT的程序呢。真没想到还能这么用！有潮水大人，咱们这坛子里
学音乐的小孩有福了。 |
|
J***A
发帖数: 1511 | 37 呵呵 也是：） 以后娃到了学琴的年龄会录下来给她自己听的， 再搞个fft 软件让她
自己看结果。
刚才还有点小后悔当初没舍得多加钱买那个能自己弹的......
录音分析的话，免费的 FFT spectrogram 足够了，其实不需要机器人钢琴那种高级货
。 |
|
|
t*******r
发帖数: 22634 | 39 这个分析器保不准做得不比专业的差，不过就是界面不好用。
说明书跟没有啥不多，基本先要去修一门数字信号处理外加
一门图像处理才能玩得比较溜。
俺这个线程是不是可以认为写了 Audacity
Spectrum/Spectrogram/Frequency Analysis for Dummies。。。
另外界面也比较差，坑爹的地方有：
（1）鼠标每次移动就动 1 Hz，其实 1 Hz 不小了，搞不好
有 8 个 cent 了。
（2）自动检测峰值的频率也就精确到 1 Hz。
（3）鼠标移动时的峰值也就精确到 1 dB。
不过好歹 Audacity 提供界面把结果送到 .txt，估计是可以
让大家把结果往到 Matlab 里插值分析。。。不过这个更坑
爹，都已经上 Matlab 了，干脆直接把原始音频送进 Matlab，
看着 Audacity 上的时间轴在 Matlab 里算 FFT，算完把结果
做一阶求导，解出所有零点，选出其中所有二阶导数为负数的
解，放入结果矩阵。然后看着结果矩阵重唱，保证唱准唱好！！！
有木有！！！有木有！！！
这个 Audacity，显然是老板强迫员工们把... 阅读全帖 |
|
t*******r
发帖数: 22634 | 40 虽然 FFT 很粗略，但是现阶段 FFT spectrogram 已经够分析音色了，
美声牛人和 classic-crossover 牛人还是牛鼻。。。 |
|
t*******r
发帖数: 22634 | 41 我已经看了很多 spectrogram 了，这个是挺有用的。的确小楼说的
正确，看 spectrogram 根本看不出声音好听与否。但 spectrogram
的价值并不在于替代耳朵，而是帮助分辨耳朵不容易听清楚的东西。
举个例子，音色亮，可能是声带支持得好，也可能是 belting，这个听得出来。
但是如果在练习中，到底更趋向于声带支持还是 belting，分辨是有难度的。
但是 FFT spectrogram 上面看这个就相对容易。
实际上这个 FFT spectrogram 已经对俺有帮助啦，俺还没有录音呢。俺看了
以前录的 spectrogram，知道有些音质相对较差，搞不清为啥。一看
spectrogram，一目了然，俺现在公鸭吼的时候又可以调整一下。嗓子更
轻松喽。。。 |
|
t*******r
发帖数: 22634 | 42 我看了看 FFT spectrogram，从 FFT 的角度看，楼主确实比以前进步不少。 |
|
t*******r
发帖数: 22634 | 43 俺实战了一下，提笑肌的确能明显提高声音的位置，使得
声音更有灵气。
提笑肌和 FEA 声带支持不矛盾。俺觉得提笑肌的作用是
改变共鸣腔。FEA 声带支持是基础。
俺做了提笑肌 vs 不提笑肌的对比测试以及 FFT spectrum
分析，提不提笑肌，谐波的分布特征大致相似。但是提笑肌
以后，总体的变化是高的地方更高，低的地方更低，同时能
量更多得移往谐波，使得声音更明亮通透。简单一句话，
提笑肌以后，FFT spectrum 上看到 B cup 变成 C cup，并
且更挺！！！
提笑肌的难点在于控制肌肉联动，不要因为提了笑肌引起不必
要的喉头紧张或上提。。。或者影响咬字。。。 |
|
t*******r
发帖数: 22634 | 44 这个是俺担心的地方。每个人的嗓子有正常的 thin 或者 thick 的特性。但如果声音
thinner than it supposed to be，可能说明声带支持和气息支持没有提供足够
的能量/能量分布。。。我的担心是基础没打好就开始花哨。。。当然我基础比较差。
。。
其实我后来一天又录过一次，声音更 open，FFT 看着也是高频更多，但声音更 thin，
有点 brittle 了。The result is consistent。。。
我对比了 FFT，声音更 thin 貌似是在每个 harmonics 上，能量更集中于峰值点，
峰值附近下降得更快。这种情况提示声带可能更加趋向线性振动而不是非线性振动（极端
一点说，更像一根弦在振动而不是一块板在振动），提示声带支持相对不足。
但是我自己并没有感觉声带支持减弱，我的猜想是更大的共鸣腔 demand 更强的声带
和气息支持，共鸣腔可能通过动态声阻反馈作用在声带和气息。
流行另说，流行或多或少总是靠些 belting 来提供能量，另外也可能靠效果器加
厚声音。 |
|
t*******r
发帖数: 22634 | 45 这个我听了并且看了 FFT，说实话声音很好，我是跟美声女高的 FFT 频谱
对比的，基本泛音列的特征啥的都有，声音也听起来很漂亮很 open。。。
老大就是恶搞的太多，不恶搞的的确非常好。。。 |
|
t*******r
发帖数: 22634 | 46 跟我自己的感觉一致，在 FFT 频谱上，出现 2kHz - 3kHz 之间的谐波列
的峰值比较平。谐波列的峰值比较平说明有些频段波形接近一点方波形状。接
近方波形状可能有 belting (yelling) 的倾向。yelling 让人感觉声音在
口腔/嗓子里。如果这个角度说得通的话，FFT 跟听感应该是一致的。 |
|
t*******r
发帖数: 22634 | 47 有谱子岂不更简单。。。
比如歌的音域也可以听，俺以前一直是听，自从会了 FFT
以后，每次都在 FFT 上看一眼。。。懒惰无止境啊。。。 |
|
t*******r
发帖数: 22634 | 48 。。。连这个都可以躺着中枪啊。。。
俺的确发现有 FFT 量着准和听着准的差别。classic crossover
可能是绝对音高上最准的（opera 颤音那么大。。。说实话我也
只好将就着量。Opera 那种颤音，我听着准，就往准的方向量，
基本上也是屁股决定脑袋。。。）。
像流行这种，在误差十几 cent 的时候，可能偏的方向比偏的
绝对值更重要。比如大三/大六往上偏，小三/小六往下偏，纯四/
纯五站中间。。。减七？。。。减七我当没看见。。。这种和弦/
爬音的色彩以及对比关系可能对主观感觉更重要。。。
不过说实话唱的时候算 FFT 估计来不及。演奏时的音准最终可能
还是靠 记忆 + 和弦/爬音色彩感觉。。。 |
|
t*******r
发帖数: 22634 | 49 老实说，我本来录完以后打算用 FFT 去量一下。录完
听了一遍，发现不用量了。。。走音太多了。。。
用 FFT 量也是很累人有木有的！！！
其实。。俺发现。。。不听自己的录音不就完事了么。。。
俺的破录音本来就是复查用的。。。 |
|
