PLOT(CURVES([[x1,y],seq([x1+l+j*d,r*sin(Pi*j/2)+y],j= 0..K),[x2,y]])): end: >
for t to N do m:=30*evalf(2+sin(4*Pi*t/N));H[t]:=pru(0,m,0):od: >
for t to N do m1:=20*evalf(2+sin(2*Pi*t/N));H2[t]:=pru(0,m1,15):od: >
a:=display(seq(H[t],t=1..N),insequence=true): >
b:=display(seq(H2[t],t=1..N),insequence=true): >
display(a,b,axes=none,scaling=constrained); 3.22-rasm
Erkinlik darajasi ikkiga teng sistemaning tebranishi Ikkita prujinaga maxkamlangan jismlarning tebranma harakatini koʻramiz A
yuk chap tomondan prujina orqali qoʻzg’almas devorga hamda oʻng tomondan
prujina orqali B yukka biriktirilgan. Bu yuklar sistemasi ikkita erkinlik darajasiga
ega boʻladi. Bu sistemani harakatini oʻrganish uchun quyidagi dasturni MAPLE
da yozamiz.
> restart; CHastotalar nisbati 3 ga teng boʻlishi uchun yuklarning massalari (kg)
quyidagicha tanlangan
> ma:=24.63:mb:=40: Prujinaning bikrligi (N/m)
312
> c:=90: A, B jismlarning kinetik energiya
Ta:=ma*vx^2/2: > Tb:=mb*(vs+vx)^2/2: Sistemaning kinetik energiyasi
> T:=Ta+Tb: Inertsiya koeffitsenti
> a11:=coeff(diff(T,vx),vx): > a12:=coeff(diff(T,vx),vs): > a22:=coeff(diff(T,vs),vs): Tebranish chastotasi omg
> B1:=c-a11*omg^2: B2:=a12*omg^2: > CHastot tenglamasi:=B1*(c-a22*omg^2)-B2^2=0: CHastota tenglamasini yechish
