% 內槽輪機構運動分析
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([-3*pi/4/dr 3*pi/4/dr -1 1])               % 定義坐標軸范圍
grid                                            % 柵格線
text(-85,-0.2,'\epsilon/\omega^{2}')            % 標注類角加速度線圖
text(10,0.55,'\omega/\omega')                   % 標注類角速度線圖
ylabel('槽輪運動線圖')                          % 定義縱坐標軸名稱
%
subplot(2,2,2);                                 % 選擇第2個子窗口
plot(c6(:,1),c6(:,3),c6(:,1),c6(:,4))           % 繪制z= 6的線圖
title('內槽輪槽數 z=6')
axis([-3*pi/4/dr 3*pi/4/dr -0.6 0.6])
grid
text(-85,-0.3,'\epsilon/\omega^{2}')
text(10,0.4,'\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/4/dr 3*pi/4/dr -0.40 0.40])
grid
text(-85,-0.3,'\epsilon/\omega^{2}')
text(10,0.32,'\omega/\omega')
ylabel('槽輪運動線圖')
%
subplot(2,2,4);                                 % 選擇第4個子窗口
plot(c10(:,1),c10(:,3),c10(:,1),c10(:,4))       % 繪制z=10的線圖
title('內槽輪槽數 z=10')
axis([-3*pi/4/dr 3*pi/4/dr -0.35 0.35])
grid
text(-80,-0.22,'\epsilon/\omega^{2}')
text(10,0.28,'\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([-3*pi/4/dr 3*pi/4/dr -0.05 0.45])
grid
text(-12,0.43,'z=4')
text(-12,0.35,'z=6')
text(-12,0.29,'z=8')
text(-15,0.21,'z=10')
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([-3*pi/4/dr 3*pi/4/dr -0.6 0.6])
grid
text(40,0.22,'z=10')
text(-60,-0.15,'z=8')
text(-110,-0.45,'z=6')
text(60,0.05,'z=4')
ylabel('槽輪類角加速度線圖')