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

@Magellan_ · Регистрация: 16.04.2011

Студворк — интернет-сервис помощи студентам

Необходимо построить амплитудно-фазовую частотную характеристику (АФЧХ) даны формулы $https://www.cyberforum.ru/cgi-bin/latex.cgi?A(\omega )=\frac{K\sqrt{({A}_{1}{\omega }^{2}+{A}_{3})^2+({A}_{0}{\omega }^{3}-{A}_{3}\omega)^2 }}{({A}_{1}{\omega }^{2}+{A}_{3})^2+({A}_{2}\omega -{A}_{0}\omega ^3)^2}$ и $https://www.cyberforum.ru/cgi-bin/latex.cgi?\varphi (\omega )=-arctg{T}_{1}\omega -arctg{T}_{2}\omega -arctg{T}_{3}\omega$ , $https://www.cyberforum.ru/cgi-bin/latex.cgi?\omega =0, 0.01, 0.05, 0.1, 0.5, 1, 2, 5, 10, 20, 50, 100, 1000, 10000, 1000000$

Visual Basic
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
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

Подскажите в том ли направлении я иду и если да, то как устранить существующие ошибки?Заранее спасибо!!!

@Magellan_ · 04.12.2011, 00:12 **[ТС]**

Так же 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

@Linoge · 05.12.2011, 11:20

Magellan_, а вас не смущает что вы говорите о афчх, а стоите годограф, который к слову, строится в полярных координатах? Или задача чисто по программированию.

Visual Basic
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
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

и кстати вместо w=array... для лучшей точности я бы поставил

Visual Basic
1
2
3
4
Dim w(1000)
For i = 0 To UBound(w)
    w(i) = 0.0001 * i ^ 2
Next

@Magellan_ · 05.12.2011, 17:40 **[ТС]**

А почему в одном случае у на sin, а в другом cos и не смогли бы в целом объяснить выполняемую роль (d * a * Cos(f), d * a * Sin(f))? Заранее спасибо!!!

@Linoge · 05.12.2011, 18:44

Это перевод из полярных в прямоугольные координаты.Годограф это траектория конца изменяющегося во времени вектора. Он задан длинной, в данном случае амплитудой, и углом - фазой.

X=A*sin(f)
Y=A*cos(f)

Коэффициент d введен здесь в качестве масштаба, т к значения x или y для данной системы меньше единицы, а функции scale принимает значение integer (получается 0) и это ей не нравится.

@Magellan_ · 05.12.2011, 19:25 **[ТС]**

Большое спасибо!!!

Новые блоги и статьи Все статьи Все блоги /
Изучаю kubernetes lagorue 13.01.2026 А пригодятся-ли мне знания kubernetes в России?	Сукцессия микоризы: основная теория в виде двух уравнений. anaschu 11.01.2026 https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/	WordPad для Windows 11 Jel 10.01.2026 WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .	Classic Notepad for Windows 11 Jel 10.01.2026 Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .
Почему дизайн решает? Neotwalker 09.01.2026 В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .	Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .

@Magellan_ 14 / 14 / 0 Регистрация: 16.04.2011 Сообщений: 236
	04.12.2011, 00:12 [ТС]
	Так же 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 0

@Magellan_ 14 / 14 / 0 Регистрация: 16.04.2011 Сообщений: 236
	05.12.2011, 17:40 [ТС]
	А почему в одном случае у на sin, а в другом cos и не смогли бы в целом объяснить выполняемую роль (d * a * Cos(f), d * a * Sin(f))? Заранее спасибо!!! 0

@Linoge 3205 / 1915 / 324 Регистрация: 25.10.2011 Сообщений: 5,564
	05.12.2011, 18:44
	Это перевод из полярных в прямоугольные координаты.Годограф это траектория конца изменяющегося во времени вектора. Он задан длинной, в данном случае амплитудой, и углом - фазой. X=Asin(f) Y=Acos(f) Коэффициент d введен здесь в качестве масштаба, т к значения x или y для данной системы меньше единицы, а функции scale принимает значение integer (получается 0) и это ей не нравится. 1

@Magellan_ 14 / 14 / 0 Регистрация: 16.04.2011 Сообщений: 236
	05.12.2011, 19:25 [ТС]
	Большое спасибо!!! 0