Нахождение периода по графику

@Myrcella · Регистрация: 21.03.2016

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

Была построена коррелограмма.
Теперь нужно найти период на этой коррелограмме,т.е. что бы программа сама его подсчитывала.
Был написан код.
Но выдает ошибку : Undefined function 'zi' for input arguments of type 'double'.
Сам код для нахождения периода:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
dt=1/f0;%период дискретизации сигнала
TN=2;
Na=0:TN;
is=0;
A=size(rnorm_2);
z=sign(A);
 
for i=2:Na;%организовываем цикл по числу элементов
if z(i)==0;  
    then z(i)=1;
if (z(i)*z(i-1))<0;
then is=is+1; 
 zi(is)=i;
    
end
end
end
T(i)=((zi(is)-zi(1))*dt)/(is-1);

Кто может помочь разобраться с ошибкой или подтолкнуть на мысль, как найти период сигнала по графику

@Centurio · 23.10.2016, 19:11

Сообщение от Myrcella

Undefined function 'zi'

Программа так и пишет: "Не определена функция zi". Откуда она взялась, если всё время использовалась z?

@Nick07 · 23.10.2016, 20:17

Ошибка выдается в 13-й или в 18-й строке?

@Centurio · 23.10.2016, 20:44

Myrcella, вот в этой строке
zi(is)=i;
может статься, что индекс окажется равным нулю (у вас is начинает отсчитываться с нуля), что приведёт к ошибке.

@Myrcella · 23.10.2016, 20:59 **[ТС]**

Сообщение от Nick07

Ошибка выдается в 13-й или в 18-й строке?

18 строка

Добавлено через 7 минут

Сообщение от Centurio

Программа так и пишет: "Не определена функция zi". Откуда она взялась, если всё время использовалась z?

Да, так и пишет.
Код еще был написан в прошлом году к диплому, но не использовался, т.к. эта же ошибка выдавалась
Поэтому точно не помню откуда брали zi.
Попробую поменять на z и is не приравнивать к 0.

@Centurio · 23.10.2016, 21:11

Myrcella, у вас ещё и цикл не будет выполняться, так как он прописан для i от 2 до 0 с положительным приращением.

@Myrcella · 23.10.2016, 21:13 **[ТС]**

Сообщение от Centurio

Myrcella, у вас ещё и цикл не будет выполняться, так как он прописан для i от 2 до 0 с положительным приращением.

Легче заново сесть и написать новый код?
Просто я даже не представляю как сделать, что бы он посчитал период.

@Centurio · 23.10.2016, 21:28

Myrcella, скорее всего, ошибка именно в задании цикла. Судя по всему, i должна принимать значение от 2 до length(A).
Предлагаю следующий код для замены строк 6-17:

Matlab M
1
2
3
z=2*(A>=0)-1;
d=diff(z);
zi=find(d)+1;

Добавлено через 3 минуты
17 строка тоже непонятная. В ней используются индексы уже за пределами цикла.

@R2D2 · 23.10.2016, 21:29

Myrcella, я, конечно, извиняюсь, но нельзя ли привести Вашу задачу к использованию данной функции?

@Myrcella · 24.10.2016, 10:46 **[ТС]**

Сообщение от R2D2

Myrcella, я, конечно, извиняюсь, но нельзя ли привести Вашу задачу к использованию данной функции?

Когда искала, как в Матлабе посчитать период сигнала, то данный функции не нашла. Смотрела в 2-х или 3-х книгах.
Сейчас попробовала, тоже не хочет считать.

@Centurio · 24.10.2016, 10:54

18 строку в дополнение к коду, написанному мной выше, можно записать в виде

Matlab M
1
T=(zi-zi(1)*dt)./(0:length(zi)-1);

@Myrcella · 24.10.2016, 11:41 **[ТС]**

Сообщение от Centurio

18 строку в дополнение к коду, написанному мной выше, можно записать в виде
Matlab MВыделить код
1
T=(zi-zi(1)*dt)./(0:length(zi)-1);

Правильно я поняла, что код должен выглядеть так:

Matlab M
1
2
3
4
5
6
7
8
9
dt=1/f0;
A=size(rnorm_2);
z=2*(A>=0)-1;
 
for i=2:length(A);
d=diff(z);
zi=find(d)+1;
end
T=(zi-zi(1)*dt)./(0:length(zi)-1);

Первый раз работаю с циклами и не очень хорошо в них разбираюсь.

