[發問] 如何判斷週期訊號?

acdx

我有一組數字fm1=32.2, fm2=44.3, fm3=63.0, fm4=79.1, fm5=95.2, fm6=111.3
我以陣列的型式將它寫入成
fm(0)=32.2
fm(1)=44.3
fm(2)=63.0
fm(3)=79.1
fm(4)=95.2
fm(5)=111.3
我要判斷這幾個數之間是否有一固定週期(間隔),並排除其中個別週期 例如:
fm(1)-fm(0)=12.1
fm(2)-fm(1)=18.7
fm(3)-fm(2)=16.1
fm(4)-fm(3)=16.1
fm(5)-fm(4)=16.1
所以固定週期(間隔)應為16.1, 而fm(1)=44.3這個數顯然有偏差故將它排除
以上的例子直接看很容易判斷 但是該如何讓成是來做判斷呢?
請板上的高手給個建議 謝謝

stillfish00

回復 1# acdx
例子怪怪的 若固定周期16.1，則 fm(1) 應為 46.9；fm(0) 應為 30.8
這樣兩個都有偏差都要排除?

acdx

回復 2# stillfish00


原始數據為量測結果 我的程式是為了協助分析
所以有些許誤差是可接受的(+/- 1.5以內可視為相同) 但是要排除個別的雜訊

Hsieh

回復 3# acdx

1. Sub nn()
   
2. Dim Ar()
   
3. k = 1.5 '容許誤差
   
4. fm = Array(32.2, 44.3, 63, 79.1, 95.2, 111.3) '常數陣列
   
5. Set d = CreateObject("Scripting.Dictionary")
   
6. For i = 0 To UBound(fm)
   
7. If i > 0 Then
   
8. ReDim Preserve Ar(s)
   
9. Ar(s) = fm(i) - fm(i - 1)
   
10. s = s + 1
    
11. End If
    
12.    d(i) = fm(i)
    
13. Next
    
14. n = Application.Mode(Ar) '最常出現差值
    
15. For i = LBound(fm) + 1 To UBound(fm) - 1
    
16.    If d(i) - d(i - 1) < n - k Or d(i + 1) - d(i) > n + k Then
    
17.    d.Remove i
    
18.    End If
    
19. Next
    
20. MsgBox Join(d.items, Chr(10))
    
21. End Sub

複製代碼

acdx

回復 4# Hsieh
我利用 n = Application.Mode(Ar) 來判斷最常出現的值
若是遇到整組數字無重複的值(無週期訊號時)程式則會出現error
這該如何修正呢?

stillfish00

回復 5# acdx
用回歸直線的方法排除差異過大的值，看看是否有幫助，
須注意這裡用的回歸直線是包含偏差值的資料所推出的直線，
所以deltaY值請自行斟酌

1. Sub Test()
   
2.     Dim d, fm
   
3.     Dim a, b, it
   
4.     Dim deltaY
   
5.    
   
6.     fm = Array(32.2, 44.3, 63, 79.1, 95.2, 111.3)
   
7.     Set d = CreateObject("scripting.dictionary")
   
8.     For i = 0 To UBound(fm)
   
9.         d(i) = fm(i)
   
10.     Next
    
11.    
    
12.     '回歸直線  fm=ax+b
    
13.     a = WorksheetFunction.Index(WorksheetFunction.LinEst(d.items, d.keys), 1)
    
14.     b = WorksheetFunction.Index(WorksheetFunction.LinEst(d.items, d.keys), 2)
    
15.    
    
16.     deltaY = 1.5  '與回歸線Y軸上的容許誤差
    
17.     For Each it In d.keys
    
18.         '排除誤差過大的值
    
19.         If Abs((it * a + b) - d(it)) > deltaY Then d.Remove it
    
20.     Next
    
21.    
    
22.     MsgBox Join(d.items, ",")
    
23.     Set d = Nothing
    
24. End Sub

複製代碼

acdx

回復 6# stillfish00
後來我把它改成這樣也可以

Sub nn()
Dim Ar()
k = 1.5 '容許誤差
fm = Array(32.2, 44.3, 63, 79.1, 95.2, 111.3) '常數陣列
Set d = CreateObject("Scripting.Dictionary")
For i = 0 To UBound(fm)
If i > 0 Then
ReDim Preserve Ar(s)
Ar(s) = fm(i) - fm(i - 1)
s = s + 1
End If
   d(i) = fm(i)
Next
n = Application.Mode(Ar) '最常出現差值
If IsNumeric(n) = False Then
    n = 0
    MsgBox "無週期訊號"
   Exit sub
End If
For i = LBound(fm) + 1 To UBound(fm) - 1
   If d(i) - d(i - 1) < n - k Or d(i + 1) - d(i) > n + k Then
   d.Remove i
   End If
Next
MsgBox Join(d.items, Chr(10))
End Sub