Фазовый портрет сигнала по точкам

@Armlet · Регистрация: 03.01.2018

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

Здравствуйте, у меня есть сигнал формата .txt (время, напряжение) и нужно построить фазовый портрет. Нашёл в интернете разные примеры, но они конкретно для дифф.ур'ов. Вот такой код получился у меня, он выглядит очень коряво и фазовый портрет тоже, всё ли правильно?

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
import numpy as np
import matplotlib.pyplot as plt
 
M1 = np.loadtxt("1.txt")
t1 = M1[:,0]
s1 = M1[:,1]
n1 = len(s1)
tAdd = t1[25201]
Fs = 202665 #частота дискретизации
t0 = np.split(t1, 2)
t0_ = t0[0];
t__ = np.append(t0_, tAdd) #время размерности графика
FFTsig1 = np.fft.rfft(s1, n1)
FFTsig1_ = np.abs(FFTsig1)
Phase0 = np.angle(FFTsig1);
F = np.arange(0, Fs/2 + Fs/(n1*2), Fs/n1);
dPhase0 = np.diff(Phase0)/np.diff(t__)
pseudoPhase = np.delete(Phase0, 25200)
plt.plot(pseudoPhase, dPhase0)

@u235 · 29.01.2020, 16:53

Может имеется ввиду АФЧХ?
Сделали fft и стройте зависимость мнимой части от действительной.

Добавлено через 10 минут
Т.е. изображаете точки на комплексной плоскости.

@Armlet · 30.01.2020, 00:35 **[ТС]**

вроде нет, зависимость производной от функции, может я просто в коде что-то не так сделал. Ещё эти цифры пришлось вбивать(частоту дискретизации, количество отсчётов/2), эти цифры приходится сначала считать из файла и каждый раз они будут разные, вбивать не самое удобное решение. Когда я пытался вектор создать от 1 до n/2 у меня писало ошибку, в общем всё диковинно и что-то может быть не правильно, в целом нужно построить, но с этим вроде относительно не возникло проблем. Раньше я пытался это сделать в матлабе, попытался импортировать код, но выдавало ошибки в элементах P*:

Python
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
#from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt
 
#from matplotlib import rcParams
#rcParams['font.sans-serif']=['Arial']
M1 = np.loadtxt("1.txt")
t1 = M1[:,0]
s1 = M1[:,1]
s = np.array(M1[:,1])
t = np.array(M1[:,0])
n1 = len(s)
FFTsig1 = np.fft(s1, n1)
def func(s1, t1):
plt.plot(t1, s1, color='blue')
        FFTsig1 = np.fft(s1, n1)
 
Phase0 = np.angle(FFTsig1);
plt.plot(t1, Phase0)
P21 = abs(FFTsig1 / n1)
P11 = P21(1, n1/2)
P11[2, n1/2] = 2*P11(2, n1/2)
P41 = np.angle(FFTsig1 / n1);
P31 = P41(1, n1/2);
P31[2, n1/2] = 2*P31(2, n1/2);  # Фазовый спектр
Phase1 = np.angle(FFTsig1); # Фазовый спектр тоже
Amp1 = abs(FFTsig1);
fs = 1 / (1 * 10^-7); # частота дискретизации
F = 0, fs/2 - fs/n1, fs/n1;
       return FFTsig1;
plt.plot(t1, FFTsig1)

@u235 · 30.01.2020, 07:33

Строка 23:
Фаза не изменится поделите вы на n или нет.
Строка 25:
Зачем умножение на 2? Для амплитуды это еще как-то объяснимо, но не для фазы.
Строка31:
Что нужно изобразить? Какую зависимость? ФЧХ, АЧХ?
Если сделали FFT, то уже не время t1, а частота.

@Armlet · 31.01.2020, 00:25 **[ТС]**

сейчас мне всё прояснили и действительно нужно было реальную по x, мнимую часть аналитического сигнала по y откладывать, через преобразование гильберта

Новые блоги и статьи Все статьи Все блоги /
Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем	10 пpимет, которые всегда сбываются Maks 31.03.2026 1. Чтобы, наконец, пришла маршрутка, надо закурить. Если сигарета последняя, маршрутка придет еще до второй затяжки даже вопреки расписанию. 2. Нaдоели зима и снег? Не надо переезжать. Достаточно. . .	Перемещение выделенных строк ТЧ из одного документа в другой Maks 31.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .
Functional First Web Framework Suave DevAlt 30.03.2026 Sauve. IO Апнулись до NET10. Из зависимостей один пакет, работает одинаково хорошо как в режиме проекта так и в интерактивном режиме. из сложностей - чисто функциональный подход. Решил. . .	Автоматическое создание документа при проведении другого документа Maks 29.03.2026 Реализация из решения ниже выполнена на нетиповых документах, разработанных в конфигурации КА2. Есть нетиповой документ "ЗаявкаНаРемонтСпецтехники" и нетиповой документ "ПланированиеСпецтехники". В. . .	Настройка движения справочника по регистру сведений Maks 29.03.2026 Решение ниже реализовано на примере нетипового справочника "ТарифыМобильнойСвязи" разработанного в конфигурации КА2, с целью учета корпоративной мобильной связи в коммерческом предприятии. . . .	Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,760
	29.01.2020, 16:53
	Может имеется ввиду АФЧХ? Сделали fft и стройте зависимость мнимой части от действительной. Добавлено через 10 минут Т.е. изображаете точки на комплексной плоскости. 0

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,760
	30.01.2020, 07:33
	Сообщение было отмечено Armlet как решение Решение Строка 23: Фаза не изменится поделите вы на n или нет. Строка 25: Зачем умножение на 2? Для амплитуды это еще как-то объяснимо, но не для фазы. Строка31: Что нужно изобразить? Какую зависимость? ФЧХ, АЧХ? Если сделали FFT, то уже не время t1, а частота. 1

@Armlet 0 / 0 / 0 Регистрация: 03.01.2018 Сообщений: 14
	31.01.2020, 00:25 [ТС]
	сейчас мне всё прояснили и действительно нужно было реальную по x, мнимую часть аналитического сигнала по y откладывать, через преобразование гильберта 0

Фазовый портрет сигнала по точкам

Решение