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

@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))

Новые блоги и статьи Все статьи Все блоги /
My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .	Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .	Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .
Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .	Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .	Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача №1: при указании работ (справочник РаботыПоРемонтуСпецтехники),. . .

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

Опции темы