/// <summary>
/// 获得射线与地形绕射棱交点、与源点距离、所在面的信息
/// </summary>
/// <param name="rays">射线</param>
/// <param name="terShadowRects">射线所经过的地形矩形</param>
/// <param name="txFrequencyBand">发射机频带信息</param>
/// <param name="multiple">绕射半径系数</param>
/// <returns>返回射线与地形绕射棱交点、与源点距离、所在面的信息</returns>
protected List <Node> GetCrossNodeWithTerDiffractionEdge(List <FrequencyBand> txFrequencyBand, double multiple)
{
List <Node> crossNodes = new List <Node>();
//获取最大波长
double maxWaveLength = 300.0 / TxFileProceed.GetMinFrequenceFromList(txFrequencyBand); //300为光速
//获取最小波长
double minWaveLength = 300.0 / TxFileProceed.GetMaxFrequenceFromList(txFrequencyBand); //300为光速
for (int i = 0; i < this.terShadowRects.Count; i++) //对经过的每个地形矩形进行遍历
{
for (int j = 0; j < this.terShadowRects[i].RectTriangles.Count; j++) //对每个矩形的三角面进行遍历
{
for (int k = 0; k < this.terShadowRects[i].RectTriangles[j].Lines.Count; k++) //对每个三角面的边进行遍历
{
if (!this.terShadowRects[i].RectTriangles[j].Lines[k].IsDiffractionEdge) //若该条边是一条绕射棱
{
continue;
}
else
{
if (this.terShadowRects[i].RectTriangles[j].Lines[k].HaveTraced)//若该条绕射棱已经被追踪
{
this.terShadowRects[i].RectTriangles[j].Lines[k].HaveTraced = false;
continue;
}
else
{
this.terShadowRects[i].RectTriangles[j].Lines[k].HaveTraced = true;
//判断凹凸性
if (this.inRay.Origin.JudgeIsConcaveOrConvexToViewPoint(this.terShadowRects[i].RectTriangles[j].Lines[k].AdjacentTriangles[0], this.terShadowRects[i].RectTriangles[j].Lines[k].AdjacentTriangles[1]))
{
Node paramNode = this.inRay.GetCrossNodeWithCylinder(this.terShadowRects[i].RectTriangles[j].Lines[k], multiple * maxWaveLength);
if (paramNode != null)//若存在交点
{
//当发射点与绕射棱在一个三角面内,或者发射点在绕射圆柱体内时，该绕射点不算
if (this.terShadowRects[i].RectTriangles[j].Lines[k].getDistanceP2Line(this.inRay.Origin) <= multiple * maxWaveLength ||
this.terShadowRects[i].RectTriangles[j].Lines[k].AdjacentTriangles[0].JudgeIfPointInFace(this.inRay.Origin) ||
this.terShadowRects[i].RectTriangles[j].Lines[k].AdjacentTriangles[1].JudgeIfPointInFace(this.inRay.Origin))
{
continue;
}
else
{
crossNodes.Add(paramNode);
}
}
}
}
}
}
}
}
this.HandleDiffractionCrossNode(crossNodes, minWaveLength);
return(crossNodes);
}
public static void Calculate(string setuppath, string terpath, string txpath, string rxpath)
{
//此判断用来将ter文件中的三角面读入Intersection.Tris中
Terrain newTer = new Terrain(terpath);
string DirectoryPath = GetDirectory(setuppath) + "result\\";
SetupObject setupOne = SetupFileProceed.GetSetupFile.GetSetup(setuppath);
int txTotal = setupOne.tr.FirstAvailableTxNumber;
string ProName = GetProName(setuppath);
List <Node> txs = TxFileProceed.GetTx(txpath, setuppath, terpath);
myLog.Info(txs.Count);
List <List <ReceiveBall> > rxs = RxFileProceed.GetRx(rxpath, setuppath, terpath);
myLog.Info(rxs.Count);
//rxs.AddRange(TxFileProceed.GetTxAsRx(txpath, setuppath, terpath));
City buildings = new City(GetDirectory(setuppath) + ProName + ".city");
buildings.RestructTerrainByBuildings(newTer.TerRect);
if (buildings.Build != null && buildings.Build.Count != 0)
{
newTer.OutPutNewTerrain(GetDirectory(setuppath));
} //输出新地形
int?[] index = { null, null };
int txindex = 0;
for (int i = 0; i < txs.Count; i++)
{
List <FrequencyBand> txFrequencyBand = TxFileProceed.GetTxFrequenceBand(txpath, setuppath, i);//获取频段信息
myLog.Debug("这是第" + (i + 1) + "个发射机的射线追踪==========================================================================");
for (int j = 0; j < rxs.Count; j++)
{
//该接收区只有一个接收机
if (rxs[j].Count == 1)
{
if (rxs[j][0].GetType().ToString().Equals("CalculateModelClasses.ReceiveArea"))
{
//电磁态势的追踪模块
ReceiveArea reArea = (ReceiveArea)rxs[j][0];
//第一步：构建态势区域
myLog.Debug("构建态势区域");
reArea.CreateAreaSituation(newTer.TerRect);
//第二步：正向追踪，获取路径
myLog.Debug("正向追踪，获取粗略路径");
RayTubeMethod areaTracing = new RayTubeMethod(txs[i], reArea, newTer, buildings, 64);
//读取发射天线的极化信息
string txPol = GetRxTxPol("tx", txs[i].TxNum);
//读取接收天线的极化信息（删除？？）
string rxPol = GetRxTxPol("areaSituation", rxs[j][0].RxNum);
//第三步：反向追踪，获取准确路径
//for (int m = 0; m < reArea.areaSituationNodes.Count; m++)
//{
// areaTracing.ReverseAreaTracingPathsAndDeleteRepeatedPaths(reArea.areaSituationNodes[m].paths);
//}
//获取发射天线的4个参数
ReadUan.GetGainPara(txs[i].UAN);
for (int m = 0; m < reArea.areaSituationNodes.Count; m++)
{
reArea.areaSituationNodes[m].classifiedFrequencyPaths = new List <List <CalculateModelClasses.Path> >();
if (reArea.areaSituationNodes[m].paths.Count != 0)
{
PathsafterPolization(reArea.areaSituationNodes[m].paths, txPol, rxPol);//极化代码
//第四步：分频段求每条路径上的场强
reArea.areaSituationNodes[m].classifiedFrequencyPaths = new List <List <CalculateModelClasses.Path> >();
//areaTracing.ScreenAreaSituationPathsByFrequencyAndCalculateEField(txFrequencyBand, reArea.areaSituationNodes[m].paths, reArea.areaSituationNodes[m].classifiedFrequencyPaths);
//计算每个频段上的各个路径叠加的总场强
for (int n = 0; n < txFrequencyBand.Count; n++)
{
EField tempEfield = new EField();
tempEfield.X = tempEfield.GetTolEx(reArea.areaSituationNodes[m].classifiedFrequencyPaths[n]);
tempEfield.Y = tempEfield.GetTolEy(reArea.areaSituationNodes[m].classifiedFrequencyPaths[n]);
tempEfield.Z = tempEfield.GetTolEz(reArea.areaSituationNodes[m].classifiedFrequencyPaths[n]);
//得到每个频段上各个路径叠加的场强
reArea.areaSituationNodes[m].totleEfields.Add(tempEfield);
//得到各个频段上各个路径后得到的总场强
reArea.areaSituationNodes[m].totleEfield.X += tempEfield.X;
reArea.areaSituationNodes[m].totleEfield.Y += tempEfield.Y;
reArea.areaSituationNodes[m].totleEfield.Z += tempEfield.Z;
}
}
}
//文件输出
string path = ".\\.\\project\\station\\areaSituationResult.txt";
if (!System.IO.File.Exists(path))
{
System.IO.File.Create(path).Close();
}
StringBuilder sb = new StringBuilder();
for (int m = 0; m < reArea.areaSituationNodes.Count; m++)
{
if (reArea.areaSituationNodes[m].paths.Count != 0)
{
//string appendText1 = m+" " + reArea.areaSituationNodes[m].paths.Count+" " ;
//sb.Append(appendText1);
//string appendText2 = "Ex " + reArea.areaSituationNodes[m].totleEfield.X.Re + "+j" + reArea.areaSituationNodes[m].totleEfield.X.Im + ",Ey " + reArea.areaSituationNodes[m].totleEfield.Y.Re + "+j" + reArea.areaSituationNodes[m].totleEfield.Y.Im + ",Ez " + reArea.areaSituationNodes[m].totleEfield.Z.Re + "+j" + reArea.areaSituationNodes[m].totleEfield.Z.Im + "\r\n";
//sb.Append(appendText2);
string appendText1 = "第" + m + "个态势点上有" + reArea.areaSituationNodes[m].paths.Count + "条路径，";
sb.Append(appendText1);
string appendText2 = "总场强为:Ex " + reArea.areaSituationNodes[m].totleEfield.X.Re + "+j" + reArea.areaSituationNodes[m].totleEfield.X.Im + ",Ey " + reArea.areaSituationNodes[m].totleEfield.Y.Re + "+j" + reArea.areaSituationNodes[m].totleEfield.Y.Im + ",Ez " + reArea.areaSituationNodes[m].totleEfield.Z.Re + "+j" + reArea.areaSituationNodes[m].totleEfield.Z.Im + "\r\n";
sb.Append(appendText2);
}
//else
//{
// string appendText1 = "没有到达第" + m + "个态势点的路径\r\n";
// sb.Append(appendText1);
//}
}
StreamWriter stw = new StreamWriter(path);
stw.Write(sb);
stw.Flush(); //清空缓冲区
stw.Close(); //关闭流
//---------------------------------
myLog.Debug("结果计算和文件输出过程结束，进入下一个接收机的追踪--------------------------------------");
}
else
{
if (("tx_" + txs[i].NodeName).Equals(rxs[j][0].RxName))
{
myLog.Debug("这是第" + (i + 1) + "个发射机第" + (j + 1) + "个接收机的射线追踪:是同一个发射机");
continue;
}
myLog.Debug("这是第" + (i + 1) + "个发射机第" + (j + 1) + "个接收机的射线追踪*****************************************************");
index[0] = j;
RayTubeMethod rayTubeMethod = new RayTubeMethod(txs[i], rxs[j][0], newTer, buildings, 32);
rayTubeMethod.ReverseTracingPathsAndDeleteRepeatedPaths();
rayTubeMethod.UpdateRayInForNodes();
// PunctiformLaunchMethod punctiformMethod = new PunctiformLaunchMethod(txs[i], rxs[j][0], newTer, buildings, 64, txFrequencyBand);
// List<CalculateModelClasses.Path> temp = punctiformMethod.GetPunctiformRxPath(tx, rxs[j][0], newTer,cityBuilding, 32, 128,TxFrequencyBand);
// if (rayTubeMethod.ReceivedPaths.Count != 0)//若存在直射射线
// {
// string txPol = GetRxTxPol("tx", txs[i].TxNum);//极化代码
// string rxPol = null;
// if (rxs[j][0].isTx == true)
// {
// rxPol = GetRxTxPol("tx", rxs[j][0].RxNum);
// }
// else
// {
// rxPol = GetRxTxPol("rx", rxs[j][0].RxNum);//
// }
// PathsafterPolization(rayTubeMethod.ReceivedPaths, txPol, rxPol);//
// }
myLog.Debug("射线追踪过程结束，进入射线筛选和计算过程-----------------------------------------");
//输出结果
ReadUan.GetGainPara(rxs[j][0].UAN);
List <List <CalculateModelClasses.Path> > ClassifiedPaths = rayTubeMethod.ScreenPunctiformPathsByFrequencyAndCalculateEField(txFrequencyBand);
for (int m = 0; m < txFrequencyBand.Count; m++)
{
P2mFileOutput.p2mfileoutput(ClassifiedPaths[m], DirectoryPath, ProName, txs[i].TxNum, txTotal, txindex, rxs[j][0], txs[i].Position, txFrequencyBand[m].MidPointFrequence, txFrequencyBand[m].FrequenceWidthMin, txFrequencyBand[m].FrequenceWidthMax);
}
List <CalculateModelClasses.Path> totalPowerPath = GetTotalPowerPath(ClassifiedPaths);
PowerOutput.GetRxTotalPower(totalPowerPath, DirectoryPath, ProName, txs[i].TxNum, txTotal, rxs[j][0], txs[i].Position);
myLog.Debug("结果计算和文件输出过程结束，进入下一个接收机的追踪---------------------------------------");
}
}
else
{
throw new Exception("没有接收机");
}
}
}
GC.Collect();
}
