Fizika pomosh Wirbelforce

[svoy]

Сос жоговурд горю!!!!

Проект на шее помошь нужна по физике!
Задание, надо описать маятник Вилберфорса (english pendulum Wilberforce, deutsch Das Pendel von Wilberforce)

Тут даже фильм есть с маятником
http://www.ap.stmarys.ca/demos/conte..._pendulum.html

Кто чем может помогите!
Кое что уже сделанно

Это началные требования!


Резулатат системы дифф уравнения


Дифференциальные уравнения движения уже решил с MAPLE и могу плот делать
На графике 2 функции одна функция (зеленая) это движение по вертикали(лонгитудиальное "longitudinal") другая (красная) вокруг своей оси (торсиональное "torsional").



А на этом графике отношение двух движений но чтото у меня не сходится при изменении т (времени ) етот график становится соовсем странным.
тут t=0..20


а тут t=0..200



t=0..2000


Или это нормально?

да до сих пор я принял что всё пропорционально то есть омега з = омега тхета
А, что произойдет если они не равны??
И еше надо мне описать движение при пропорциональном затухании колебаний.

Полный код с графикамин здесь
http://www.grigz.com/kmp/wirbelforce.html

Заранее благодарен!
Gev

[Ektich]

Не знаю как на счет "bug in Maple", то что ты называешь "численное решение" никакого физического смысла не имеет, в этом ты прав...
А что за численный метод? Я бы в первую очередь постарался бы понять чего конкретно он считает, и почему не сходится...

Кстати, что такое этот start_bed?

[svoy]

Chislenniy method eto protivopolozhnoe k analisis ili symbolic method!

start_bed eto list nachalniy usloviy. NAprimer z(0)= 0.5 Wilberforce podnyali na 0.5 ili theta (0=) = -0.5 povernuli na 0.5 a ostal'nie 2 eto nachal'naya skorost' etoi dannie nuzhni dlya resheniya sistemi Diff uravneniy.

EKTICH mozhesh' pozhaluysta poprobovat' vischitat' eti uravneniya smosh'yu tvoe progi lili kto nibud' drugo nuzhno eti uravneniya soschitat' lyuboy programmoy krome maple9.

Glavnoe pritom ispol'zovat' dsolve no numeric a ne symbolic. Spasibo zaranee!

[Ektich]

Svoy, ya znaju chto takoe "chislennij metod", ya sprashival kakoj konkretno metod Maple ispolzuet.

"svoej" programmki u menya net, to chem ya grafiki stroil eto tolko grafopostroitel. No esli poluchitsa ya mozhet napishu programku na FORTRAN-e chtobi chislenno re**** eti uravneniya. Xotya ne obeshaju, u menya v poslednee vremya napryag so svobodnim vremenem.

Ektich

[svoy]

Originally Posted by Ektich

"svoej" programmki u menya net, to chem ya grafiki stroil eto tolko grafopostroitel. No esli poluchitsa ya mozhet napishu programku na FORTRAN-e chtobi chislenno re**** eti uravneniya. Xotya ne obeshaju, u menya v poslednee vremya napryag so svobodnim vremenem.
Ektich

Vot vot Ektich Jan mne i nuzhno prosto postroit' plot etix funkciy no s nachal'nimi usloviyami z ( 0 ) = - theta ( 0 ) No prosto ya neznayu kak tvoya grafik proga vsyo eto ocenivayet....

Merci zaranee!

[svoy]

Думал, думал и наконец-то придумал numpoints=600 и всё проблем больше не будет!

большое спасибо всем за внимание и за коментарии осбенно Ektich!!!
Без тебя Ектич я бы это сам не на"плотил" бы !!!! Биг танкс!
последняя версия тут solved

[svoy]

Привет всем сострадальцам этого проэкта!

В пятницу сдали то-есть сделали презантацию получили оценку 5.5 (из 6) это по сов. стандартам 5- радости полные штаны!

Спасибо всем за хелп! Особенно Ektich! А вот и пикча сокурсник снял моей мобилой