Matlab课程设计
李俊机自091
设计题目一：凸轮机构设计
已知轮廓为圆形的凸轮（圆的半径为100mm、偏心距为20mm），推杆与凸轮运动中心的距离20mm，滚子半径为10mm，请利用matlab仿真出凸轮推杆的运动轨迹和运动特性（速度，加速度），并利用动画演示出相关轨迹和运动特性。
总程序代码
clcclfclearpfigurepositio
1001001200600fori1360
画圆形凸轮
R100
圆形凸轮半径
A000062pi
Bipi180
e20
偏心距
aecosBbesi
BxRcosAayRsi
Absubplot121plotxybLi
eWidth3
填充fillxyyaxisReReReRe100setgcaXlimReResetgcaYlimReRe100axisequalaxisma
ualaxisoffholdo
plotabogplote0orplot00orLi
eWidth3
画滚子gcx0r10
滚子中心X坐标滚子半径
fgcysqrtRr2a2b滚子中心Y坐标
gxrcosAgcx
滚子X坐标
gyrsi
Agcy
滚子Y坐标
plotgxgybLi
eWidth2
画其它部分plot0a0bkLi
eWidth4plot33170190mLi
eWidth4plot33170190mLi
eWidth4画顶杆
gc120dgx00dgygcygcygcplotdgxdgyLi
eWidth4holdoff
画位移图
sxiBsyigcysubplot322plotsxsybLi
eWidth3title位移线图
grido
holdoff
画速度图
vxiBvyi20cosB40cosBsi
B1214cosB212subplot324plotvxvygLi
eWidth3title速度线图
grido
holdoff
画加速度图
axiBayi40cosB21214cosB21220si
B40si
B21214cosB212160cosB2si
B21214cosB232subplot326plotaxayrLi
eWidth3xlabelBtitle加速度线图
grido
holdoff
Mgetframe
fe
d
截图
附：通过求导求速度和加速度
求速度symsBaecosBbesi
BssqrtRr2a2bvdiffs结果：v20cosB40cosBsi
B1214cosB212求加速度symsBv20cosB40cosBsi
B1214cosB212adiffv结果：a40cosB21214cosB21220si
B40si
B21214cosB212160cosB2si
B21214cosB232
ffr