-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathPartial3.m
35 lines (31 loc) · 1.26 KB
/
Partial3.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
%% 逕ィ莠手ョ。邂怜刈蜈・蝓コ蝮先�ウサ蛻ーCMM蝮先�ウサ霓ャ謐「遏ゥ髦オ譬。蜃�噪讓。蝙句¥蠕ョ蛻?
function [dC] = Partial2(SphereTrans)
%% 隶。邂苓スャ謐「遏ゥ髦オ逧�¥蠕ョ蛻�洸髦オ
% 霎灘�荳コ霓ャ謐「遏ゥ髦オ蜈ウ莠主推荳ェ蜈ウ闃ら噪蝗帑クェ蜿よ焚逧�¥蠕ョ蛻
%% 蛛丞セョ蛻�ョ。邂鈴。コ蠎丈クコp,a,t,d,蜀嶺ス吝盾謨ー荳コd2�慧6蜥荊6
syms a1 p1 d1 t1;
syms a2 p2 d2 t2;
syms a3 p3 d3 t3;
syms a4 p4 d4 t4;
syms a5 p5 d5 t5;
syms a6 p6 d6 t6;
syms x y z a b g;
% 蛻�悪荳コa�径lpha�慧�荊heta
%縲?桷蟒コ莉擦MM蝮先�ウサ蛻ー蝓コ蝮先�ウサ逧�スャ謐「遏ゥ髦?
T0 = Rotz(a)*Roty(b)*Rotx(g)+transl(x,y,z)-...
[1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1];
T1 = transl(0,0,d1)*Rotz(t1)*transl(a1,0,0)*Rotx(p1);
T2 = transl(0,0,d2)*Rotz(t2)*transl(a2,0,0)*Rotx(p2);
T3 = transl(0,0,d3)*Rotz(t3)*transl(a3,0,0)*Rotx(p3);
T4 = transl(0,0,d4)*Rotz(t4)*transl(a4,0,0)*Rotx(p4);
T5 = transl(0,0,d5)*Rotz(t5)*transl(a5,0,0)*Rotx(p5);
T6 = transl(0,0,d6)*Rotz(t6)*transl(a6,0,0)*Rotx(p6);
T7 = SphereTrans;
% 隶。邂佑MM蝮先�ウサ蛻ー蟾・蜈キ逅�攝譬�ウサ逧�スャ謐「遏ゥ髦?
T = T0*T1*T2*T3*T4*T5*T6*T7;
% 霓ャ蛹紋クコ��,y,z,psi,theta,phi�牙錘蛻�悪豎ょ¥蟇?
[C] = CalPos(T,0);
% 隶。邂怜ッケ蜷�クェ蜿よ焚逧�¥蟇シ
% mテ溶,m譏ッ蝮先�盾謨ー�系譏ッ豈丈クェ霑樊揉逧ДH蜿よ焚
dC = jacobian(C,[p1;a1;t1;d1;p2;a2;t2;p3;a3;t3;d3;p4;a4;t4;d4;p5;a5;t5;d5;p6;a6;x;y;z;a;b;g]);
end