Обратная матрица

@by-ga-ga · Регистрация: 16.05.2013

Студворк — интернет-сервис помощи студентам

1 элемент считает а дальше нули(проверял по калькулятору 1 элемент сходиться) вот код

C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
 
namespace Obr
{
    class Program
    {
        static void Vxod(out double[,] Mas, int n)
        {
            
            Mas = new double[n, n];
            Random rnd = new Random();
            for (int i = 0; i < n; i++)
                for (int j = 0; j < n; j++)
                    Mas[i, j] = rnd.Next(1, 9);
        }
        static void Print(double[,] Mas, int n)
        {
            for (int i = 0; i < Mas.GetLength(0); i++)
            {
                for (int j = 0; j < Mas.GetLength(1); j++)
                    Console.Write("{0,10:N}", Mas[i, j]);
                Console.WriteLine();
            }
        }
       static Double AlgDop(int size, double[,] m, int row, int col )
        {
            int i, j;  
           double[,] temp=new double[m.GetLength(0),m.GetLength(0)];
             for(i=0; i<row-1;i++)
           {
               for(j=0;j<col-1;j++)
                   temp[i,j]=m[i,j];
                   for(j=col+1;j<size;j++)
                       temp[i,j-1]=m[i,j];
           }
           for(i=row+1;i<size;i++)
           {
               for(j=0;j<col-1;j++)
                   temp[i-1,j]=m[i,j];
               for(j=col+1;j<size;j++)
                   temp[i-1,j-1]=m[i,j];
           }
           if ((row+col) % 2 == 0)
               return  Det(size-1,ref temp);
           else return  (-1)*Det(size-1,ref temp);
       }
 
        static Double Det(int size ,ref double[,] m)
       {
           int i; double r;
           if (size == 1) return  m[0, 0];
           else
               if (size== 2) return  (m[0, 0] * m[1, 1] - m[0, 1] * m[1, 0]);
               else
               {
                   r = 0;
                   for (i = 0; i < size; i++)
                   {
                       r = r + m[0, i] * AlgDop(size,m,0,i);
                   }
                   return  r;
               }
 
      }
        static void ObrMatr(double[,] a, double[,] am, int n)
        {    
            int i, j; double dt;
            dt = Det(n,ref a);
            for (i = 0; i < n; i++)
            {
                for (j = 0; j < n; j++)
                    am[i, j] = AlgDop(n,a,j,i) / dt;
            }
        }
        static void Main(string[] args)
        {
            Console.Write("Введите размерность квадратной матрицы m=");
            int n = int.Parse(Console.ReadLine());
            double[,] a;
            double[,] Obr= new double[n,n];
            Vxod(out a, n);
            Print(a, n);
            ObrMatr(a,Obr, n);
            Console.WriteLine();
            Print(Obr, n);
        }
    }
}

@Cheshire94 · 03.06.2013, 02:28

Сообщение от by-ga-ga

1 элемент считает а дальше нули(проверял по калькулятору 1 элемент сходиться) вот код

C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
 
namespace Obr
{
    class Program
    {
        static void Vxod(out double[,] Mas, int n)
        {
            
            Mas = new double[n, n];
            Random rnd = new Random();
            for (int i = 0; i < n; i++)
                for (int j = 0; j < n; j++)
                    Mas[i, j] = rnd.Next(1, 9);
        }
        static void Print(double[,] Mas, int n)
        {
            for (int i = 0; i < Mas.GetLength(0); i++)
            {
                for (int j = 0; j < Mas.GetLength(1); j++)
                    Console.Write("{0,10:N}", Mas[i, j]);
                Console.WriteLine();
            }
        }
       static Double AlgDop(int size, double[,] m, int row, int col )
        {
            int i, j;  
           double[,] temp=new double[m.GetLength(0),m.GetLength(0)];
             for(i=0; i<row-1;i++)
           {
               for(j=0;j<col-1;j++)
                   temp[i,j]=m[i,j];
                   for(j=col+1;j<size;j++)
                       temp[i,j-1]=m[i,j];
           }
           for(i=row+1;i<size;i++)
           {
               for(j=0;j<col-1;j++)
                   temp[i-1,j]=m[i,j];
               for(j=col+1;j<size;j++)
                   temp[i-1,j-1]=m[i,j];
           }
           if ((row+col) % 2 == 0)
               return  Det(size-1,ref temp);
           else return  (-1)*Det(size-1,ref temp);
       }
 
        static Double Det(int size ,ref double[,] m)
       {
           int i; double r;
           if (size == 1) return  m[0, 0];
           else
               if (size== 2) return  (m[0, 0] * m[1, 1] - m[0, 1] * m[1, 0]);
               else
               {
                   r = 0;
                   for (i = 0; i < size; i++)
                   {
                       r = r + m[0, i] * AlgDop(size,m,0,i);
                   }
                   return  r;
               }
 
      }
        static void ObrMatr(double[,] a, double[,] am, int n)
        {    
            int i, j; double dt;
            dt = Det(n,ref a);
            for (i = 0; i < n; i++)
            {
                for (j = 0; j < n; j++)
                    am[i, j] = AlgDop(n,a,j,i) / dt;
            }
        }
        static void Main(string[] args)
        {
            Console.Write("Введите размерность квадратной матрицы m=");
            int n = int.Parse(Console.ReadLine());
            double[,] a;
            double[,] Obr= new double[n,n];
            Vxod(out a, n);
            Print(a, n);
            ObrMatr(a,Obr, n);
            Console.WriteLine();
            Print(Obr, n);
        }
    }
}

Мой код: Но я его по другому писал и с косиками , мб что полезным будет.

C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
using System;
using System.IO;
 
namespace ConsoleApplication72
{
    static class Program
    {
        public static void Main(string[] args)
        {
            int createType = 0;
            const int dimension = 3;
            Console.Title = "Решение систем ленейных алгебраических уравнений";
            Console.WriteLine("Введите метод создания матрицы.\n1:Ввод с клавиатуры.\n2:Ввод с файла" +
                "\n3:Инициализация(матрица будет задана программой)");
            createType = Convert.ToInt32(Console.ReadLine());
            MyMatrix m = new MyMatrix(createType, dimension, dimension);
            MyMatrix1 qq = new MyMatrix1(m);
            m.Opredelitel();
            m.Transponirovanie();
            m.Soyznaya(); m.Obratnaya();
            m.rezmt();
            qq.prov();
            Console.ReadKey();
        }
    }
    public class MyMatrix
    {
        public MyMatrix()
        {
            //myArray = new double[3,3] {{0,0,0},{0,0,0},{0,0,0}};
        }
        public MyMatrix(int rows, double dat, int columns, int Rows) { rows = 0; dat = 0; columns = 0; Rows = 0; }
        int rows;
        public double[,] myArray;
        double[,] Y;
        public double[] C = new double[3];
        public double dat;
        int columns;
        public int Rows
        {
            get { return rows; }
            set { rows = value; }
        }
 
        public int Columns
        {
            get { return columns; }
            set { columns = value; }
        }
        public MyMatrix(int createType, int rows, int columns)
        {
            Rows = rows;
            Columns = columns;
            dat = 0;
            myArray = new double[rows, columns];
            Y = new double[rows, columns];
            switch (createType)
            {
                case 1:
                    getMatrixFromKeyboard(rows, columns);
                    printMyMatrix();
                    break;
                case 2:
                    getMatrixFromFile(rows, columns);
                    printMyMatrix();
                    break;
                case 3:
                    getMatrixFromInitialization(rows, columns);
                    printMyMatrix();
                    break;
                default:
                    Console.WriteLine("Вы указали не действительный метод создания матрицы");
                    break;
            }
        }
 
        private void getMatrixFromInitialization(int rows, int columns)
        {
            myArray = new double[3, 3] { { 7, 4, 1 }, { 2, 8, 7 }, { 3, 9, 7 } };
        }
 
        private void getMatrixFromFile(int rows, int columns)
        {
            string sp = "";
            string line = "";
            string myLine = "";
            do
            {
 
                Console.WriteLine("Введите путь к файлу");
                sp = Console.ReadLine();
            }
            while (!File.Exists(sp));
 
 
            StreamReader file = new StreamReader(sp);
 
            while ((line = file.ReadLine()) != null)
            {
                myLine += line;
            }
            file.Close();
 
            for (int i = 0; i < columns; i++)
            {
                for (int j = 0; j < rows; j++)
                {
                    string s = myLine.Split(new char[] { ' ', ';' })[i * rows + j];
                    myArray[i, j] = Convert.ToDouble(s);
                }
            }
        }
        private void getMatrixFromKeyboard(int rows, int columns)
        {
 
            int r = 0;
            string line = "";
            string myLine = "";
            do
            {
                Console.WriteLine("Введите {0}-ую строку матрицы:", r + 1);
                line = Console.ReadLine() + ";";
                myLine += line;
                r++;
            } while (r < rows);
 
 
            for (int i = 0; i < columns; i++)
                for (int j = 0; j < rows; j++)
                {
                    myArray[i, j] = Convert.ToDouble(myLine.Split(new char[] { ' ', ';' })[i * rows + j]);
                }
        }
        private void printMyMatrix()
        {
            Console.WriteLine("Имеем матрицу 3x3");
            for (int i = 0; i < columns; i++)
            {
                for (int j = 0; j < rows; j++)
                {
                    Console.Write("{0:F1} ", myArray[i, j]);
 
                }
                Console.WriteLine();
            }
        }
        public double Opredelitel()
        {
            dat = myArray[0, 0] * myArray[1, 1] * myArray[2, 2] - myArray[0, 0] * myArray[1, 2] * myArray[2, 1]
                - myArray[0, 1] * myArray[1, 0] * myArray[2, 2] + myArray[0, 1] * myArray[1, 2] * myArray[2, 0] +
                    myArray[0, 2] * myArray[1, 0] * myArray[2, 1] - myArray[0, 2] * myArray[1, 1] * myArray[2, 0];
            if (dat == 0)
            {
                Console.WriteLine("dat=0 => Матрица не имеет единственного решения.Матрица является выражденной.");
                Console.ReadKey();
            }
            else
                Console.WriteLine("\nОпределитель матрицы={0:f2}", dat);
            Console.Write("----------------\n");
            return dat;
        }
        public double[,] Transponirovanie()
        {
            Console.WriteLine("Транспонированая матрица:\n");
 
            for (int i = 0; i < 3; i++)
            {
                for (int j = 0; j < 3; j++)
                {
                    Console.Write("{0:f1} ", myArray[j, i]);
                }
                Console.WriteLine();
            }
            return myArray;
        }
        public double[,] Soyznaya()
        {
            for (int i = 0; i < 3; i++)
            {
                for (int j = 0; j < 3; j++)
                {
                    Y[0, 0] = myArray[1, 1] * myArray[2, 2] - myArray[2, 1] * myArray[1, 2];
 
                    Y[0, 1] = myArray[1, 2] * myArray[2, 0] - myArray[1, 0] * myArray[2, 2];
 
                    Y[0, 2] = myArray[1, 0] * myArray[2, 1] - myArray[1, 1] * myArray[2, 0];
 
                    Y[1, 0] = myArray[0, 2] * myArray[2, 1] - myArray[0, 1] * myArray[2, 2];
 
                    Y[1, 1] = myArray[0, 0] * myArray[2, 2] - myArray[2, 0] * myArray[0, 2];
 
                    Y[1, 2] = myArray[2, 0] * myArray[0, 1] - myArray[0, 0] * myArray[2, 1];
 
                    Y[2, 0] = myArray[0, 1] * myArray[1, 2] - myArray[1, 1] * myArray[0, 2];
 
                    Y[2, 1] = myArray[0, 2] * myArray[1, 0] - myArray[0, 0] * myArray[1, 2];
 
                    Y[2, 2] = myArray[0, 0] * myArray[1, 1] - myArray[0, 1] * myArray[1, 0];
                }
            }
            Console.WriteLine("\nСоюзная матрица");
            for (int i = 0; i < 3; i++)
            {
                Console.WriteLine();
                for (int j = 0; j < 3; j++)
                {
                    Console.Write("{0:f1} ", Y[i, j]);
                }
            }
            return Y;
        }
        public double[,] Obratnaya()
        {
            //double[,] obr = new double[3, 3];
            for (int i = 0; i < 3; i++)
            {
                for (int j = 0; j < 3; j++)
                {
                    Y[i, j] = Y[i, j] / dat;
                }
            }
            Console.WriteLine("\nОбратная матрица");
            for (int i = 0; i < 3; i++)
            {
                Console.WriteLine();
                for (int j = 0; j < 3; j++)
                {
                    Console.Write(" {0:f3}", Y[j, i]);
                }
            }
 
            return Y;
        }
        public double[] rezmt()
        {
            Console.WriteLine("\nВведите вектор свободных членов (3 числа)");
            string inputData = Console.ReadLine();
            string[] myarr = inputData.Split(new char[] { ' ' });
            //Вектор на который нужно умножать матрицу
            double[] B = new double[3];
 
            for (int k = 0; k < 3; k++)
            {
                B[k] = Convert.ToDouble(myarr[k]);
            }
            //Результат умножения матрицы на вектор
            double[,] obr = new double[3, 3];
            for (int i = 0; i < 3; i++)
                for (int j = 0; j < 3; j++)
                    obr[i, j] = Y[j, i];
 
            C[0] = obr[0, 0] * B[0] + obr[0, 1] * B[1] + obr[0, 2] * B[2];
            C[1] = obr[1, 0] * B[0] + obr[1, 1] * B[1] + obr[1, 2] * B[2];
            C[2] = obr[2, 0] * B[0] + obr[2, 1] * B[1] + obr[2, 2] * B[2];
            Console.WriteLine("\nРешение системы 3х3");
            Console.WriteLine("x1={0:f1}  x2={1:f1}  x3={2:f1} ", C[0], C[1], C[2]);
            return C;
        }
    }
    public class MyMatrix1 : MyMatrix
    {
 
