Модераторы: Daevaorn
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> Помогите с алгоритмом, реш. нелин.ур-я sin(x)=0методом итераций 
:(
    Опции темы
Georgiii
Дата 21.12.2005, 02:23 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 21
Регистрация: 19.12.2005

Репутация: нет
Всего: нет



smile
Дано f(x)=0 (1)
Заменим (1) x=w(x) (2)
Выберем грубое,приближенное значение х0,принадлежащее [a,b],подставим его в правую часть ур-я (2),получим:
х1=w(x0) (3)
далее подставим х1 в правую часть ур-я (3) получим:
х2=w(x1) (4)
x3=w(x2) (5)
Проделаем данный процесс n раз получим х(n)=w(x(n-1))
Если эта посл-ть является сходящейся,т.е сущ-ет предел x*=lim(x(n)), то данный алгоритм позволяет опр-ть искомый корень. Вфражение (5) запишем как x*=w(x*) (6)
Выражение (6) явл-ся решением выражения (2),теперь необходимо рассмотреть в каких случаях посл-ть х1....х(n) явл-ся сходящейся. Условием сходимости явл-ся если во всех точках х принадл.[a,b] выполняется условие: |w'(x)|<1 (производная <1)
Вопрос в том что не получается написать прогу решения sin(x)=0 на С++ т.к непонятен алгоритм,непонятно какое х выводить.
Если кто может напишите пожалуйста или помогите понять суть метода....лучше конечно первое smile
PM MAIL   Вверх
Georgiii
Дата 21.12.2005, 02:56 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 21
Регистрация: 19.12.2005

Репутация: нет
Всего: нет



smile smile
PM MAIL   Вверх
nworm
Дата 21.12.2005, 23:32 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 502
Регистрация: 22.10.2005

Репутация: 1
Всего: 8



А чего неясного? Надо ввести с клавиатуры x[0] и присваивать в цикле x[j+1]=sin(x[j]) - x[j] до момента |x[j] - x[j-1]|/|x[j]| < погрешность.
PM MAIL WWW   Вверх
Georgiii
Дата 22.12.2005, 00:34 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 21
Регистрация: 19.12.2005

Репутация: нет
Всего: нет



А как же с производной?
PM MAIL   Вверх
nworm
Дата 22.12.2005, 19:36 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 502
Регистрация: 22.10.2005

Репутация: 1
Всего: 8



1. Опечатался x[j+1]=x[j] - sin(x[j]).
2. А что? Требуется исследование? Если надо просто метод реализовать, делай примерно как я предлагал.
Можно поисследовать сходимость (для этого производная нужна), например, воспользовавшись теоремой
теорема

Это сообщение отредактировал(а) nworm - 22.12.2005, 19:37
PM MAIL WWW   Вверх
Partizan
Дата 23.12.2005, 02:51 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Let's do some .NET
****


Профиль
Группа: Модератор
Сообщений: 2828
Регистрация: 19.12.2005
Где: Санкт-Петербург

Репутация: 4
Всего: 67



да не парься...возьми приближённое значение производной в точке...возьми какое-то Х....и ООООО чень мальенькое dX...и вычисли по формуле y'=(F(X+dX)-F(X))/dX ...


--------------------
СУВ,
       Partizan.
PM MAIL WWW ICQ Skype GTalk Jabber   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
Правила форума "С++:Общие вопросы"
Earnest Daevaorn

Добро пожаловать!

  • Черновик стандарта C++ (за октябрь 2005) можно скачать с этого сайта. Прямая ссылка на файл черновика(4.4мб).
  • Черновик стандарта C (за сентябрь 2005) можно скачать с этого сайта. Прямая ссылка на файл черновика (3.4мб).
  • Прежде чем задать вопрос, прочтите это и/или это!
  • Здесь хранится весь мировой запас ссылок на документы, связанные с C++ :)
  • Не брезгуйте пользоваться тегами [code=cpp][/code].
  • Пожалуйста, не просите написать за вас программы в этом разделе - для этого существует "Центр Помощи".
  • C++ FAQ

Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, Earnest Daevaorn

 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | C/C++: Общие вопросы | Следующая тема »


 




[ Время генерации скрипта: 0.1489 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.