Нелинейная регрессия (оценка частоты синуса)

@jetairnav · Регистрация: 25.07.2016

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

Час добрый, не работает (не оценивается частота) следующий код:

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
32
33
34
35
36
37
38
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import least_squares
 
 
def generate_data(t, A0, A, omega, phase=0, noise=0, n_outliers=0, random_state=0):
    y = A0 + A * np.sin(omega * t + phase)
    rnd = np.random.RandomState(random_state)
    error = noise * rnd.randn(t.size)
    outliers = rnd.randint(0, t.size, n_outliers)
    error[outliers] *= 30
    return y + error
 
 
#  Инициализация значений модельной функции
A = 7.744444
A0 = 0
omega = 2*np.pi*1 #2*np.pi*0.2
noise = 0.0
phase = 0
n_outliers = 0
t_min = 0
t_max = 100
 
t_model = np.linspace(t_min, t_max, 100)
y_model = generate_data(t_model, A0= A0, A=A, omega=omega, phase=phase, noise=noise, n_outliers=n_outliers, random_state=1)
plt.plot(t_model,y_model)
 
 
def fun(x, t, y):
    return x[0]*np.sin(x[1] * t ) - y
 
 
x0 = np.ones(4)
res_lsq = least_squares(fun, x0, args=(t_model, y_model))
print('res_lsq.x = ', res_lsq.x)
print(omega)
plt.show()

Не понимаю, в чем причина..

@jetairnav · 21.11.2018, 14:02 **[ТС]**

Нашел ошибку... Конечно, t_max должно быть равно 1))

Новые блоги и статьи Все статьи Все блоги /
Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .	Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .
модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .	Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка. Hrethgir 08.04.2026 Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .	Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ

@jetairnav 2 / 1 / 1 Регистрация: 25.07.2016 Сообщений: 3
	21.11.2018, 14:02 [ТС]
	Нашел ошибку... Конечно, t_max должно быть равно 1)) 0