Форум программистов, компьютерный форум, киберфорум
Наши страницы
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.57/7: Рейтинг темы: голосов - 7, средняя оценка - 4.57
Nikita174rus
0 / 0 / 0
Регистрация: 28.03.2016
Сообщений: 8
1

Работа с матрицами

20.05.2016, 00:03. Просмотров 1316. Ответов 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
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
#include<iostream>
#include "math.h"
 
using namespace std;
 
double eps;//погрешность
 
class Matrix{
    double **a;//матрица это двумерный массив, нам будет удобней работать с квадратными матрицами, поэтому делаем их такими
    int n;//размер матрицы
public:
    Matrix(int k){
        n=k;
        a=new double *[n];
        for(int i=0;i<n;i++){
            a[i]=new double [n];
            for(int j=0;j<n;j++){
                a[i][j]=i+(j*(i/2));
            }
        }
    };//конструктор
    ~Matrix(){//деконструктор
    };
    friend const Matrix operator+(const Matrix &x, const Matrix &y);//перегрузка для сложения двух матриц
    friend const Matrix operator*( const Matrix& x, const double z);//перегрузка для умножения матрицы на число
    friend const Matrix operator*(const Matrix &x, const Matrix &y);//перегрузка для умножения матрицы на матрицу
    void transp();//функция транспонирования
    double d();//функция которая считает дискрименант
    void obr();//функция обращающая матрицу
    void Show();//функция выводит матрицу на экран
};
void Matrix::Show(){
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cout<<a[i][j]<<" ";
        }
        cout<<"\n";
    }
    cout<<"\n";
}
void Matrix::transp(){
    double **b;//создаем вспомогающую матрицу
    b=new double *[n];
    for(int i=0;i<n;i++){
        b[i]=new double [n];
        for(int j=0;j<n;j++){
            b[i][j]=a[j][i];//транспонируем матрицу
        }
    }
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            a[i][j]=b[i][j];//меняем элементы главную матрицу на транспонированую
        }
    }
    for(int i=0;i<n;i++){
        delete []b[i];
    }
    delete []b;
}
 
double Matrix::d(){
    double D=0;//сюда записываем результат
    double **b;//создаем вспомагательный массив
    b=new double *[n];
    for(int i=0;i<n;i++){
        b[i]=new double [n+n-1];//создаем удобную для просчета определителя матрицу
        for(int j=0;j<n+n-1;j++){
            b[i][j]=a[i][j%n];//заполняем матрицу b
        }
    }
    for(int i=0;i<n;i++){
        double temp=1;
        for(int j=0;j<n;j++){
            temp=temp*b[j%n][j+i];//просчитываем главные диагонали
        }
        D=D+temp;
    }
    for(int i=0;i<n;i++){
        double temp=1;
        for(int j=0;j<n;j++){
            temp=temp*b[n-1-(j%n)][j+i];//просчитываем побочные диагонали
        }
        D=D-temp;
    }
    for(int i=0;i<n;i++){
        delete []b[i];
    }
    delete []b;
    return D;
}
void Matrix::obr(){
    double r=d();
    if(abs(r)<eps){
        cout<<"obratnoi matrici net!\n";
    }else{
        Matrix M(n);//минор матрицы
        if(n>2){//если размер матрицы больше 2 то, чтобы найти минор нужно считать определители, для этого создаем новую вспомогательную матрицу
            Matrix temp(n-1);
            int ik=0,jk=0;
            for(int i=0;i<n;i++){//циклы для 
                for(int j=0;j<n;j++){//просчета минора
                    for(int k=0;k<n;k++){//циклы для
                        for(int l=0;l<n;l++){//занесения чисел в temp
                            if(k!=i&&l!=j){//если элемент не из вычеркнутого строки или столбца
                                temp.a[ik][jk]=a[k][l];//занести в temp
                                if(jk==n-2){//перенос места под элемент в temp-е
                                    ik++;
                                    jk=0;
                                }else jk++;
                                if(ik==n-1)break;
                            }
                        }
                    }
                    M.a[i][j]=temp.d()*((-1)*(i+j+2));//записываем минор для элемента [i][j]
                }
            }
        }else{
            M.a[0][0]=a[1][1];//иначе матрица
            M.a[1][0]=a[0][1]*(-1);//2Х2 и ее
            M.a[0][1]=a[1][0]*(-1);//можно сделать
            M.a[1][1]=a[0][0];//вручную
        }
        M.transp();//транспонируем
        Matrix v(n);
        v=M*(double(1)/r);
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++){
                a[i][j]=v.a[i][j];
            }
        }
    }
}
 
