pos=cell（130）;
for?i=1:30
????pos{i}=[rand（11）*10rand（11）*10rand（11）*10]
end
x=zeros（130）;
y=zeros（130）;
z=zeros（130）;
v=0.1;
direct1=cell（130）;
direct2=cell（130）;
direct3=cell（130）;
direct4=cell（130）;
direct5=cell（130）;
direct6=cell（130）;
direct7=cell（130）;
direct8=cell（130）;
cons=zeros（130）;
for?i=1:30
????direct1{i}=[rand（11）*10rand（11）*10rand（11）*10]-pos{i};
????direct2{i}=zeros（13）;
????direct3{i}=zeros（13）;
????direct4{i}=zeros（13）;
end
for?t=1:1000
????for?i=1:30
????????m=0;
????????n=0;
????????k=0;
????????l=[000];
????????direct7{i}=zeros（13）;
????????direct8{i}=zeros（13）;
????????for?j=1:30
????????????if?j==i
????????????????continue
????????????else
????????????????for?r=1:30
????????????????????direct7{i}=（direct7{i}+pos{r}-pos{i}）./r;
????????????????????direct8{i}=（direct8{i}+direct1{r}）./r;
????????????????end
????????????????if?（norm（pos{i}-pos{j}）<=5）&&（norm（pos{i}-pos{j}）>=0.1）
????????????????????m=m+1;
????????????????????k=k+norm（direct1{j}）;
????????????????????l=l+direct1{j};
????????????????????direct2{i}=（direct2{i}+（pos{j}-pos{i}））./m;
????????????????????direct3{i}=（direct3{i}+direct1{j}）./m;
????????????????????cons（i）=norm（l）/k;
????????????????elseif?norm（pos{i}-pos{j}）<0.1;
????????????????????n=n+1;
????????????????????direct4{i}=（direct2{i}+（pos{i}-pos{j}））./n;
????????????????end
????????????end
????????end
????????if?（pos{i}（1）>=10）||（pos{i}（2）>=10）||（pos{i}（1）<=0）||（pos{i}（2）<=0）||（pos{i}（3）>=10）||（pos{i}（3）<=0）
????????????direct1{i}=direct1{i}*（-1）;
????????else
?????????????direct1{i}=0.1*direct1{i}+0.1*direct2{i}+0.1*direct3{i}+0.1*direct4{i}+0.3*direct7{i}+0.3*direct8{i};
????????end
????????pos{i}=pos{i}+direct1{i}/norm（direct1{i}）*v*exp（5*（cons（i）-1））;
????end
????for?i=1:30
????????x（i）=pos{i}（1）;
????????y（i）=pos{i}（2）;
????????z（i）=pos{i}（3）;
????end
????plot3（xyz‘^‘）
????axis（[0?10?0?10?0?10]‘square‘‘manual‘）
????grid?on;
????pause（0.01）
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-16?18:28??Boid模型代碼\
?????文件????????2028??2018-04-16?18:26??Boid模型代碼\Untitled.m
?????文件????????3590??2018-04-16?18:27??Boid模型代碼\Untitled2.m