Esempio n. 1
0
        public void DrawIntegr(ref Grid gr, Fraction a, Fraction b, int res = 100)
        {
            int      por    = this.NCol - 2;
            int      row    = 1;
            Fraction xfrMin = new Fraction();
            Fraction yfrMin = new Fraction();
            Fraction xfrMax = new Fraction();
            Fraction yfrMax = new Fraction();

            //Fraction xfr0 = new Fraction((int)((gr.Width - 20) / 2), 1);
            //Fraction yfr0 = new Fraction((int)((gr.Height - 20) / 2), 1);
            xfrMin = this.matr.A[1, 0].CopyFr;
            //yfrMin.FractTxt = this.txtBxs[1, 1].Text;
            //double yMin = yfrMin.ToDouble();
            Fraction h = new Fraction();

            xfrMax = this.matr.A[this.NRow - 1, 0].CopyFr;
            //double xMax = xfrMax.ToDouble();
            yfrMax = this.matr.A[this.NRow - 1, 1].CopyFr;
            yfrMin = this.matr.A[1, 1].CopyFr;
            for (Fraction i = xfrMin; i <= xfrMax; i += (xfrMax - xfrMin) / 50)
            {
                if (GetPolyNewtonVal(i, por, row) >= yfrMax)
                {
                    yfrMax = GetPolyNewtonVal(i, por, row);
                }
            }
            for (Fraction i = xfrMin; i <= xfrMax; i += (xfrMax - xfrMin) / 50)
            {
                if (GetPolyNewtonVal(i, por, row) < yfrMin)
                {
                    yfrMin = GetPolyNewtonVal(i, por, row);
                }
            }
            CoordAxis axis = new CoordAxis(ref gr, xfrMin.ToDouble(), xfrMax.ToDouble(), yfrMin.ToDouble(), yfrMax.ToDouble(), 10);

            //Fraction yMult = new Fraction((int)((gr.Height - 20) * 1000), 1000);
            h = (b - a) / res;
            PointCollection pc = new PointCollection();

            for (Fraction i = a; i <= b; i += h)
            {
                //pc.Add(new System.Windows.Point((xfr0 + i*xMult).ToDouble(), (yfr0 - GetPolyNewtonVal(i, por, row)*yMult).ToDouble()));
                pc.Add(axis.ConvertCoord(i.ToDouble(), GetPolyNewtonVal(i, por, row).ToDouble()));
            }

            //double wdt = gr.Width - 20;
            //double hgt = gr.Height - 20;

            PointCollection pgPt = new PointCollection();

            pgPt.Add(axis.ConvertCoord(a.ToDouble(), 0));
            foreach (Point pt in pc)
            {
                pgPt.Add(pt);
            }
            pgPt.Add(axis.ConvertCoord(b.ToDouble(), 0));
            Polygon pg = new Polygon();

            pg.Points = pgPt;
            pg.Fill   = Brushes.Pink;
            gr.Children.Add(pg);
        }
