IPB
ЛогинПароль:

> Внимание!

1. Пользуйтесь тегами кода. - [code] ... [/code]
2. Точно указывайте язык, название и версию компилятора (интерпретатора).
3. Название темы должно быть информативным. В описании темы указываем язык!!!

> Вычисление относительной погрешности для функции.
Krjuger
сообщение 12.04.2011 17:19
Сообщение #1


Профи
****

Группа: Пользователи
Сообщений: 652
Пол: Мужской
Реальное имя: Алексей

Репутация: -  20  +


В общем это некое продолжение моей прошлой темы,только обрастающее новыми подробностями.
У меня дана функция x(expx-1).Я эту функцию раскладываю я ряд Тейлора и получаю сумма от 1 до N от xn+1/n!.
Чтобы найти относительную погрешность мне надо а -n ый член разделить на сумму ряда,все это по модулю.Ну в общем то ,я думаю, вы и так понимаете как это делать.Суть заключается в том,у меня опять есть ограниченная разрядность мантисы и мне надо посмотреть как она будет влиять на результат.И как будет меняться N при которых мы будем выходить за граници возможностей нашей машины.

Я немного абстрагировался от этой задачи.что я сделал, при помощи маткада я посчитал,при каком N будет достигаться относительная погрешность на грани машинного эпсилон 10-16,это N =24.Да забыл сказать,что точка в которой мы раскладываем ряд у меня дана.Это -2.3.Затем я высчитал значение этой погрешности она составила 8.602898672363349*10^-16.Дальше я посомтрел как оно себя будет вести при разрядах мантисы от 10 до 25,при фиксированных исходных данных.Для рязрядности с 10 до 16 я получил ожидаемый результат,но дальше начало твориться что то вообще непонятное.Ну или точнее я не могу понять, как это интерпретировать.
ИзображениеФайл с самой программой тоже прикрепляю.


Прикрепленные файлы
Прикрепленный файл  ChMet_Laba_1.cpp ( 1.24 килобайт ) Кол-во скачиваний: 250
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
 
 Ответить  Открыть новую тему 
Ответов
Krjuger
сообщение 12.04.2011 19:56
Сообщение #2


Профи
****

Группа: Пользователи
Сообщений: 652
Пол: Мужской
Реальное имя: Алексей

Репутация: -  20  +


Я не много изменил основную программу,она стала выглядить так

void main()
{
double num=0;
double res=0;
double x1,x2;
int m=0;
int N;
x1=-2.3;
x2=-17;
for (int i=5;i<=20;i++)
{
cout << "Razryadnost* manticy:" << i << endl;
N=0;
for (int s=1;s<25;s++)
{
res=DeltaPogr(x1,s,i);
if (res!=0)
{
cout << setprecision (10) << res << endl;
N++;
}
}
cout << "Kol-vo chlenov ryada:" << N << endl;
}
}


Чтобы стало немного нагляднее.Принципи,большинство цифр получились вполне ожидаемыми,но возникли вопросы,почему,при разрядности мантисы 19 получается самое точное число ,а именно 8.603*10^-16,но дальше начинается уже какой то бред,так же непонятно,потому при 19,по идеи это должно было произойти при 16,либо я что то сильно недопонимаю.И я не понял почему при 20 я не получил на 1 значащую цифру более точный результат?В общем вопросов много....
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 

Сообщений в этой теме


 Ответить  Открыть новую тему 
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 



- Текстовая версия 28.07.2025 14:35
Хостинг предоставлен компанией "Веб Сервис Центр" при поддержке компании "ДокЛаб"