const Matrix operator+(const Matrix &x, const Matrix &y){
    if(x.n!=y.n)cout<<"matrici ne ravni!\n";
    else{
        Matrix b(x.n);
        for(int i=0;i<x.n;i++){
            for(int j=0;j<x.n;j++){
                b.a[i][j]=x.a[i][j]+y.a[i][j];
            }
        }
        return b;
    }
}
const Matrix operator*( const Matrix& x, const double z){
    Matrix b(x.n);
    for(int i=0;i<x.n;i++){
        for(int j=0;j<x.n;j++){
            b.a[i][j]=x.a[i][j]*z;
        }
    }
    return b;
}
const Matrix operator*(const Matrix &x, const Matrix &y){
    Matrix b(x.n);
    if(x.n!=y.n)cout<<"matrici ne ravni!\n";
    else{
        for(int i=0;i<x.n;i++){
            for(int j=0;j<x.n;j++){
                b.a[i][j]=0;
                for(int k=0;k<x.n;k++){
                    b.a[i][j]=b.a[i][j]+x.a[i][k]*y.a[k][j];
                }
            }
        }
        return b;
    }
}
int main(){
    int n1, n2;
    cout<<"vvedite razmer pervoi matrici: ";
    cin>>n1;
    cout<<"vvedite pogreshnost: ";
    cin>>eps;
    Matrix Mat1(n1);
    cout<<"vvedite razmer vtoroi matrici: ";
    cin>>n2;
    Matrix Mat2(n2);
    Mat1.Show();
    Mat2.Show();
    Mat2=Mat1;
    Mat1.obr();
    Mat1.Show();
    Mat1.transp();
    Mat1.Show();
    Matrix v(n1);
    v=Mat1*Mat2;
    v.Show();
    v=Mat1+Mat2;
    v.Show();
    v=Mat1*0.1;
    v.Show();
    return 0;
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.05.2016, 00:03
Ответы с готовыми решениями:

Работа с матрицами
Мне необходимо заполнить матрицу размера 9 на 9 случайными числами и повернуть...

Работа с for и матрицами
Помогите пожалуйста, есть задача: Дана матрица 8х8, нужно поделить каждый...

Работа с матрицами
Всем привет. У меня возникла не большая проблема. Никак не могу разобраться со...

Работа с двухмерными матрицами
Дана квадратная матрица целых чисел размером NxM. N - нечетное число от 5 до...

Работа с файлами и матрицами
Всем здравствуйте! У меня есть проблема. Имеется такое задание: В заданной...

1
Ilot
Эксперт С++
1832 / 1190 / 342
Регистрация: 16.05.2013
Сообщений: 3,139
Записей в блоге: 5
Завершенные тесты: 1
20.05.2016, 12:17 2
В чем вопрос?
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.05.2016, 12:17

Матрицы. Работа с матрицами
Дана матрица А размера 8*6. Найти элемент, котоpый встречается в массиве чаще...

работа с матрицами в файлах
Люди, помогите сделать задачи, очень прошу на языке С Задача 1. Дан файл...

Работа с квадратными матрицами
Помогите решить пожалуйста! 1)Дана матрица размера M ×N (N —четное...


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

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru