public void BeeDown(double min, double max, int count, int t, int maxmax = 100) { Vectors old = Vectors.Union2(cnew.SubVector(t), d.SubVector(t)); //slau.GaussSelection(t); //cnew= var r = BeeHiveAlgorithm.GetGlobalMin(this.F, (t) * 2, min, max, 1e-15, count, 60, old, maxmax); //$"Hive method ---> {r.Item2}".Show(); VecToCD(r.Item1); ErrorsMasL[t - 1] = r.Item2; Functional ff = (Point x) => { double s = 0; for (int i = 0; i < t; i++) { s += cnew[i] * alpha(x, i) + d[i] * beta(x, i); } return(Math.Abs(s - KursMethods.U(x))); }; ErrorsMasQ[t - 1] = IntegralClass.Integral(ff, CIRCLE - 1); }
private Tuple <Vectors, double, double> HalfSolve(int len) { SqMatrix iD1 = D1.SubMatrix(len), iD2 = D2.SubMatrix(len), iS = S.SubMatrix(len); Vectors iR1 = R1.SubVector(len), iR2 = R2.SubVector(len); SqMatrix iD1m = D1m.SubMatrix(len), iD2m = D2m.SubMatrix(len), idl; Vectors idr; idl = (D2 - S * (D1m * S)).SubMatrix(len); idr = (R2 + S * D1m * R1).SubVector(len); Vectors id = idl.Solve(idr), icnew = (D1m * (R1 + S * id)).SubVector(len); Vectors res = Vectors.Union2(icnew, id); double p1 = IntegralClass.IntegralCurve((Point y) => (TestFuncAndCurve.Grads[GF - 1](y) * TestFuncAndCurve.Norm[CIRCLE - 1](y)).Sqr(), CIRCLE - 1); double p2 = IntegralClass.IntegralCurve((Point y) => (U(y)).Sqr(), CIRCLE - 1); Functional fif = (Point x) => { double s = 0; for (int i = 0; i < id.Deg; i++) { s += icnew[i] * alpha(x, i) + id[i] * beta(x, i); } return(Math.Abs(s - KursMethods.U(x))); }; double L = IntegralClass.Integral(fif, CIRCLE - 1); $"len = {len}, L = {this.F(res)}".Show(); return(new Tuple <Vectors, double, double>(res, F(res), L)); }
private void button2_Click(object sender, EventArgs e) { textBox1.Text += "Процесс запущен" + Environment.NewLine; int k = Convert.ToInt32(Program.Form1.textBox7.Text); if (Program.Form1.radioButton1.Checked) { int n = Int32.Parse(Program.Form1.textBox3.Text); int cu = Int32.Parse(Program.Form1.textBox1.Text); int gf = Int32.Parse(Program.Form1.textBox2.Text); KursMethods.Desigion(n, gf, cu); textBox1.Text += "Данные считаны" + Environment.NewLine; if (Program.Form1.checkBox1.Checked) { KursMethods.bstr = ""; KursMethods.sl = ""; KursMethods.Illustrating(); Program.FORM.chart1.Refresh(); } if (Program.Form1.checkBox2.Checked) { KursMethods.bstr = ""; KursMethods.sl = ""; KursMethods.Fixity(); Program.FORM.chart1.Refresh(); } if (Program.Form1.checkBox3.Checked) { KursMethods.bstr = ""; KursMethods.sl = ""; KursMethods.Quality(n, k, gf, cu); Program.FORM.chart1.Refresh(); } } if (Program.Form1.radioButton2.Checked) { int a = Int32.Parse(Program.Form1.textBox4.Text); int b = Int32.Parse(Program.Form1.textBox5.Text); int h = Int32.Parse(Program.Form1.textBox6.Text); textBox1.Text += "Данные считаны" + Environment.NewLine; if (Program.Form1.checkBox6.Checked) { KursMethods.Pictures_ill(a, b, h); //графики приближения для 4-10 функций с шагом 3 } if (Program.Form1.checkBox5.Checked) { KursMethods.Pictures_fix(a, b, h); //картинки зависимости погрешности аппроксимации для от 30 до 40 функций, шаг 30 } if (Program.Form1.checkBox4.Checked) { KursMethods.Pictures_qua(k, a, b, h); //картинки зависимости погрешности от кривой для 20 функций с кривыми от 40 до 100 и шагом 20 } } }
/// <summary> /// Вычисление криволинейного интеграла первого рода по этой кривой от функции BasisFuncPow(int i,int j,Point z) /// </summary> /// <param name="i"></param> /// <param name="j"></param> /// <returns></returns> public double Firstkind(int i, int j) { //return this.Firstkind((Point x)=>KursMethods.BasisFuncPow(i,j,(Point)x)); Functional f = (Point x) => KursMethods.BasisFuncPow(i, j, (Point)x); Func <double, double> func = (double t) => f(this.Transfer(t)); return(IntegralClass.Integral(func, this.a, this.b)); }
public void UltraHybrid(int t) { //UltraCount.Show(); if (UltraCount == 0)//если вообще не решалось { //"Вошло".Show(); x[0] = b[0] / A[0, 0]; //x[0].Show(); VALUE_FOR_ULTRA = Error(1); ErrorsMas[0] = VALUE_FOR_ULTRA; Functional f = (Point x) => { double sum = this.x[0] * KursMethods.masPoints[0].PotentialF((Point)x); double s = sum - KursMethods.fig(x); return(s * s); }; Functional fM1 = (Point x) => { //это я ещё не использовал!!! //Functional M1ApproxD = (Point u) => //{ // Functional ro = (Point y) => // { // double sss = this.x[0] * KursMethods.masPoints[0].PotentialF(u); // return sss * KursMethods.Exy(u, y); // }; // return DoubleIntegral(ro, curve, curve.S, FuncMethods.DefInteg.Method.GaussKronrod15, 0.001, FuncMethods.DefInteg.countY) /*IntegralClass.Integral(ro, KursMethods.CIRCLE - 1)*/ /*- right(x)*/; //}; double su = KursMethods.M1Approx(x); double s = su - KursMethods.TestFuncAndCurve.DFunctions[KursMethods.GF - 1](x); return(s * s); }; ErrorMasP[0] = Math.Sqrt(DoubleIntegral(f, curve, curve.S, FuncMethods.DefInteg.Method.GaussKronrod15, 0.001, FuncMethods.DefInteg.countY)); ErrorsM1[0] = Math.Sqrt(IntegralClass.Integral2(fM1, KursMethods.CIRCLE - 1)); $"{1} = {ErrorsM1[0]}".Show(); UltraCount++; for (int i = UltraCount + 1; i <= t; i++) { UltraHybridLast(i); } UltraCount = t; } else if (UltraCount == t - 1)//если надо решить только по последней координате { UltraHybridLast(t); UltraCount++; } else { for (int i = UltraCount + 1; i <= t; i++) { UltraHybridLast(i); } UltraCount = t; } }
public void UltraHybrid(int t) { Functional f = (Point x) => { double sum = 0; for (int i = 0; i < dim; i++) { sum += c[i] * alpha(x, i) + d[i] * beta(x, i); } return(Math.Abs(sum - KursMethods.U(x))); }; if (UltraCount == 0)//если вообще не решалось { ErrorsMasQ = new double[dim]; ErrorsMasL = new double[dim]; //"Вошло".Show(); d[0] = dr[0] / dl[0, 0]; //x[0].Show(); VALUE_FOR_ULTRA = Error; ErrorsMasL[0] = VALUE_FOR_ULTRA; ErrorsMasQ[0] = IntegralClass.Integral(f, CIRCLE - 1); //BeeDown(-20, 20, 1000, 1); UltraCount++; for (int i = UltraCount + 1; i <= t; i++) { UltraHybridLast(i); } UltraCount = t; } else if (UltraCount == t - 1)//если надо решить только по последней координате { UltraHybridLast(t); UltraCount++; } else { for (int i = UltraCount + 1; i <= t; i++) { UltraHybridLast(i); } UltraCount = t; } }
public static Tuple <double[], double[], int[], Func <Point, double> > LastMethod2(int n, int g, int cu, SLAUpok SYSTEM = null, SLAUpok SYSTEMQ = null) { ForDesigion.Building(n, g, cu, SYSTEM, SYSTEMQ); double[] ErrorsMasL = new double[n], ErrorsMasQ = new double[n]; int[] Times = new int[n]; var tmp = new BiharmonicEquation(); double range = 200; Vectors old = Vectors.Union2(tmp.cnew, tmp.d); var r = BeeHiveAlgorithm.GetGlobalMin(tmp.F, 2 * n, -range, range, 1e-15, 1200, 60, old, -1); for (int t = 0; t < n; t++) { Functional fif = (Point x) => { double s = 0; for (int i = 0; i < t; i++) { s += tmp.cnew[i] * alpha(x, i) + tmp.d[i] * beta(x, i); } return(Math.Abs(s - KursMethods.U(x))); }; ErrorsMasQ[t] = IntegralClass.Integral(fif, CIRCLE - 1); ErrorsMasL[t] = tmp.F(Vectors.Union2(tmp.cnew.SubVector(t), tmp.d.SubVector(t))); } Func <Point, double> ff = (Point x) => { double s = 0; for (int i = 0; i < tmp.dim; i++) { s += tmp.cnew[i] * alpha(x, i) + tmp.d[i] * beta(x, i); } return(s); }; return(new Tuple <double[], double[], int[], Func <Point, double> >(ErrorsMasL, ErrorsMasQ, Times, ff)); }
/// <summary> /// Ультра-гибридный метод суперского решения по последней координате /// </summary> /// <param name="t"></param> public void UltraHybridLast(int t) //гибридный с координатной минимизацией по последней координате { double[] c = new double[t], cn = new double[t]; //cnew = this.c.dup; for (int i = 0; i < t - 1; i++) { c[i] = slau.x[i]; cn[i] = cnew[i]; } Vectors mk1 = new Vectors(c), mk2 = new Vectors(c); double sum = 0; slau.GaussSpeedy(t); cnew = this.c.dup; double tmp = Error; if (VALUE_FOR_ULTRA < tmp) //если погрешность выросла - исправить это, потому что новое решение не годится { $"{VALUE_FOR_ULTRA} < {tmp} при t = {t} (до покоординатной минимизации результата СПИДГАУССА)".Show(); //покоординатная минимизация результата СПИДГАУССА MINIMAKA(t); tmp = Error; if (VALUE_FOR_ULTRA < tmp) { $"{VALUE_FOR_ULTRA} < {tmp} при t = {t} (до полной покоординатной минимизации вектора с1 с2 ... 0)".Show(); for (int i = 0; i < t; i++)//исправили, теперь пробуем новый метод { slau.x[i] = c[i]; cnew[i] = cn[i]; } //покоординатная минимизация MINIMAKA(t); double tmp1 = Error; if (VALUE_FOR_ULTRA < tmp1) //погрешность опять выросла - тогда просто оставляем 0 на конце { $"{VALUE_FOR_ULTRA} < {tmp1} при t = {t} (до полной покоординатной минимизации на конце)".Show(); for (int i = 0; i < t; i++)//исправили, теперь пробуем новый метод { slau.x[i] = c[i]; cnew[i] = cn[i]; } for (int j = 0; j < t - 1; j++) { sum += c[j] * S[t - 1, j]; sum += d[j] * D2[t - 1, j]; } sum += c[t - 1] * S[t - 1, t - 1]; slau.x[t - 1] = (R2[t - 1] - sum) / slau.A[t - 1, t - 1]; sum = 0; tmp = Error; if (VALUE_FOR_ULTRA < tmp) { for (int i = 0; i < t - 1; i++)//исправили, теперь пробуем новый метод { slau.x[i] = c[i]; cnew[i] = cn[i]; } slau.x[t - 1] = 0; cnew[t - 1] = 0; } else { $"Погрешность уменьшена МИНИМАКОЙ НА КОНЦЕ на {(VALUE_FOR_ULTRA - tmp) / VALUE_FOR_ULTRA * 100} %".Show(); VALUE_FOR_ULTRA = tmp; } } else { $"Погрешность уменьшена ПОЛНОЙ МИНИМАКОЙ на {(VALUE_FOR_ULTRA - tmp1) / VALUE_FOR_ULTRA * 100} %".Show(); VALUE_FOR_ULTRA = tmp1; } } else { $"Погрешность уменьшена МИНИМАКОЙ СПИДГАУССА на {(VALUE_FOR_ULTRA - tmp) / VALUE_FOR_ULTRA * 100} %".Show(); VALUE_FOR_ULTRA = tmp; } } else { $"Погрешность уменьшена СПИДГАУССОМ на {(VALUE_FOR_ULTRA - tmp) / VALUE_FOR_ULTRA * 100} %".Show(); VALUE_FOR_ULTRA = tmp; } ErrorsMasL[t - 1] = VALUE_FOR_ULTRA; Functional f = (Point x) => { double s = 0; for (int i = 0; i < t; i++) { s += c[i] * alpha(x, i) + d[i] * beta(x, i); } return(Math.Abs(s - KursMethods.U(x))); }; ErrorsMasQ[t - 1] = IntegralClass.Integral(f, CIRCLE - 1); }
private async void button2_Click(object sender, EventArgs e) { textBox1.Text += "Процесс запущен" + Environment.NewLine; int k = Convert.ToInt32(Program.Form1.textBox7.Text); //BiharmonicEquation.Test(); //Tables(); if (Program.Form1.radioButton1.Checked) { int n = Int32.Parse(Program.Form1.textBox3.Text); int cu = Int32.Parse(Program.Form1.textBox1.Text); int gf = Int32.Parse(Program.Form1.textBox2.Text); //FonPotok(new ThreadStart(() => KursMethods.Desigion(n, gf, cu))); textBox1.Text += "Данные считаны" + Environment.NewLine; if (Program.Form1.checkBox1.Checked) { KursMethods.Desigion(n, gf, cu); KursMethods.bstr = ""; KursMethods.sl = ""; /*await Task.Run(()=>*/ KursMethods.Illustrating() /*)*/; Program.FORM.chart1.Refresh(); } if (Program.Form1.checkBox2.Checked) { KursMethods.Desigion(n, gf, cu); KursMethods.bstr = ""; KursMethods.sl = ""; /*await Task.Run(()=>*/ KursMethods.Fixity() /*)*/; Program.FORM.chart1.Refresh(); } if (Program.Form1.checkBox3.Checked) { KursMethods.Desigion(n, gf, cu); KursMethods.bstr = ""; KursMethods.sl = ""; /*await Task.Run(()=> */ KursMethods.Quality(n, k, gf, cu) /*)*/; Program.FORM.chart1.Refresh(); } if (Program.Form1.checkBox7.Checked) { KursMethods.Desigion(n, gf, cu); KursMethods.bstr = ""; KursMethods.sl = ""; KursMethods.LnGraf(n, gf, cu); } if (Program.Form1.checkBox9.Checked) { KursMethods.bstr = ""; KursMethods.sl = ""; double[] Ls = Program.Form1.textBox12.Text.ToDoubleMas(); BigarmRes(n, gf, cu, Ls); } } if (Program.Form1.radioButton2.Checked) { int a = Int32.Parse(Program.Form1.textBox4.Text); int b = Int32.Parse(Program.Form1.textBox5.Text); int h = Int32.Parse(Program.Form1.textBox6.Text); textBox1.Text += "Данные считаны" + Environment.NewLine; if (Program.Form1.checkBox6.Checked) /*await Task.Run(() => */ KursMethods {
/// <summary> /// Ультра-гибридный метод суперского решения по последней координате /// </summary> /// <param name="t"></param> public void UltraHybridLast(int t) //гибридный с координатной минимизацией по последней координате { double[] c = new double[t]; for (int i = 0; i < t - 1; i++) { c[i] = x[i]; } Vectors mk1 = new Vectors(c), mk2 = new Vectors(c); bool existres = true; double sum = 0; GaussSpeedy(t); double tmp = Error(t); if (VALUE_FOR_ULTRA < tmp) //если погрешность выросла - исправить это, потому что новое решение не годится { $"{VALUE_FOR_ULTRA} < {tmp} при t = {t} (до покоординатной минимизации результата СПИДГАУССА)".Show(); //покоординатная минимизация результата СПИДГАУССА for (int k = 0; k <= t - 1; k++) { for (int j = 0; j < k; j++) { sum += x[j] * A[k, j]; } for (int j = k + 1; j < t - 1; j++) { sum += x[j] * A[k, j]; } x[k] = (b[k] - sum) / A[k, k]; sum = 0; } tmp = Error(t); if (VALUE_FOR_ULTRA < tmp) { $"{VALUE_FOR_ULTRA} < {tmp} при t = {t} (до полно покоординатной минимизации вектора с1 с2 ... 0)".Show(); for (int i = 0; i < t; i++)//исправили, теперь пробуем новый метод { x[i] = c[i]; } //покоординатная минимизация for (int k = 0; k <= t - 1; k++) { for (int j = 0; j < k; j++) { sum += x[j] * A[k, j]; } for (int j = k + 1; j < t - 1; j++) { sum += x[j] * A[k, j]; } x[k] = (b[k] - sum) / A[k, k]; sum = 0; } double tmp1 = Error(t); if (VALUE_FOR_ULTRA < tmp1) //погрешность опять выросла - тогда просто оставляем 0 на конце { $"{VALUE_FOR_ULTRA} < {tmp1} при t = {t} (до полной покоординатной минимизации на конце)".Show(); for (int i = 0; i < t; i++)//исправили, теперь пробуем новый метод { x[i] = c[i]; } for (int j = 0; j < t - 1; j++) { sum += x[j] * A[t - 1, j]; } x[t - 1] = (b[t - 1] - sum) / A[t - 1, t - 1]; sum = 0; tmp = Error(t); if (VALUE_FOR_ULTRA < tmp) { for (int i = 0; i < t - 1; i++) { x[i] = c[i]; } x[t - 1] = 0; existres = false; } else { $"Погрешность уменьшена МИНИМАКОЙ НА КОНЦЕ на {(VALUE_FOR_ULTRA - tmp) / VALUE_FOR_ULTRA * 100} % (now {tmp})".Show(); VALUE_FOR_ULTRA = tmp; } } else { $"Погрешность уменьшена ПОЛНОЙ МИНИМАКОЙ на {(VALUE_FOR_ULTRA - tmp1) / VALUE_FOR_ULTRA * 100}% (now {tmp1})".Show(); VALUE_FOR_ULTRA = tmp1; } } else { $"Погрешность уменьшена МИНИМАКОЙ СПИДГАУССА на {(VALUE_FOR_ULTRA - tmp) / VALUE_FOR_ULTRA * 100} % (now {tmp})".Show(); VALUE_FOR_ULTRA = tmp; } } else { $"Погрешность уменьшена СПИДГАУССОМ на {(VALUE_FOR_ULTRA - tmp) / VALUE_FOR_ULTRA * 100} % (now {tmp})".Show(); VALUE_FOR_ULTRA = tmp; } ErrorsMas[t - 1] = VALUE_FOR_ULTRA; Functional f = (Point x) => { sum = 0; for (int ii = 1; ii <= t; ii++) { sum += this.x[ii - 1] * KursMethods.masPoints[ii - 1].PotentialF(x); } double s = sum - KursMethods.fig(x); return(s * s); }; ErrorMasP[t - 1] = Math.Sqrt(DoubleIntegral(f, curve, curve.S, FuncMethods.DefInteg.Method.GaussKronrod15, 0.001, FuncMethods.DefInteg.countY)); //UltraCount = t; NEVA = Nev(A, x, b, t); Functional fM1 = (Point x) => { // Functional M1ApproxD = (Point u) => //{ // Functional ro = (Point y) => // { // double sss = 0; for (int i = 0; i < t; i++) sss += this.x[i] * KursMethods.masPoints[i].PotentialF(u); // return sss * KursMethods.Exy(u, y); // }; // //return DoubleIntegral(ro, curve, curve.S, FuncMethods.DefInteg.Method.GaussKronrod15, 0.001, FuncMethods.DefInteg.countY) /*IntegralClass.Integral(ro, KursMethods.CIRCLE - 1)*/ /*- right(x)*/; // return IntegralClass.Integral(ro, KursMethods.CIRCLE - 1); //}; double su = KursMethods.M1Approx(x); double s = su - KursMethods.U(x); return(s * s); }; if (existres) { $"Считается погрешность метода 1".Show(); ErrorsM1[t - 1] = Math.Sqrt(IntegralClass.Integral2(fM1, KursMethods.CIRCLE - 1)); } else { ErrorsM1[t - 1] = ErrorsM1[t - 2]; } $"{t} = {ErrorsM1[t - 1]}".Show(); }