Амплитудно - фазовая частотная характеристика

Function Amp(ByVal A As Variant)
A = K * Sqr((A1 * x(i) ^ 2 + A3) ^ 2 + (A0 * x(i) ^ 3 - A3 * x(i)) ^ 2) / ((A1 * x(i) ^ 2 + A3) ^ 2 + (A2 * x(i) - A0 * x(i) ^ 3) ^ 2)
Amp = A
End Function
Function Faz(ByVal F As Variant)
F = -Atn(T1 * x(i)) - Atn(T2 * x(i)) - Atn(T3 * x(i))
Faz = F
End Function
 
Private Sub Command1_Click()
Picture1.Scale (-100, 100)-(100, -100)
DrawWidth = 2
Line (-100, 0)-(100, 0)
Line (0, 100)-(-100, 0)
DrawWith = 1
x = Array(0, 0.01, 0.05, 0.1, 0.5, 1, 2, 5, 10, 20, 50, 100, 1000, 10000, 1000000)
For i = LBound(x) To UBound(x)
Line -(Amp(A), Faz(F))
Next i
End Sub

Dim A0, A1, A2, A3, T1, T2, T3, K, d
 
Function Amp(ByVal w As Double)
Amp = K * Sqr( _
((A1 * w ^ 2 + A3) ^ 2 + (A0 * w ^ 3 - A3 * w) ^ 2)) _
/ ((A1 * w ^ 2 + A3) ^ 2 + (A2 * w - A0 * w ^ 3) ^ 2)
End Function
Function Faz(ByVal w As Double)
Faz = -Atn(T1 * w) - Atn(T2 * w) - Atn(T3 * w)
End Function
 
Private Sub Form_Load()
'd введен для того чтоб scale не ругался на слишком малые (нулевые) значения
d = 100
K = 1.8
T1 = 10
T2 = 3
T3 = 1
A0 = T1 * T2 * T3
A1 = T1 * T2 + T2 * T3 + T1 * T3
A2 = T1 + T2 + T3
A3 = 1
End Sub
 
Private Sub Command1_Click()
 
w = Array(0, 0.01, 0.05, 0.1, 0.5, 1, 2, 5, 10, 20, 50, 100, 1000, 10000, 1000000)
'нахождение максимальных размеров графика
min_xx = 0
max_xx = 0
min_yy = 0
max_yy = 0
For i = LBound(w) To UBound(w)
    a = Amp(w(i))
    f = Faz(w(i))
    If d * a * Cos(f) < min_xx Then min_xx = d * a * Cos(f)
    If d * a * Cos(f) > max_xx Then max_xx = d * a * Cos(f)
    If d * a * Sin(f) < min_yy Then min_yy = d * a * Sin(f)
    If d * a * Sin(f) > max_yy Then max_yy = d * a * Sin(f)
    Print a & " " & f
Next
 
If Abs(min_xx) > max_xx Then
    xx = Int(Abs(min_xx) * 1.2)
Else
    xx = Int(max_xx * 1.2)
End If
 
If Abs(min_yy) > max_yy Then
    yy = Int(Abs(min_yy) * 1.2)
Else
    yy = Int(max_yy * 1.2)
End If
'задание параметров и координат
Me.Scale (-xx, yy)-(xx, -yy)
DrawWidth = 2
Line (-xx, 0)-(xx, 0)
Line (0, yy)-(0, -yy)
DrawWith = 1
'построение графика
a = Amp(w(0))
f = Faz(w(0))
Line (d * a * Cos(f), d * a * Sin(f))-(d * a * Cos(f), d * a * Sin(f))
For i = LBound(w) To UBound(w)
    a = Amp(w(i))
    f = Faz(w(i))
    Line -(d * a * Cos(f), d * a * Sin(f)), vbRed
Next
End Sub

Dim w(1000)
For i = 0 To UBound(w)
    w(i) = 0.0001 * i ^ 2
Next