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

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [Mathematica] Выборка корней системы уравнений, Как выбрать из 8ми корней один нужный?! 
:(
    Опции темы
factorina
Дата 23.12.2011, 23:34 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет
Помогите, пожалуйста! 
Дана система уравнений, я её решаю, получаю 8 корней. Мне нужен один! 
А конкретнее этот: 
{x -> ArcCos[1/Sqrt[1 + t^2]], V -> 1/(1 + t^2)^(1/4)}}
Как сделать так, чтобы программа выдавала только его??? Какие функции использовать,  что делать? Я совершенно не представляю, а зачет в понедельник...
f[
Код

in:
f1[x_,V_]:=-Sin[x] - t*V^2
f2[x_,V_]:=(V^2 - Cos[x])/V
m=Solve[{f1[x,V]==0,f2[x,V]==0},{x,V}]

out:
{{x -> -ArcCos[1/Sqrt[1 + t^2]], V -> -(1/(1 + t^2)^(1/4))}, 
  {x -> ArcCos[1/Sqrt[1 + t^2]], V -> -(1/(1 + t^2)^(1/4))}, 
  {x -> -ArcCos[-(1/Sqrt[1 + t^2])], V -> -(I/(1 + t^2)^(1/4))},
  {x -> ArcCos[-(1/Sqrt[1 + t^2])], V -> -(I/(1 + t^2)^(1/4))},
  {x -> -ArcCos[-(1/Sqrt[1 + t^2])], V -> I/(1 + t^2)^(1/4)},
  {x -> ArcCos[-(1/Sqrt[1 + t^2])], V -> I/(1 + t^2)^(1/4)},
  {x -> -ArcCos[1/Sqrt[1 + t^2]], V -> 1/(1 + t^2)^(1/4)},
  {x -> ArcCos[1/Sqrt[1 + t^2]], V -> 1/(1 + t^2)^(1/4)}}
PM MAIL   Вверх
Фантом
Дата 24.12.2011, 00:06 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


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

Репутация: 7
Всего: 49
Что значит "мне нужен один"?  smile 

Это, собственно, намек. Вам надо сформулировать какой-то критерий, который отличает этот корень от остальных, и добавить этот критерий (надо полагать, это будет какое-то неравенство или несколько неравенств) в список уравнений решаемой системы.
PM   Вверх
factorina
Дата 24.12.2011, 21:25 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет
Да, я понимаю. Но  функция Solve у меня не решает систему уравнений и неравенств. А вообще должна? Я пишу просто к примеру: 
Код

Solve[{f1[x,V]==0,f2[x,V]==0,V>0},{x,V}]

Может по-другому нужно?

Решает Reduce, но здесь другая проблема: 
Мне этот корень потом нужно подставить в частную производную функции в точке m, я пишу:
Код

Derivative[1,0][f1[x,V]][x/.m,V/.m]

Для Solve считает правильно, для Reduce ругается на ReplaceAll:reps
PM MAIL   Вверх
Фантом
Дата 24.12.2011, 21:44 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


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

Репутация: 7
Всего: 49
Цитата(factorina @  24.12.2011,  22:25 Найти цитируемый пост)
Но  функция Solve у меня не решает систему уравнений и неравенств. А вообще должна?

Вообще-то она для этого и предназначена.  smile 
PM   Вверх
factorina
Дата 25.12.2011, 13:37 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет
Тогда что я делаю не так? Ваши ответы пользу приносят нулевую)
Можно пример?

Я пишу:
Код

Solve[{f1[x,V]==0,f2[x,V]==0,V>0},{x,V}]

А в ответ:
Код

Solve::nsmet: This system cannot be solved with the methods available to Solve

м?
PM MAIL   Вверх
Фантом
Дата 25.12.2011, 14:42 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


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

Репутация: 7
Всего: 49
Цитата(factorina @  25.12.2011,  14:37 Найти цитируемый пост)
Тогда что я делаю не так?

По-видимому, не осознаете математический смысл (буквально математический) того, что пытаетесь скормить Mathematica.

Цитата(factorina @  25.12.2011,  14:37 Найти цитируемый пост)
Я пишу:
1: Solve[{f1[x,V]==0,f2[x,V]==0,V>0},{x,V}]
А в ответ:
1: Solve::nsmet: This system cannot be solved with the methods available to Solve

Естественно. Посмотрите внимательнее на получавшиеся изначально корни. Про какой из них можно точно сказать, что V>0, при условии, что t неизвестно? Ни про какой. Сформулированное Вами дополнительное условие недостаточно.
PM   Вверх
factorina
  Дата 25.12.2011, 16:35 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет
Ладно другой вопрос. Как подставить полученный Reduce-ом корень в выражение для производной. Мне нужно найти производную в полученной точке.
Код

a=D[f1[x,V],V]
PM MAIL   Вверх
  
 Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman
 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Центр помощи | Следующая тема »

Подписаться на тему | Подписка на этот форум | Скачать/Распечатать тему


 




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


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

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