基于MATLAB曲柄滑块机构运动仿真报告..doc

???????????????????????? 计算机仿真技术matlab 报告 ???????????????????????? 曲柄滑块机构 目录 一、基于GUI的曲柄滑块机构运动仿真 二、基于simulink的曲柄滑块机构运动仿真 曲柄滑块机构 1.题目描述 题目：对如图1所示的曲柄滑块机构的运动过程进行仿真，并用动画的方式显示曲柄滑块机构的运动过程，位移曲线、速度曲线和加速度曲线。图中，AB长R2，BC长R3mm，A点为坐标原点。 图1 曲柄滑块机构示意图 2.实现方法 利用GUI界面设计来对曲柄滑块机构的运动过程进行仿真，并用动画的方式显示曲柄滑块机构的运动过程。 3.界面设计 1. Gui 设计 1）新建GUI：菜单-新建-gui，并保存为test5 2）界面设计：拖拽左侧图标到绘图区，创建GUI界面 拖拽左侧图标值绘图区 设置如下的按钮 最终的仿真界面如图所示 3）代码添加： 进入代码界面 4.代码编程 %模型求解 a1=str2double(get(handles.edit1,String)); a2=str2double(get(handles.edit2,String)); a3=str2double(get(handles.edit3,String)); a4=str2double(get(handles.edit4,String)); a5=str2double(get(handles.edit5,String)); a=a1*((1-cos(a4*a5))+0.25*(a1/a2)*(1-cos(2*a4*a5))); set(handles.edit6,String,a); a0=(a4*a1)*(sin(a4*a5)+0.5*(a1/a2)*sin(2*a4*a5)); set(handles.edit7,String,a0); a6=(a4*a4*a1)*(cos(a4*a5)+(a1/a2)*cos(a4*a5)); set(handles.edit8,String,a6); %绘制位移、速度、加速度曲线 axes(handles.axes3); r1=str2double(get(handles.edit1,String)); r2=str2double(get(handles.edit2,String)); omiga1=str2double(get(handles.edit4,String)); x11=1:720; for i=1:720 x1(i)=i*pi/180; %sin(x2(i)=r1/r2*sin(x1(i)); x2(i)=asin(-r1/r2*sin(x1(i))); x22(i)=x2(i)*180/pi; r3(i)=r1*cos(x1(i))+r2*cos(x2(i)); B=[-r1*omiga1*sin(x1(i));r1*omiga1*cos(x1(i))]; A=[r2*sin(x2(i)) 1;-r2*cos(x2(i)) 0]; X=inv(A)*B; omiga2(i)=X(1,1); v3(i)=X(2,1); end plot(x11/60,0.5*r1*sin(x1)); xlabel(ê±???á t/sec) ylabel(á????êD??úY?áé?μ?????/mm) figure(2) plot(x11/60,r3); title(??ò???í?) grid on hold off; xlabel(ê±??t/sec) ylabel(???é??ò? r3/mm) figure(3) plot(x11/60,omiga2); title(á??????ù?è) grid on hold off; xlabel(ê±?? t/sec) ylabel(á??????ù?è omiga2/rad/sec) figure(4) plot(x11/60,v3*pi/180); title(???é?ù?è) grid on hold off; xlabel(ê±?? t/sec) ylabel(???é?ù?è v3/mm/sec) %绘制表格 axes(handles.axes3); grid on axes(handles.axes1); grid on %制作动画 axes(handles.axes1); hf=figure(name,?ú±ú