Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Vx = j(nrow(jari),1); /*matrix initialization*/
- Vy = j(nrow(jari),1); /*matrix initialization*/
- do i = 1 to nrow(jari);
- /* velocity vector (vx,vy)*/
- Vx[i]= velo[i] * cosm[i];
- Vy[i]= velo[i] * sinm[i];
- end;
- print Vx, Vy;
- /* intersection point*/
- Xpcij = j(nrow(m), nrow(m));
- Ypcij = j(nrow(m), nrow(m));
- do i = 1 to nrow(m);
- do j=1 to nrow(m);
- Xpcij[i,j] = (koordinatY[j,]-((Vy[j,]/Vx[j,])*koordinatX[j,])-koordinatY[i,]-((Vy[i,]/Vx[i,])*koordinatX[i,]))/((Vy[i,]/Vx[i,])-(Vy[j,]/Vx[j,]));
- Ypcij[i,j] = ((Vy[i,]/Vx[i,])*(koordinatY[j,]-((Vy[j,]/Vx[j,])*koordinatX[j,]))-(Vy[j,]/Vx[j,])*(koordinatY[i,]-((Vy[i,]/Vx[i,])*koordinatX[i,])))/((Vy[i,]/Vx[i,])-(Vy[j,]/Vx[j,]));
- end;
- end;
- print Xpcij;
- print Ypcij;
- /* false conflict criteria*/
- d1ij = j(nrow(m), nrow(m));
- d1ji = j(nrow(m), nrow(m));
- d2ij = j(nrow(m), nrow(m));
- d2ji = j(nrow(m), nrow(m));
- do i = 1 to nrow(m);
- do j=1 to nrow(m);
- d1ij[i,j] = sqrt(((koordinatX[i,]-Xpcij[i,j])**2)+((koordinatY[i,]-Ypcij[i,j])**2));
- d1ji[i,j] = sqrt(((koordinatX[j,]-Xpcij[i,j])**2)+((koordinatY[j,]-Ypcij[i,j])**2));
- d2ij[i,j] = sqrt((((koordinatX[i,]+cosm[i,])-Xpcij[i,j])**2)+(((koordinatY[i,]+sinm[i,])-Ypcij[i,j])**2));
- d2ji[i,j] = sqrt((((koordinatX[j,]+cosm[j,])-Xpcij[i,j])**2)+(((koordinatY[j,]+sinm[j,])-Ypcij[i,j])**2));
- end;
- end;
- print d1ij;
- print d1ji;
- print d2ij;
- print d2ji;
- di= d2ij-d1ij;
- dj= d2ji-d1ji;
- print di,dj;
- /*False Conflict*/
- FC = j(nrow(m), nrow(m));
- do i = 1 to nrow(m);
- do j=1 to nrow(m);
- if di[i,j] ^=missing(di[i,j]) | dj[i,j] ^= missing(dj[i,j]) then
- do;
- if di[i,j]>0 & dj[i,j]>0 then FC[i,j]= 1;
- else FC[i,j] =0;
- end;
- else
- do;
- FC[i,j] =0;
- end;
- end;
- end;
- print FC;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement