[問題] 求救~請教如何讓迴圈快一點 謝謝~已回收
小弟想用數值方法(循環座標下降法CCD)來解一組方程式
請教各位先進 如何讓我的迴圈執行快一點
有哪些地方可以改的 會不會用太多變數?
懇請指教 感謝^^
syms a1 a2 a3 b1 b2 b3 r1 r2 r3 d1x d2x d1y d2y d1z d2z fx fy fz ux uy uz l2 l3 t2 t3 c3 c2 o3 o2 vx vy vz kx ky
syms kz epx epy epz dx dy dz
rft2=[fx^2*(1-cos(t2))+cos(t2),fx*fy*(1-cos(t2))-fz*sin(t2),fx*fz*(1-cos(t2))+fy*sin(t2);fx*fy*(1-cos(t2))+fz*sin(t2),...
fy^2*(1-cos(t2))+cos(t2),fy*fz*(1-cos(t2))-fx*sin(t2);fx*fz*(1-cos(t2))-fy*sin(t2),fy*fz*(1-cos(t2))+fx*sin(t2),fz^2*(1-cos(t2))+cos(t2)];
rft3=[fx^2*(1-cos(t3))+cos(t3),fx*fy*(1-cos(t3))-fz*sin(t3),fx*fz*(1-cos(t3))+fy*sin(t3);fx*fy*(1-cos(t3))+fz*sin(t3),...
fy^2*(1-cos(t3))+cos(t3),fy*fz*(1-cos(t3))-fx*sin(t3);fx*fz*(1-cos(t3))-fy*sin(t3),fy*fz*(1-cos(t3))+fx*sin(t3),fz^2*(1-cos(t3))+cos(t3)];
rzc2=[0^2*(1-cos(c2))+cos(c2),0*0*(1-cos(c2))-1*sin(c2),0*1*(1-cos(c2))+0*sin(c2);0*0*(1-cos(c2))+1*sin(c2),...
0^2*(1-cos(c2))+cos(c2),0*1*(1-cos(c2))-0*sin(c2);0*1*(1-cos(c2))-0*sin(c2),0*1*(1-cos(c2))+0*sin(c2),1^2*(1-cos(c2))+cos(c2)];
rzc3=[0^2*(1-cos(c3))+cos(c3),0*0*(1-cos(c3))-1*sin(c3),0*1*(1-cos(c3))+0*sin(c3);0*0*(1-cos(c3))+1*sin(c3),...
0^2*(1-cos(c3))+cos(c3),0*1*(1-cos(c3))-0*sin(c3);0*1*(1-cos(c3))-0*sin(c3),0*1*(1-cos(c3))+0*sin(c3),1^2*(1-cos(c3))+cos(c3)];
epx2=(rzc2*[epx;epy;epz]);
epx2=epx2(1);
epy2=(rzc2*[epx;epy;epz]);
epy2=epy2(2);
epz2=(rzc2*[epx;epy;epz]);
epz2=epz2(3);
epx3=(rzc3*[epx;epy;epz]);
epx3=epx3(1);
epy3=(rzc3*[epx;epy;epz]);
epy3=epy3(2);
epz3=(rzc3*[epx;epy;epz]);
epz3=epz3(3);
rep2o2=[epx2^2*(1-cos(o2))+cos(o2),epx2*epy2*(1-cos(o2))-epz2*sin(o2),epx2*epz2*(1-cos(o2))+epy2*sin(o2);epx2*epy2*(1-cos(o2))+epz2*sin(o2),...
epy2^2*(1-cos(o2))+cos(o2),epy2*epz2*(1-cos(o2))-epx2*sin(o2);epx2*epz2*(1-cos(o2))-epy2*sin(o2),epy2*epz2*(1-cos(o2))+epx2*sin(o2),epz2^2*(1-cos(o2))+cos(o2)];
rep3o3=[epx3^2*(1-cos(o3))+cos(o3),epx3*epy3*(1-cos(o3))-epz3*sin(o3),epx3*epz3*(1-cos(o3))+epy3*sin(o3);epx3*epy3*(1-cos(o3))+epz3*sin(o3),...
epy3^2*(1-cos(o3))+cos(o3),epy3*epz3*(1-cos(o3))-epx3*sin(o3);epx3*epz3*(1-cos(o3))-epy3*sin(o3),epy3*epz3*(1-cos(o3))+epx3*sin(o3),epz3^2*(1-cos(o3))+cos(o3)];
rfl2=[fx^2*(1-cos(l2))+cos(l2),fx*fy*(1-cos(l2))-fz*sin(l2),fx*fz*(1-cos(l2))+fy*sin(l2);fx*fy*(1-cos(l2))+fz*sin(l2),...
fy^2*(1-cos(l2))+cos(l2),fy*fz*(1-cos(l2))-fx*sin(l2);fx*fz*(1-cos(l2))-fy*sin(l2),fy*fz*(1-cos(l2))+fx*sin(l2),fz^2*(1-cos(l2))+cos(l2)];
rfl3=[fx^2*(1-cos(l3))+cos(l3),fx*fy*(1-cos(l3))-fz*sin(l3),fx*fz*(1-cos(l3))+fy*sin(l3);fx*fy*(1-cos(l3))+fz*sin(l3),...
fy^2*(1-cos(l3))+cos(l3),fy*fz*(1-cos(l3))-fx*sin(l3);fx*fz*(1-cos(l3))-fy*sin(l3),fy*fz*(1-cos(l3))+fx*sin(l3),fz^2*(1-cos(l3))+cos(l3)];
rot1=[cos(r1)*cos(b1),-sin(r1),cos(r1)*sin(b1);cos(a1)*sin(r1)*cos(b1)+sin(a1)*sin(b1),cos(a1)*cos(r1),cos(a1)*sin(r1)*sin(b1)-sin(a1)*cos(b1);...
sin(a1)*sin(r1)*cos(b1)-cos(a1)*sin(b1),sin(a1)*cos(r1),sin(a1)*sin(r1)*sin(b1)+cos(a1)*cos(b1)];
rot2=[cos(r2)*cos(b2),-sin(r2),cos(r2)*sin(b2);cos(a2)*sin(r2)*cos(b2)+sin(a2)*sin(b2),cos(a2)*cos(r2),cos(a2)*sin(r2)*sin(b2)-sin(a2)*cos(b2);...
sin(a2)*sin(r2)*cos(b2)-cos(a2)*sin(b2),sin(a2)*cos(r2),sin(a2)*sin(r2)*sin(b2)+cos(a2)*cos(b2)];
rot3=[cos(r3)*cos(b3),-sin(r3),cos(r3)*sin(b3);cos(a3)*sin(r3)*cos(b3)+sin(a3)*sin(b3),cos(a3)*cos(r3),cos(a3)*sin(r3)*sin(b3)-sin(a3)*cos(b3);...
sin(a3)*sin(r3)*cos(b3)-cos(a3)*sin(b3),sin(a3)*cos(r3),sin(a3)*sin(r3)*sin(b3)+cos(a3)*cos(b3)];
I=eye(3);
multi1=(rzc2-I)*[vx;vy;vz]+(rep2o2*rzc2-I)*[kx;ky;kz]+(rot2-rot1)*[dx;dy;dz];
multi2=(rzc3-I)*[vx;vy;vz]+(rep3o3*rzc3-I)*[kx;ky;kz]+(rot3-rot1)*[dx;dy;dz];
f1=multi1(1)-d1x;f2=multi1(2)-d1y;f3=multi1(3)-d1z;
f4=multi2(1)-d2x;f5=multi2(2)-d2y;f6=multi2(3)-d2z;
f7=(rot1*[dx;dy;dz])'*[ux;uy;uz]-(rot2*[dx;dy;dz])'*(rfl2*rft2*[ux;uy;uz]);
f8=(rot1*[dx;dy;dz])'*[ux;uy;uz]-(rot3*[dx;dy;dz])'*(rfl3*rft3*[ux;uy;uz]);
f9=epx^2+epy^2+epz^2-1;
F=f1^2+f2^2+f3^2+f4^2+f5^2+f6^2+f7^2+f8^2+f9^2;
FF=subs(F,[d1x,d1y,d1z,d2x,d2y,d2z,a1,b1,r1,a2,b2,r2,a3,b3,r3],[-10,-6.5,20,-12,-7.5,26,pi/6,0,0,pi/4,0,pi/12,pi/3,0,pi/6]);
FFF=subs(FF,[fx,fy,fz,l2,l3,t2,t3,ux,uy,uz,o2,vz,kz,c2,c3,vy,epz],[0,0,1,6.33,5.87,pi/6,pi/3,-10,124.36,125.9,pi/12,40,70,pi/6,pi/3,0,0.5]);
x=[vx kx ky dx dy dz epx epy o3];
y=[10 10 10 10 10 10 10 10 10];
val=subs(FFF,x,y)
while val>0.1
ans1=subs(FFF,[x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9)],[y(2),y(3),y(4),y(5),y(6),y(7),y(8),y(9)]);
ans11=@(vx) eval(ans1);
[y(1),val]=fminbnd(ans11,-200000,200000);
ans2=subs(FFF,[x(1),x(3),x(4),x(5),x(6),x(7),x(8),x(9)],[y(1),y(3),y(4),y(5),y(6),y(7),y(8),y(9)]);
ans22=@(kx) eval(ans2);
[y(2),val]=fminbnd(ans22,-2000,2000);
ans3=subs(FFF,[x(1),x(2),x(4),x(5),x(6),x(7),x(8),x(9)],[y(1),y(2),y(4),y(5),y(6),y(7),y(8),y(9)]);
ans33=@(ky) eval(ans3);
[y(3),val]=fminbnd(ans33,-2000,2000);
ans4=subs(FFF,[x(1),x(2),x(3),x(5),x(6),x(7),x(8),x(9)],[y(1),y(2),y(3),y(5),y(6),y(7),y(8),y(9)]);
ans44=@(dx) eval(ans4);
[y(4),val]=fminbnd(ans44,-2000,2000);
ans5=subs(FFF,[x(1),x(2),x(3),x(4),x(6),x(7),x(8),x(9)],[y(1),y(2),y(3),y(4),y(6),y(7),y(8),y(9)]);
ans55=@(dy) eval(ans5);
[y(5),val]=fminbnd(ans55,-2000,2000);
ans6=subs(FFF,[x(1),x(2),x(3),x(4),x(5),x(7),x(8),x(9)],[y(1),y(2),y(3),y(4),y(5),y(7),y(8),y(9)]);
ans66=@(dz) eval(ans6);
[y(6),val]=fminbnd(ans66,-2000,2000);
ans7=subs(FFF,[x(1),x(2),x(3),x(4),x(5),x(6),x(8),x(9)],[y(1),y(2),y(3),y(4),y(5),y(6),y(8),y(9)]);
ans77=@(epx) eval(ans7);
[y(7),val]=fminbnd(ans77,-2000,2000);
ans8=subs(FFF,[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(9)],[y(1),y(2),y(3),y(4),y(5),y(6),y(7),y(9)]);
ans88=@(epy) eval(ans8);
[y(8),val]=fminbnd(ans88,-2000,2000);
ans9=subs(FFF,[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8)],[y(1),y(2),y(3),y(4),y(5),y(6),y(7),y(8)]);
ans99=@(o3) eval(ans9);
[y(9),val]=fminbnd(ans99,0,2*pi);
val
y
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.118.197.101