Программа 17. Задача 26
> | restart;read "C:\\Podpr1.m"; with(plots): |
Угловая скорость
> | omega[1]:=4: epsilon[1]:=0: |
Координаты
> | x[1]:=0: y[1]:=0: |
> | x[2]:=4: y[2]:=3: |
> | x[3]:=4: y[3]:=-1: |
> | x[4]:=8: y[4]:=-1: |
Уравнения 3х угловых скоростей
> | eq1:=sum(omega[i]*(x[i]-x[i+1]),i=1..3): |
> | eq2:=sum(omega[i]*(y[i]-y[i+1]),i=1..3): |
> | R:=solve({eq1,eq2},{omega[2],omega[3]}):assign(R): |
> | omega2=omega[2],omega3=omega[3]; |
Уравнения 3х угловых ускорений
> | eq3:=sum(epsilon[i]*(x[i]-x[i+1]),i=1..3)-sum(omega[i]^2*(y[i]-y[i+1]),i=1..3): |
> | eq4:=sum(epsilon[i]*(y[i]-y[i+1]),i=1..3)+sum(omega[i]^2*(x[i]-x[i+1]),i=1..3): |
> | solve({eq3,eq4},{epsilon[2],epsilon[3]}); |
> | display(ЛиниЯ(1,2),ЛиниЯ(2,3),ЛиниЯ(3,4),Опора(1,0.1),Опора(4,0.1),axes=NONE,scaling=constrained, title=`Четырехзвенник`); |
> |