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();
}
