2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 RootOf, allvalues не срабатывает
Сообщение30.06.2011, 14:56 


30/06/11
5
вот кусочек программы
restart;
> Do:=E*h^3/12/(1-mu^2);
> E:= 2*10^11; mu:=0.3; ro:=7800; R:=0.15; h:=0.001; l:=0.4; l1:=0.07; R1:=0.03; teta1:=arctg(0.12/0.15); teta2:=arctg(0.18/0.15); x:=r*cos(teta); y:=r*sin(teta); ksi:=x/R; phi:=y/R;
> k:=lambda/n; LL:=(1-mu^2)/E*ro*R^2*omega^2; Omega:=LL*omega^2; beta:=h^2/12/R^2;
> solve(Omega^3-Omega^2*(1+(3-mu)/2*n^2*(1-k^2)+beta*n^4*(1-k^2)^2)+Omega*((3-mu)/2*n^2*(1-k^2)-n^2*(1-mu^2*k^2)+(1-mu)/2*n^4*(1-k^2)^3+(3-mu)/2*beta*n^6*(1-k^2)^3)-(1-mu)/2*(1-mu^2)*k^2*n^4-(1-mu)/2*beta*n^8*(1-k^2)^4=0,lambda);
выдает ответ в таком вот виде
RootOf(_Z^2-RootOf(3240740740000000000000000000000000000*n^2*_Z^4+
(698709375000000000000000000000000*omega^4+9981562498003687500000000000*omega^4*
n^2-12962962960000000000000000000000000000*n^4)*_Z^3+
(-2096128125000000000000000000000000*omega^4*n^2+
19444444440000000000000000000000000000*n^6-29944687494011062500000000000*
omega^4*n^4+5904094217396740617*omega^8*n^2)*_Z^2+
(-12962962960000000000000000000000000000*n^8+2515353750000000000000000000000000*
omega^4*n^2+2096128125000000000000000000000000*omega^4*n^4+
29944687494011062500000000000*omega^4*n^6+
796250000000000000000000000000000000000000*n^4-
11808188434793481234*omega^8*n^4)*_Z-1272933045750000*omega^12*n^2+
1594105439000000000000000*omega^8*n^2+2152042342650000000000000*omega^8*n^4+
5904094217396740617*omega^8*n^6-698709375000000000000000000000000*omega^4*n^4-
698709375000000000000000000000000*omega^4*n^6-9981562498003687500000000000*omega^4*
n^8+3240740740000000000000000000000000000*n^10))
далее пытаюсь с помощью команды allvalues(%) перевести в человеческий вид. пробывала считать в версиях
с 5 по 15, в одних вылетает (выдает сообщение, что мапл нашел у себя ошибку) или зависает с концами. бьюсь
над проблемой больше месяца

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 11:07 
Заслуженный участник


27/06/08
4063
Волгоград
В версиях, начиная с 6-й в "человеческий вид" приводят не командой allvallues, а присваивая переменной окружения _EnvExplicit значение true. Это позволяет представить корни уравнений в 4-й степени в радикалах.
Во в Вашем случае это не поможет. Ответ станет еще страшнее.Причем намного.
Можно посчитать округленно, но и это не избавит Вас от громоздкости. Ведь искомое lambda зависит от параметров n и omega.

