发帖数: 1 | 1 我觉得你说的有问题,虽然弃婴扯到泛函我也不懂,但是卷积在频率域上就是两者fft
变换的乘积,如果是低通滤波器的话那就是把高频直接截掉了,想完全恢复过去是不可
能的
:高端数学象牙塔把数字卷积叫泛函分析。
: |
|
发帖数: 1 | 2 数值信号还是可以当做原信号的数字化采样嘛,原信号f*g卷积做完fft就是 F·G,理论
上就是g这个卷积核在频率域上滤波,能不能恢复基本要看g的性质了。
不知道理解的对不对
:数字卷积器本身并不低通滤波。。。如果卷积器加频域大幅低通滤波器,那确实理论
上信息量损失太大而还原不出来。
: |
|
t******l 发帖数: 10908 | 3 理想低通滤波器是不可能实现,但马工坑爹非理想数字频率大剪刀是可以实现的,也就
是先 FFT,然后写个 if (freq > 1Hz) return 0。。。虽然老板会很恼怒,但编译器
还是过了不是?
: 你不可能截去高频的部分,因为理想低通滤波器是无法实现的。
: 你可以压制高频的部分。但是只要是线性处理的话,总可以恢复出来,畸变可以
保持很
: 小。
|
|
t******l 发帖数: 10908 | 4 数字频率截断是很常见,只不过不是虎肉老师说的理想低通滤波器(这个行话里面指模
拟滤波器)。
这个最关键的原因,是在于数字图像是有限的离散信号。所以我们不严格的说 FFT 的
时候,实际上是在说 DCT discrete cosine transform。。。其结果在数学上不是无限
区间(实数区间)上的 Fourier Transform 连续函数,而是有限长度 Fourier Series
有限长度离散序列。。。而这样使得数字频率截断可以实现,而理想模拟低通滤波器
不可实现。
这也就是我前面构造那道 AMC12 题目的用意之一,有限离散序列而不是无限连续函数。
弃婴这次在完全不是其领域的地方发力,丢大发了。。。 |
|
发帖数: 1 | 5 你连什么是采样都不明白就不要在这现什么数字图像了。
还奈奎斯特,离散,FFT,你先把冲击函数搞明白是正理。 |
|
|
|
w***u 发帖数: 17713 | 8 人类的知识体系,不都是搭积木嘛,能把别人的building block用好就是技术进步,我
们搞冶金搞材料的还觉得搭飞船搭摩天大楼的没啥呢,可我们被放到鄙视链的最低端了
。用C能把fft,eigenvalue写得又快又好的有几个?用好matlab做下一步工作,不该被
鄙视吧。 |
|
e***d 发帖数: 8248 | 9 比如说,相关运算,在信号似稳时,相关最大等效于差模最小,
这样模数转换进来的信号只有整数加运算和符号位置位运算,
所有的乘法都变成加法,数百倍的效率提升。
再比如说,FFT的Divide and Conquer.
码农只能做些特定条件下的改进,比如例子一。
革命性的东西还要靠数学家,比如例子二。 |
|
w***u 发帖数: 17713 | 10 计算机的普及岂只是消灭了低端人口的工作机会?大量中高端人口的机会也减少很多。
看看工学院学生就知道了。原来会个微积分、机原机零、土力学、传热传质、电路分析
傅立叶变换,就算可以凭脑袋和手艺吃饭的中高端人口了。现在计算机暴力破解一切,
无论是工程设计图还是有限元,无论FFT还是VLSI设计,从原子弹爆炸仿真到大尺度天
气预报模型,普通计算机几秒钟超越成千上万工程师的工作量,只剩下一点现场监控工
程师和给程序提供计算框架的理论分析人员。 |
|
f***y 发帖数: 4447 | 11 4月23日,中国电科38所在福州举行的首届数字中国建设峰会上发布了实际运算性能业
界同类产品最强的数字信号处理器——“魂芯二号A”。该芯片由38所完全自主设计,
在一秒钟内能完成千亿次浮点操作运算,单核性能超过当前国际市场上同类芯片性能4
倍。
高性能芯片被誉为“工业粮草”,代表了一个国家信息技术水平。一直以来,我国在高
性能数字信号处理器(DSP)方面始终依赖进口。 12年前,38所就开始进入数字信号处理
器芯片领域。2012年,该所推出我国自主研发的首款实用型高性能浮点通用DSP芯片—
—“魂芯一号”,性能高于同期市场同类DSP芯片4~6倍,并成功应用在我国空警-500预
警机雷达等多个国防科技装备上,成为我国首款广泛应用于国防科技装备的高端自主数
字信号处理器。
“魂芯二号A”采用全自主体系架构,研发历时6年,突破了控制器设计等多个技术难题
,获得国家技术发明专利、软件著作权等科技成果30余项;拥有当前业界性能最强的
DSP核,实现了对国内外同类产品性能指标的超越。相对于“魂芯一号”,“魂芯二号A
”性能提升了6倍,通过单核变多核、扩展运算部件、升级指令系统等手段,使器件性
能千... 阅读全帖 |
|
s*x 发帖数: 8041 | 12 据安徽日报报道,4月23日,中国电科38所在福州举行的首届数字中国建设峰会上发布
了实际运算性能业界同类产品最强的数字信号处理器——“魂芯二号A”。该芯片由38
所完全自主设计,在一秒钟内能完成千亿次浮点操作运算,单核性能超过当前国际市场
上同类芯片性能4倍。
资料图
高性能芯片被誉为 “工业粮草”,代表了一个国家信息技术水平。一直以来,我国在
高性能数字信号处理器(DSP)方面始终依赖进口。 12年前,38所就开始进入数字信号
处理器芯片领域。2012年,该所推出我国自主研发的首款实用型高性能浮点通用DSP芯
片——“魂芯一号”,性能高于同期市场同类DSP芯片4~6倍,并成功应用在我国空警-
500预警机雷达等多个国防科技装备上,成为我国首款广泛应用于国防科技装备的高端
自主数字信号处理器。
“魂芯二号A”采用全自主体系架构,研发历时6年,突破了控制器设计等多个技术难题
,获得国家技术发明专利、软件著作权等科技成果30余项;拥有当前业界性能最强的
DSP核,实现了对国内外同类产品性能指标的超越。相对于“魂芯一号”,“魂芯二号A
”性能提升了6倍,通过单核变多核、扩展运算部件、升级指令... 阅读全帖 |
|
s*x 发帖数: 8041 | 13 吗?
中兴将重获生机?我国5G芯片获重大突破,量产在即
陷入困境的中兴将要迎来新的生机?中科院自动化研究所耕耘9年,具有完全自主
知识产权的5G芯片获得重大突破,预计明年即将量产。这一成果犹如一场及时雨,很可
能会使中国在即将到来的5G时代实现继续领跑世界,令欲打算借助“禁令”阻断中国5G
领跑势头的美国计划落空。具体怎么样?我们一起来看看。
首先我们应当知道,本次美国对中兴的封杀绝非是一个中国企业和一个最发达国
家之间的恩怨,而是在即将到来的5G时代,中国和美国谁更有话语权、谁能够成为5G行
业标准的制定者之间的较量。所以,美国人的封杀令实质上是要封住中国率先跨入5G大
门的势头,从而为美国超越中国继续成为新的行业标准制定者争取时间。总之,在嗜血
成性的美国人心目中,中国人只能永远是他们榨取血汗的奴隶,而不能是挺直腰杆成为
西方世界的对等者!
那么,在美国人阴毒的迟滞中国的计划里,中国真的就没希望了吗?不过我中华
威武、天佑中华!日前,中科院特聘研究员、自动化研究所原所长王东琳对媒体宣布,
具有完全自主知识产权的5G芯片获得重大突破,已经具备量产条件。... 阅读全帖 |
|
f***y 发帖数: 4447 | 14 http://news.sciencenet.cn/htmlnews/2018/4/410090.shtm
曾经如日中天的中兴通讯或许不会想到,一枚小小的芯片按下了企业快速发展的“暂停
键”:美国商务部日前宣布,7年内禁止美国企业向中兴通讯出售任何电子技术或通讯
元件。
突如其来的中兴危机折射出我国通信产业“缺芯少魂”的尴尬处境。未来“中兴们”何
去何从?“中国芯”弯道超车的集结号已经吹响。中国科学院自动化研究所王东琳团队
日前正式对外宣布,其历时九年研制出了具有自主知识产权的5G通信设备强大“心脏”
——UCP(通用通信芯片)内核,在计算能力和性能功耗比方面具有国际领先水平。
基站、光通信、手机,这三项是中兴的主营业务。美国的一纸禁令究竟意味着什么?透
过一组数据或可管窥一斑。根据有关报道,仅中兴智能手机产品中就有一半以上使用的
是高通芯片。
而这只是冰山一角。“国产芯片目前主要应用于消费领域,在通信、军事等行业应用中
,与国外产品差距明显,诸多短板造就了我国高端芯片受制于人的局面。”中科院西安
光机所副研究员、中科创星创始人米磊说。
米磊的担忧并非多余。《2017年中国集成电路产业... 阅读全帖 |
|
发帖数: 1 | 15 Fft
1235680
: 网上都不会告诉你三角函数怎么算出来的
: 当年我国内老板主持祖国第一个函数计算器的算法,这么小的rom,这么弱的算
力,需
: 要在一眨眼功夫算出三角函数,用的是什么算法?
: 后来这台函数计算器作为礼物送给第一次回共产中国的羊82,羊82知道算法后大
为惊异
|
|
z****g 发帖数: 3509 | 16 中电科报告发明世界最厉害芯片魂芯二号A 但没说能否救中兴
作者 小山 发表时间 30-04-2018 更改时间 30-04-2018 发表时间 13:34
多家媒体透露中国研制发明成功一项芯片技术,该芯片称为“魂芯二号A”,据信是世
界实际运算性能业界同类产品最强的数字信号处理芯片。这一成果在中国电科38所近日
在福州举行的首届数字中国建设峰会上发布,但没有引起中国官方媒体高调轰动性报道。
据美国媒体表明引自多维消息报道,中国电科38所近日在福州举行的首届数字中国建设
峰会上发布了实际运算性能业界同类产品最强的数字信号处理器 “魂芯二号A”。该
芯片完全由电科38所自主设计,一秒内能完成千亿次浮点操作运算,单核性能超过当前
国际市场上同类芯片性能4倍。
报道说,中国电科首席科学家、“魂芯二号A”总设计师洪一介绍,中国人可以做出自
己的芯片,魂芯一号、魂芯二号A的陆续推出,充分验证了国产处理器的可用性和可靠
性。
据介绍,美国德州仪器公司生产的TMS320C6678芯片,是目前国际市场上通用的DSP产品
。与国际市场同类产品相比,“魂芯二号A”同样表现出色,并且实现性能的全面超越。
报... 阅读全帖 |
|
p***n 发帖数: 17190 | 17 用聽不到的低頻
靠在牆上震盪
低頻聽不到但可以從人體傳到腦部去
不用幾天就會有影響
可以用個小壓電器
接FFT 就能抓到
以前我做STM時就這麼抓雜訊 |
|
发帖数: 1 | 18 太湖之光的CPU sw26010在设计取向上是针对“暴力浮点”优化的,出现上面这个结果
(FLOPS强,某些其他测试不太好)意料之中。
比较显著的一个指标是这个CPU的内存带宽/浮点比率达到令人难以置信的1/22.4=0.
045 Byte/FLOP。作为对比(数据来源Yokota,J Algo Comp Tech, 7,3,2013):
INTEL Xeon E5 2690: 0.211
AMD Opteron 6284 SE: 0.235
NVIDIA Fermi GF110: 0.266
IBM PowerPC A2(BG/Q):0.208
Fujitsu Sparc64(FX10):0.359
相比较而言sw26010的内存带宽太小,很多操作都会卡在内存上。实际应用上,优化好
的FMM或者DGEMM大概需要0.2,使用特殊向量指令还可以更低,3D FFT大概需要0.6~0.
9, stencil大概2~3,spmv大概5。在x86上FMM/DGEMM是compute bound,后面几个都
是memory bound,而在sw26010上这些操作几乎全都是memory boun... 阅读全帖 |
|
发帖数: 1 | 19 全球一盤棋:漫談2018年的中美博弈(下)
2018/06/14 16:00:13
中國的半導體產業
中美貿易戰最引發人們的關注點是中興通訊的休克事件,所以半導體技術是一個重要話
題。
中國的半導體產業的確落後美國,但是沒有一般人想像的落後這麼多。美國對中興通訊
禁售高通的晶片和谷歌的安卓操作系統,中興立刻倒地宣布休克。這是非常戲劇性和吸
引眼球的事件,自然成為最佳話題和群體討論的焦點。吹捧美國的人立刻拿這個大作文
章,說什麼中國科技如何差勁、如何沒有根基、落後美國多少年...等。其實中興休
克倒地是事實,但是它再度爬起來並不難,因為中國的科技是有根基的。當中興爬起來
的時候,前景會更好。
首先,美國禁運晶片,導致中興休克,是一個犯規的出擊,並不是美國有多厲害,用拳
擊的術語來說,這是一個“cheap low blow”。拳擊賽中,腰部以下是不能打的。所謂
"low blow"就是打肚子,在拳擊賽中是犯規的。“cheap low blow” 直接白話翻譯就
是「犯規的賤招」。
在WTO的組織架構下,產業鏈各有各的分工,上游產業的國家禁售產品給下游產業的國
家是不被允許的,是非常廉價和惡... 阅读全帖 |
|
p***n 发帖数: 17190 | 20 能不能寫個程式
每天自動掃描 人民日報頭版的標題有幾個包帝
然後當做y 值 對日期x 畫成曲線圖
再取FFT |
|
p***n 发帖数: 17190 | 21 你要放長了看
過幾天打臉
過幾周就是持續上揚
也就是說
取個FFT會有幾天的peak
但取個幾周的線性回歸
還是持續上升的
你看看老將前幾周還在炒 6.6
現在GDP跟人民幣匯率都已經確定華麗嬌插了 |
|
m***y 发帖数: 14763 | 22 就算从用FFT开始,OFDM也是七十年代的技术了。华为中研某高级干部就敢冒充人爹,
卡桑你可千万别回去,人家那“一将成名万骨枯”,你就是那个万分之一。 |
|
p***n 发帖数: 17190 | 23 那華北跟華中
以後還要種麥嗎
還是要改種米?
還是種甘藷確保大家能吃上飯?
我覺得要看土地的不同
氣候暖化也是一年一年過去
所以要對各地的土壤跟氣候
濕度與溫度
降雨與積水
做長時期大數據的分析
加上AI與FFT等工具
配合最新育種的成效
來決定要不要改麥為米
快寫成梁家河大學問了
你們給補上吧 |
|
s***h 发帖数: 487 | 24 内置音乐遥控情趣跳蛋不是简单的节拍器,是 FFT spectrogram 频域时域分析仪。
: 都是专业的,再用节拍器让人笑话
: :我就想知道三个女生有没有被内置音乐遥控情趣跳蛋,比如 Lovense Lush 之
类的 。
: :。。合着音乐的节拍跳舞震动 。。。
|
|
|
s*****V 发帖数: 21731 | 26 稍微强度高的科学计算要是用python做慢死,比如FFT, linear algebra, geometry。
估计能有1000倍的差距,1000X对很多应用就是存在和不存在的差距了。
C |
|
s***h 发帖数: 487 | 27 科学计算大不了还可以等一年,Robotics 才更坑爹。
比如花了一秒才算完,但就这一秒钟机器人已经跌在阴沟里了。这种算还不如不算,至
少省点电。
: 稍微强度高的科学计算要是用python做慢死,比如FFT, linear algebra,
geometry。
: 估计能有1000倍的差距,1000X对很多应用就是存在和不存在的差距了。
: C
|
|
|
p***n 发帖数: 17190 | 29 這個以前在做STM
用FFT 看電信號
發現整個建築物在抖動
把光學桌浮起來
所謂量子力學
其實是
整個宇宙 都在抖動
你以為的局域測量 是廣域抖動的投影 |
|
p***n 发帖数: 17190 | 30 共產黨可以盡量把各個景氣循環的相位往後推遲
最後結果就是相位重合
就好比 step function 的 FFT 之後的
重合 |
|
h*********n 发帖数: 11319 | 31 laf小kernel,不具有通用性
人类大部分科学计算程序都有规整的数学表达式,规整的数据结构和规整的“热点”循
环。特别是底层的数学库,往往非常适合用gpu加速
看看blas, fft, acml oski这些典型数学api的加速比再来说话吧。别告诉我他们的cpu
版本都是“不怎么样的baseline”。在上cuda showcase瞧瞧有多少商用软件开源软件有
了gpu版本。 |
|
m*********e 发帖数: 533 | 32 全球一盤棋:漫談2018年的中美博弈(下)
中國的半導體產業
中美貿易戰最引發人們的關注點是中興通訊的休克事件,所以半導體技術是一個重要話
題。
中國的半導體產業的確落後美國,但是沒有一般人想像的落後這麼多。美國對中興通訊
禁售高通的晶片和谷歌的安卓操作系統,中興立刻倒地宣布休克。這是非常戲劇性和吸
引眼球的事件,自然成為最佳話題和群體討論的焦點。吹捧美國的人立刻拿這個大作文
章,說什麼中國科技如何差勁、如何沒有根基、落後美國多少年...等。其實中興休
克倒地是事實,但是它再度爬起來並不難,因為中國的科技是有根基的。當中興爬起來
的時候,前景會更好。
首先,美國禁呔瑢е轮信d休克,是一個犯規的出擊,並不是美國有多厲害,用拳
擊的術語來說,這是一個“cheap low blow”。拳擊賽中,腰部以下是不能打的。所謂
"low blow"就是打肚子,在拳擊賽中是犯規的。“cheap low blow” 直接白話翻譯就
是「犯規的賤招」。
在WTO的組織架構下,產業鏈各有各的分工,上游產業的國家禁售產品給下游產業的國
家是不被允許的,是非常廉價和惡劣的傷害行為,表現的是愚蠢的衝動,並不能顯示... 阅读全帖 |
|
|
d*****w 发帖数: 684 | 34 另外就是从来没有考虑过做时序序列的那种分析,像你说的自相关那种。一次我无聊登
陆美国的一个网站,发现有人用FFT变换来做分析呢,觉得有点走火入魔了。
实际上,TA是易学易用难以盈利,所以最终大部分高手们都转到了难学难用但是容易盈
利的FA去了。当然了,坚持TA到底的高手人物也有,就是太稀少了一点。 |
|
s****i 发帖数: 1593 | 35 FFT我也想过,哈哈!主要最近就在考虑谱分析的东西,所以联想到股票上了。
其实好多多东西和其他领域都是通的,关于均线,在我这块算是条件平均的一种。
还有你说的胳膊拗不过大腿,在我们这块就是尺度分析,大波大尺度,小波小尺度。
哈哈,现在工作不好找,为转行打基础了。以后就拿操盘记录去应聘算了。
国内好像有个牛人就是拿高三一年的操盘记录,拿到牛津还是剑桥的offer的吗? |
|
G**Y 发帖数: 33224 | 36 今天读paper,讲一算法,读着读着,人家作者说我们这个是专利的。这算法看上去还
不错,作者似乎有点资历的,不像乱来的。
可这是啥意思呀,别人再也不能用这个idea了?
我以前只知道FFT是专利的。
这要是大家想出点新主意就专利一下,还怎么玩呀。 |
|
n*******l 发帖数: 2911 | 37 似乎看不到哪个版面更合适,就问问众位经验丰富的牛牛们吧。
现在我在做一个耗时间非常长的计算,用谱方法解微分方程,主要就是
重复的FFT。现在我用Matlab,在我的iMac上算一次要用好几个礼拜。
这要是增加一下规模的话,算一次就是好几个月了,受不了。而且学校
太烂,也没有像样的有足够内存的并行机可用。
据说c++要比Matlab快一点儿,但是我想问的是,按说Matlab的向量运
算和并行处理已经优化得很好了,反正是要调用一些已经优化好了的库函
数,C++真的会比Matlab强一点儿么?强在哪儿呢? |
|
a**********y 发帖数: 2367 | 38 并行计算?are you sure 你使用了matlab的并行计算吗,这个可不是自动的,是需要
编程解决的,matlab的并行很呆瓜,但对于一些固定的计算也可以将就一下
大型的并行计算fft很成熟了吧,并行效率很高,这个据说是最快的,但是你得会用http://www.fftw.org/ |
|
l********n 发帖数: 804 | 39 matlab 最新的版本应该fft是自动多核运算的
不知道你说的并行运算是不是这个
matlab慢的是for loop
你做time domain pseudospectral method?
估计也是要用for loop吧
C应该会快至少2-3倍
但是你这个要算几个星期太夸张了
可以考虑加GPU
我试过可以再快7-8倍 |
|
j***o 发帖数: 5096 | 40 没听说过严肃的计算数学问题用matlab搞的,不就是fft吗,又不是什么特殊的toolbox
,推翻用fortran或者c++重写吧,速度差了三到10倍,而且可以传到cluster上算,学
校里难道没有HPC? |
|
c*******h 发帖数: 1096 | 41 想要数量级上的速度提升,申请一个computing cluster的帐号是正道
不过fft其实不是特别scabable,因为communication开销太大,
所以你还得找一个网络特别快的cluster
换个数值算法是根本解决方法。。。要么就心安理得不要算那么大的n |
|
n*******l 发帖数: 2911 | 42 你是说FFT?估计很多人根本不在意那个复杂度问题吧,就知道很快就行了。 |
|
p*e 发帖数: 308 | 43 你问的是快速傅立叶变换FFT吧,跟傅立叶变换不是同一件事。 |
|
t*******t 发帖数: 1067 | 44 对,FFT,在不在意是一回事,但这应该知道的,如果要用的话。 |
|
f**********y 发帖数: 164 | 45 FFT有蝶形算法,可以并行计算,时间复杂度log(N) |
|
n*******l 发帖数: 2911 | 46 不知道你的新算法比现有的能有多大提高。我觉得,算法除非是像FFT那种开天辟地的,
小打小闹的提高个20-30%的效率,也有点意思但是不会让人惊艳。 |
|
a********a 发帖数: 61 | 47 如果楼主程序优化水准达到了 打败(比如大部分条件下性能高5%以上吧) FFTW 和
MKL FFT 的程度,找计算机这边的工作应该不是问题。 |
|
i*******g 发帖数: 168 | 48 到底谁逻辑有问题呀?
USC满校园的台湾香港人,一大半是在那些花钱就能上的 PROGRAM里面。
USC这样的水货PROGRAM,估计比整个CHICAGO的系加一起都多。
跟UCHICAGO比谁是精英,这不是搞笑吗? 只建设TOP 10的学科,这是UCHICAGO的
办学方针。所以UC没有EE,没有工学院,但是UC的经济系学生对FFT的算法发明
起了决定性的作用。USC工学院大的吓人,整个学校有几个TOP10的PROGRAM呢?
USC要是想和CHICAGO比,先和那些烂亲戚脱离关系,再把自己的门槛抬高点在说吧。
40%多的录取率,那也叫精英不是让人笑掉大牙吗?
一个还不够?一年毕业9个PHD,还有40多个在读的,比CHICAGO 的CS整个系的
博士生都多了,这样的例子,你只能到中国才能找到第二个。 |
|
Z*****o 发帖数: 4090 | 49 http://www.stereophile.com/features/368/index.html
By Robert Harley • Posted: Nov 1, 1993
Not that long ago, digital audio was considered perfect if all the bits
could be stored and retrieved without data errors. If the data coming off
the disc were the same as what went on the disc, how could there be a sound-
quality difference with the same digital/analog converter? This "bits is
bits" mentality scoffs at sonic differences between CD transports, digital
interfaces, and CD tweaks. Becaus... 阅读全帖 |
|
D**********s 发帖数: 3139 | 50 睁着眼睛说瞎话,fft频响相似的很,但是时间响度图明显polk糊得厉害。
如果要画waterfall decay的话,polk就是一直在谐振。
另外是用condenser mic对着woofer和tweeter中轴开外1m测的嘛。 |
|