// linearni regrese
// y=ax+b
public static void LinearRegresion(SLValueList list, ref double a, ref double b)
{
SLValueList newlist = new SLValueList();
double sumax = 0;
double suma2x = 0;
double sumaxy = 0;
double sumay = 0;
double suma2y = 0;
double prumx = 0;
double prumy = 0;
for (int i = 0; i < list.Count; i++)
{
SLValueXY val = list[i];
sumax += val.X;
suma2x += val.X * val.X;
sumay += val.Y;
suma2y += val.Y * val.Y;
sumaxy += val.X * val.Y;
}
prumx = sumax / list.Count;
prumy = sumay / list.Count;
double sx = (1 / (float)list.Count * suma2x) - prumx * prumx;
double sy = (1 / (float)list.Count * suma2y) - prumy * prumy;
double sxy = (1 / (float)list.Count * sumaxy) - prumx * prumy;
b = sxy / sx;
a = prumy - b * prumx;
}
public int Add(double x, double y)
{
SLValueXY val = new SLValueXY();
val.X = x;
val.Y = y;
return(Add(val));
}
public bool SaveToFile(string fname)
{
using (System.IO.StreamWriter sw = new System.IO.StreamWriter(fname))
{
for (int i = 0; i < this.Count; i++)
{
SLValueXY val = (SLValueXY)this[i];
sw.WriteLine("{0};{1}", val.X, val.Y);
}
}
return(true);
}
//lepsi nez EveluateFce
// lze parametrizovat
// --------------------
public SLValueList EvalFunction(string expr, double xfrom, double xto, double step)
{
if (xfrom >= xto)
{
return(null);
}
if ((xto - xfrom) < step)
{
return(null);
}
SLValueList newdata = new SLValueList();
newdata.Title = expr;
math.Expresion = expr;
MathVar var;
var = math.Get("x");
if (var == null)
{
var = math.AddVariable("x", 0);
}
for (double i = xfrom; i < xto + step; i += (double)step)
{
if (Math.Abs(i) < 0.00001)
{
i = 0;
}
i = Math.Round(i, 5);
var.value = i;
double x = i;
double y = math.Evaluate();
SLValueXY val = new SLValueXY();
val.X = x;
val.Y = y;
newdata.Add(val);
}
//Add(newdata);
return(newdata);
}
public static SLValueList EvalPolarFce(string xexpr, string yexpr, double from, double to, double step)
{
if (from >= to)
{
return(null);
}
if ((to - from) < step)
{
return(null);
}
SLValueList newdata = new SLValueList();
SLMath.SLMathExp mathx = new SLMathExp();
SLMath.SLMathExp mathy = new SLMathExp();
mathx.Expresion = xexpr;
mathy.Expresion = yexpr;
MathVar var1 = mathx.AddVariable("x", 0);
MathVar var2 = mathy.AddVariable("x", 0);
for (double i = from; i < to + step; i += (double)step)
{
if (Math.Abs(i) < 0.00001)
{
i = 0;
}
i = Math.Round(i, 5);
var1.value = i;
var2.value = i;
double x = mathx.Evaluate();
double y = mathy.Evaluate();
SLValueXY val = new SLValueXY();
val.X = x;
val.Y = y;
newdata.Add(val);
}
//Add(newdata);
return(newdata);
}
// smaze polozku
public void Remove(SLValueXY data)
{
base.Remove(data);
}
// prida polozku
public void Insert(int index, SLValueXY data)
{
base.Insert(index, data);
}
// prida dalsi polozku k datum
public int Add(SLValueXY data)
{
return(base.Add(data));
}
//*********** kresleni dat **********
protected void PaintDataList(PaintEventArgs e)
{
for (int i = 0; i < DataList.Count; i++)
{
SLMath.SLValueList ValueList = DataList[i];
int navic = 3 - (ValueList.Count - 4) % 3;
Point[] points = new Point[ValueList.Count + navic];
Pen barva = new Pen(ValueList.ColorLine);
Brush BrBarva = new SolidBrush(ValueList.ColorLine);
if (ValueList.ShowPoints == false && ValueList.ShowLine == false)
{
continue;
}
for (int j = 0; j < ValueList.Count; j++)
{
//int pocetF = ValueList.Count;
SLMath.SLValueXY valXY = (SLMath.SLValueXY)ValueList[j];
if (valXY.isValid() == false)
{
continue;
}
int bodX1 = (int)(x_abs((float)valXY.X));
int bodY1 = (int)(y_abs((float)valXY.Y));
points[j] = new Point(bodX1, bodY1);
if (j == ValueList.Count - 1)
{
for (int c = 0; c < navic; c++)
{
points[ValueList.Count + c] = new Point(bodX1, bodY1);
}
}
// kresli kolecka ... pokud najedeme mysi na bod
/*
* if ((RealMouseX >= (float)valXY.X-0.1 && RealMouseX <= (float)valXY.X+0.1) && (RealMouseY >= (float)valXY.Y-0.1 && RealMouseY <= (float)valXY.Y+0.1))
* e.Graphics.FillEllipse(BrBarva,bodX1-ValueList.PointSize-1,bodY1-ValueList.PointSize-1,2*ValueList.PointSize+2,2*ValueList.PointSize+2);*/
//
if (ValueList.ShowPoints == true)
{
if (ValueList.PointStyle == SLMath.SLValueList.PointStyles.Cross)
{
e.Graphics.DrawLine(barva, bodX1 - ValueList.PointSize, bodY1, bodX1 + ValueList.PointSize, bodY1);
e.Graphics.DrawLine(barva, bodX1, bodY1 - ValueList.PointSize, bodX1, bodY1 + ValueList.PointSize);
}
if (ValueList.PointStyle == SLMath.SLValueList.PointStyles.Circle)
{
e.Graphics.FillEllipse(BrBarva, bodX1 - ValueList.PointSize, bodY1 - ValueList.PointSize, 2 * ValueList.PointSize + 1, 2 * ValueList.PointSize + 1);
}
if (ValueList.PointStyle == SLMath.SLValueList.PointStyles.Cube)
{
e.Graphics.FillRectangle(BrBarva, bodX1 - ValueList.PointSize, bodY1 - ValueList.PointSize, ValueList.PointSize * 2 + 1, ValueList.PointSize * 2 + 1);
}
if (ValueList.PointStyle == SLMath.SLValueList.PointStyles.Diamond)
{
// Create points that define polygon.
Point point1 = new Point(bodX1 - ValueList.PointSize, bodY1);
Point point2 = new Point(bodX1, bodY1 - ValueList.PointSize);
Point point3 = new Point(bodX1 + ValueList.PointSize, bodY1);
Point point4 = new Point(bodX1, bodY1 + ValueList.PointSize);
Point[] diamPoints = { point1, point2, point3, point4 };
// Draw polygon to screen.
e.Graphics.FillPolygon(BrBarva, diamPoints);
}
if (ValueList.PointStyle == SLMath.SLValueList.PointStyles.Triangle)
{
// Create points that define polygon.
Point point1 = new Point(bodX1 - (2 * ValueList.PointSize / 3) - 1, (bodY1 + ValueList.PointSize / 3) + 1);
Point point2 = new Point(bodX1 + (2 * ValueList.PointSize / 3) + 1, (bodY1 + ValueList.PointSize / 3) + 1);
Point point3 = new Point(bodX1, (bodY1 - 2 * ValueList.PointSize / 3) - 1);
Point point4 = new Point(bodX1 - (2 * ValueList.PointSize / 3) - 1, (bodY1 + ValueList.PointSize / 3) - 1);
Point[] tryPoints = { point1, point2, point3, point4 };
// Draw polygon to screen.
e.Graphics.FillPolygon(BrBarva, tryPoints);
}
}
/*
* if (ValueList.ShowLine == true)
* {
* if (ValueList.AproxStyle == SLMath.SLValueList.AproximationStyles.aLine)
* {
* e.Graphics.DrawLine(barva,bodX1,bodY1,bodX2,bodY2);
* }
* if (ValueList.AproxStyle == SLMath.SLValueList.AproximationStyles.aBezier)
* {
*
* }
* }*/
}
if (ValueList.ShowLine == true)
{
if (ValueList.AproxStyle == SLMath.SLValueList.AproximationStyles.aLine)
{
e.Graphics.DrawLines(barva, points);
}
if (ValueList.AproxStyle == SLMath.SLValueList.AproximationStyles.aBezier)
{
e.Graphics.DrawBeziers(barva, points);
}
if (ValueList.AproxStyle == SLMath.SLValueList.AproximationStyles.aSpline)
{
e.Graphics.DrawCurve(barva, points);
}
}
}
}
public static SLValueList DerivateFce(SLValueList Fx)
{
SLValueList newdata = new SLValueList();
for (int k = 0; k < Fx.Count - 1; k++)
{
SLValueXY val1 = (SLMath.SLValueXY)Fx[k];
SLValueXY val2 = (SLMath.SLValueXY)Fx[k + 1];
SLValueXY newval = new SLValueXY();
float krok = (float)val2.X - (float)val1.X;
newval.X = val1.X;
newval.Y = ((float)val2.Y - (float)val1.Y) / krok; //!!!!!!!!!!
newdata.Add(newval);
}
return(newdata);
/*double point = System.Convert.ToDouble(bod.Text);
* double krok1 = System.Convert.ToDouble(k1.Text);
* double krok2 = System.Convert.ToDouble(k2.Text);
* double krok3 = System.Convert.ToDouble(k3.Text);
* double krok4 = System.Convert.ToDouble(k4.Text);
* double krok5 = System.Convert.ToDouble(k5.Text);
* double[] prok1 = new double[2];
* double[] prok2 = new double[2];
* double[] prok3 = new double[2];
* double[] prok4 = new double[2];
* double[] prok5 = new double[2];
*
*
* if (puntik1.Checked == true) // dvoubodovka prvni derivace
* {
* prok1[0] = (Fce(point+krok1) - Fce(point))/krok1;
* prok2[0] = (Fce(point+krok2) - Fce(point))/krok2;
* prok3[0] = (Fce(point+krok3) - Fce(point))/krok3;
* prok4[0] = (Fce(point+krok4) - Fce(point))/krok4;
* prok5[0] = (Fce(point+krok5) - Fce(point))/krok5;
* }
*
* if (puntik2.Checked == true) // tribodovka prvni derivace
* {
* prok1[0] = (-3*Fce(point) + 4*Fce(point+krok1) - Fce(point+2*krok1))/2*krok1;
* prok2[0] = (-3*Fce(point) + 4*Fce(point+krok2) - Fce(point+2*krok2))/2*krok2;
* prok3[0] = (-3*Fce(point) + 4*Fce(point+krok3) - Fce(point+2*krok3))/2*krok3;
* prok4[0] = (-3*Fce(point) + 4*Fce(point+krok4) - Fce(point+2*krok4))/2*krok4;
* prok5[0] = (-3*Fce(point) + 4*Fce(point+krok5) - Fce(point+2*krok5))/2*krok5;
* }
*
* if (puntik3.Checked == true) // ctyrbodovka prvni derivace
* {
* prok1[0] = (-11*Fce(point) + 18*Fce(point+krok1) - 9*Fce(point+2*krok1) + 2*Fce(point+3*krok1))/6*krok1;
* prok2[0] = (-11*Fce(point) + 18*Fce(point+krok2) - 9*Fce(point+2*krok2) + 2*Fce(point+3*krok2))/6*krok2;
* prok3[0] = (-11*Fce(point) + 18*Fce(point+krok3) - 9*Fce(point+2*krok3) + 2*Fce(point+3*krok3))/6*krok3;
* prok4[0] = (-11*Fce(point) + 18*Fce(point+krok4) - 9*Fce(point+2*krok4) + 2*Fce(point+3*krok4))/6*krok4;
* prok5[0] = (-11*Fce(point) + 18*Fce(point+krok5) - 9*Fce(point+2*krok5) + 2*Fce(point+3*krok5))/6*krok5;
* }
*
* if (puntik4.Checked == true) // petibodovka prvni derivace
* {
* prok1[0] = (-25*Fce(point) + 48*Fce(point+krok1) - 36*Fce(point+2*krok1) + 16*Fce(point+3*krok1) - 3*Fce(point+4*krok1))/12*krok1;
* prok2[0] = (-25*Fce(point) + 48*Fce(point+krok2) - 36*Fce(point+2*krok2) + 16*Fce(point+3*krok2) - 3*Fce(point+4*krok2))/12*krok2;
* prok3[0] = (-25*Fce(point) + 48*Fce(point+krok3) - 36*Fce(point+2*krok3) + 16*Fce(point+3*krok3) - 3*Fce(point+4*krok3))/12*krok3;
* prok4[0] = (-25*Fce(point) + 48*Fce(point+krok4) - 36*Fce(point+2*krok4) + 16*Fce(point+3*krok4) - 3*Fce(point+4*krok4))/12*krok4;
* prok5[0] = (-25*Fce(point) + 48*Fce(point+krok5) - 36*Fce(point+2*krok5) + 16*Fce(point+3*krok5) - 3*Fce(point+4*krok5))/12*krok5;
* }
*
* if (puntik5.Checked == true) // tribodovka druha derivace
* {
* prok1[1] = (Fce(point) - 2*Fce(point+krok1) + Fce(point+2*krok1))/krok1*krok1;
* prok2[1] = (Fce(point) - 2*Fce(point+krok2) + Fce(point+2*krok2))/krok2*krok2;
* prok3[1] = (Fce(point) - 2*Fce(point+krok3) + Fce(point+2*krok3))/krok3*krok3;
* prok4[1] = (Fce(point) - 2*Fce(point+krok4) + Fce(point+2*krok4))/krok4*krok4;
* prok5[1] = (Fce(point) - 2*Fce(point+krok5) + Fce(point+2*krok5))/krok5*krok5;
* }
*
* if (puntik6.Checked == true) // ctyrbodovka druha derivace
* {
* prok1[1] = (2*Fce(point) - 5*Fce(point+krok1) + 4*Fce(point+2*krok1) + Fce(point+3*krok1))/krok1*krok1;
* prok2[1] = (2*Fce(point) - 5*Fce(point+krok2) + 4*Fce(point+2*krok2) + Fce(point+3*krok2))/krok2*krok2;
* prok3[1] = (2*Fce(point) - 5*Fce(point+krok3) + 4*Fce(point+2*krok3) + Fce(point+3*krok3))/krok3*krok3;
* prok4[1] = (2*Fce(point) - 5*Fce(point+krok4) + 4*Fce(point+2*krok4) + Fce(point+3*krok4))/krok4*krok4;
* prok5[1] = (2*Fce(point) - 5*Fce(point+krok5) + 4*Fce(point+2*krok5) + Fce(point+3*krok5))/krok5*krok5;
* }
*
*
*
* v1.Text = prok1[0].ToString();
* v2.Text = prok2[0].ToString();
* v3.Text = prok3[0].ToString();
* v4.Text = prok4[0].ToString();
* v5.Text = prok5[0].ToString();
* v6.Text = prok1[1].ToString();
* v7.Text = prok2[1].ToString();
* v8.Text = prok3[1].ToString();
* v9.Text = prok4[1].ToString();
* v10.Text = prok5[1].ToString();*/
}