Esempio n. 2
0
        public void DrawPolynom(ref Grid gr, int por, Fraction x, int row = 1, int res = 100)
        {
            Fraction xfrMin = new Fraction();
            Fraction yfrMin = new Fraction();
            Fraction xfrMax = new Fraction();
            Fraction yfrMax = new Fraction();

            xfrMin = this.matr.A[1, 0].CopyFr;
            yfrMin = this.matr.A[1, 1].CopyFr;
            //double yMin = yfrMin.ToDouble();
            Fraction h = new Fraction();

            xfrMax = this.matr.A[this.NRow - 1, 0].CopyFr;
            //double xMax = xfrMax.ToDouble();
            yfrMax = this.matr.A[this.NRow - 1, 1].CopyFr;
            for (Fraction i = xfrMin; i <= xfrMax; i += (xfrMax - xfrMin) / res)
            {
                if (GetPolyNewtonVal(i, por, row) >= yfrMax)
                {
                    yfrMax = GetPolyNewtonVal(i, por, row);
                }
            }
            for (Fraction i = xfrMin; i <= xfrMax; i += (xfrMax - xfrMin) / res)
            {
                if (GetPolyNewtonVal(i, por, row) < yfrMin)
                {
                    yfrMin = GetPolyNewtonVal(i, por, row);
                }
            }
            CoordAxis axis = new CoordAxis(ref gr, xfrMin.ToDouble(), xfrMax.ToDouble(), yfrMin.ToDouble(), yfrMax.ToDouble(), 10);

            axis.DrawAxis();
            //Fraction xMult = new Fraction((int)((gr.Width - 20)*1000), 1000);
            //xMult /= xfrMax*2;
            //Fraction yMult = new Fraction((int)((gr.Height - 20)*1000), 1000);
            //yMult /= yfrMax*2;
            h = (xfrMax - xfrMin) / res;

            Polyline pl = new Polyline();

            Ellipse el = new Ellipse();

            el.HorizontalAlignment = HorizontalAlignment.Left;
            el.VerticalAlignment   = VerticalAlignment.Top;
            el.Width  = 8;
            el.Height = 8;
            Point valXY = axis.ConvertCoord(x.ToDouble(), (GetPolyNewtonVal(x, por, row).ToDouble()));

            System.Windows.Thickness th = new System.Windows.Thickness(valXY.X - el.Width / 2, valXY.Y - el.Height / 2, 0, 0);
            el.Margin = th;

            el.Stroke = Brushes.Red;
            el.Fill   = Brushes.Red;

            PointCollection pc = new PointCollection();

            for (Fraction i = xfrMin; i <= xfrMax; i += h)
            {
                //pc.Add(new System.Windows.Point((xfr0 + i*xMult).ToDouble(), (yfr0 - GetPolyNewtonVal(i, por, row)*yMult).ToDouble()));
                pc.Add(axis.ConvertCoord(i.ToDouble(), GetPolyNewtonVal(i, por, row).ToDouble()));
            }
            pl.Points          = pc;
            pl.Stroke          = Brushes.Blue;
            pl.StrokeThickness = 4;
            double wdt  = gr.Width - 20;
            double hgt  = gr.Height - 20;
            Line   xVal = new Line();

            //xVal.X1 = (xfr0 + x * xMult).ToDouble();
            xVal.X1              = axis.ConvertCoord(x.ToDouble(), GetPolyNewtonVal(x, por, row).ToDouble()).X;
            xVal.Y1              = axis.ConvertCoord(x.ToDouble(), GetPolyNewtonVal(x, por, row).ToDouble()).Y;
            xVal.X2              = xVal.X1;
            xVal.Y2              = axis.crossAxis.Y;
            xVal.Stroke          = Brushes.Red;
            xVal.StrokeDashArray = new DoubleCollection()
            {
                8, 4
            };
            Line yVal = new Line();

            yVal.X1              = axis.ConvertCoord(x.ToDouble(), GetPolyNewtonVal(x, por, row).ToDouble()).X;
            yVal.Y1              = axis.ConvertCoord(x.ToDouble(), GetPolyNewtonVal(x, por, row).ToDouble()).Y;
            yVal.X2              = axis.crossAxis.X;
            yVal.Y2              = yVal.Y1;
            yVal.Stroke          = Brushes.Red;
            yVal.StrokeDashArray = new DoubleCollection()
            {
                8, 4
            };
            gr.Children.Add(pl);
            gr.Children.Add(xVal);
            gr.Children.Add(yVal);


            TextBox valX = new TextBox();

            valX.Text   = "" + x.ToDouble();
            valX.Width  = 100;
            valX.Height = 25;
            valX.HorizontalAlignment = HorizontalAlignment.Left;
            valX.VerticalAlignment   = VerticalAlignment.Top;
            TextBox valY = new TextBox();

            valY.Text   = "" + GetPolyNewtonVal(x, por, row).ToDouble();
            valY.Width  = 200;
            valY.Height = 25;
            valY.HorizontalAlignment = HorizontalAlignment.Left;
            valY.VerticalAlignment   = VerticalAlignment.Top;
            TextBox polyErr  = new TextBox();
            TextBox diff1    = new TextBox();
            TextBox diff2    = new TextBox();
            TextBox diff1Err = new TextBox();
            TextBox diff2Err = new TextBox();

            polyErr.Text    = "погр. интерполирования равна " + GetPolyNewtonErr(x, por, row).ToDouble();
            diff1.Text      = "первая производная равна " + GetPolyNewtonDiff1Val(x, por, row).ToDouble();
            diff1Err.Text   = "ошибка первой производной " + GetPolyNewtonDiff1Err(x, por, row).ToDouble();
            diff2.Text      = "вторая производная равна " + GetPolyNewtonDiff2Val(x, por, row).ToDouble();
            diff2Err.Text   = "ошибка второй производной " + GetPolyNewtonDiff2Err(x, por, row).ToDouble();
            polyErr.Width   = 300;
            polyErr.Height  = 25;
            diff1.Width     = 300;
            diff1.Height    = 25;
            diff2.Width     = 300;
            diff2.Height    = 25;
            diff1Err.Width  = 300;
            diff1Err.Height = 25;
            diff2Err.Width  = 300;
            diff2Err.Height = 25;
            SolidColorBrush br = new SolidColorBrush(Color.FromArgb(0, 0, 0, 0));

            valY.BorderBrush                 = br;
            valX.BorderBrush                 = br;
            polyErr.BorderBrush              = br;
            diff1.BorderBrush                = br;
            diff2.BorderBrush                = br;
            diff1Err.BorderBrush             = br;
            diff2Err.BorderBrush             = br;
            valX.Background                  = br;
            valY.Background                  = br;
            polyErr.Background               = br;
            diff1.Background                 = br;
            diff2.Background                 = br;
            diff1Err.Background              = br;
            diff2Err.Background              = br;
            valX.HorizontalContentAlignment  = System.Windows.HorizontalAlignment.Center;
            diff1.HorizontalContentAlignment = System.Windows.HorizontalAlignment.Left;
            diff2.HorizontalContentAlignment = System.Windows.HorizontalAlignment.Left;
            polyErr.Margin  = new System.Windows.Thickness(0, 0, wdt - diff1.Width, hgt - diff1.Height);
            diff1.Margin    = new System.Windows.Thickness(0, 0, wdt - diff1.Width, hgt - 3 * diff1.Height);
            diff1Err.Margin = new System.Windows.Thickness(0, 0, wdt - diff1.Width, hgt - 5 * diff1.Height);
            diff2.Margin    = new System.Windows.Thickness(0, 0, wdt - diff2.Width, hgt - 7 * diff2.Height);
            diff2Err.Margin = new System.Windows.Thickness(0, 0, wdt - diff2.Width, hgt - 9 * diff2.Height);
            if (valXY.Y <= axis.crossAxis.Y)
            {
                valX.Margin = new System.Windows.Thickness(valXY.X - valX.Width / 2, axis.crossAxis.Y - valY.Height, 0, 0);
            }
            else
            {
                valX.Margin = new System.Windows.Thickness(valXY.X - valX.Width / 2, axis.crossAxis.Y + valY.Height, 0, 0);
            }
            if (valXY.X <= axis.crossAxis.X)
            {
                valY.Margin        = new System.Windows.Thickness(axis.crossAxis.X - axis.bord - valY.Width, valXY.Y - valY.Height, 0, 0);
                valY.TextAlignment = axis.taY;
            }
            else
            {
                valY.Margin = new System.Windows.Thickness(axis.crossAxis.X + axis.bord, valXY.Y - valY.Height, 0, 0);
                //valY.HorizontalContentAlignment = System.Windows.HorizontalAlignment.Left;
                valY.TextAlignment = axis.taY;
            }

            gr.Children.Add(el);
            gr.Children.Add(valX);
            gr.Children.Add(valY);
            gr.Children.Add(polyErr);
            gr.Children.Add(diff1);
            gr.Children.Add(diff2);
            gr.Children.Add(diff1Err);
            gr.Children.Add(diff2Err);
        }