Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- N = 3;
- numericUpDown1.Value = N;
- UstawTablice();
- int blad;
- AZ = new Complex[N + 1, N + 1];
- BZ = new Complex[N + 1];
- XZ = new Complex[N + 1];
- Complex e1 = 15;
- Complex e2 = 10;
- double r1 = 5, r2 = 10, r3 = 15;
- double xl1 = 15, xl2=20, xc = 10;
- AZ[1, 1] = 1; AZ[1, 2] = -1; AZ[1, 3] = -1;
- AZ[2, 1] = r1; AZ[2, 2] = new Complex(r2, xl1); AZ[2, 3] = 0;
- AZ[3, 1] = 0; AZ[3, 2] = new Complex(-r2, -xl1); AZ[3, 3] = new Complex(r3, xl2-xc);
- BZ[1] = Complex.Zero;
- BZ[2] = e1;
- BZ[3] = e2;
- for (int i = 1; i <= N; i++)
- {
- for (int j = 1; j <= N; j++)
- {
- MacierzA.Rows[i - 1].Cells[j - 1].Value = AZ[i, j].ToString();
- }
- WektorB.Rows[i - 1].Cells[0].Value = BZ[i].ToString();
- }
- if (radioButton1.Checked)
- {
- blad = MetodaGaussa.RozRowMacGaussa(AZ, BZ, XZ, 1e-30);
- }
- else
- {
- blad = MetodaRozkladuLU.RozRowMacCroutDoolitta(AZ, BZ, XZ, 1e-30);
- }
- if (blad == 0)
- for (int i = 1; i <= N; i++)
- WektorX[0, i - 1].Value = XZ[i].ToString("0.000000");
- textBox1.Text = "IpK bilans: " + (XZ[1] - XZ[2] - XZ[3]).ToString();
- textBox2.Text = "IIpK bilans: " + (e1-XZ[1]*r1-XZ[2]*(r2+Complex.j*xl1)).ToString();
- double Pz = (e1 * XZ[1].Conjugate + e2 * XZ[3].Conjugate).Re;
- double Podb = (XZ[1] * XZ[1].Conjugate * r1).Re + (XZ[2] * XZ[2].Conjugate * (r2 + Complex.j * xl1)).Re + (XZ[3] * XZ[3].Conjugate * (r3 + Complex.j * xl2 - Complex.j * xc)).Re;
- textBox3.Text = " Pz = " + Pz.ToString();
- textBox4.Text = "Podb = " + Podb.ToString();
- Complex Sz = new Complex(0, 0);
- Complex Sodb = new Complex(0, 0);
- Sz = e1 * XZ[1].Conjugate + e2 * XZ[3].Conjugate;
- textBox5.Text = " Sz= " + Sz.ToString();
- Sodb = (XZ[1] * XZ[1].Conjugate * r1) + (XZ[2]*XZ[2].Conjugate*(r2 + Complex.j*xl1)) + (XZ[3]*XZ[3].Conjugate*(r3+Complex.j*xl2-Complex.j*xc));
- textBox6.Text = "Sodb = " + Sodb.ToString();
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement