?? 04-槽輪機構的運動分析-1.m__.htm
字號:
<html><head><title>MATLAB.rar 04-槽輪機構的運動分析-1.m</title>
<LINK href="/inc/read_style.css" type=text/css rel=stylesheet></head>
<body>
<p><a href=http://www.pudn.com>www.pudn.com</a> > <a href="http://www.pudn.com/downloads111/sourcecode/graph/detail461284.html">MATLAB.rar</a> > 04-槽輪機構的運動分析-1.m</p><!-- saved from http://www.pudn.com -->
<script src="/inc/gg_read1.js"></script><BR>
<pre name="code" class="m">
% 外槽輪機構運動分析
dr=pi/180.0; % 角度與弧度的轉換系數
% 銷輪2轉角范圍:-f20<f2<f20,步長為bc度,計算運動參數
for z=4:2:10 % 設定槽輪槽數
f30=pi/z; % 計算槽輪槽間半角
f20=pi/2-f30; % 計算銷輪運動半角
lmd=sin(pi/z); % 計算曲柄2與機架1的長度比
bc=10; % 循環步長
cz=-f20/dr; % 循環初值
zz=f20/dr; % 循環終值
i=1; % 根據步長變化的運動參數矩陣cs行數計數器
for f2=cz:bc:zz % 計算槽輪角位移、類角速度、類角加速度
wy=atan(lmd*sin(f2*dr)/(1-lmd*cos(f2*dr)));
sd=lmd*(cos(f2*dr)-lmd)/(1-2*lmd*cos(f2*dr)+lmd^2);
jsd=-lmd*sin(f2*dr)*(1-lmd^2)/(1-2*lmd*cos(f2*dr)+lmd^2)^2;
switch z % 矩陣c(i,:)表示第i行的各列元素
case 4,c4(i,:)=[f2 wy/dr sd jsd];
case 6,c6(i,:)=[f2 wy/dr sd jsd];
case 8,c8(i,:)=[f2 wy/dr sd jsd];
case 10,c10(i,:)=[f2 wy/dr sd jsd];
end
i=i+1;
end
end
% 輸出外槽輪機構運動參數
['輪槽數 z=4']
[' 銷輪轉角',' 槽輪角位移',' 角速度',' 角加速度']
% 矩陣c(:,j)表示第j列的各行元素
[c4(:,1),c4(:,2),c4(:,3),c4(:,4)]
['輪槽數 z=6']
[' 銷輪轉角',' 槽輪角位移',' 角速度',' 角加速度']
% 矩陣c(:,j)表示第j列的各行元素
[c6(:,1),c6(:,2),c6(:,3),c6(:,4)]
['輪槽數 z=8']
[' 銷輪轉角',' 槽輪角位移',' 角速度',' 角加速度']
% 矩陣c(:,j)表示第j列的各行元素
[c8(:,1),c8(:,2),c8(:,3),c8(:,4)]
['輪槽數 z=10']
[' 銷輪轉角',' 槽輪角位移',' 角速度',' 角加速度']
% 矩陣c(:,j)表示第j列的各行元素
[c10(:,1),c10(:,2),c10(:,3),c10(:,4)]
%
% 繪制槽輪機構運動參數曲線
figure(1); % 生成槽輪運動線圖窗口
subplot(2,2,1); % 選擇第1個子窗口
plot(c4(:,1),c4(:,3),c4(:,1),c4(:,4)) % 繪制z= 4的線圖
title('外槽輪槽數 z=4') % 標注子窗口名稱
axis([-pi/4/dr pi/4/dr -6 6]) % 定義坐標軸范圍
grid % 柵格線
text(-2,4.2,'\epsilon/\omega^{2}') % 標注類角加速度線圖
text(20,1.6,'\omega/\omega') % 標注類角速度線圖
ylabel('槽輪運動線圖') % 定義縱坐標軸名稱
%
subplot(2,2,2); % 選擇第2個子窗口
plot(c6(:,1),c6(:,3),c6(:,1),c6(:,4)) % 繪制z= 6的線圖
title('外槽輪槽數 z=6')
axis([-pi/3/dr pi/3/dr -1.5 1.5])
grid
text(10,-0.7,'\epsilon/\omega^{2}')
text(30,0.7,'\omega/\omega')
ylabel('槽輪運動線圖')
%
subplot(2,2,3); % 選擇第3個子窗口
plot(c8(:,1),c8(:,3),c8(:,1),c8(:,4)) % 繪制z= 8的線圖
title('外槽輪槽數 z=8')
axis([-3*pi/8/dr 3*pi/8/dr -0.8 0.8])
grid
text(12,-0.3,'\epsilon/\omega^{2}')
text(40,0.4,'\omega/\omega')
ylabel('槽輪運動線圖')
%
subplot(2,2,4); % 選擇第4個子窗口
plot(c10(:,1),c10(:,3),c10(:,1),c10(:,4)) % 繪制z=10的線圖
title('外槽輪槽數 z=10')
axis([-2*pi/5/dr 2*pi/5/dr -0.5 0.5])
grid
text(15,-0.2,'\epsilon/\omega^{2}')
text(40,0.3,'\omega/\omega')
ylabel('槽輪運動線圖')
%
figure(2); % 生成類線圖窗口
subplot(1,2,1); % 選擇第1個子窗口
plot(c4(:,1),c4(:,3),c6(:,1),c6(:,3),c8(:,1),c8(:,3),c10(:,1),c10(:,3))
title('\omega/\omega')
axis([-f20/dr f20/dr -0.1 2.5])
grid
text(-10,0.35,'z=10')
text(-8,0.7,'z=8')
text(-8,1.1,'z=6')
text(10,2.1,'z=4')
ylabel('槽輪類角速度線圖')
%
subplot(1,2,2); % 選擇第2個子窗口
plot(c4(:,1),c4(:,4),c6(:,1),c6(:,4),c8(:,1),c8(:,4),c10(:,1),c10(:,4))
title('\epsilon/\omega^{2}')
axis([-f20/dr f20/dr -5.5 5.5])
grid
text(-50,0.2,'z=10')
text(-30,0.9,'z=8')
text(-25,1.6,'z=6')
text(0,3.5,'z=4')
ylabel('槽輪類角加速度線圖')
</pre>
<script src="/inc/gg_read2.js"></script><BR>
<script src="http://s117.cnzz.com/stat.php?id=1236358&web_id=1236358&show=pic" language="JavaScript" charset="gb2312"></script>
</body></html>
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -