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
| inf2 = function (matrix)
{
n = dim(matrix)[2]
angle = NULL
for(i in 2: (n-1))
{
aa = paste("a",i,sep="")
y = assign(aa,matrix[1,i],pos=1)
bb = paste("b",i,sep="")
x = assign(bb,i,pos=1)
cc = paste("c",i-1,sep="")
yy = assign(cc,matrix[1,i-1],pos=1)
dd = paste("d",i-1,sep="")
xx = assign(dd,i-1,pos=1)
ee = paste("e",i+1,sep="")
yyy= assign(ee,matrix[1,i+1],pos=1)
ff = paste("f",i+1,sep="")
xxx= assign(ff,i+1,pos=1)
angle = c(angle,assign(aa,((atan2(y-yy,x-xx)+atan2(yyy-y,xxx-x))/2),pos=1))
}
return(angle)
} |
Partager