реклама на сайте
подробности

 
 
> метод наименьших квадратов, как оценить необходимую разрядность представления чисел
net
сообщение Apr 19 2009, 13:15
Сообщение #1


Знающий
****

Группа: Свой
Сообщений: 858
Регистрация: 9-08-04
Пользователь №: 473



есть набор данных эксперимента N*{x,y}
нужно подобрать коэффициенты полинома степени M для аппроксимации эксперимента

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


ВОПРОС
как оценить необходимую разрядную сетку представления чисел в ЭВМ необходимую для обеспечения данного метода


в литературе как то смутно написано что типа работает и при больших степенях полинома
но хотелось бы
1 посмотреть на вывод данного утверждения
2 посмотреть формулы для разрядности представления чисел не влияющей на точность аппроксимации

или дайте ссылку с данными выводами
спасибо
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Xenia
сообщение Nov 29 2011, 14:47
Сообщение #2


Гуру
******

Группа: Модератор FTP
Сообщений: 4 479
Регистрация: 20-02-08
Из: Москва
Пользователь №: 35 237



Кто подскажет, как из коэффициентов разложения в ряд Чебышева, можно получить коэффициенты обычного степенного разложения типа ax5+bx4+cx3+dx2+ex+f ?

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

Несколько раз бросалась искать ту рекомендацию в Google, но из-за обилия паразитных ссылок так и не нашла. А в этой теме, похоже, как раз намекается на тот метод. Подскажите, где искать про это, если кто помнит!

P.S. Среди студенческих рефератов (большего я так и не нашла) бродит такой абзац:

Цитата
Преобразование коэффициентов полинома Чебышева в коэффициенты традиционного многочлена
Вводим коэффициенты a0, a1, …, an многочлена T(x) и образуем массив ai
Для j = 2, 3, …, n и k = n, n-1, …, j в первом случае поднимаясь, а во втором спускаясь, осуществляем преобразование коэффициентов по следующим формулам:
а) ak-1 = ak-2 - ak
б) ak = 2ak
В результате получаем коэффициенты полинома Pn(x)
http://100balov.com/data/rus/Drygoe/26/POL...ChEBIrShEVA.rtf


Однако ж ничего из этого объяснения не поняла sm.gif, а большего там и не было.
Go to the top of the page
 
+Quote Post
iiv
сообщение Nov 29 2011, 17:33
Сообщение #3


вопрошающий
*****

Группа: Свой
Сообщений: 1 726
Регистрация: 24-01-11
Пользователь №: 62 436



Цитата(Xenia @ Nov 29 2011, 20:47) *
Когда-то даже видела рекомендацию, не считайте, мол, степенные коэффициенты в лоб, т.к. при решении нормальных уравнений матрица окажется плохо обусловлена. А считайте сперва полиномы Чебышева, а степенные коэффициенты получите из них.

У Вас координаты сетки интерполяции на чебышевских точках, или на равномерных заданы? Если на равномерных, то при конвертации из Чебышевской в степенное представление, Вы и получите опять эту же ошибку. Если у Вас сетка на Чебышевских точках, то матрица интерполяционных коэффициентов будет хорошо обусловлена.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- net   метод наименьших квадратов   Apr 19 2009, 13:15
- - Oldring   Цитата(net @ Apr 19 2009, 17:15) для реше...   Apr 20 2009, 09:04
|- - net   Цитата(Oldring @ Apr 20 2009, 13:04) . Ка...   Apr 20 2009, 11:39
|- - Oldring   Цитата(net @ Apr 20 2009, 15:39) все дело...   Apr 20 2009, 12:16
|- - net   Цитата(Oldring @ Apr 20 2009, 16:16) Все ...   Apr 20 2009, 15:16
|- - Oldring   Цитата(net @ Apr 20 2009, 19:16) но еще и...   Apr 20 2009, 16:20
|- - net   Цитата(Oldring @ Apr 20 2009, 20:20) Орто...   Apr 21 2009, 05:17
|- - Oldring   Цитата(net @ Apr 21 2009, 09:17) вы по де...   Apr 21 2009, 08:28
- - 314   В подобных ситуациях иногда можно выкрутится с пом...   May 3 2009, 16:10
|- - net   Цитата(314 @ May 3 2009, 20:10) В подобны...   May 12 2009, 10:31
- - 314   Прошу прощения, если не очень понятно выразился. П...   May 12 2009, 19:32
|- - net   Цитата(314 @ May 12 2009, 23:32) Предлага...   May 13 2009, 08:15
- - Xenia   Цитата(iiv @ Nov 29 2011, 21:33) У Вас ко...   Nov 29 2011, 19:52
- - iiv   Цитата(Xenia @ Nov 30 2011, 01:52) А я со...   Nov 29 2011, 20:42
- - _Pasha   Цитата(iiv @ Nov 30 2011, 00:42) Коэффици...   Nov 29 2011, 22:33
- - iiv   Цитата(_Pasha @ Nov 30 2011, 03:33) Для п...   Nov 30 2011, 00:36


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 19th July 2025 - 23:39
Рейтинг@Mail.ru


Страница сгенерированна за 0.01373 секунд с 7
ELECTRONIX ©2004-2016