예제 #1
0
파일: Form1.cs 프로젝트: MMShatalin/Differ
        private void пересчитатьРеактивностиToolStripMenuItem_Click(object sender, EventArgs e)
        {
            Reactivity mynewR = new Reactivity();

            mynewR.AddFirstData(mydata.Data[0][0], mydata.Data[MyConst.tok1][0]);
            for (int i = 1; i < mydata.Data[MyConst.tok1].Count; i++)
            {
                mynewR.AddData(mydata.Data[0][i], mydata.Data[MyConst.tok1][i]);
                chart1.Series[1].Points.AddXY(mydata.Data[0][i], mynewR.Ro);
            }
        }
예제 #2
0
        /// <summary>
        /// КНОПКА РАСЧЕТА ТЕКУЩЕГО ВОЗМУЩЕНИЯ - КНОПКА "МНК"
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void button5_Click(object sender, EventArgs e)
        {
            //ЭТА КНОПКА НАЧАЛА РАСЧЕТА
            //ИТАК НУЖНО ДЛЯ НАЧАЛА СЧЕТА ПОДГОТОВИТЬ ВСЕ ДАННЫЕ
            ///НА ПОНАДОБЯТЬСЯ:
            ///

            //Пускай это будет у нас для удобства в отдельных массивах

            List <double> time = new List <double>();
            List <double> I    = new List <double>();
            List <double> P    = new List <double>();
            List <double> pdH  = new List <double>();

            /*Итак у нас есть:
             * 1) исходный двумерный массив входных данных MyInputData
             * 2) индекс когда начались все перемещения  myOneVozmuchenie.Ir или что тоже самое myOneVozmuchenie.PerIndex[0, 0];
             * 3) кол-во перемещений myOneVozmuchenie.NPer;
             * 4) тепловая мощность в момент начала перемещений myOneVozmuchenie.Nt
             * 5) индекс окончания последнего перемещения  myOneVozmuchenie.PerIndex[myOneVozmuchenie.NPer - 1, 1];
             */


            //обозначили границы нашего перемещения
            myOneVozmuchenie.Ir = myOneVozmuchenie.PerIndex[0, 0];

            myOneVozmuchenie.Iend = myOneVozmuchenie.PerIndex[myOneVozmuchenie.NPer - 1, 1] + MyConst.delayTime; //при условии что частота 40 Гц добавим 5 секунд



            //Reactivity mynewR = new Reactivity();
            //mynewR.AddFirstData(mydata.Data[0][0], mydata.Data[MyConst.tok1][0]);
            //for (int i = 1; i < mydata.Data[MyConst.tok1].Count; i++)
            //{
            //    mynewR.AddData(mydata.Data[0][i], mydata.Data[MyConst.tok1][i]);
            //    chart1.Series[1].Points.AddXY(mydata.Data[0][i], mynewR.Ro);
            //}

            Reactivity mynewR = new Reactivity();

            for (int i = myOneVozmuchenie.Ir; i < myOneVozmuchenie.Iend; i++)
            {
                time.Add(myInputData.Data[0][i]);
                //ну так у Боева в методе GetJ
                I.Add(myInputData.Data[MyConst.tok1][i] / myInputData.Data[MyConst.tok1][myOneVozmuchenie.Ir]);


                if (i == myOneVozmuchenie.Ir)
                {
                    mynewR.AddFirstData(myInputData.Data[0][i], I[0]);
                }
                else
                {
                    mynewR.AddData(myInputData.Data[0][i], I[I.Count - 1]);
                }
                P.Add(mynewR.Ro);

                this.form1Chart1.Series[3].Points.AddXY(myInputData.Data[0][i], mynewR.Ro);
                //this.form1Chart1.Series[4].Points.AddXY(myInputData.Data[0][i], myInputData.Data[MyConst.H10][i]);


                //myOneVozmuchenie.PerIndex[myOneVozmuchenie.NPer - 1, 1]
            }



            ///нужно раздробить перемещения
            double ddH = 0;

            for (int i = 0; i < myOneVozmuchenie.NPer; i++)
            {
                //ведь каждый шаг должен в сумме давать 2
                for (int j = myOneVozmuchenie.PerIndex[i, 0] + 1; j <= myOneVozmuchenie.PerIndex[i, 1]; j++)
                {
                    ddH -= (2 / (myOneVozmuchenie.Per[i, 1] - myOneVozmuchenie.Per[i, 0])) * (myInputData.Data[0][j] - myInputData.Data[0][j - 1]);
                    pdH.Add(ddH);
                }

                //ДОБАВИТЬ КУСОК КОГДА НЕ ПЕРЕМЕЩАЛОСЬ

                if (i != myOneVozmuchenie.NPer - 1)
                {
                    for (int jj = myOneVozmuchenie.PerIndex[i, 1] + 1; jj <= myOneVozmuchenie.PerIndex[i + 1, 0]; jj++)
                    {
                        pdH.Add(ddH);
                    }
                }
            } // for i по всем перемещениям



            for (int i = 0; i < MyConst.delayTime; i++)
            {
                pdH.Add(ddH);
            }

            // MessageBox.Show(pdH.Count.ToString());



            this.form1Chart1.Series[4].Points.Clear();



            //Ну все якобы ВСЕ ГОТОВО К СЧЕТУ


            //List<double> time = new List<double>();
            //List<double> I = new List<double>();
            //List<double> P = new List<double>();
            //List<double> pdH = new List<double>();
            //Сверка всех длин массива прошла успешно
            // MessageBox.Show(time.Count + "  " + I.Count + " " + P.Count + " " + pdH.Count);



            double       Ss    = 1000;
            pertubResult tempR = new pertubResult();

            for (int i = 0; i < 200; i++)
            {
                tempR = Calc(3 + i / 200, time, I, P, pdH);
                if (tempR.SS > Ss)
                {
                    tempR = Calc(3 + (i - 1) / 200, time, I, P, pdH);
                    break;
                }
                Ss = tempR.SS;
            }

            tempR.Ro = tempR.Ro * MyConst.Rect.Beff;
            tempR.aH = tempR.aH * MyConst.Rect.Beff;
            tempR.b  = tempR.b * MyConst.Rect.Beff;
            tempR.SS = tempR.SS * MyConst.Rect.Beff;


            // this.mydataTableResult.Rows.Add(80, "DR/DH", "DR/DN", "tau");


            //  MessageBox.Show(tempR.FF.Count.ToString());


            this.mydataTableResult.Rows.Add(myInputData.Data[MyConst.H10][myOneVozmuchenie.Ir], tempR.aH, tempR.b / myOneVozmuchenie.Nt, tempR.tau);


            StreamWriter sr = new StreamWriter("Возмущение " + myInputData.Data[MyConst.H10][myOneVozmuchenie.Ir] + "_" + myInputData.Data[MyConst.H9][myOneVozmuchenie.Ir] + "_" + myInputData.Data[MyConst.H8][myOneVozmuchenie.Ir] + " .csv");

            sr.WriteLine("h10 s; h10 e; h10 mid; h9 s; h9 e; h9 mid; h8 s; h8 e; h8 mid; N; b; aH; aN; tau ; Ro;  'PO PICHKAM'; 'Rm = myresult.Ro + myresult.aH * dH[i] + myresult.b * F;'");

            sr.Write(myInputData.Data[MyConst.H10][myOneVozmuchenie.Ir] + ";" + myInputData.Data[MyConst.H10][myOneVozmuchenie.Iend] + ";" + (myInputData.Data[MyConst.H10][myOneVozmuchenie.Ir] + myInputData.Data[MyConst.H10][myOneVozmuchenie.Iend]) / 2 + ";");
            sr.Write(myInputData.Data[MyConst.H9][myOneVozmuchenie.Ir] + ";" + myInputData.Data[MyConst.H9][myOneVozmuchenie.Iend] + ";" + (myInputData.Data[MyConst.H9][myOneVozmuchenie.Ir] + myInputData.Data[MyConst.H9][myOneVozmuchenie.Iend]) / 2 + ";");
            sr.Write(myInputData.Data[MyConst.H8][myOneVozmuchenie.Ir] + ";" + myInputData.Data[MyConst.H8][myOneVozmuchenie.Iend] + ";" + (myInputData.Data[MyConst.H8][myOneVozmuchenie.Ir] + myInputData.Data[MyConst.H8][myOneVozmuchenie.Iend]) / 2 + ";");

            sr.Write(myOneVozmuchenie.Nt);

            sr.Write(";" + tempR.b);

            sr.Write(";" + tempR.aH);
            sr.Write(";" + tempR.b / myOneVozmuchenie.Nt);
            sr.Write(";" + tempR.tau);
            sr.Write(";" + tempR.Ro);



            //посчитаем aH по пичкам

            double DroAver = 0;

            for (int i = 0; i < myOneVozmuchenie.NPer; i++)
            {
                MessageBox.Show(myInputData.Data[MyConst.r1][myOneVozmuchenie.PerIndex[i, 1]].ToString() + " " + myInputData.Data[MyConst.r1][myOneVozmuchenie.PerIndex[i, 0]]);
                DroAver = DroAver + (myInputData.Data[MyConst.r1][myOneVozmuchenie.PerIndex[i, 1]] - myInputData.Data[MyConst.r1][myOneVozmuchenie.PerIndex[i, 0]]);


                //ведь каждый шаг должен в сумме давать 2
                //for (int j = myOneVozmuchenie.PerIndex[i, 0] + 1; j <= myOneVozmuchenie.PerIndex[i, 1]; j++)
                //{
                //    ddH -= (2 / (myOneVozmuchenie.Per[i, 1] - myOneVozmuchenie.Per[i, 0])) * (myInputData.Data[0][j] - myInputData.Data[0][j - 1]);
                //    pdH.Add(ddH);
                //}
            }
            MessageBox.Show(myOneVozmuchenie.NPer.ToString());
            MessageBox.Show(DroAver.ToString());
            DroAver = DroAver / myOneVozmuchenie.NPer;
            MessageBox.Show(DroAver.ToString());
            DroAver = DroAver / -2;
            MessageBox.Show(DroAver.ToString());
            DroAver = DroAver * MyConst.Rect.Beff;
            MessageBox.Show(DroAver.ToString());
            sr.Write(";" + DroAver);

            sr.WriteLine();

            // sr.WriteLine("начало -" + myOneVozmuchenie.Ir + "; конец - " + myOneVozmuchenie.Iend);

            sr.WriteLine("Time; Current; React; pdH; F");

            for (int i = 0; i < pdH.Count; i++)
            {
                this.form1Chart1.Series[4].Points.AddXY(myInputData.Data[0][myOneVozmuchenie.Ir + i], pdH[i]);
                sr.WriteLine(time[i] + ";" + I[i] + ";" + P[i] + ";" + pdH[i] + " ; " + tempR.FF[i]);
            }



            sr.Close();



            //for (int i = myOneVozmuchenie.Ir; i < myOneVozmuchenie.Iend; i++)
            //{
            //    if ((i > myOneVozmuchenie.PerIndex[currPerem, 0]) && (i < myOneVozmuchenie.PerIndex[currPerem, 1]))
            //    {

            //        ddH+=(2/(myOneVozmuchenie.Per[currPerem,1]-myOneVozmuchenie.Per[currPerem,0]))*(myInputData.Data[0][i])

            //    }
            //    else
            //    {



            //        if (true)
            //        {

            //        }


            //    }



            //}



            //Занесем пока все что есть, а все что просто. А дальше будем думать



            //СКОРЕЕ ВСЕГО ЭТО ВСЕ ЧТО НАМ ПОНАДОБИТЬСЯ:
            //double[] I -  МАССИВ ТОКОВ С ВЫСОКОЙ ЧАСТОТОЙ
            ///double[] pm - МАССИВ РЕАКТИВНОСТЕЙ РАСЧИТАННЫХ ПО ТОКАМ с такой де высокой частосто
            ///double[] H - МАССИВ ПОЛОЖЕНИЯ ГРУПП ДЛЯ КАЖДОГО ОТЧЕТА (перессчитанная также с высокой частостой для большой частоты регистрации)


            //this.mydataTableResult.Rows.Add(80, 1E-3, 1E-4);
            //this.mydataTableResult.Rows.Add(80, 1E-3, 1E-4);
            //MessageBox.Show(this.form1Chart1.ChartAreas[0].CursorX.Position.ToString());
        }