func (dbs *DBS) optimizePowerAllocation() {
var meanPtot, meanPd, meanPtotPd, meanPePd, meanPr, meanPe float64
//var meanBER float64
for e := dbs.Connec.Front(); e != nil; e = e.Next() {
c := e.Value.(*Connection)
M := c.E
meanPtotPd += M.BERT() / M.Req()
meanPePd += c.BER / M.Req()
meanPtot += M.BERT()
meanPe += c.BER
meanPd += M.Req()
meanPr += c.Pr
}
nbconnec := float64(dbs.Connec.Len())
meanPtot /= nbconnec
meanPtotPd /= nbconnec
meanPePd /= nbconnec
meanPe /= nbconnec
meanPd /= nbconnec
meanPr /= nbconnec
for e := dbs.Connec.Front(); e != nil; e = e.Next() {
c := e.Value.(*Connection)
M := c.E
if c.Status == 0 && c.E.GetCh() != 0 { // if master connection
var b, need, delta, alpha float64
b = M.BERT() / M.Req()
//b = M.BERT() / M.Req()/meanPtotPd;
b = b * 1.5
alpha = 1.0
need = 2.0*math.Exp(-b)*(b+1.0) - 1.0
delta = math.Pow(geom.Abs(need), 1) *
math.Pow(M.GetPower(), 1) *
geom.Sign(need-M.GetPower()) * alpha *
math.Pow(geom.Abs(need-M.GetPower()), 1.5)
if delta > 0 {
v := (1.0 - M.GetPower()) / 2.0
if delta > v {
delta = v
}
} else {
v := -M.GetPower() / 2.0
if delta < v {
delta = v
}
}
if delta > 1 || delta < -1 {
fmt.Println("Power Error ", delta)
}
M.PowerDelta(delta)
//M.SetPower(1.0/math.Pow(2.0+dbs.R.Pos.Distance(M.GetPos()),4 ) )
}
}
}
func DrawReceptionField(dbs []synstation.DBS, name string) {
var data [Size][Size]float64
//dbs:= make([]synstation.DBS,1)
//dbs[0].R.SetPos(geom.Pos{6000,6000})
e := new(synstation.Emitter)
e.Power = 1
var TMax, TMin float64
TMax = -300
TMin = +300
ch := synstation.NChRes
sums := float64(0)
for x := 0; x < Size; x++ {
e.Pos.X = inField(x)
sumy := float64(0)
for y := 0; y < Size; y++ {
e.Pos.Y = inField(y)
Pr := float64(0.0)
for k := 0; k < len(dbs); k++ {
/* p:= e.GetPos().Minus(*dbs[k].R.GetPos())
theta := math.Atan2(p.Y,p.X) * 180/ math.Pi
if theta<0 {theta=theta+360}
if theta <0 {fmt.Println(" theta <0!!!!")}
dbs[k].R.Orientation[ch]=theta //+ (dbs.Rgen.Float64()*30-15)
*/
if geom.Abs(inField(x)-dbs[k].R.GetPos().X) < 1500 &&
geom.Abs(inField(y)-dbs[k].R.GetPos().Y) < 1500 {
p, _ := dbs[k].R.EvalSignalPr(e, ch)
if p > Pr {
Pr = p
}
}
}
Pr = 10 * math.Log10(Pr)
data[x][y] = Pr
sumy += Pr
if Pr > TMax {
TMax = Pr
}
if Pr < TMin {
TMin = Pr
}
}
sums += sumy
}
fmt.Println("mean ", sums/600/600)
fmt.Println("TMax ", TMax)
fmt.Println("TMin ", TMin)
if TMin < -150 {
TMin = -150
}
// -15 -125
// -5 -115
// -9 -117
im := image.NewNRGBA(Size, Size)
for x := 0; x < Size; x++ {
for y := 0; y < Size; y++ {
var r, g, b float64
v := (data[x][y] - TMin) / (TMax - TMin) * -8.0
switch {
case v > -2.0:
b = 2.0 + v
g = -v
case v > -4:
g = 1
r = -v/2.0 - 1.0
default:
r = 1.0
b = -v/2.0 - 2.0
g = 1.0 - b
}
r = r * 255
g = g * 255
b = b * 255
if r < 0 {
r = 0
}
if r > 255 {
r = 255
}
if g < 0 {
g = 0
}
if g > 255 {
g = 255
}
if b < 0 {
b = 0
}
if b > 255 {
b = 255
}
im.Pix[y*im.Stride+x] = uint8(r) //image.NRGBAColor{uint8(r), uint8(g), uint8(b), uint8(255)}
//im.Pix[y*im.Stride+x] = image.NRGBAColor{uint8(v*255),uint8(v*255),uint8(v*255),255}
}
}
f, err := os.Create(name)
if err = png.Encode(f, im); err != nil {
os.Exit(1)
}
}