У меня получилось:
Код:
[.138873e-8*(69282.*((.162732e18*((.410936e71*n^8+.177376e50*omega^12*n^2+.777490e58*omega^8*n^4-.573470e48*omega^8*n^6+.259630e63*omega^4*n^6+.155647e41*omega^16+.774090e48*omega^8*n^8-.216818e45*omega^12*n^4+.145326e40*omega^12*n^6+.375528e36*omega^16*n^2+.300177e36*omega^16*n^4-.392010e21*omega^20*n^4-.124288e32*omega^20-.295924e27*omega^20*n^2+304899.*n^4*omega^24+.100000e12*(.102046e103*omega^12*n^8+.168868e120*n^16+.104808e99*omega^12*n^10+.710972e106*omega^8*n^12+.651420e103*omega^8*n^14+.213382e112*omega^4*n^14+.105374e112*omega^4*n^16+.967085e94*omega^16*n^6+.604490e94*omega^16*n^8+.305502e86*omega^20*n^4+.978110e94*n^12*omega^12+.307008e90*n^10*omega^16+.193076e86*n^6*omega^20+.338975e76*n^4*omega^24+.242259e59*omega^32-.386900e50*omega^36-.651660e115*n^20+.154474e41*omega^40+.330060e72*omega^24*n^6+.518910e71*omega^24*n^8+.125733e63*omega^28*n^6+.265528e59*omega^32*n^2-.946350e45*omega^36*n^4-.201483e50*omega^36*n^2+.926445e31*omega^40*n^4+.188546e77*omega^20*n^10-.274802e67*omega^24*n^10-.194587e59*omega^28*n^8+.326775e50*omega^32*n^6-.184288e41*omega^36*n^6+.201075e26*omega^40*n^6-.392178e85*omega^16*n^12-.478455e71*omega^20*n^12+.901055e49*omega^32*n^8-.242662e35*omega^36*n^8-.284790e94*n^14*omega^12-.285122e80*n^14*omega^16+.281428e53*n^10*omega^28-.146034e39*n^10*omega^32-.216544e11*omega^44*n^6+.735592e36*omega^40*n^2-.757905e15*omega^44*n^4-.298870e98*n^16*omega^8-.295495e107*n^18*omega^4+.419228e62*n^12*omega^24+.285495e20*n^8*omega^40+.158461e66*n^14*omega^20+.511728e23*n^10*omega^36-.943175e88*n^16*omega^12+.379138e74*n^16*omega^16-.819250e47*n^12*omega^28+.224294e32*n^12*omega^32-35857.*n^8*omega^44+.321695e77*omega^24*n^2+.303580e67*omega^28*n^4-.273530e68*omega^28*n^2+.107009e55*omega^32*n^4+.967215e81*omega^20*n^8)^(1/2))*n^2)^(1/3)*omega^8+.464950e13*((.410936e71*n^8+.177376e50*omega^12*n^2+.777490e58*omega^8*n^4-.573470e48*omega^8*n^6+.259630e63*omega^4*n^6+.155647e41*omega^16+.774090e48*omega^8*n^8-.216818e45*omega^12*n^4+.145326e40*omega^12*n^6+.375528e36*omega^16*n^2+.300177e36*omega^16*n^4-.392010e21*omega^20*n^4-.124288e32*omega^20-.295924e27*omega^20*n^2+304899.*n^4*omega^24+.100000e12*(.102046e103*omega^12*n^8+.168868e120*n^16+.104808e99*omega^12*n^10+.710972e106*omega^8*n^12+.651420e103*omega^8*n^14+.213382e112*omega^4*n^14+.105374e112*omega^4*n^16+.967085e94*omega^16*n^6+.604490e94*omega^16*n^8+.305502e86*omega^20*n^4+.978110e94*n^12*omega^12+.307008e90*n^10*omega^16+.193076e86*n^6*omega^20+.338975e76*n^4*omega^24+.242259e59*omega^32-.386900e50*omega^36-.651660e115*n^20+.154474e41*omega^40+.330060e72*omega^24*n^6+.518910e71*omega^24*n^8+.125733e63*omega^28*n^6+.265528e59*omega^32*n^2-.946350e45*omega^36*n^4-.201483e50*omega^36*n^2+.926445e31*omega^40*n^4+.188546e77*omega^20*n^10-.274802e67*omega^24*n^10-.194587e59*omega^28*n^8+.326775e50*omega^32*n^6-.184288e41*omega^36*n^6+.201075e26*omega^40*n^6-.392178e85*omega^16*n^12-.478455e71*omega^20*n^12+.901055e49*omega^32*n^8-.242662e35*omega^36*n^8-.284790e94*n^14*omega^12-.285122e80*n^14*omega^16+.281428e53*n^10*omega^28-.146034e39*n^10*omega^32-.216544e11*omega^44*n^6+.735592e36*omega^40*n^2-.757905e15*omega^44*n^4-.298870e98*n^16*omega^8-.295495e107*n^18*omega^4+.419228e62*n^12*omega^24+.285495e20*n^8*omega^40+.158461e66*n^14*omega^20+.511728e23*n^10*omega^36-.943175e88*n^16*omega^12+.379138e74*n^16*omega^16-.819250e47*n^12*omega^28+.224294e32*n^12*omega^32-35857.*n^8*omega^44+.321695e77*omega^24*n^2+.303580e67*omega^28*n^4-.273530e68*omega^28*n^2+.107009e55*omega^32*n^4+.967215e81*omega^20*n^8)^(1/2))*n^2)^(1/3)*omega^8*n^2+.162028e8*omega^8*n^4*((.410936e71*n^8+.177376e50*omega^12*n^2+.777490e58*omega^8*n^4-.573470e48*omega^8*n^6+.259630e63*omega^4*n^6+.155647e41*omega^16+.774090e48*omega^8*n^8-.216818e45*omega^12*n^4+.145326e40*omega^12*n^6+.375528e36*omega^16*n^2+.300177e36*omega^16*n^4-.392010e21*omega^20*n^4-.124288e32*omega^20-.295924e27*omega^20*n^2+304899.*n^4*omega^24+.100000e12*(.102046e103*omega^12*n^8+.168868e120*n^16+.104808e99*omega^12*n^10+.710972e106*omega^8*n^12+.651420e103*omega^8*n^14+.213382e112*omega^4*n^14+.105374e112*omega^4*n^16+.967085e94*omega^16*n^6+.604490e94*omega^16*n^8+.305502e86*omega^20*n^4+.978110e94*n^12*omega^12+.307008e90*n^10*omega^16+.193076e86*n^6*omega^20+.338975e76*n^4*omega^24+.242259e59*omega^32-.386900e50*omega^36-.651660e115*n^20+.154474e41*omega^40+.330060e72*omega^24*n^6+.518910e71*omega^24*n^8+.125733e63*omega^28*n^6+.265528e59*omega^32*n^2-.946350e45*omega^36*n^4-.201483e50*omega^36*n^2+.926445e31*omega^40*n^4+.188546e77*omega^20*n^10-.274802e67*omega^24*n^10-.194587e59*omega^28*n^8+.326775e50*omega^32*n^6-.184288e41*omega^36*n^6+.201075e26*omega^40*n^6-.392178e85*omega^16*n^12-.478455e71*omega^20*n^12+.901055e49*omega^32*n^8-.242662e35*omega^36*n^8-.284790e94*n^14*omega^12-.285122e80*n^14*omega^16+.281428e53*n^10*omega^28-.146034e39*n^10*omega^32-.216544e11*omega^44*n^6+.735592e36*omega^40*n^2-.757905e15*omega^44*n^4-.298870e98*n^16*omega^8-.295495e107*n^18*omega^4+.419228e62*n^12*omega^24+.285495e20*n^8*omega^40+.158461e66*n^14*omega^20+.511728e23*n^10*omega^36-.943175e88*n^16*omega^12+.379138e74*n^16*omega^16-.819250e47*n^12*omega^28+.224294e32*n^12*omega^32-35857.*n^8*omega^44+.321695e77*omega^24*n^2+.303580e67*omega^28*n^4-.273530e68*omega^28*n^2+.107009e55*omega^32*n^4+.967215e81*omega^20*n^8)^(1/2))*n^2)^(1/3)+126347.*((.410936e71*n^8+.177376e50*omega^12*n^2+.777490e58*omega^8*n^4

Не уверен, что такой ответ Вас обрадует :)

А откуда уверенность, что ответ должен быть компактным?

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 11:32 


30/06/11
5
спасибо вам огромное . громоздкость не имеет значения

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 13:47 


30/06/11
5
нет, к сожалению, этот вариант не подходит. почему то не учитываются комплексные корни

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 16:35 
Заслуженный участник


27/06/08
4063
Волгоград
Unona в сообщении #463958 писал(а):
нет, к сожалению, этот вариант не подходит. почему то не учитываются комплексные корни
Почему не учитываются?
Например, при $n=1, \ -1\le \omega \le 1$ два из восьми значений $\lambda$ оказываются вещестеннымим, а остальные комплексны.

Есть подозрение, что в прошлый раз не все скопировал (больно коротенькое выражение получилось :-) Похоже, буфера обмена не хватает).
Поэтому просто приведу команды, которую я вводил:
Код:
Do:=E*h^3/12/(1-mu^2);
> E:= 2*10^11; mu:=0.3; ro:=7800; R:=0.15; h:=0.001; l:=0.4; l1:=0.07; R1:=0.03; teta1:=arctg(0.12/0.15); teta2:=arctg(0.18/0.15); x:=r*cos(teta); y:=r*sin(teta); ksi:=x/R; phi:=y/R;
> k:=lambda/n; LL:=(1-mu^2)/E*ro*R^2*omega^2; Omega:=LL*omega^2; beta:=h^2/12/R^2;
> _EnvExplicit:=true:
> s:=evalf(simplify([solve(Omega^3-Omega^2*(1+(3-mu)/2*n^2*(1-k^2)+beta*n^4*(1-k^2)^2)+Omega*((3-mu)/2*n^2*(1-k^2)-n^2*(1-mu^2*k^2)+(1-mu)/2*n^4*(1-k^2)^3+(3-mu)/2*beta*n^6*(1-k^2)^3)-(1-mu)/2*(1-mu^2)*k^2*n^4-(1-mu)/2*beta*n^8*(1-k^2)^4=0,lambda)]),6);

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 23:12 


30/06/11
5
никогда бы сама не додумалась. ансис в этом плане проще. спасибо вам еще раз. а что в вашем расчете шестерка означает, та что в самом конце, после лямбды?

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение02.07.2011, 00:58 
Заслуженный участник


27/06/08
4063
Волгоград
Unona в сообщении #464126 писал(а):
а что в вашем расчете шестерка означает, та что в самом конце, после лямбды?
Точность округления.

 Профиль  
                  
 
 Re: RootOf, allvalues не срабатывает
Сообщение02.07.2011, 22:47 


30/06/11
5
спасибо вам еще раз, вы очень выручили

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 8 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



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

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


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

Найти:
cron
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group