Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Vitalyalya
3 / 3 / 0
Регистрация: 07.04.2013
Сообщений: 13
#1

Матрицы не хотят умножаться (внутри практически рабочий код!) - C++

07.04.2013, 21:25. Просмотров 304. Ответов 2
Метки нет (Все метки)

Добрый день! Есть задание:
Задан массив C(n) . Сформировать матрицу A(n,n),как произведение массивов C и СT и матрицу B(n,n), элементы которой вычислить по формуле:
Aij=Bij/max|A|. Решить матричное уравнение XA=3B-E, где E – единичная max[A].
Массив нужно вводить из файла.
Ввод из файла работает, массивы перемножаются, а вот при создании матрицы B появляются неприятности. В общем код рабочий, но в отдельных аспектах появляются проблемы. Я не прошу решить (ибо уже практически все решено ), но прошу помочь с решением.

Вот сам код:
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
#include <iostream>
#include <math.h>
#include <fstream>
using namespace std;
 
double determinant(double **matrica_a, int n)
{ int i,j,k,r; double c,M,max,det=1, **a;
a=new double *[n];
for(i=0;i<n;i++) a[i]=new double[n];
for(i=0;i<n;i++)
 for(j=0;j<n;j++)
a[i][j]=matrica_a[i][j];
for(k=0;k<n;k++)
{
max=fabs(a[k][k]); r=k;
for(i=k+1;i<n;i++)
if (fabs(a[i][k])>max)
{ max=fabs(a[i][k]); r=i; }
if (r!=k) det=-det;
for(j=0;j<n;j++)
{c=a[k][j]; a[k][j]=a[r][j]; a[r][j]=c; }
for(i=k+1;i<n;i++)
for(M=a[i][k]/a[k][k],j=k;j<n;j++)
a[i][j]-=M*a[k][j];
}
for(i=0;i<n;i++)
det*=a[i][i];
return det;
}
 
int main()
{
    int n,i,j,max=0;
    cout<<"Ввод количества элементов массива \n";
    cin>>n;
    cout<<endl;
    double q;
    double **a;
    double **b;
    int *c=new int [n];
    a=new double *[n];
    for (i=0;i<n;i++)
    a[i]=new double (n);
    fstream fi;
    fi.open("/home/vitalya/abcd.txt", ios::in);
    if (!fi)
    cout<< "ошибка!\n"; 
    for(i=0;i<n;i++)
    fi>>c[i];
    cout<<"Чудо-массив \n";
    for(i=0;i<n;i++)
    cout<<c[i]<<' ';
    cout<<endl;
    fi.close();
    cout<<"Транспонированный массив \n";
    for (i=0;i<n;i++)
    cout<<c[i]<<endl;
    cout<<endl;
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    a[i][j]=c[i]*c[j];
    cout<<"Произведение (Матрица A[i][j]) \n";
    for (i=0;i<n;cout<<endl,i++)
    for (j=0;j<n;cout<<"\t",j++)
    cout<<a[i][j];
    cout<<endl;
    max=a[0][0];
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    if (max<a[i][j])
    max=a[i][j];
    b=new double *[n];
    for (i=0;i<n;i++)
    b[i]=new double [n];
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    b[i][j]=max/a[i][j];
    cout<<"Матрица, состоящая из A(i,j)/max|A| (Матрица B[i][j]) \n";
    for (i=0;i<n;cout<<endl,i++)
    for (j=0;j<n;j++)
    cout<<b[i][j]<<"\t";
    cout<<endl;
    for (i=0;i<n;i++)
    for (j=0;j<n;j++)
    if (i==j)
    b[i][j]=3*b[i][j]-1;
    else 
    b[i][j]=3*b[i][j];
    cout<<"3*В-Е \n";
    for (i=0;i<n;cout<<endl,i++)
    for (j=0;j<n;cout<<"\t",j++)
    cout<<b[i][j];
    cout<<endl;
    q=determinant(a,n);
    cout<<"Определитель равен "<<q<<", то есть нулю";
}
Добавлено через 2 часа 26 минут
Никто не поможет..?
Я же не прошу решить за меня, а просто помочь
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
07.04.2013, 21:25
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Матрицы не хотят умножаться (внутри практически рабочий код!) (C++):

У меня есть файл и внутри этого файла есть код, я хочу его открыть в IDLE питон так, чтобы внутри код вывести - Python
У меня есть файл и внутри этого файла есть код, я хочу его открыть в IDLE питоне так, чтобы код можно было выполнить в idle: Файл с таким...

DirectX 9.0 Ищу нормальный рабочий код для вывода текста. В книгах код не работает! - DirectX
Подскажите нормальный рабочий код для вывода текста на экран. Сложность в том, что при разных версиях DirectX разные коды. В интернете...

Вирус устанавливает программы, добавляет ссылки на рабочий стол (логи внутри) - Удаление вирусов
Здравствуйте! Хватанула вирус. При запуске системы открывался браузер с workno.ru, поиск перенаправлялся на рамблер через какой-то...

Вычислить N степень каждого элемента этой матрицы( не могу понять что от меня хотят ) - Pascal ABC
Даны натуральное N и элементы квадратной вещественной матрицы А пятого порядка. Вычислить N-ю степень каждого элемента этой матрицы...

.NET 4.x Внутри программы на C# выполнить JS код и передать результаты в основной код C# - C#
Добрый день. Требуется внутри программы на C# выполнить JS код и передать результаты в основной код C# , где с ними будут производиться...

Код читаю, а сам пишу с трудом. Код внутри - Web
Много учил, читал. На бесплатных обучалках проходил курсы. В общем и целом, читать код я могу, понимаю(я не хвастаюсь, я многого не знаю)....

2
Vitalyalya
3 / 3 / 0
Регистрация: 07.04.2013
Сообщений: 13
09.04.2013, 17:08  [ТС] #2
Все еще актуально.
Это я так, на всякий случай.
0
-=ЮрА=-
Заблокирован
Автор FAQ
09.04.2013, 17:55 #3
Vitalyalya, касательно детерминанта
http://www.cyberforum.ru/faq/thread436065-page2.html#post2471796 и все посты ниже.
На счёт умножения матриц сюда
Как умножить матрицы
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
09.04.2013, 17:55
Привет! Вот еще темы с ответами:

Рабочий код - jQuery
Здравствуйте, не работает элементарная вещь, в чем может быть проблема?на jsfiddle проверял работает,тут нет &lt;head &gt; &lt;link...

Проверьте ,рабочий ли код - Turbo Pascal
uses crt; var f, g : text; a, b : array of real; i, m1, m2 : byte; begin

Переделать рабочий код - C++
Нужно изменить код программы. Она рабочая, но нужно, чтобы была прозрачной и отличалась от оригинала: #include &lt;iostream.h&gt; #include...

Допольнить рабочий код - Pascal
Найти количество четных элементов одномерного массива целых до первого встреченного числа, равного наперед заданному числу а. Размерность...


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

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

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