作者:
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如果可以用工具我還談什麼五萬筆無法解傅立葉?所以我很前面就說了：數學和程式的差別，在於程式必需思