Операции над статическими матрицами

@ik4 · Регистрация: 12.12.2020

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

Добрый вечер! Разработала код для решения выражения прикрепленного ниже (фото). Работаю со статическими матрицами. Объявила их одномерными, так как при умножении на V (матрица строка) то, что в скобках должно являться вектором столбцом, соответственно, одномерным массивом. Но не работает. (Результат работы программы прикреплю. Подскажите, где ошибаюсь. Некоторые функции отличаются от классического языка С по той причине, что работаю в среде разработки VS 2019 и там свои правила. Заранее спасибо!

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
#include "iostream"
#include "locale.h"
#include "windows.h"
#include "math.h"
#include <conio.h>
int main()
{
    setlocale(LC_ALL, "RUS");
    double A[100], B[100], A_t[100], V[100];
    int i, j, x, y, q;
    /*Размер матрицы А*/
    int n, m;
    printf("Введите размер матрицы A: ");
    scanf_s("%d", &n);
    printf("Матрица A:\n");
    for (i = 0; i < n; i++) {
            A[i] = -10 + 20.0 * rand() / RAND_MAX;
            printf("%8.2f", A[i]);
        printf("\n");
    }
    /*Транспонирование матрицы A*/
    double tr_A[100];
    printf("Транспон. матрица A:\n");
    for (i = 0; i < n; i++) {
            tr_A[i] = A[i];
            printf("%8.2f", tr_A[i]);
        printf("\n");
    }
    /*Размер матрицы B*/
    printf("Введите размер матрицы B: ");
    scanf_s("%d", &x);
    printf("Матрица B:\n");
    for (i = 0; i < x; i++) {
            B[i] = -10 + 20.0 * rand() / RAND_MAX;
            printf("%8.2f", B[i]);
        printf("\n");
    }
    /*Размер матрицы V*/
    int k;
    printf("Введите размер матрицы C: ");
    scanf_s("%d", &k);
    printf("Матрица C:\n");
    for (i = 0; i < k; i++) {
            V[i] = -10 + 20.0 * rand() / RAND_MAX;
            printf("%8.2f", V[i]);
        printf("\n");
    }
    /*Вычисление выражения*/
    double rez=0.0;
    printf("(A_t-B)*V:\n");
    for (i = 0; i < n; i++) {
            rez = (A_t[i] - B[i]) * V[i];
            printf("%8.2f", rez);
        printf("\n");
    }
    _getch();
    return 0;
}

@Annemesski · 14.07.2021, 12:48

Catstail, зачем все эти рассуждения про вектора и матрицы, есть же простое правило о совместимости матриц: считая вектор-строку как матрицу 1xN, а вектор-столбец как матрицу Nx1 и проверяя матрицы на совместимость выдавать либо результат, либо сообщение о невозможности выполнения умножения.

@Volga_ · 14.07.2021, 13:17

Annemesski, согласен с вами. И вот код может быть так:

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
#include <locale.h>
#include <windows.h>
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <time.h>
 
int main()
{
    setlocale(LC_ALL, "RUS");
    srand(time(0));
 
    double A[100], B[100], V[100], rez[100][100];
    int i, j;
 
    /*Вектор строки V*/
    int n;
    printf("Введите количество элементов вектора строки V: ");
    scanf_s("%d", &n);
    if (n > 100)
    {
        printf("Ошибка ! Нужен не больше 100...");
        return 0;
    }
    printf("Вектор строки V (размер 1x%d):\n", n);
    for (i = 0; i < n; i++) {
        V[i] = -10 + 20.0 * rand() / RAND_MAX;
        printf("%8.2f\t", V[i]);
    }
    printf("\n\n");
 
    /*Матрица А*/
    int rowA, colA;
    printf("Введите размер матрицы A:\n");
    printf("\trowA = ");
    scanf_s("%d", &rowA);
    printf("\tcolA = ");
    scanf_s("%d", &colA);
    if (rowA > 100 || colA > 100 || rowA != 1 || colA != n)
    {
        printf("Ошибка ! Нужны: rowA<=100 && colA<=100 && rowA=1 && colA=%d", n);
        return 0;
    }
    printf("Матрица A (размер %dx%d):\n", rowA, colA);
    for (i = 0; i < n; i++) {
        A[i] = -10 + 20.0 * rand() / RAND_MAX;
        printf("%8.2f\t", A[i]);
    }
    printf("\n\n");
 
    /*Транспонирование матрицы A*/
    printf("Транспон. матрица A (размер %dx%d):\n", colA, rowA);
    for (i = 0; i < n; i++) {
        printf("%8.2f\n", A[i]);
    }
    printf("\n");
 
    /*Матрица B*/
    int rowB, colB;
    printf("Введите размер матрицы B:\n");
    printf("\trowB = ");
    scanf_s("%d", &rowB);
    printf("\tcolB = ");
    scanf_s("%d", &colB);
    if (rowB > 100 || colB > 100 || rowB != n || colB != 1)
    {
        printf("Ошибка ! Нужны: rowB<=100 && colB<=100 && rowB=%d && colB=1", n);
        return 0;
    }
    printf("Матрица B (размер %dx%d):\n",rowB,colB);
    for (i = 0; i < n; i++) {
        B[i] = -10 + 20.0 * rand() / RAND_MAX;
        printf("%8.2f\n", B[i]);
    }
    printf("\n");
 
    /*Вычисление выражения*/
    printf("Результат (A_t-B)*V с размером %dx%d:\n", n, n);
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            rez[i][j] = (A[i] - B[i]) * V[j];
            printf("%8.2f\t", rez[i][j]);
        }
        printf("\n");
    }
 
    _getch();
    return 0;
}

