Форум программистов, компьютерный форум, киберфорум
Python: Научные вычисления
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.57/7: Рейтинг темы: голосов - 7, средняя оценка - 4.57
7 / 7 / 1
Регистрация: 24.01.2017
Сообщений: 229

Полиномиальная регрессия. Не удаётся подобрать коэффициенты

04.05.2019, 12:28. Показов 1729. Ответов 1

Студворк — интернет-сервис помощи студентам
Пробовал для полиномов меньшей степени. Код работает хорошо, всё подбирается. А вот для полинома 9 степени перестаёт работать, после 1000 эпох все коэффициенты равны nan.

Код для jupyter notebook, а так же сам ноутбук прилагаю:

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
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
import tensorflow as tf
import numpy as np
from tqdm import tqdm_notebook as tqdm
from matplotlib import pyplot as plt
 
learning_rate = 0.01
training_epoches = 500
 
y_coeffs = [12, 5, 5, 245, 3, 59, -190, 99, 8]
num_coeffs = len(y_coeffs)
 
x_train = np.linspace(-16, 11, 1000)
 
y_train = 0
for i in range(num_coeffs):
    y_train += y_coeffs[i] * np.power(x_train, i)
 
y_train += np.random.randn(*x_train.shape) * 50000000.4
 
fig = plt.figure(figsize=(10, 10))
plt.scatter(x_train, y_train)
plt.show()
 
X = tf.placeholder(tf.float32)
Y = tf.placeholder(tf.float32)
w = tf.Variable([0.0] * num_coeffs, name='weights')
 
def F(x, w):
    terms = []
    for i in range(num_coeffs):
        term = tf.multiply(w[i], tf.pow(x, i))
        terms.append(term)
        
    return tf.add_n(terms)
 
model = F(X, w)
cost = tf.pow(Y - model, 2)
 
optimizer = tf.train.AdadeltaOptimizer(learning_rate)
train_operation = optimizer.minimize(cost)
 
sess = tf.Session()
init = tf.global_variables_initializer()
sess.run(init)
 
for epoche in tqdm(range(training_epoches), desc='epoches: '):
    for (x, y) in zip(x_train, y_train):
        sess.run(train_operation, feed_dict={X : x, Y : y})
        
w_end = sess.run(w)
sess.close()
 
fig = plt.figure(figsize=(10, 10))
plt.scatter(x_train, y_train)
 
trY2 = 0
for i in range(num_coeffs):
    trY2 += w_end[i] * np.power(x_train, i)
    
plt.plot(x_train, trY2, linewidth=5, color='green')
plt.show()
 
print(w_end)
В данном коде использовал только 500 эпох. При этом левая часть графика никак не хочет аппроксимировать.





Подскажите, где ошибка?
Вложения
Тип файла: zip PolynomialRegression.ipynb.zip (87.9 Кб, 0 просмотров)
0
Programming
Эксперт
39485 / 9562 / 3019
Регистрация: 12.04.2006
Сообщений: 41,671
Блог
04.05.2019, 12:28
Ответы с готовыми решениями:

Полиномиальная регрессия используя Градиентный Спуск и Матричный Способ
Здравствуйте! Мне нужно было применить линейную и полиномиальную регрессии (без использования Scikit-Learn и специальных функций, например...

Полиномиальная регрессия
Здравствуйте, уважаемые знатоки математики и программирования!! Необходимо Разработать программу построения графика экспериментальных...

Полиномиальная регрессия
Помогите найти ошибку. Чет не пойму где ошибся.

1
7 / 7 / 1
Регистрация: 24.01.2017
Сообщений: 229
06.05.2019, 11:03  [ТС]
Всем спасибо.
Просто нужно было нормализовать данные и применить регуляризацию.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
inter-admin
Эксперт
29715 / 6470 / 2152
Регистрация: 06.03.2009
Сообщений: 28,500
Блог
06.05.2019, 11:03
Помогаю со студенческими работами здесь

Полиномиальная регрессия в Python
Всем привет! Может кто-нибудь помочь с написанием кода? Нужно написать реализацию класса для обучения полиномиальной регрессии, то...

Полиномиальная логистическая регрессия из документации sklearn
Здравствуйте, уважаемые форумчане! Изучаю полиномиальную логистическую регрессию и только нашёл только 1 код по ней, а именно - из...

Подобрать коэффициенты в выражении
Добрый день! Подскажите, пожалуйста, можно ли средствами MATLAB решить следующую задачу. Необходимо сделать так, чтобы система...

Подобрать коэффициенты для описания функции с помощью гиперболического синуса
Здравствуйте! Есть такая небольшая проблема: необходимо подобрать коэффициенты для описания функции с помощью гиперболического синуса, то...

Как подобрать коэффициенты для многочлена, чтобы разложить его на множители?
Всем привет. Подскажите как подобрать коэфф. чтобы разложить данный многочлен: {x}^{2}+x-\sqrt{2x}+2 И вообще, можно ли проверить,...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Новые блоги и статьи
модель ЗдравоСохранения 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
Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики.
Hrethgir 07.04.2026
Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .
Горе от ума
kumehtar 07.04.2026
Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .
Использование значений реквизитов справочника в документе, с определенными условиями и правами
Maks 07.04.2026
1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .
Доступность команды формы по условию
Maks 07.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru