2014 dxdy logo

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

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




 
 RootOf, allvalues не срабатывает
Сообщение30.06.2011, 14:56 
вот кусочек программы
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 
В версиях, начиная с 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 
спасибо вам огромное . громоздкость не имеет значения

 
 
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 13:47 
нет, к сожалению, этот вариант не подходит. почему то не учитываются комплексные корни

 
 
 
 Re: RootOf, allvalues не срабатывает
Сообщение01.07.2011, 16:35 
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 
никогда бы сама не додумалась. ансис в этом плане проще. спасибо вам еще раз. а что в вашем расчете шестерка означает, та что в самом конце, после лямбды?

 
 
 
 Re: RootOf, allvalues не срабатывает
Сообщение02.07.2011, 00:58 
Unona в сообщении #464126 писал(а):
а что в вашем расчете шестерка означает, та что в самом конце, после лямбды?
Точность округления.

 
 
 
 Re: RootOf, allvalues не срабатывает
Сообщение02.07.2011, 22:47 
спасибо вам еще раз, вы очень выручили

 
 
 [ Сообщений: 8 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group