и результат моделирования:

Code
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
Введите количество элементов вектора строки V: 4
Вектор строки V (размер 1x4):
    7,74            7,44           -8,99            6,85
 
Введите размер матрицы A:
        rowA = 1
        colA = 4
Матрица A (размер 1x4):
    6,00           -0,38            2,45            2,31
 
Транспон. матрица A (размер 4x1):
    6,00
   -0,38
    2,45
    2,31
 
Введите размер матрицы B:
        rowB = 4
        colB = 1
Матрица B (размер 4x1):
    9,45
    2,51
    5,93
   -0,42
 
Результат (A_t-B)*V с размером 4x4:
  -26,71          -25,68           31,00          -23,64
  -22,44          -21,57           26,04          -19,85
  -26,99          -25,94           31,32          -23,88
   21,12           20,30          -24,51           18,69

@Catstail · 14.07.2021, 13:41

Сообщение от Annemesski

зачем все эти рассуждения про вектора и матрицы, есть же простое правило о совместимости матриц

- я это и имел в виду. Я в чем-то ошибся?

@Volga_ · 14.07.2021, 14:19

Сообщение от Catstail

Я в чем-то ошибся?

По-моему у вас делать V*(A^T-B), но в условии задачи нужно делать (A^T-B)*V.

@Annemesski · 14.07.2021, 14:49

Сообщение от Catstail

Я в чем-то ошибся?

нет, у Вас все верно, мой комментарий к посту №14 (надо было цитату сделать)

@Catstail · 14.07.2021, 15:12

Вот вполне универсальный код. Он работает с векторами как с матрицами:

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
#include <stdio.h>
#include <stdlib.h>
 
// Дать индекс элемента [i,j] в прямоугольной 
// матрице с количеством столбцов = cols
 
int ind(int i, int j, int cols)
{
    return i*cols+j;
}
 
// Дать элемент по индексу
 
double get_el(double *matr, int cols, int i, int j)
{
    int k=ind(i,j,cols);
    return matr[k];
}
 
// Занести элемент по индексу
 
void put_el(double *matr, int cols, int i, int j, double v)
{
    int k=ind(i,j,cols);
    matr[k]=v;
}
 
// Умножение  матриц
 
double *mult_matr(double *a, int ca, int ra, double *b, int cb, int rb)
{
    if (ca != rb) return NULL;
    
    double *res=(double *)calloc(sizeof(double),ra*cb);
    
    int i,j,k,ik,kj;
    
    double s;
    
    for (i=0; i<ra; i++)
        for (j=0; j<cb; j++)
        {
            s=0;
            for (k=0; k<ca; k++)
                s+=get_el(a,ca,i,k)*get_el(b,cb,k,j);
            put_el(res,cb,i,j,s);
        }
    return res;
}
 
// Транспонирование матрицы
 
double *transp_matr(double *a, int cols, int rows)
{
    double *res=(double *)calloc(sizeof(double),cols*rows);
    int i,j,k;
    double v;
 
    for (i=0; i<rows; i++)
        for (j=0; j<cols; j++)
        {
            v=get_el(a,cols,i,j);
            put_el(res,rows,j,i,v);
        }
    return res;
}
 
// Печать матрицы
 
void print_matr(double *a, int cols, int rows)
{
    int i,j;
    for (i=0; i<rows; i++)
    {
        for (j=0; j<cols; j++) printf("%10.5lf ",get_el(a,cols,i,j));
        printf("\n");
    }   
}
 
// Ввод матрицы
 
void input_matr(char nam,double *a,int cols, int rows)
{
    int i,j;
    for (i=0; i<rows; i++)
        for (j=0; j<cols; j++)
        {
            printf("%c[%d,%d]=",nam,i,j);
            scanf("%lf",&a[ind(i,j,cols)]);
        }
}
 
int main()
{
    int i,n;
    
    printf("n=");
    scanf("%d",&n);
    
    // Ввод вектора-строки A 
    
    double *a=(double *)calloc(sizeof(double),n); 
    input_matr('A',a,n,1);
    
    // Ввод вектора-столбца B 
    
    double *b=(double *)calloc(sizeof(double),n); 
    input_matr('B',b,1,n);
 
    // Ввод вектора-строки V 
    
    double *v=(double *)calloc(sizeof(double),n); 
    input_matr('V',v,n,1);
 
    // Транспонируем A, получаем вектор-столбец
 
    double *at=transp_matr(a,n,1);
 
    // Вычисляем At-B
    
    for (i=0; i<n; i++) at[i]=at[i]-b[i];
    
    // Умножаем результат на вектор строку V
    
    double *r=mult_matr(at,1,n,v,n,1);
 
    print_matr(r,n,n);
 
    free(a);
    free(at);
    free(b);
    free(v);
    free(r);
 
    return 0;
}

Новые блоги и статьи Все статьи Все блоги /
GitHub Copilot в JetBrains: модель по умолчанию, Fix и контроль MCP and_y87 12.10.2026 10 октября 2026 года GitHub обновил Copilot для JetBrains: добавил настройки стартовой модели, исправления из диагностики и управление автозапуском MCP. Поддержка IDE 2025. 1 прекращена — требуется. . .	День 5 из 28: Codex подсказывает следующий запрос and_y87 12.10.2026 9 октября 2026 года OpenAI запустила бету Composer predictions — подсказок следующего сообщения в текущем диалоге Codex. Это помощь с формулировкой запроса, а не автоматическое продолжение задачи. . . .	Запустил конкурс "тем и промптов для текстовых квестов созданных почти чисто ИИ" Adler 06.10.2026 Всем привет! За последние три-четыре дня я создал более 16 текстовых квестовых игр используя преимущественно по одному запросу к ИИ на игру. Мне так понравилось смотреть все ветки/ сцены во всех. . .	ИИ не может найти нужный язык в списке Supersumestria 05.10.2026 Я ему даю вот такое изображение и прошу найти и подчеркнуть немецкий язык. Возвращает он вот это: https:/ / i. **********/ vqBWLe2. png Нужную строчку в 3й колонке просто выдумал. . Это. . .
Новая последняя моя музыка в SUNO zorxor 05.10.2026 Здравствуйте, дорогие мои друзья! С большой радостью я хотел бы представить вам свою новую последнею музыку, которую сгенерировала мне по моей просьбе нейросеть SUNO. С уважением, zorxor. Это. . .	Программный домашний кинотеатр russiannick 27.09.2026 Сподобился на программный домашний кинотеатр. В качестве ЯВУ по традиции выбрал js. В помощники взял Яндекс-Алису. Было создано три зала на разные интересы. исторические и ретро сериал Хичкок. . .	Беседа с ИИ о программистах, недопускающих к созданию и правке кода генеративные ИИ и причины этого zorxor 21.09.2026 Раньше я радовался или получал некоторые эмоции, пусть небольшие, но всё же, от самого процесса написания кода, рекомпиляции и запуска, видя постепенное развитие программы и прочее. А теперь лень. . .	Мобильное приложение ColorStep pavlinmavlin 17.09.2026 Реализовал приложение Красный, Зеленый, Синий в Unity3d + c#. Название изменил на ColorStep. Приложение прошло модерацию и теперь доступно для скачивания. Делал его сам, шаг за шагом — и вот,. . .

@Annemesski 2692 / 1349 / 485 Регистрация: 08.11.2016 Сообщений: 3,739
	14.07.2021, 12:48
	Catstail, зачем все эти рассуждения про вектора и матрицы, есть же простое правило о совместимости матриц: считая вектор-строку как матрицу 1xN, а вектор-столбец как матрицу Nx1 и проверяя матрицы на совместимость выдавать либо результат, либо сообщение о невозможности выполнения умножения. 1