Работа с динамической памятью

@TrIX · Регистрация: 09.12.2007

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

люди выручите меня, я только недавно начал изучать язык СИ, и еще не знаю как делать задачи вот такого типа:

Дана квадратная матрица A порядка n. Получить матрицу bT(A-E)b, где E ¬– единичная матрица порядка n; элементы вектора b вычисляются по формуле:

bi=1/i*i*i+2
bT - транспонированная матрица.
E-единичная матрица.Т.е сначала считается A-E, затем b*(A-E) и потом полученный результат умножают на транпонированную матрицу b

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

помогите пожалуйста

@dioxidin · 10.12.2007, 11:01

Про умножение матриц где-то в этом форуме было. Выражения надеюсь знаете... Скоро со своей прогой разделаюсь и выложу функцию умножения матриц. Если дело очень-очень срочное, то пишите...
Про динамическое выделение памяти (а не выделение динамической памяти

) читайте тут: Передача массива в функцию

PS: Люди пользуйтесь поиском, пожалуйста...

@TrIX · 10.12.2007, 12:43 **[ТС]**

вообщем я попыхтел тут на паре и вот что у меня получилось на данный момент:

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
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include<windows.h>
char* rus(const char* text);
 char bufrus[256];
        char* rus(const char* text){
        CharToOem(text, bufrus);
        return bufrus;
        }
 
FILE *in;
 
int main()
{
 int *n=new int();
 int temp;
        int kol;
 char name[256];
  do{
  clrscr();
 printf(rus("Введите имя файла:"));
 scanf("%s",name);
 do{
 if ((in = fopen(name,"rt")) == NULL)
 {
  printf(rus("Невозможно открыть файл !!!\n"));
  continue;
 } else
 {
  kol=0;
  float t;
    do{
  while(fscanf(in,"%f",&t) == 1)
  kol++;
    }while(!feof(in));
  rewind(in);
 }
 *n=(sqrt(kol));
 if (kol!=(*n)*(*n))
 {
        printf(rus("Матрица не квадратная!!!\n"));
        continue;
 }
 int i=0,j=0;
 float *b=new float[(*n)*(*n)];
 float *c=new float[(*n)*(*n)];
 float *a=new float[(*n)*(*n)];
 float *d=new float[(*n)*(*n)];
 float *x=new float[(*n)*(*n)];
       
 int st=0;
 do{
  j=0;
  while(fscanf(in,"%f",(a+i*(*n)+j)) == 1)
  {
   j++;
   if (j==*n)
   {
    i++;
    j=0;
    fscanf(in,"\n");
   }
  }
 }while (!feof(in));
 if (st!=*n && j!=0)
  {
   printf(rus("Файл задан неверно!!!\n"));
   continue;
  }
                if(*n==0)
                { printf(rus("Файл пуст.\n"));
                continue;}
                clrscr();
        printf(rus("Чтение файла %s прошло успешно.\n\n"),name);
 printf(rus(" --===Матрица размера %d x %d===--\n\n"),*n,*n);
 printf(rus("-Матрица А:\n"));
        printf("------------------\n");
 for (i=0;i<(*n);i++)
 {
  for (j=0;j<(*n);j++)
  printf("    %.2lf",*(a+i*(*n)+j));
  printf("\n");
 }
 fclose(in);
//-----------------------------------------
for (i=0;i<(*n);i++)
for (j=0;j<(*n);j++)
if (i==j)
*(b+i*(*n)+j)=1; else
if (i<j)
*(b+i*(*n)+j)=0;
printf(rus("\n-Единичная матрица:\n"));
        printf("------------------\n");
 for (i=0;i<(*n);i++)
 {
  for (j=0;j<(*n);j++)
  printf("    %.2lf",*(b+i*(*n)+j));
  printf("\n");
 }
//---------------------------------------------
int z;
for (i=0;i<(*n);i++)
{
z=i;
*(d+i)=1.0/((z+1)*(z+1)*(z+1)+2); }
printf(rus("\n-Вектор b:\n"));
printf("------------------\n");
for (i=0;i<(*n);i++)
printf("    %.2lf",*(d+i));
printf("\n");
//------------------------------------------------
int k;
for (i=0;i<(*n);i++)
for (j=0;j<(*n);j++)
*(c+i*(*n)+j)=((*(a+i*(*n)+j))-(*(b+i*(*n)+j)));
 printf(rus("\n-Матрица A-E:\n"));
        printf("------------------\n");
 for (i=0;i<(*n);i++)
 {
  for (j=0;j<(*n);j++)
  printf("    %.2lf",*(c+i*(*n)+j));
  printf("\n");
 }
//------------------------------------------
for (i=0;i<(*n);i++)
for (j=0;j<(*n);j++){
*(x+i*(*n)+j)=0.0;
 
*(x+i*(*n)+j)=*(x+i*(*n)+j)+(*(c+i*(*n)+j) * (*(d+i)));}
printf(rus("-Перемножение b*(A-E):\n"));
        printf("------------------\n");
 for (i=0;i<(*n);i++)
 {
  for (j=0;j<(*n);j++)
  printf("    %.2lf",*(x+i*(*n)+j));
  printf("\n");
 }
//-------------------------------------------
}
        while(!temp);
     printf(rus("\nВы хотите продолжить? 1=Да 2=Нет\n"));
     scanf("%i",&temp);
    }
    while(temp==1);
 return 0;
        }

