Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ Сообщений: 2 ] 
Автор Сообщение
 Заголовок сообщения: Система уравнений
СообщениеДобавлено: 24 май 2021, 07:43 
Не в сети
Продвинутый
Зарегистрирован:
24 дек 2019, 00:16
Сообщений: 58
Cпасибо сказано: 18
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Доброго времени суток) Нужна помощь в решении системы.
restart;
with(plots):
tau1:=mu:
tau2:=mu/5:
H:=beta0:
mu:=1:
r0:=1:
#Q(lambda,psi,beta0)
R2_Q:=-(2*lambda)/(beta0^2*psi^2)*ln(1-(beta0*psi*r)/(lambda))-(2*r)/(beta0*psi):
dRdr2_Q:=r^2/(R2_Q*(lambda-psi*r*H)^2):
I1_Q:=dRdr2_Q+R2_Q/r^2+H^2+psi^2*R2_Q+lambda^2:
I2_Q:=(lambda-psi*r*H)^2*R2_Q/r^2+dRdr2_Q*(R2_Q/r^2+H^2+psi^2*R2_Q+lambda^2):
Q:=unapply(Pi*(int((tau1*diff(I1_Q,lambda)+tau2*diff(I2_Q,lambda))*r, r=0..r0)), lambda, psi, beta0):
#M(lambda,psi,beta0)
R2_M:=-(2*lambda)/(beta0^2*psi^2)*ln(1-(beta0*psi*r)/(lambda))-(2*r)/(beta0*psi):
dRdr2_M:=r^2/(R2_M*(lambda-psi*r*H)^2):
I1_M:=dRdr2_M+R2_M/r^2+H^2+psi^2*R2_M+lambda^2:
I2_M:=(lambda-psi*r*H)^2*R2_M/r^2+dRdr2_M*(R2_M/r^2+H^2+psi^2*R2_M+lambda^2):
M:=unapply(Pi*(int((tau1*diff(I1_M,psi)+tau2*diff(I2_M,psi))*r, r=0..r0)), lambda, psi, beta0):

В виде:
solve({Q(lambda, psi, beta0)=0, M(lambda, psi, beta0)=0},{lambda, psi});

при любом значении beta0, кроме нуля
например:
solve({Q(lambda, psi, 0.5)=0, M(lambda, psi, 0.5)=0},{lambda, psi});

получаю сообщение:
Warning, solutions may have been lost

Как иначе можно решить систему или точно узнать, что она не совместна?

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Численное решение с fsolve
СообщениеДобавлено: 27 май 2021, 22:33 
Не в сети
Продвинутый
Зарегистрирован:
14 фев 2020, 10:46
Сообщений: 92
Cпасибо сказано: 12
Спасибо получено:
42 раз в 38 сообщениях
Очков репутации: 26

Добавить очки репутацииУменьшить очки репутации
Я уже писал ранее, что подобные интегралы от достаточно сложных трансцендентных функций могут быть вычислены только численно. Соответственно, систему следует решать, используя fsolve вместо solve :
restart;
with(plots):
tau1:=mu:
tau2:=mu/5:
H:=beta0:
mu:=1:
r0:=1:
#Q(lambda,psi,beta0)
R2_Q:=-(2*lambda)/(beta0^2*psi^2)*ln(1-(beta0*psi*r)/(lambda))-(2*r)/(beta0*psi):
dRdr2_Q:=r^2/(R2_Q*(lambda-psi*r*H)^2):
I1_Q:=dRdr2_Q+R2_Q/r^2+H^2+psi^2*R2_Q+lambda^2:
I2_Q:=(lambda-psi*r*H)^2*R2_Q/r^2+dRdr2_Q*(R2_Q/r^2+H^2+psi^2*R2_Q+lambda^2):
Q:=(lambda, psi, beta0)->Pi*(Int((tau1*diff(I1_Q,lambda)+tau2*diff(I2_Q,lambda))*r, r=0..r0)):
R2_M:=-(2*lambda)/(beta0^2*psi^2)*ln(1-(beta0*psi*r)/(lambda))-(2*r)/(beta0*psi):
dRdr2_M:=r^2/(R2_M*(lambda-psi*r*H)^2):
I1_M:=dRdr2_M+R2_M/r^2+H^2+psi^2*R2_M+lambda^2:
I2_M:=(lambda-psi*r*H)^2*R2_M/r^2+dRdr2_M*(R2_M/r^2+H^2+psi^2*R2_M+lambda^2):
M:=(lambda, psi, beta0)->Pi*(Int((tau1*diff(I1_M,psi)+tau2*diff(I2_M,psi))*r, r=0..r0)):
fsolve({evalf(eval(Q(lambda,psi, beta0),beta0=0.5)),evalf(eval(M(lambda,psi, beta0),beta0=0.5))},{lambda=-5..5, psi=-5..5});

{lambda = 0.9096005059, psi = -0.5877529067}

Вернуться к началу
 Профиль  
Cпасибо сказано 
За это сообщение пользователю Kitonum "Спасибо" сказали:
Susanna Gaybaryan
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему      Страница 1 из 1 [ Сообщений: 2 ]

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Система уравнений

в форуме Алгебра

[Alexa]

7

263

09 янв 2022, 19:15

Система уравнений

в форуме Алгебра

Tenken

13

1263

02 авг 2016, 21:40

Система уравнений

в форуме Алгебра

vitlik2409

9

763

08 окт 2014, 22:13

Система уравнений

в форуме Алгебра

ivashenko

20

886

07 май 2016, 00:00

Система уравнений

в форуме Алгебра

neeara

10

485

08 июн 2018, 08:06

Система уравнений

в форуме Алгебра

uiiiiiii

5

259

22 апр 2020, 17:21

Система уравнений

в форуме Алгебра

DeD

4

638

24 авг 2016, 22:05

Система уравнений

в форуме Алгебра

qwer

1

233

16 янв 2016, 21:52

Система уравнений

в форуме Алгебра

[Alexa]

5

274

29 дек 2021, 20:33

Система уравнений

в форуме Тригонометрия

Musk

16

825

23 апр 2018, 13:15


Часовой пояс: UTC + 3 часа [ Летнее время ]



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 6


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved