Advertisement
pVinc

kozik_walotek

Oct 22nd, 2021
1,013
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C# 2.48 KB | None | 0 0
  1.             N = 3;
  2.             numericUpDown1.Value = N;
  3.             UstawTablice();
  4.             int blad;
  5.             AZ = new Complex[N + 1, N + 1];
  6.             BZ = new Complex[N + 1];
  7.             XZ = new Complex[N + 1];
  8.  
  9.             double R1 = 10;
  10.             double R2 = 10;
  11.             double R3 = 10;
  12.             double XC2 = 10;
  13.             double XL1 = 10;
  14.             double XL3 = 10;
  15.             Complex E1 = 10;
  16.             Complex E3 = 10;
  17.  
  18.  
  19.             AZ[1, 1] = 1; AZ[1, 2] = -1; AZ[1, 3] = -1;
  20.             AZ[2, 1] = new Complex(R1, XL1); AZ[2, 2] = new Complex(R2, -XC2); AZ[2, 3] = 0;
  21.             AZ[3, 1] = 0; AZ[3, 2] = new Complex(R2, -XC2); AZ[3, 3] = new Complex(R3, XL3);
  22.  
  23.             BZ[1] = Complex.Zero;
  24.             BZ[2] = E1;
  25.             BZ[3] = E3;
  26.  
  27.             for (int i = 1; i <= N; i++)
  28.             {
  29.                 for (int j = 1; j <= N; j++)
  30.                 {
  31.                     MacierzA.Rows[i-1].Cells[j-1].Value = AZ[i,j].ToString();
  32.                 }
  33.                 WektorB.Rows[i - 1].Cells[0].Value = BZ[i].ToString();
  34.             }
  35.             if (radioButton1.Checked)
  36.             {
  37.                 blad = MetodaGaussa.RozRowMacGaussa(AZ, BZ, XZ, 1e-30);
  38.             }
  39.             else
  40.             {
  41.                 blad = MetodaRozkladuLU.RozRowMacCroutDoolitta(AZ, BZ, XZ, 1e-30);
  42.             }
  43.             if (blad == 0)
  44.                 for (int i = 1; i <= N; i++)
  45.                     WektorX[0, i - 1].Value = XZ[i].ToString("0.000000");
  46.  
  47.            
  48.             textBox1.Text = "IpK bilans: " + (XZ[1] - XZ[2] - XZ[3]).ToString();
  49.             textBox2.Text = "IIpK bilans: " + (E1 - XZ[1] * (R1 + Complex.j * XL1) - XZ[2] * (R2 - Complex.j * XC2)).ToString();
  50.  
  51.             double Pz = (E1 * XZ[1].Conjugate + E3 * XZ[3].Conjugate).Re;
  52.             textBox3.Text = "  Pz = " + Pz.ToString();
  53.  
  54.             double Podb = (XZ[1].Re * XZ[1].Re * (R1 + XL1) + XZ[2].Re * XZ[2].Re * (R2 + XC2) + XZ[3].Re * XZ[3].Re * (R3 + XL3));
  55.             textBox4.Text = "Podb = " + Podb.ToString();
  56.            
  57.             Complex Sz = new Complex(0, 0);
  58.             Sz = E1 * XZ[1].Conjugate + E3 * XZ[3].Conjugate;
  59.             textBox5.Text = "  Sz= " + Sz.ToString();
  60.  
  61.             Complex Sodb = new Complex(0, 0);
  62.             Sodb = (XZ[1].Abs * XZ[1].Abs * (R1 + Complex.j * XL1)) + (XZ[2].Abs * XZ[2].Abs * (R2 - Complex.j * XC2)) + (XZ[3].Abs * XZ[3].Abs * (R3 + Complex.j * XL3));
  63.             textBox6.Text = "Sodb = " + Sodb.ToString();
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement