Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.55/11: Рейтинг темы: голосов - 11, средняя оценка - 4.55
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93

Функция time для вычисления скорости нахождения Определителя.Dev-C++7.3.1.3

27.03.2011, 14:16. Показов 2269. Ответов 17
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Добрый день.
Помогите немного подкорректировать код.
В main'е switch'ом выбирается 3 пользовательских варианта:
1.Вычисление определителя матрицы;
2.Вычисление корней системы линейных уравнений;
3.0-для выхода из этой безнадежной программы =)

С вычислением корней системы линейных уравнений проблема.
Определитель вроде сносно отрабатывает.
Пытаюсь "засечь" скорость выполнения алгоритма нахождения определителя.
На выходе - стабильно "0".

Код:
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
#include <cstdlib>
#include <iostream>
#include <time.h>
#include <stdio.h>
#include <vector>
using namespace std;
void showMenu();
void showInputMatrix(std::vector< std::vector<double> > sq_matrix);
void showInputMatrix_Vector(std::vector< std::vector<double> > sl_matrix, std::vector<double> lf);
                    
double determinant(std::vector< std::vector<double> >& matrix);
void slau_funct(std::vector< std::vector<double> > sl_matrix, std::vector<double> lf, std::vector<double> lx);
int main(int argc, char *argv[])
{
    
    int choice;
    double dtm;
    //clock_t time;
         
   showMenu();
   cin>>choice;
   
   while(choice != 0){
   switch(choice)
   {
     case 1:{
          vector<vector<double> > square_matrix; 
          vector<double> line;
          int n; 
          cout<<"Vvedite razmernost' matrix:\n";
          cin>>n;
          double temp;
          cout<<"Vvedite znachenie:\n";
          for(int i=0;i<n;i++)
          {
            line.clear();
             for(int j=0;j<n;j++)
             {
               cin>>temp;
               line.push_back(temp);                  
                cout<<"Vvedite znachenie:\n";
             }
      
            square_matrix.push_back(line);
          }
    
          showInputMatrix(square_matrix);
          //time = clock();
          dtm = determinant(square_matrix);
          //time = clock() - time;
          cout<<"Determinant = ";
          dec(cout);
          cout<<dtm<<endl;
          //cout<<"time: "<<time/CLOCKS_PER_SEC<<endl;
     break;}//end case 1
    
     case 2:{
          vector<vector<double> > slau_matrix; 
          vector<double> line2;
          vector<double> line_x;
          vector<double> line_f;    // *****slau_matrix * line_x = line_f*****
    
        int m,i2,j2;
        double temp2;
        double temp_f;
        cout<<"Vvedite razmernost' matrix:\n";
        cin>>m;
        cout<<"Vvedite znacheniya v matrix:\n";
           for( i2=0;i2<m;i2++)
           {
             line2.clear();             //ochishaem vector line
             line_x.clear();           //ochishaem vector line_x
             for(j2=0;j2<m;j2++)
             {
                cin>>temp2;
                line2.push_back(temp2);                  
                cout<<"Vvedite znacheniya v matrix:\n";
             }
      
            slau_matrix.push_back(line2);    //sozdaem 2merny massiv
           }
    
            cout<<"Vvedite znacheniya v vector_f (znachenia dlya funczcii_f):\n";
            for(i2=0;i2<m;i2++)
            {
               cin>>temp_f;
               line_f.push_back(temp_f);      //zanosim znachenia v vector line_f
               cout<<"Vvedite znacheniya v vector_f (znachenia dlya funczcii_f):\n";
            }
    
         showInputMatrix_Vector(slau_matrix, line_f);
         slau_funct(slau_matrix, line_f, line_x);
         
        /* cout<<"Naidennie znachenia x:\n";
         for(i2=0;i2<line_x.size();i2++)
         {
                cout<<line_x[i2]<<"|";
         }*/
         
     break;}//end case 2
     
     default :
            cout<<"exit...\n";
     break;   
   }//end switch
  
   showMenu();
  cin>>choice;
}//end while 
 
 
    system("PAUSE");
    return EXIT_SUCCESS;
}
 
//################opredelenie Function's##################
void showInputMatrix(std::vector< std::vector<double> > matrix)
{    int i,j;
     cout<<"Matrix\n";
       for(i=0;i<matrix.size();i++)
       {
         for(j=0;j<matrix.size();j++)
         { 
           dec(cout);
           cout.width(11);
           cout.precision(10);
           cout<<matrix[i][j]<<"|";
         }
        cout<<endl;
       }
}
 
 
void showInputMatrix_Vector(std::vector< std::vector<double> > sl_matrix, std::vector<double> lf)
{       
      int i,j;
        cout<<"Sistema lineynih uravneniy\n";
          for(i=0;i<sl_matrix.size();i++)
          {
            for(j=0;j<sl_matrix.size();j++)
            { 
              dec(cout);
              cout.width(11);
              cout.precision(10);
              cout<<sl_matrix[i][j]<<"|";
            }
           cout<<endl;
          } 
           
           for(j=0;j<lf.size();j++)   //vivod vectora line_f
            { 
              dec(cout);
              cout.width(11);
              cout.precision(10);
              cout<<lf[j]<<"\n";
            } 
            
}
 
 
void showMenu()
 {
    cout<<"vvedite 1,2 ili 0:\n"
    "1 - Vichislyaem opredelitel matrizcy\n"
    "2 - Reshaem sistemu lineinih uravneny\n"
    "0 - Vihod\n";
 }
 
double determinant(std::vector< std::vector<double> >& matrix)  //nahodim opredelitel 
{   //function find determinant
    //
    int i=0, j=0;
    double koef = 1;
    double lastKoef = 0;
    bool checkKoefIsNull = false;
    int length = matrix.size();
    int fullCol = length;
    
    clock_t time_on;
    time_on = clock();
     for(i=0; i<length && koef; i++)
    {
        if(matrix[i][j] != 0)
        {
            if(matrix[i][j] < 0)
            {
                for(int j=0; j<fullCol; j++)
                    matrix[i][j] *= -1;
                koef *= -1;
            }
            double lastKoef = matrix[i][j];
            for(int k=0; k<fullCol; k++)
                matrix[i][k] /= lastKoef;
            koef *= lastKoef;
            for(int nextRow=i+1; nextRow<length; nextRow++)
            {
                lastKoef = matrix[nextRow][j];
                for(int nextCol=j; nextCol<fullCol; nextCol++)
                    matrix[nextRow][nextCol] -= lastKoef*matrix[i][nextCol];
            }
            ++j;
            checkKoefIsNull=false;
        }
        else
        {
            for(int k=0; k<fullCol; k++)
                std::swap(matrix[i][k],matrix[length-1][k]);
            --i;
            koef *= -1;
            if(checkKoefIsNull)
            {
                ++i;
                koef = 0;
            }
            checkKoefIsNull = true;
        }
     }//end for
     
  
     cout<<"time: "<<((double)clock() - time_on)/CLOCKS_PER_SEC<<endl;
        return koef;
}
 
 
void slau_funct(std::vector< std::vector<double> > sl_matrix, std::vector<double> lf, std::vector<double> lx)
{   //nahodim korni x1,x2,...xn  v  s.l.a.u. 
    int n = lx.size();
    int l = n - 1;
    int j;
    int m;
    int i2;
    double sum;
    lx[l] = lf[l] / sl_matrix[l][l];     // *****slau_matrix * line_x = line_f*****
    for(m = l-1; m > -1; m--)    
    {
      sum = 0;
         for(j = m+1; j < n; j++) 
         {  
          sum = sum + (sl_matrix[m][j] * lx[j]);   
         }
    
      lx[m] = ( 1 /  sl_matrix[m][m]  ) * ( lf[m]  -  sum );
      cout<<endl;
      cout<<lx[m]<<"|";
    }
    
}
Перечисление небольшой суммы на wmz гарантирую.
Вложения
Тип файла: rar opredelitel.rar (197.1 Кб, 15 просмотров)
0
Programming
Эксперт
39485 / 9562 / 3019
Регистрация: 12.04.2006
Сообщений: 41,671
Блог
27.03.2011, 14:16
Ответы с готовыми решениями:

Функция для нахождения определителя по Гауссу
Я написал функцию для нахождения определителя квадратных матриц по Гауссу, но она у меня не работает и выдает ошибку помогите исправить код...

И снова класс Матрицы с функция-член класса для нахождения определителя
Проблема в следующем испробовал все варианты, которые закомментированы, все рабочие ну некоторые с отрицательными числами неправильно...

Используя рекурсию и рекуррентную формулу для вычисления определителя матрицы написать алгоритм вычисления определителя матрицы
Помогите, пожалуйста.

17
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 14:26
Размерность матрицы какую задаёшь?
Чтобы успеть "засечь" хотя бы 1000х1000 задай.
1
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 14:37  [ТС]
Цитата Сообщение от Deviaphan Посмотреть сообщение
Размерность матрицы какую задаёшь?
Чтобы успеть "засечь" хотя бы 1000х1000 задай.
Я в ручную заношу данные)
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 14:42
Можешь на пальцах прикинуть, сколько операций выполняется для поиска определителя?
И сколько тактов процессора выполняется за 1/1000 секунды?
И как замерить промежуток времени, который меньше 1/1000 при помощи time?
И главное, какая разница, сколько это займёт времени, если временной промежуток столь маг, что его нельзя измерить.)
Если хочешь реально проверить изменение производительности используй матрицу большой размерности. На самом деле и 1000х1000 это маленькая матрица.
1
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 14:54  [ТС]
Попробовал.
Затянулся процесс...
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 15:06
100х100
0
27.03.2011, 15:07

Не по теме:

Цитата Сообщение от USHЁL Посмотреть сообщение
Я в ручную заношу данные)
И какой же способ от этого отойти? Эм... Может быть, заполнять матрицу автоматически, с использованием генератора случайных чисел?

0
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 15:25  [ТС]
Дело в том,что изначально я задумывал програму как вычисление определителя для,введенных пользователем матриц.
Со временем добавился вопрос с нахождением времени выполнения алгоритма.
Щас подправлю

Добавлено через 4 минуты
Цитата Сообщение от Deviaphan Посмотреть сообщение
100х100
100*100 задал.На выходе всё по нулям.
Code
1
2
3
4
5
6
7
8
9
10
11
12
 for(int i=0;i<n;i++)
          {
            line.clear();
             for(int j=0;j<n;j++)
             {
               //cin>>temp;
               temp = rand();
               line.push_back(temp);                  
                //cout<<"Vvedite znachenie:\n";
             }
      
            square_matrix.push_back(line);

Вывод для 100 на 100 матрицы:
time:0.031
Determinant = -inf <-----wtf?
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 15:27
Это лаба или для себя?
Если для себя, то в сети есть куча библиотек для работы с матрицами. Там в стопицот раз оптимальнее реализовано, чем ты сможешь сам написать.)

Добавлено через 29 секунд
Цитата Сообщение от USHЁL Посмотреть сообщение
100*100 задал.На выходе всё по нулям.
500х500
.)
0
27.03.2011, 15:32

Не по теме:

Цитата Сообщение от USHЁL Посмотреть сообщение
Дело в том,что изначально я задумывал програму как вычисление определителя для,введенных пользователем матриц
Да, этот факт кардинально всё меняет... Ведь компилятор не допускает больше одного изменения кода в неделю...
Это я к тому, что можно сначала сменить ввод на автоматически, а потом обратно на ручной.

1
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 15:34  [ТС]
Это лаба.
Нужно как бы и самому немного подразобраться.
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 15:35
А считает определитель правильно?
Это должен быть первый вопрос.)
0
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 15:37  [ТС]
=)
А вот мне то как раз не до шуток.
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 15:39
Цитата Сообщение от USHЁL Посмотреть сообщение
Определитель вроде сносно отрабатывае
ТАк какая тебе разница, сколь мало времени на это уходит.)
Или по заданию засечь надо? Тогда размер матрицы увеличивай. Постепенно.
0
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 15:49  [ТС]
мне потом нужен письменный отчёт составлять по этому методу.

И график чертить зависимости размерности от времени,хотябы для 3-4 случаев.
Это реально как Вы считаете или надомной поглумился инопланетный препод(прошу прошения за ненорматив)?=)
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 16:00
Чтобы засечь время матрица должна быть относительно большая.
График чертить смысл очень есть - зависимость нелинейная.
Размерность матрицы увеличивай, чтобы и измерить успеть и не сидеть пол дня.
1
 Аватар для USHЁL
0 / 0 / 1
Регистрация: 15.11.2009
Сообщений: 93
27.03.2011, 16:11  [ТС]
Благодаря Вам понял,что дело в размерности.
Не понятен вывод после ввода матрицы 1000 на 1000
На выходе Detrminant = -inf
Что это означает?

Для матрицы 1000*1000
time = 31.797
Для матрицы 1111*1111
time = 43.547
Determinent = inf

Как Вы считаете код сносно реализован или нужно менять?

main:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
for(int i=0;i<n;i++)
          {
            line.clear();
             for(int j=0;j<n;j++)
             {
               //cin>>temp;
               temp = rand();
               line.push_back(temp);                  
                //cout<<"Vvedite znachenie:\n";
             }
      
            square_matrix.push_back(line);
          }
determinant:
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
double determinant(std::vector< std::vector<double> >& matrix)  //nahodim opredelitel 
{   //function find determinant
    //
    int i=0, j=0;
    double koef = 1;
    double lastKoef = 0;
    bool checkKoefIsNull = false;
    int length = matrix.size();
    int fullCol = length;
    
    clock_t time_on;
    time_on = clock();
     for(i=0; i<length && koef; i++)
    {
        if(matrix[i][j] != 0)
        {
            if(matrix[i][j] < 0)
            {
                for(int j=0; j<fullCol; j++)
                    matrix[i][j] *= -1;
                koef *= -1;
            }
            double lastKoef = matrix[i][j];
            for(int k=0; k<fullCol; k++)
                matrix[i][k] /= lastKoef;
            koef *= lastKoef;
            for(int nextRow=i+1; nextRow<length; nextRow++)
            {
                lastKoef = matrix[nextRow][j];
                for(int nextCol=j; nextCol<fullCol; nextCol++)
                    matrix[nextRow][nextCol] -= lastKoef*matrix[i][nextCol];
            }
            ++j;
            checkKoefIsNull=false;
        }
        else
        {
            for(int k=0; k<fullCol; k++)
                std::swap(matrix[i][k],matrix[length-1][k]);
            --i;
            koef *= -1;
            if(checkKoefIsNull)
            {
                ++i;
                koef = 0;
            }
            checkKoefIsNull = true;
        }
     }//end for
     
  
     cout<<"time: "<<((double)clock() - time_on)/CLOCKS_PER_SEC<<endl;
        return koef;
}
0
Делаю внезапно и красиво
Эксперт С++
 Аватар для Deviaphan
1313 / 1228 / 72
Регистрация: 22.03.2011
Сообщений: 3,744
27.03.2011, 17:10
Это либо переполнение, либо были не инициализированные данные. Или деление на ноль было.
1
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
inter-admin
Эксперт
29715 / 6470 / 2152
Регистрация: 06.03.2009
Сообщений: 28,500
Блог
27.03.2011, 17:10
Помогаю со студенческими работами здесь

неправильно работает программа для вычисления определителя определителя (c++)
работает вроде нормально, но иногда пишет неправильный знак определителя, причём ошибается в одном и том же...

Функция нахождения определителя матрицы с вычислением времени
Ребят, помогите доделать программу. Программа ищет определитель матрицы, данные заполняются рандомно либо от руки. Надо сделать еще функцию...

Рекурсивная функция вычисления определителя квадратной матрицы
Написать рекурсивную функцию вычисления определителя квадратной матрицы (использовать разложение по строке/столбцу).

Создать класс «матрица» с членами определитель и функция вычисления определителя
Доброго времени суток! Помогите реализовать задачу: Класс «матрица», включающий матрицу, ее определитель и функцию вычисления ...

Формула для нахождения определителя матрицы
Здравствуйте! У меня имеется небольшой вопрос: как записать формулу для нахождения определителя матрицы(ранг&gt;2) на с++? Считать...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
18
Ответ Создать тему
Новые блоги и статьи
Конвертировать закладки radiotray-ng в m3u-плейлист
damix 19.02.2026
Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI
Eddy_Em 18.02.2026
Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .
Камера Toupcam IUA500KMA
Eddy_Em 12.02.2026
Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу
zbw 12.02.2026
И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.
«Знание-Сила»
zbw 12.02.2026
«Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров
8Observer8 12.02.2026
Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image)
8Observer8 11.02.2026
Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image
8Observer8 10.02.2026
Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru