計算機仿真技術matlab報告曲柄滑塊機構目錄一、基于GUI的曲柄滑塊機構運動仿真二、基于simulink的曲柄滑塊機構運動仿真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;fori=1:720x1(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);endplot(x11/60,0.5*r1*sin(x1));xlabel('ê±???át/sec')ylabel('á????êD??úY?áé?μ?????/mm')figure(2)plot(x11/60,r3);title('??ò???í?')gridonholdoff;xlabel('ê±??t/sec')ylabel('???é??ò?r3/mm')figure(3)plot(x11/60,omiga2);title('á??????ù?è')gridonholdoff;xlabel('ê±??t/sec')ylabel('á??????ù?èomiga2/rad/sec')figure(4)plot(x11/60,v3*pi/180);title('???é?ù?è')gridonholdoff;xlabel('ê±??t/sec')ylabel('???é?ù?èv3/mm/sec')%繪制表格axes(handles.axes3);gridonaxes(handles.axes1);gridon%制作動畫axes(handles.axes1);hf=figure('name','?ú±ú???é?ú11');set(hf,'color','r');holdonaxis([-6,6,-4,4]);gridonaxis('off');xa0=-5;%??è?×ó?￥μ?×?±êxa1=-2.5;%??è?óò?￥μ?×?±êxb0=-2.5;%á???×ó?￥μ?×?±êxb1=2.2;%á???óò?￥μ?×?±êx3=3.5;%×a??×?±êy3=0;%×a??×?±êx4=xb1;%éè??á???í·μ?3?ê?????oá×?±êy4=0;%éè??á???í·μ?3?ê?????×Y×?±êx5=xa1;y5=0;x6=x3;%éè??á??á3?ê?oá×?±êy6=0;%éè??á??á3?ê?×Y×?±êa=0.7;b=0.7c=0.7a1=line([xa0;xa1],[0;0],'color','b','linestyle','-','linewidth',40);%éè????è?a3=line(x3,y3,'color',[0.50.60.3],'linestyle','.','markersize',300);%éè??×a??a2=line([xb0;xb1],[0;0],'color','black','linewidth',10);%éè??á???a5=line(x5,y5,'color','black','linestyle','.','markersize',40);%éè??á?????è?á??óí·a4=line(x4,y4,'color','black','linestyle','.','markersize',50);%éè??á???á??óí·a6=line([xb1;x3],[0;0],'color','black','linestyle','-','linewidth',10);a7=line(x3,0,'color','black','linestyle','.','markersize',50);%éè?????ˉ?DD?a8=line([-5.1;-0.2],[0.7;0.7],'color','y','linestyle','-','linewidth',5);%éè?????×±úa9=line([-5.1;-0.2],[-0.72;-0.72],'color','y','linestyle','-','linewidth',5);%éè?????×±úa10=line([-5.1;-5.1],[-0.8;0.75],'color','y','linestyle','-','linewidth',5);%éè?????×±úa11=fill([-5,-5,-5,-5],[0.61,0.61,-0.61,-0.61],[a,b,c]);%éè?????×??ì?len1=4.8;%á???3¤len2=2.5;%??è?3¤r=1.3;%???ˉ°???dt=0.015*pi;t=0;while1t=t+dt;ift>2*pit=0;endlena1=sqrt((len1)^2-(r*sin(t))^2);%á????ú???ˉ1y3ì?Doá?áé?μ?óDD§3¤?èrr1=r*cos(t);%°????ú???ˉ1y3ì?Doá?áé?μ?óDD§3¤?èxaa1=x3-sqrt(len1^2-(sin(t)*r)^2)-(r*cos(t));%??è??ú???ˉ1y3ì?Dμ?óò?￥μ?×?±ê????xaa0=xaa1-2.5;%%??è??ú???ˉ1y3ì?Dμ?×ó?￥μ?×?±ê????x55=x3-cos(t)*r;%á????ú???ˉ1y3ì?Doá×?±ê????y55=y3-sin(t)*r;%á????ú???ˉ1y3ì?D×Y×?±ê????set(a4,'xdata',x55,'ydata',y55);%éè??á????￥μ????ˉset(a1,'xdata',[xaa1-2.5;xaa1],'ydata',[0;0]);%éè????è????ˉset(a2,'xdata',[xaa1;x55],'ydata',[0;y55]);set(a5,'xdata',xaa1);%éè????è?ó?á???á??óí·μ????ˉset(a6,'xdata',[x55;x3],'ydata',[y55;0]);set(a11,'xdata',[-5,xaa0,xaa0,-5]);%éè????ì?μ?ì?3?set(gcf,'doublebuffer','on');%??3y?e?ˉdra