【正文】 *L2*cos(th2(i))]。w=inv(A)*B。w3(i)=w(1)。w4(i)=w(2)。endfor i=1:length(th2)C=[L3*sin(th34(i,1)) L4*sin(th34(i,2))。L3*cos(th34(i,1)) L4*cos(th34(i,2))]。D=[w2^2*L2*cos(th2(i))+w3(i)^2*L3*cos(th34(i,1))w4(i)^2*L4*cos(th34(i,2))。... w2^2*L2*sin(th2(i))+w3(i)^2*L3*sin(th34(i,1))w4(i)^2*L4*sin(th34(i,2))]。a=inv(C)*D。a3(i)=a(1)。a4(i)=a(2)。endval=get(hObject,39。Value39。)。str=get(hObject,39。String39。)。switch str{val}case 39。連桿3的幾個位置點39。th2=[0:1/6:2]*pi。 th34=zeros(length(th2),2)。 options=optimset(39。display39。,39。off39。)。 for m=1:length(th2) th34(m,:)=fsolve(39。fourbarposition39。,[1 1],options,th2(m),L2,L3,L4,L1)。 end y=L2*sin(th2)+L3*sin(th34(:,1)39。)。 x=L2*cos(th2)+L3*cos(th34(:,1)39。)。 xx=[L2*cos(th2)]。 yy=[L2*sin(th2)]。 plot([x。xx],[y。yy],39。k39。,[0 L1],[0 0],39。k^39。,x,y,39。ko39。,xx,yy,39。ks39。) title(39。連桿3的幾個位置點39。) xlabel(39。水平方向39。) ylabel(39。垂直方向39。) axis equal grid oncase 39。角位移線圖39。 plot(th2*180/pi,th34(:,1)*180/pi,th2*180/pi,th34(:,2)*180/pi) axis([0 360 0 170]) grid on xlabel(39。主動件轉角\theta_2(度)39。) ylabel(39。從動件角位移(度)39。) title(39。角位移線圖39。) text(120,120,39。搖桿4角位移39。) text(150,40,39。連桿3角位移39。)case 39。角速度線圖39。 plot(th2*180/pi,w3,th2*180/pi,w4)。 axis([0 360 175 200]) text(50,160,39。搖桿4角速度(\omega_4)39。) text(220,130,39。連桿3角速度(\omega_3)39。) grid on xlabel(39。主動件轉角\theta_2(度)39。) ylabel(39。從動件角速度(rad\cdot s^{1})39。) title(39。角速度線圖39。)case 39。角加速度線圖39。plot(th2*180/pi,a3,th2*180/pi,a4)。 axis([0 360 50000 65000]) text(50,50000,39。搖桿4角加速度(\alpha_4)39。) text(220,12000,39。連桿3角加速度(\alpha_3)39。) grid on xlabel(39。主動件轉角\theta_2(度)39。) ylabel(39。從動件角加速度(rad\cdot s^{2})39。) title(39。角加速度線圖39。)endguidata(hObject,handles)%其中，guidata(hObject,handles)命令用于更新句柄，當輸入不同參數(shù)是，程序能夠做出相應的相應。(3)打開M文件編輯器(Mfile Edit),點擊 向下的三角圖標，在Open菜單回調(diào)函數(shù)下添加以下代碼：file = uigetfile(39。*.fig39。)。if ~isequal(file, 0) open(file)。end此菜單用以打開fig文件。其相應界面如圖9所示： 圖9(4)打開M文件編輯器(Mfile Edit),點擊 向下的三角圖標，在Print菜單回調(diào)函數(shù)下添加以下代碼：printdlg() %用于圖形輸出。(5)打開M文件編輯器(Mfile Edit),點擊 向下的三角圖標，在Close菜單回調(diào)函數(shù)下添加以下代碼：selection = questdlg([39。Close 39。 get(,39。Name39。) 39。?39。],... [39。Close 39。 get(,39。Name39。) 39。...39。],39。Yes39。,39。No39。,39。Yes39。)。if strcmp(selection,39。No39。) return。enddelete()此菜單用于關閉界面，其響應界面如圖10所示： 圖10(6)完成M文件編寫后，運行程序進行檢驗，單擊 圖標或M文件工具欄的 圖標，在編輯欄輸入個構件參數(shù)，運行結果如圖11： 圖11 界面運行效果圖 分別點擊下拉菜單 下各選項，則在繪圖區(qū)域分別顯示圖2~圖5。4 小結在本次基于MATLAB平面四桿機構機構的運動分析課程設計中，不僅用到了MATLAB的m語言編程，還有GUI圖形用戶界面設計。對于課程的設計來說，m語言本身來說，功能相當強大，但是由于其復雜的編程方法，讓大多數(shù)初學者望而卻步；而GUI圖形界面則正好彌補了它的不足，它采用的是所見即所得的編程方式，用它來做軟件的界面就如圖制作網(wǎng)頁一樣簡單明了，用它制作出來的軟件不需要太多的編程知識都可以輕松搞定；在此次課程設計中，本人盡量揚長避短，把這些工具的優(yōu)點結合到一起，發(fā)揮其最大的作用。通過此次課程設計，不僅使我們對使我們對四桿機構有了進一步認識，并且使我們認識到計算機技術對工程應用的重要性。本此設計凝結了團隊合作的結晶，是我們利用三周的時間努力學習、設計的成果，同時，在設計過成中得到了鄭XX老師的悉心指導，在此表示衷心的感謝！參考文獻 [1] 孫桓,[M].:高等教育出版社,2006.[2] .[M].1版. 長沙:中南工業(yè)大學出版社,1995.[3] MATLAB原理與工程應用[M].:電子工業(yè)出版社,2002.[4] [5] [6] 22