        public MyMatrix1(MyMatrix m)
        {
            C = m.C;
            myArray = m.myArray;
        }
        public double[] prov()
        {
            double[] x = new double[3];
            Console.WriteLine("\nПодстановка в исходную матрицу для проверки решения");
            x[0] = myArray[0, 0] * C[0] + myArray[0, 1] * C[1] + myArray[0, 2] * C[2];
            x[1] = myArray[1, 0] * C[0] + myArray[1, 1] * C[1] + myArray[1, 2] * C[2];
            x[2] = myArray[2, 0] * C[0] + myArray[2, 1] * C[1] + myArray[2, 2] * C[2];
            Console.WriteLine("1 строка = {0:f1}  1 строка = {1:f1}     3 строка = {2:f1} ", x[0], x[1], x[2]);
            return x;
        }
    }
}

Psilon · 03.06.2013, 07:58

Уберите этот ужас... Интересно, что бы вы сделали, если бы матрица была не 3х3 а 10х10

@by-ga-ga · 03.06.2013, 10:31 **[ТС]**

мне нужно рассчитать матрицу любого порядка или хотя бы до 10 а не 3....

Psilon · 03.06.2013, 13:01

Отладку запустите: правильно ли считаются определители и тд
И есть подозрения на дублированный код

Новые блоги и статьи Все статьи Все блоги /
Программный контроль заполнения реквизита табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать контроль заполнения реквизита "ПричинаСписания". . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .
Программное заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем

Psilon Master of Orion 6102 / 4958 / 905 Регистрация: 10.07.2011 Сообщений: 14,522 Записей в блоге: 5
	03.06.2013, 07:58
	Уберите этот ужас... Интересно, что бы вы сделали, если бы матрица была не 3х3 а 10х10 0

@by-ga-ga 0 / 0 / 0 Регистрация: 16.05.2013 Сообщений: 13
	03.06.2013, 10:31 [ТС]
	мне нужно рассчитать матрицу любого порядка или хотя бы до 10 а не 3.... 0

Psilon Master of Orion 6102 / 4958 / 905 Регистрация: 10.07.2011 Сообщений: 14,522 Записей в блоге: 5
	03.06.2013, 13:01
	Отладку запустите: правильно ли считаются определители и тд И есть подозрения на дублированный код 0