%?function?four_pole?
%?四連桿機構?
%close?all;clc;clear;?
%?xm=[50.048.542342930344248555651];
%?ym=[91111107906745281712142452];
%?plot（xmym‘ro‘）
%?hold?on;
global?l;
global?beta;
global?gama1;
global?XA;
global?YA;
%?l=[45.0157171.1939113.1270108.943857.66472];
%?beta=0.5695;
%?gama1=0.22643;
%?XA=67;
%?YA=10;
%?l1=l;
%?beta1=beta;
%?gama2=gama1;
%?XA1=XA;
%?YA1=YA;
t1=sin（0:.1:2*pi）;?%?生成固定點的橫軸數據
t2=cos（0:.1:2*pi）;?%?生成固定點的縱軸數據?
r=2;?%?固定點的半徑
fill（XA+r*t1YA+r*t2‘r‘）;set（gcf‘doublebuffer‘‘on‘）;?%?繪制左側固定點
hold?on;?
Cx=XA+real（l（4）*exp（i*beta））;%計算右端點的橫坐標
Cy=YA+imag（l（4）*exp（i*beta））;%計算右端點的縱坐標
fill（Cx+r*t1Cy+r*t2‘r‘）;?%?繪制右側固定點?
axis（[-100200-100150]）;?%?設置坐標軸范圍及屬性
t=1;?%?記錄時間的參數
A=t;?%?旋轉的角度
z=l（1）*exp（i*A）;
Ax=XA+real（l（1）*exp（i*A））;?%?圓周運動的端點的橫坐標
Ay=YA+imag（l（1）*exp（i*A））;%圓周遠動的端點的縱坐標
m=l（4）*exp（i*beta）;
H1=plot（[XAAx][YAAy]‘linewidth‘4）;?%?畫出作圓周運動的連桿OA
AC=sqrt（（Cx-Ax）^2+（Cy-Ay）^2）;?%?計算當前線段AC的長度
AB=l（2）;?%?設置AB的長度，其為固定值
BC=l（3）;?%?設置BC的長度，其為固定值
JA=acos（（AC^2+AB^2-BC^2）/（2*AC*AB））?;%?計算AB和AC夾角大小
Bx=Ax+real（AB*exp（（（angle（m-z））+JA）*i））;%?獲得B點當前的橫坐標
By=Ay+imag（AB*exp（（（angle（m-z））+JA）*i））;%獲得B點當前的縱坐標
%B=z+AB*exp（（（angle（m-z））+JA）*i）;?%?獲得B點當前的坐標
Dx=Ax+real（l（5）*exp（（angle（m-z）+JA+gama1）*i））;%獲得D點當前的橫坐標
Dy=Ay+imag（l（5）*exp（（angle（m-z）+JA+gama1）*i））;%獲得D點當前的縱坐標
%D=z+57.66472*exp（（angle（m-z）+JA+0.2264）*i）;
k=1;?%?控制循環是否執行的參數?
H2=plot（[AxBx][AyBy]‘g‘‘linewidth‘4）;?%?繪制AB連桿
H3=plot（[BxCx][ByCy]‘linewidth‘4‘color‘[0.40.20.5]）;??%?繪制BC連桿
H4=plot（[AxDx][AyDy]‘y?‘‘linewidth‘4）;?%繪制AD連桿
%Ga=plot（AxAy）;?%?繪制A點當前的軌跡曲線
%Gb=plot（BxBy）;?%?繪制B點當前的軌跡曲線?
Gd=plot（DxDy）;%?繪制D點軌跡
%set（[GaGb]‘color‘‘r‘）;?%?設置軌跡曲線的顏色
%zga=z;?%?定義記錄A點軌跡點的變量
%zgb=B;?%?定義記錄B點軌跡點的變量?
zgdx=Dx;
zgdy=Dy;
xlabel（‘Please?press?“space“?key?and?stop?this?program!‘...?
???‘fontsize‘14‘color‘‘r‘）;?%?加注標注
T=title（[‘time=??‘num2str（t）]）;?%?實時顯示時間
text（0100‘鉸鏈四連桿機構‘）;?%?加注文字
while?k;?%?模擬四連桿的運動過程?
???s=get（gcf‘currentkey‘）;?%?檢測鍵入按鍵的名稱
???if?strcmp（s‘space‘）;?%?判斷鍵入按鍵是否為空格鍵
???????clc;k=0;?%?設置參數使循環結束
???end?????
???pause（0.01）;?%?暫停一下，顯示動畫效果
???t=t+0.02;?%?更新時間?
???z=l（1）*exp（i*t）;?
???Ax=XA+real（l（1）*exp（i*t））;?%?圓周運動的端點的橫坐標
???Ay=YA+imag（l（1）*exp（i*t））;%圓周遠動的端點的縱坐標
???AC=sqrt（（Cx-Ax）^2+（Cy-Ay）^2）;?%?計算當前線段AC的長度
???%AC=abs（m-z）;?%?計算當前線段AC的長度
???JA=acos（（AC^2+AB^2-BC^2）/（2*AC*AB））;?%?計算當前線段AB和BC的夾角
???%B=z+AB*exp（（（angle（m-z））+JA）*i）;?%?獲得B點當前的坐標
???Bx=Ax+real（AB*exp（（（angle（m-z））+JA）*i））;%?獲得B點當前的橫坐標
???By=Ay+imag（AB*exp（（（angle（m-z））+JA）*i））;%獲得B點當前的縱坐標
???Dx=Ax+real（l（5）*exp（（angle（m-z）+JA+gama1）*i））;%獲得D點當前的橫坐標
???Dy=Ay+imag（l（5）*exp（（angle（m-z）+JA+gama1）*i））;%獲得D點當前的縱坐標
???%D=z+57.66472*exp（（angle（m-z）+JA+0.2264）*i）;
???set（H1‘xdata‘[XAAx]‘ydata‘[YAAy]）;?%?更新線段OA的位置?
???set（H2‘xdata‘[AxBx]‘ydata‘[AyBy]）;?%?更新線段AB的位置
???set（H3‘xdata‘[BxCx]‘ydata‘[ByCy]）;?%?更新線段BC的位置
???set（H4‘xdata‘[DxAx]‘ydata‘[DyAy]）;?%更新線段AD的位

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3886??2009-12-07?10:57??Mechanisms5\fourbarmotion2.asv

?????文件???????3968??2009-12-07?15:05??Mechanisms5\four_pole.asv

?????文件???????3973??2009-12-23?17:18??Mechanisms5\four_pole.m

?????文件??????42034??2009-11-04?16:44??Mechanisms5\fuorbarGUI.asv

?????文件???????7274??2009-12-07?15:15??Mechanisms5\fuorbarGUI.fig

?????文件??????42011??2009-12-07?11:17??Mechanisms5\fuorbarGUI.m

?????文件????????721??2009-10-31?21:32??Mechanisms5\my.asv

?????文件????????702??2009-12-23?17:17??Mechanisms5\my.m

?????文件????????771??2009-10-30?19:21??Mechanisms5\myfun.asv

?????文件????????741??2009-11-01?17:34??Mechanisms5\myfun.m

?????文件?????586752??2009-11-12?20:06??Mechanisms5\圖片.doc

?????目錄??????????0??2009-12-23?17:20??Mechanisms5

-----------?---------??----------?-----??----

???????????????692833????????????????????12