huatu.asv

这个一个遗传算法的实例

clear
D=[50,7,1000];%依次表示行,列,页
N=3;n=8;%依次表示轴数和插值点数
%以P表示位置,V表示速度,A表示加速度,J表示加加速度
%T表示时间,h表示时间间隔
VC=[100,95,1000];%依次设置第一轴,第二轴,第三轴的速度约束
AC=[70,75,75];%依次设置第一轴,第二轴,第三轴的加速度约束
JC=[70,85,70];%依次设置第一轴,第二轴,第三轴的冲击约束
V=zeros(N,n);%设置速度边界
A=zeros(N,n);%设置加速度边界

t1=0;%起始时间

P(1,1)=120;P(1,3)=100;P(1,4)=60;P(1,5)=30;P(1,6)=10;P(1,8)=0;
P(2,1)=-10;P(2,3)=20;P(2,4)=0;P(2,5)=30;P(2,6)=60;P(2,8)=80;
P(3,1)=0;P(3,3)=20;P(3,4)=50;P(3,5)=70;P(3,6)=50;P(3,6)=30;

H_ZY=[1.2733    1.1170    1.4075    1.3836    1.1348    0.6830    0.9450];%输入量
H_ZY_2=[1.2516    1.4499    1.8011    1.6056    1.5306    0.8843    1.0704];%输入量

T_HT=SOV_T(H_ZY,n,t1);
T_HT_2=SOV_T(H_ZY_2,n,t1);

QA_HT=SOL_QA(H_ZY,V,A,P,N,n);
QA_HT_2=SOL_QA(H_ZY_2,V,A,P,N,n);

P_HT=SOV_P(P,QA_HT,H_ZY,V,N,n);
P_HT_2=SOV_P(P,QA_HT_2,H_ZY_2,V,N,n);

HTBL=3;%可选择1,2,3表示第一,二,三轴
%for HTBL=1:3;
PHTL=P_HT(HTBL,:);%画图的位置
PHTL_2=P_HT_2(HTBL,:);

QAHTL=QA_HT(HTBL,:);%画图用的加速度
QAHTL_2=QA_HT_2(HTBL,:);

VL=V(HTBL,:);%zeros(1,n);%画图用的速度

HL=H_ZY;%画图用的时间间隔
HL_2=H_ZY_2;

TL=T_HT;%画图用的插值点时间
TL_2=T_HT_2;

[QJ_MAX,QJ_ZH]=SOV_QJ(QA_HT,HL,N,n)%求出最大的冲击
[QJ_MAX_2,QJ_ZH_2]=SOV_QJ(QA_HT_2,HL_2,N,n)

T1=TL(1):0.01:TL(2);
T2=TL(2):0.01:TL(3);
T3=TL(3):0.01:TL(4);
T4=TL(4):0.01:TL(5);
T5=TL(5):0.01:TL(6);
T6=TL(6):0.01:TL(7);
T5=TL(7):0.01:TL(8);

T1_2=TL_2(1):0.01:TL_2(2);
T2_2=TL_2(2):0.01:TL_2(3);
T3_2=TL_2(3):0.01:TL_2(4);
T4_2=TL_2(4):0.01:TL_2(5);
T5_2=TL_2(5):0.01:TL_2(6);


J1=(T1*(QAHTL(2)-QAHTL(1))/HL(1))./(T1+eps);
J2=T2*(QAHTL(3)-QAHTL(2))/HL(2)./T2;
J3=T3*(QAHTL(4)-QAHTL(3))/HL(3)./T3;
J4=T4*(QAHTL(5)-QAHTL(4))/HL(4)./T4;
J5=T5*(QAHTL(6)-QAHTL(5))/HL(4)./T5;

J1_2=(T1_2*(QAHTL_2(2)-QAHTL_2(1))/HL_2(1))./(T1_2+eps);
J2_2=T2_2*(QAHTL_2(3)-QAHTL_2(2))/HL_2(2)./T2_2;
J3_2=T3_2*(QAHTL_2(4)-QAHTL_2(3))/HL_2(3)./T3_2;
J4_2=T4_2*(QAHTL_2(5)-QAHTL_2(4))/HL_2(4)./T4_2;
J5_2=T5_2*(QAHTL_2(6)-QAHTL_2(5))/HL_2(4)./T5_2;
%subplot(3,1,HTBL)
plot(T1,J1,'k-',T2,J2,'k-',T3,J3,'k-',T4,J4,'k-',T5,J5,'k-')
%plot(T1,A1,'k-',T2,A2,'k-',T3,A3,'k-',T4,A4,'k-',T5,A5,'k-')
%plot(T1,V1,'k-',T2,V2,'k-',T3,V3,'k-',T4,V4,'k-',T5,V5,'k-')
%plot(T1,Q1,'k-',T2,Q2,'k-',T3,Q3,'k-',T4,Q4,'k-',T5,Q5,'k-'),%legend('位移','速度','加速度','冲击')
hold on
%plot(T1_2,Q1_2,'k:',T2_2,Q2_2,'k:',T3_2,Q3_2,'k:',T4_2,Q4_2,'k:',T5_2,Q5_2,'k:'),%legend('位移','速度','加速度','冲击')
%plot(T1_2,V1_2,'k:',T2_2,V2_2,'k:',T3_2,V3_2,'k:',T4_2,V4_2,'k:',T5_2,V5_2,'k:')
%plot(T1_2,A1_2,'k:',T2_2,A2_2,'k:',T3_2,A3_2,'k:',T4_2,A4_2,'k:',T5_2,A5_2,'k:')
plot(T1_2,J1_2,'k:',T2_2,J2_2,'k:',T3_2,J3_2,'k:',T4_2,J4_2,'k:',T5_2,J5_2,'k:')
%axis(0,11,-23,23)
%xlabel('t/s');
%axis tight
%ylabel('angle');
hold off
%end;
%title('第三关节的位置、速度、加速度、冲击曲线')
%grid