осталось только перемножить (A-E) c b и полученный результат умножить на транспонированный вектор b.
а как его транспонировать?если бы было например M[i][j] то транспонированная - M[j][i].А если есть b[i] то как?не могу понять

@dioxidin · 10.12.2007, 13:34

В прикрепленном файле TMatrix находится описание класса для динамического создания двумерных массивов. Посмотри - может это тебе лучше подойдет...
Можешь просто скопировать файлы в папку проекта и подключать к заголовкам модулей заголовок модуля в прикрепленном архиве...

Про транспонирование... Все просто...
Объявляешь двумерный массив - напр., 4 строки и 1 столб - и создаешь двойной цикл. В одном из циклов индекс просто будет от 0 до 0... и делаешь все также как и с обычными двумер. массивами.

Касаемо умножения матриц - вот моя функция, использующая класс TMatrix:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
void KalmanThread::MMulty(TMatrix *Mat1, TMatrix *Mat2, TMatrix *REZ)
{
    int k = Mat1->GetRowCount()-1;
    int l = Mat1->GetColCount()-1;
    int m = Mat2->GetColCount()-1;
 
    for (int i1 = 0; i1 <= k; i1++) {
        for (int i2 = 0; i2 <= m; i2++) {
            for (int i3 = 0; i3 <= l; i3++) {
                if ((abs(Mat1->GetM(i1,i3))>=1E-5)||(abs(Mat2->GetM(i3,i2))>=1E-5)) {
                    float t;
                    t = REZ->GetM(i1,i2) + Mat1->GetM(i1,i3)*Mat2->GetM(i3,i2);
                    REZ->SetM(t,i1,i2);
                }
            }
        }
    }
}

@TrIX · 10.12.2007, 19:29 **[ТС]**

большое спасибо.Код дописал, все работает отлично

Новые блоги и статьи Все статьи Все блоги /
Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .	Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .
Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .	Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизитов табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/

@TrIX 2 / 2 / 0 Регистрация: 09.12.2007 Сообщений: 26

	Работа с динамической памятью 09.12.2007, 20:01. Показов 4658. Ответов 4 Метки нет (Все метки) люди выручите меня, я только недавно начал изучать язык СИ, и еще не знаю как делать задачи вот такого типа: Дана квадратная матрица A порядка n. Получить матрицу bT(A-E)b, где E ¬– единичная матрица порядка n; элементы вектора b вычисляются по формуле: bi=1/iii+2 bT - транспонированная матрица. E-единичная матрица.Т.е сначала считается A-E, затем b*(A-E) и потом полученный результат умножают на транпонированную матрицу b Все используемые массивы и переменные должны быть реализованы как динамические структуры. Считывание исходных данных должно осуществляться из файла. Доступ к элементам массива должен осуществляться через указатели. Программа сама должна определять размеры матрицы, исходя из содержимого файла, и выдавать сообщение, если матрица не удовлетворяет условиям задачи. помогите пожалуйста 0

@dioxidin 132 / 99 / 11 Регистрация: 21.11.2007 Сообщений: 544
	10.12.2007, 11:01
	Про умножение матриц где-то в этом форуме было. Выражения надеюсь знаете... Скоро со своей прогой разделаюсь и выложу функцию умножения матриц. Если дело очень-очень срочное, то пишите... Про динамическое выделение памяти (а не выделение динамической памяти ) читайте тут: Передача массива в функцию PS: Люди пользуйтесь поиском, пожалуйста... 0

@TrIX 2 / 2 / 0 Регистрация: 09.12.2007 Сообщений: 26
	10.12.2007, 19:29 [ТС]
	большое спасибо.Код дописал, все работает отлично 0