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

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

20.05.2016, 00:03. Показов 2146. Ответов 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

Помощь в написании контрольных, курсовых и дипломных работ здесь.

Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
20.05.2016, 00:03
Ответы с готовыми решениями:

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

Работа с матрицами. С+
Здравствуйте! Помогите,пожалуйста,выполнить задание

Работа с матрицами
Помогите написать программу что бы можно было выбирать размерность матрицы(квадратная) и...

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

1
Эксперт по математике/физикеЭксперт С++
1932 / 1270 / 365
Регистрация: 16.05.2013
Сообщений: 3,317
Записей в блоге: 6
20.05.2016, 12:17 2
В чем вопрос?
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
20.05.2016, 12:17

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

Работа с матрицами
Не могу понять где ошибка, помогите пожалуйста!!package matrix; /** * * @author kat */...

Работа с матрицами
Доброго времени суток! Может кто-то сможет подсказать по поводу следующего: Задана матрица...

работа с матрицами
здравствуйте! решите примеры пожалуйста. А то я совсем не разбираюсь в delphi. 1.Даны...


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

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

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