public IIteration DoIteration(IIteration prevIter) { it = (HJIteration)prevIter; double f = ip.GetFuncValue(it.arrX); while (it.arrE == null) { f = ip.GetFuncValue(it.arrX); double[] arrENext, arrXNext; double fNext = Research(it.arrX, out arrENext, out arrXNext); if (Math.Abs(fNext - f) < ip.fEpsilon) { for (int i = 0; i < it.arrXDelta.Length; i++) { it.arrXDelta[i] /= 2; } for (int i = 0; i < it.arrXDelta.Length; i++) { if (it.arrXDelta[i] < ip.xEpsilon) { it.CalcResult(ip); return(null); } } } else { it.arrE = arrENext; } } double[] arrXSmp = new double[it.arrX.Length]; for (int i = 0; i < arrXSmp.Length; i++) { arrXSmp[i] = it.arrX[i] + it.arrE[i] * it.mult; } double fSmp = ip.GetFuncValue(arrXSmp); if (Math.Abs(fSmp - f) < ip.fEpsilon || fSmp > f) { it.fRes = f; return(new HJIteration(it.arrX, it.arrXDelta, null, 1)); } double[] arrXRes, arrERes; double fRes = Research(arrXSmp, out arrERes, out arrXRes); HJIteration itNext; if (Math.Abs(fRes - f) < ip.fEpsilon) { it.fRes = f; itNext = new HJIteration(it.arrX, it.arrXDelta, null, 1); } else { it.fRes = fSmp; itNext = new HJIteration(arrXSmp, it.arrXDelta, it.arrE, it.mult + 1); } return(itNext); }
void hJToolStripMenuItem_Click(object sender, EventArgs e) { try { HJIteration it; ReadHJParams(out ipHJ, out it); HJOptimizer opt = new HJOptimizer(); opt.Initialize(ipHJ); listHJ = new List <HJIteration>(); it.CalcResult(ipHJ); listHJ.Add(it); it = new HJIteration(it.arrX, it.arrXDelta); it.CalcResult(ipHJ); int iterNum = 0; do { listHJ.Add(it); it = (HJIteration)opt.DoIteration(it); iterNum++; }while (it != null && iterNum < maxIt); lbIterHJ.Items.Clear(); lbIterHJ.Items.AddRange(listHJ.ToArray()); psGraph.Clear(); psItGraph.Clear(); if (listHJ[0].arrX.Length != 2) { return; } LinePlot[] arrLp = GetPotentials(ipHJ); foreach (LinePlot lp in arrLp) { psGraph.Add(lp); psItGraph.Add(lp); } try { DrawHJTrajectory(); } catch { } } catch (Exception ex) { MessageBox.Show("Ошибка: " + ex.Message, "Расчеты прерваны", MessageBoxButtons.OK, MessageBoxIcon.Error); } }
static void Main() { Application.EnableVisualStyles(); Application.SetCompatibleTextRenderingDefault(false); Application.Run(new MainForm()); return; NMOptimizer opt = new NMOptimizer(); NMInitialParams ip = new NMInitialParams(1, 0.5, 2, 0.001, new string[] { "x", "y" }, new string[] { "x", "y" }, "x ^ 2 + y ^ 2"); opt.Initialize(ip); List <NMIteration> listIter = new List <NMIteration>(); NMIteration it = new NMIteration(new double[][] { new double[] { 1, 2 }, new double[] { 5, 6 }, new double[] { 7, 10 } }); it.CalcFuncAndResult(ip); do { listIter.Add(it); it = (NMIteration)opt.DoIteration(it); }while (it != null); HJOptimizer opt2 = new HJOptimizer(); HJInitialParams ip2 = new HJInitialParams(0.01, 0.1, new string[] { "x", "y" }, new string[] { "x", "y" }, "x * x + 2 * x + 1 + y * y + 2 * y + 1"); opt2.Initialize(ip2); List <HJIteration> listIter2 = new List <HJIteration>(); HJIteration it2 = new HJIteration(new double[] { 20, 20 }, new double[] { 0.1, 0.1 }, null, 1); it2.CalcResult(ip2); do { listIter2.Add(it2); it2 = (HJIteration)opt2.DoIteration(it2); }while (it2 != null); }