开个板很好啊 不管是 ptt2 的个板 还是脸书的私人页面 自己的空间总是很好的 当我说到个板时，想到的是 po 自己出去旅游的照片 想到的是 po 自己炒的菜 我还追星，我有很多本黄妃相簿 我没有把每一本黄妃相簿 po 在熟板 如果我这麽做，基於狂热的程度 搞不好也会有人叫我回自己个板 po 那就回去 po 吧.. 板上没有黄妃歌迷，在这里 po 是不会有共鸣的 我还开了黄妃社团，我当管理员，目前会员约三千人 po 在那边就不会被嫌弃，就会有共鸣 点阅率高的文章通常都会有一两千次的点阅 个板没什麽不好 问题不在个板啊 有的人和别人相处时，分享开心 独处时，也是自在开心 ---------------- 那如果我要 po 日本的感谢祭呢？ 我可以 po 在自己个板 也可以去 po 性板 八卦板口味重，大概也可以 问题出在特别的文章应该在对的地方 po 这是整个 ptt 都很常提的事 跑错板也是违反板规，所有板友都会摃上来，板主也得处理 我们不会说板友有问题，还罗织罪名套在板友的头上 对了，这次真的可以说'大家' 回想一下发生这种事时，会跳出来的人并不少 不是谁澎风，也不是分身 可不可以 po，制定板规的是板主，板友是检举，判断准则是拿捏在板主 因此有问题，摃上的是板主！想清楚喔！ 我从来不认为板主好当，有时甚至觉得，当不好就辞，溜了也好 （八卦板的板主就是高折损率） --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 111.82.25.183 (台湾) ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/RIPE_gender/M.1572212553.A.02F.html ※ 编辑: HuangJC (111.82.25.183 台湾), 10/28/2019 05:45:40 512 笔的傅立叶阵列完成了 我当初的问题是：输入 512 笔，输出也是同一个阵列 512 笔 可是输出怎麽看？ 有人回一句：时域转频域 这句我懂，如果我找一位老师也只教我这个，那不够细节 可是学问本来就要自己做，所以我努力很久终於搞懂了 https://webptt.com/cn.aspx?n=bbs/Math/M.1571655907.A.7B7.html 这边，发表在数学板，有程式有结果 然後凑答案发现，s[0] 要除以 4，才变成直流成份 1 而 s[1] 要除以 2，才变成振幅 1 为什麽这麽奇怪？但能回答我的人非常少 因为他们都只注重频率响应，而不需要关注直流成份 更甚且，除以 2 也很少人在做，因为他们要看的常是'特徵频率在哪里' 而不是振幅有多大 只有我在用的东西当然只有我研究 https://www.itread01.com/content/1544320336.html 不过这个网址里有和我相同的结论 他也没有推导，而是直接给结论，不知是否和我一样凑答案 反正就执行看看，然後看着结论凑一下就懂了 我接下来的挑战是要做五万笔阵列的傅立叶 同事马上叫我不要做，他说没人有这个需求，我一定是搞错需求了 但他後来还是给了我方向： 你的应用只用到FFT. 如果资料太多可以切段作转换 再将结果加权後相加。 FFT的长度看你需要的解析度决定 因为 你只需要显示频谱的大小而已，可以用上述方法进行。 以上教完了。 他搞错了，谁说我只要大小... 加权要加多少权，他也没说（应该是他不知道） 这边我还是得自己研究 但有人给我一个网页 http://www.fftw.org/ 说是最权威的傅立叶函式库 我同事给我的傅立叶明明很简短，翻页两个 page down 就没了 怎麽有人大费周章做函式库？！ 也许这种专业的有本事处理五万笔数据 行的话当然用，反正又不用钱 （软体工程师有很强烈的分享性格啊，一堆不用钱又强大的东西） 数学和程式的差别，在於程式必需再考虑速度和溢位 65536笔很有可能是过大的阵列 如果引起溢位或堆叠不足，那我这程式也不用跑了 直流是 a[0] 没错，我程式也付上了，跑了就明白 你是说，对称的部份先叠加在一起，然後除一样多的数？ 因为对称的关系，我是有一半的结果不看的 ※ 编辑: HuangJC (111.82.25.183 台湾), 10/28/2019 15:08:12 而且是以 s[256] 为中央向两侧对称，这个中央我还不太敢用 不要讲 256 点那麽多，我们来讲 16 点 （16hz 取样） 输入 s[i] = 1 + 1 * sin(x + 30.0 / 180 * M_PI) + 2 * sin(2 * x + 15.0 / 180 * M_PI) + 3 * sin(3 * x) + 4 * sin(4 * x) + 5 * sin(5 * x) + 6 * sin(6 * x); 因为我会让 x 走完一个周期, 所以 sin(x) 就是 1hz, sin(6x)就是 6hz 我後来才搞懂，所谓几 hz 那是我自己的解读 因为程式里根本没含入这件事 但我把 sin(x)解读为 1hz，这件事定义下来 也就会把输出时的 s[1] 定义为 1hz 的系数项 上面 30/180*PI 这部份算式，都只是在基频上转个角度 这样在观察共轭上会更方便 来看输出 after FFT s[0]=(15.999999,0.000000) = 15.999999 deg 0.000000 s[1]=(4.000000,-6.928203) = 8.000000 deg -60.000003 s[2]=(4.141104,-15.454813) = 15.999999 deg -75.000001 s[3]=(-0.000000,-24.000000) = 24.000000 deg 89.999999 s[4]=(0.000000,-32.000000) = 32.000000 deg -89.999999 s[5]=(0.000001,-39.999999) = 39.999999 deg -89.999999 s[6]=(0.000001,-47.999999) = 47.999999 deg -89.999999 s[7]=(0.000001,-0.000000) = 0.000001 deg -6.930822 s[8]=(0.000001,0.000000) = 0.000001 deg 0.000000 s[9]=(0.000001,0.000000) = 0.000001 deg 6.930820 s[10]=(0.000001,47.999999) = 47.999999 deg 89.999999 s[11]=(0.000001,39.999999) = 39.999999 deg 89.999999 s[12]=(0.000000,32.000000) = 32.000000 deg 89.999999 s[13]=(-0.000000,24.000000) = 24.000000 deg -89.999999 s[14]=(4.141104,15.454813) = 15.999999 deg 75.000001 s[15]=(4.000000,6.928203) = 8.000000 deg 60.000003 s[0] 是 0hz s[1] & s[15] 是 1hz s[2] & s[14] 2hz s[3] & s[13] 3hz . . . s[7] & s[9] 7hz s[8]........ 最好不要用 这就我说不敢用中央项的原因 ※ 编辑: HuangJC (111.82.25.183 台湾), 10/28/2019 15:39:55 我一开始就开了取样定理一个玩笑 当 s= sin(x); // 1hz 那我用 2hz 去取样，2hz 有达到两倍嘛，符合取样定理 结果取到全是 0... 所以我想是要'高於两倍'才安全，刚好两倍用起来还是毛毛的 又如上面 16 项的输出，中央项是 s[8], 按顺序推算是 8hz 的系数 那我们也可以实验一下含 8hz 的输入 s[i] = 1 + 1 * sin(x) + 2 * sin(2 * x) + 3 * sin(3 * x) + 4 * sin(4 * x) + 5 * sin(5 * x) + 6 * sin(6 * x) + 7 * sin(7 * x) + 8 * sin(8 * x); 这次不含角度，不观察共轭，只查看中央项的实用性 s[0]=(16.000000,0.000000) = 16.000000 deg 0.000000 s[1]=(0.000000,-7.999999) = 7.999999 deg -89.999999 s[2]=(0.000000,-15.999999) = 15.999999 deg -90.000000 s[3]=(-0.000000,-23.999999) = 23.999999 deg 89.999999 s[4]=(0.000000,-32.000000) = 32.000000 deg -90.000000 s[5]=(0.000000,-40.000000) = 40.000000 deg -90.000000 s[6]=(-0.000000,-48.000001) = 48.000001 deg 90.000000 s[7]=(0.000000,-56.000000) = 56.000000 deg -90.000000 s[8]=(-0.000000,0.000000) = 0.000000 deg 0.000000 <= 这项， 0 s[9]=(0.000000,56.000000) = 56.000000 deg 90.000000 s[10]=(-0.000000,48.000001) = 48.000001 deg -90.000000 s[11]=(0.000000,40.000000) = 40.000000 deg 90.000000 s[12]=(0.000000,32.000000) = 32.000000 deg 90.000000 s[13]=(-0.000000,23.999999) = 23.999999 deg -89.999999 s[14]=(0.000000,15.999999) = 15.999999 deg 90.000000 s[15]=(0.000000,7.999999) = 7.999999 deg 89.999999 8hz 的振幅应该是 1hz 的 8倍，这里却是 0 所以中央项我不敢用 ※ 编辑: HuangJC (111.82.25.183 台湾), 10/28/2019 15:57:47 取样 32hz 时，中央项已经变成 16 了 我不是在谈 8hz 解不解得出来，我是在谈 '中央项不要用' 所以你还是没用到中央项啊，不是嘛... 而且别项都有共轭两组，中央项没有 真要用你把它当什麽？怎麽解读呢？ 会不会它就是叠加，两组叠在一起 那既然是共轭，虚数项应该会消失，实数项变两倍？ 喔～ 你有在处理时序资料？ 那太好了，就是有在实战的战友 我的问题你应该都可以解了 目前我还在写 UI，暂时把数学搁置 因为 UI 不好，做对做错我也不知道.. 要再等等才有下一次的回馈.. ※ 编辑: HuangJC (111.82.25.183 台湾), 10/28/2019 17:16:12 我的问题解了啊，就是 s[0] 是什麽，s[1] 是什麽，s[2]是什麽... 之前只知道是频域，现在知道是直流成份，1hz, 2hz, 3hz...的参数 而且也知道了必需除以 N 倍 --- 下一个问题是我要写出高通，低通数位滤波器 不过这得等我把 UI 完成了 直接来的，这题完全不要去想前面的傅立叶 因为我要做五万笔的资料 五万笔变成频域，删掉不要的频率 再重新从频域转回时域 这种做法会疯掉 XD 理论是理论，我要的是五万笔刷～ for loop 跑完就做完高通，低通 程式我已经有了，拷来的 但跑起来不太对劲 所以等 UI 弄好再看问题 我看不懂你在打什麽，句子必需完整点 是的 因为没拿到程式 XD 我也看不懂你的句子 s[time] = v; //仪器所读到的电压，time 间隔单位是千分之一秒或百分之一秒 阵列从 0~5万 这好像你讲的两个都是嘛.. 如果是五万个阵列，这我不会提出来 因为五万个之间没有交互关系 我自己 for loop 解决就好 我的问题对懂的人来说会非常简单，不过一次谈一件就好 DaimlerHuang八成就是混着看文章，所以回了我看不懂的句子 频谱分析是频谱分析 带通是带通 两件我都要做，现在只谈频谱分析，而且我已经看懂怎麽做了 带通等我 UI 完成，试过有想法再来改 https://imgur.com/a/JTorrBn 程式介面就是这样 老板不会管我用高通低通合并起来做带通 反正做出来就好 这年头，写程式当然是先求有再求好 XD 如果带通其实就是高通加低通 那万一执行效率也一样呢？ 直接跑一个带通花的时间和跑一次高通再跑一次低通花的时间一样 如果是，我没必要管它们是不是一起的.. 所以搞懂弄出来就好 https://www.facebook.com/notes/%E9%BB%83%E7%91%9E%E6%98%8C/%E9%97%9C%E6%96%BC%E5%BB%BA%E7%AB%8B%E5%B7%A2%E7%8B%80-csplitterwnd/2696877440372540/ https://tinyurl.com/yymrh8jn 我举另一例子好了 这是我写的笔记 这例子我 google 到好几个地方有，他们都错成一样的 照着写程式都不会过 阿共那边现在是不是很多内容农场啊？都抄来抄去的.. 手上的程式有做到啊 其实我是竞争者，我不做的话目前已经有程式了 只需要对其中一个周期做FFT转换?(其它周期视做相同 交流电应该是你举的例，前面我举 1+sin(x) 也是我为了了解傅利叶而做的测试 我实际的资料不是周期波 因此我认为要做五万个资料做 FFT 转换 我同事认为不可能，是因为他觉得还没听过这种需求 关於需求的部份我现在也不是很清楚，我只知道有一串数字是五万多笔 这五万笔是十秒内还是一秒内，我也还在问 不过这不影响傅立叶不是嘛 反正这五万笔要喂进去 反正输出阵列的频率解析度，是根据输入阵列的频率来定义 现在或许可以赌：还真没有五万笔一次塞进去的需求 （这就是在猜老板要什麽了） 其实或许真的没需求吧！ 因为塞进去愈多笔的差别是可以分析出愈低频的资料 了不起就少了些超低频的频谱值 但我已经不想管老板要不要了 身为工程师，我突然任性的想：为什麽我拿到的程式不能塞五万笔？ 而我同事则是说：别闹了，快点赚钱，五万笔没人在要求的 我是一个工程师，我的工作是产出一支程式能做 而不是教客户：你去用别人的东西，把钱付给别人就好 就算 octave 能做，我也得挖出 octave 的程式包起来（它是 open source 吧） 变成执行档交差 而不能教客户说：来，装套 octave,我教你，它很简单 那我不成了 octave 的推广业务？ 不过如果你说的是： 必需要有开发工具，大家都推 matlab，我看你用 octave 当工具好了 我同意，有人说要教我 UI 完成了阶段任务，来做带通吧～ 上面有误会澄清的东西我消一消，然後开新的一篇回文 ※ 编辑: HuangJC (111.83.189.192 台湾), 10/29/2019 18:04:40