@Centurio · 24.10.2016, 12:38

Myrcella, здесь цикл не нужен. Кроме того, период - это расстояние между двумя точками, в которых функция одинаково меняет знак, т.е. и в первой, и во второй точке либо с плюса на минус, либо с минуса на плюс. Поэтому код будет такой.:

Matlab M
1
2
3
4
5
dt=1/f0;
A=size(rnorm_2);
d=diff(z<0);
zi=find(d)+1;
T=(zi-zi(1)*dt)./(0:length(zi)-1);

@Myrcella · 24.10.2016, 13:42 **[ТС]**

Centurio,
теперь выдает ошибку:
Attempted to access zi(1); index out of bounds because numel(zi)=0.
Как я понимаю, это уже на массив он ругается

@Centurio · 24.10.2016, 14:18

Myrcella, А, ну да. Там происходит деление на ноль. Скорее всего, последняя строчка должна быть

Matlab M
1
T=(zi-zi(1)*dt)./(1:length(zi));

Вообще, неплохо было бы, если бы предоставили данные коррелограммы и более подробно объяснили, каким способом вы хотите найти период.

@Myrcella · 24.10.2016, 17:24 **[ТС]**

Centurio,
Задание было: Найти средний период колебаний П, который вычисляется путем измерения периодов колебаний на коррелограмме за определенный временной отрезок и последующего усреднения полученных периодов.

Cначала хотела вообще делать через fft, но сказали, что надо находить период не пользуясь преобразованием Фурье.
Поэтому думала находить через пики.
Но могу что-нибудь перепутать, т.к. это было еще в мае, а бумаги с описанием этого кода куда-то дела.

@Centurio · 24.10.2016, 17:50

Myrcella, а данные, по которым строилась коррелограмма, есть?

@Myrcella · 24.10.2016, 17:55 **[ТС]**

Centurio, да, есть, но они в формате "edf".

@Centurio · 24.10.2016, 19:31

Myrcella, выложите.

@Myrcella · 24.10.2016, 19:50 **[ТС]**

Centurio,

и параметры, которые были взяты с сайта

Matlab M
1
2
3
4
5
%Задаем частоту дискретизации
f0=165; %Гц
 
%Задаем время наблюдения
Tn=20; %c

Новые блоги и статьи Все статьи Все блоги /
Программный контроль заполнения реквизита табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать контроль заполнения реквизита табличной части. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .
Функция заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	23.10.2016, 20:17
	Ошибка выдается в 13-й или в 18-й строке? 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	23.10.2016, 20:44
	Myrcella, вот в этой строке zi(is)=i; может статься, что индекс окажется равным нулю (у вас is начинает отсчитываться с нуля), что приведёт к ошибке. 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	23.10.2016, 21:11
	Myrcella, у вас ещё и цикл не будет выполняться, так как он прописан для i от 2 до 0 с положительным приращением. 0

@R2D2 939 / 839 / 121 Регистрация: 23.11.2012 Сообщений: 2,488
	23.10.2016, 21:29
	Myrcella, я, конечно, извиняюсь, но нельзя ли привести Вашу задачу к использованию данной функции? 0

@Myrcella 0 / 0 / 0 Регистрация: 21.03.2016 Сообщений: 9
	24.10.2016, 13:42 [ТС]
	Centurio, теперь выдает ошибку: Attempted to access zi(1); index out of bounds because numel(zi)=0. Как я понимаю, это уже на массив он ругается 0

@Myrcella 0 / 0 / 0 Регистрация: 21.03.2016 Сообщений: 9
	24.10.2016, 17:24 [ТС]
	Centurio, Задание было: Найти средний период колебаний П, который вычисляется путем измерения периодов колебаний на коррелограмме за определенный временной отрезок и последующего усреднения полученных периодов. Cначала хотела вообще делать через fft, но сказали, что надо находить период не пользуясь преобразованием Фурье. Поэтому думала находить через пики. Но могу что-нибудь перепутать, т.к. это было еще в мае, а бумаги с описанием этого кода куда-то дела. Миниатюры 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	24.10.2016, 17:50
	Myrcella, а данные, по которым строилась коррелограмма, есть? 0

@Myrcella 0 / 0 / 0 Регистрация: 21.03.2016 Сообщений: 9
	24.10.2016, 17:55 [ТС]
	Centurio, да, есть, но они в формате "edf". 0

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	24.10.2016, 19:31
	Myrcella, выложите. 0