楼主:
HuangJC (吹笛牧童)
2019-10-28 05:42:29开个板很好啊
不管是 ptt2 的个板
还是脸书的私人页面
自己的空间总是很好的
当我说到个板时，想到的是 po 自己出去旅游的照片
想到的是 po 自己炒的菜
我还追星，我有很多本黄妃相簿
我没有把每一本黄妃相簿 po 在熟板
如果我这么做，基于狂热的程度
搞不好也会有人叫我回自己个板 po
那就回去 po 吧..
板上没有黄妃歌迷，在这里 po 是不会有共鸣的
我还开了黄妃社团，我当管理员，目前会员约三千人
po 在那边就不会被嫌弃，就会有共鸣
点阅率高的文章通常都会有一两千次的点阅
个板没什么不好
问题不在个板啊
有的人和别人相处时，分享开心
独处时，也是自在开心
楼主:
HuangJC (吹笛牧童)
2019-10-28 05:50:00上一个会在标题打我 ID的，我记得是 D姐..两个果然半斤八两
作者:
joned (极度肥胖)
2019-10-28 07:38:00上一个会叫我去别的板不要在这个是D姐 两人要在一起?
作者:
st3336 (Strance)
2019-10-28 10:19:00傅立叶玩出心得了吗?
这里已经有点像个板了....题外话,你用的傅立叶是快速复利叶,资料点数量必须是2的次方,近6万笔数字,使用65536输入输出实部虚部各65536笔,至于直流部分,我不懂理论上有频率就不算是直流,直流大概是a[0]
作者:
st3336 (Strance)
2019-10-28 14:27:00其实应该都要除4 因N=4 叠加四倍 1与3共轭 振福2叠加/4=1你的取样点就是512 实际上1~256 对称 257-512光看FFTW的feature基本上就很猛了 快速转换/多工是重点你看s[1]跟s[3]是共轭的 实际上能量属于同一个频率总之在公式推导的结论上 s[1]+s[3] 代表真正的震幅而且转换是一个叠加过程 你有n个点就会sigma n次所以最后要除n你math版有一个测试512的 你如果会画图 你画1-512 就知道更精确地说 s[0]是第一点 2-512震幅谱会对称s[0]就是你所谓的DC level并不是中央向不要用 而是不要用到nyquist frequency如果你硬要他取样16的话 8当然解不出来 取32就解决了以我处理时序资料为例子 不会遇到这么低取样率的主要还是端看你要处理的资料类型我是没自己写拉 现成很多东西都写好了 主要玩的还在后面
楼主:
HuangJC (吹笛牧童)
2019-10-28 17:28:00我也无法再笔算积分，只想'叫用'已写好的快速傅立叶但光是使用，这些讨论就要懂啊，不然不会解读传回值请问你直流成份，各频率系数是多少？我也只要这程度
作者:
st3336 (Strance)
2019-10-28 18:07:00你的问题是什么 我用matlab带你的case 结果是一样的ˊ顶多差在截断误差而已 比如说7.9999 跟 8 的差异没什么参数吧 也不一定是所谓 1 2 3HZ1 2 3端看你的点数/频率取样绿高通低通就是写transfer function 看你要几个pole这样看你是要频谱相乘 还是直接转时域做convolution
删掉不要的频率?五万多笔数据怎么能只宣告[512]???讨论提到取样,Nyquist却违背取样定理?
作者:
st3336 (Strance)
2019-10-29 11:15:00实际上你把滤波器转时域做convolution会比转五万比傅立叶来的快你又要高通又要低通 干嘛不直接带通我这边不太懂你是五万笔时间序列 还是一笔五万点的时间序五万笔 N点未知的时间序列 vs 一笔 五万点的时间序列这里这样讲 就是一笔 五秒的资料 取样率1ms 这样那就简单多了 你也只是要做频谱分析而已吧维基跟我想的差不多 带通就是频域高通*低通我说维基百科的 你之前问得我在维基百科都有看到
楼主:
HuangJC (吹笛牧童)
2019-10-29 13:42:00我说我有拿到程式，这些 google 不难
作者:
st3336 (Strance)
2019-10-29 13:42:00我只是单纯觉得带通 应该就是频率域的高通*低通
作者:
st3336 (Strance)
2019-10-29 13:44:00频率滤波器应该没有太多分歧吧 我看都是transfer func
楼主:
HuangJC (吹笛牧童)
2019-10-29 13:45:00等等我附上我研究到一半的东西你就知道我说哪里没解决了就像一个傅立叶，你说程式做的和 matlab 一样一样没告诉我必需除以 N 啊，我必需除以 N 才能用而直流成份必需除以 2N，特别不一样咧所以大部份谈数位滤波器的，都是附一个多项式公式然后跑了 matlab 说：看，这式子有用，它就是一阶滤波
作者:
st3336 (Strance)
2019-10-29 13:49:00所以说 你的程式只做到一半 还没有达到理论推导的结果
楼主:
HuangJC (吹笛牧童)
2019-10-29 13:49:00但他们附的程式其实怪怪的 XDDD我手上的程式输出会去修改输入阵列，我看不懂这什么逻辑
作者:
st3336 (Strance)
2019-10-29 13:51:00阵列大小吧 要砍一半
楼主:
HuangJC (吹笛牧童)
2019-10-29 13:51:00如果他想表达这叫 IIR, 程式设计技巧不该是这样的
楼主:
HuangJC (吹笛牧童)
2019-10-29 13:52:00因这种改变输入阵列的做法催毁原值，那我要另外保留原值?低通我做了有用，高通没用，我准备留三套阵列了所以那程式有必要修正..
作者:
st3336 (Strance)
2019-10-29 14:01:00其实5万点 应该也只是一秒内的事情
ms间隔取样5万笔,所以收了50秒?1秒1000个资料台湾交流电60Hz,一个周期有1000/60个资料你需要对完整5万个资料做FFT转换? 还是
作者:
st3336 (Strance)
2019-10-29 14:18:00我是觉得用octave就好了
楼主:
HuangJC (吹笛牧童)
2019-10-29 16:21:00所以我的工作其实就是读档，数学，绘图这三个动作去操作 octave 是可以完成，但档案有它的格式这个读档就是我程式的价值；本来就很简单，已完成我是在补数学这一块；对懂数学的人来说，它很简单
作者:
st3336 (Strance)
2019-10-29 17:17:00如果是要制作产品那就另当别论了 我一直以为你只是要研究其实也不用纠结五万笔数值 一笔时间序列本来就是可变只要dimention开够大 就不会出错你的结果只要比对现有程式的结果 答案一样就能用了你能用C写GUI 我还蛮佩服的
楼主:
HuangJC (吹笛牧童)
2019-10-29 18:05:00如果可以用工具我还谈什么五万笔无法解傅立叶?所以我很前面就说了：数学和程式的差别，在于程式必需思