public_system.bas

来自「系统辨识」· BAS 代码 · 共 38 行

Attribute VB_Name = "Public_system"
Option Explicit
Public old_seita(10) As Double
Public Current_Y(-Order To 0) As Double
Public Current_U(-Order To 0) As Double
Sub system_ini()
    Dim k As Integer
    For k = -Order To 0
    Current_Y(k) = 0
    Current_U(k) = 0
    Next k
End Sub

Sub system()
    Dim k As Integer
    '注释 下列模拟系统生成数据的方法是正确的
    k = 0
    Current_U(k) = Rnd()
    '无噪声模型
    Current_Y(k) = 2.718 ^ (-10) * Current_Y(k - 1) + 22 * (1 - 2.71828 ^ (-10)) * Current_U(k - 6)
   ' Current_Y(k) = -0.3 * Current_Y(k - 1) - 0.15 * Current_Y(k - 2) + 0.2 * Current_U(k - 1) + 0.25 * Current_U(k - 2)
    '白噪声模型
    'Current_Y(k) = -0.1 * Current_Y(k - 1) - 0.15 * Current_Y(k - 2) + 0.2 * Current_U(k - 1) + 0.25 * Current_U(k - 2) + 0.01 * (Rnd() - 0.5)
    
    '有色噪声模型
    'Current_Y(k) = -0.1 * Current_Y(k - 1) - 0.15 * Current_Y(k - 2) + 0.2 * Current_U(k - 1) + 0.25 * Current_U(k - 2) + 0.01 * Rnd()
    
    
    For k = -Order To -1
    Current_Y(k) = Current_Y(k + 1)
    Current_U(k) = Current_U(k + 1)
    Next k
    
    
    
End Sub