Степенной метод нахождения наибольшего по модулю собственного значения и соответствующего ему собственного вектора

@Alexandr1966 · Регистрация: 17.12.2012

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

Привет всем. Есть код программы, но почему то он не работает, подскажите в чем дело, может ввожу неправильно.

Тестовый пример:

Матрица А:
4 2 2
2 5 1
2 1 6

Матрица x1:
0,8077c1
0,7720c1
c1

Ответ:
L1=8,3874

Код программы:

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
// ВМ6.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include "math.h"
#include "conio.h"
 
void Input(int n,int A[][100])
{
int i,j;
for (i=0;i<n;i++)
for (j=0;j<n;j++)
scanf("%d",&A[i][j]);
}
 
void Matr(int A[][100],int n)
{
float w0[100],w[100],summ=0,w0norm[100],e,d,d0;
int i,j,k;
for (i=0;i<n;i++)
w0[i]=0;
w0[0]=1;
do
{
for (i=0;i<n;i++)
summ=summ+w0[i]*w0[i];
d0=sqrt(summ);
for (i=0;i<n;i++)
w0norm[i]=w0[i]/d0;
for (i=0;i<n;i++)
    {
    w[i]=0;
    for (j=0;j<n;j++)
    w[i]=w[i]+A[i][j]*w0norm[j];
    }
summ=0;
for (i=0;i<n;i++)
summ=summ+w[i]*w[i];
d=sqrt(summ);
e=fabs(d-d0);
for (i=0;i<n;i++)
w0[i]=w[i];
summ=0;
}
while(e>0.001);
printf("%f\n",d);
for (i=0;i<n;i++)
printf("%f\n",w0norm[i]);
}
 
void main()
{
int n,A[100][100];
scanf("%d",&n);
Input(n,A);
Matr(A,n);
}

Добавлено через 12 часов 47 минут
есть предложения?

Добавлено через 1 час 47 минут
кто-нибудь в курсе как исправить?

@Alexandr1966 · 23.03.2014, 09:28 **[ТС]**

ребята отзовитесь)

@recoder · 23.03.2014, 09:41

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
#include <math.h>
#include <conio.h>
#include <stdio.h>
 
void Input(int n,int A[][100])
{
    int i,j;
    for (i=0;i<n;i++)
        for (j=0;j<n;j++){
            printf("A[%d][%d]: ", i, j);
            scanf("%d",&A[i][j]);
        }
}
 
void Matr(int A[][100],int n)
{
    float w0[100],w[100],summ=0,w0norm[100],e,d,d0;
    int i,j,k;
    for (i=0;i<n;i++)
    w0[i]=0;    
    w0[0]=1;
    do
    {
        for (i=0;i<n;i++)
            summ=summ+w0[i]*w0[i];
        d0=sqrt(summ);
        for (i=0;i<n;i++)
            w0norm[i]=w0[i]/d0;
        for (i=0;i<n;i++)
        {
            w[i]=0;
            for (j=0;j<n;j++)
                w[i]=w[i]+A[i][j]*w0norm[j];
        }
        summ=0;
        for (i=0;i<n;i++)
            summ=summ+w[i]*w[i];
        d=sqrt(summ);
        e=fabs(d-d0);
        for (i=0;i<n;i++)
            w0[i]=w[i];
        summ=0;
    } while(e>0.001);
    for (i=0;i<n;i++){
        for (j=0;j<n;j++)
            printf("%d ", A[i][j]);
        printf("\n");
    }
    printf("L = %f\n",d);   
}
 
int main()
{
    int n,A[100][100];
    printf("N = ");
    scanf("%d",&n);
    Input(n,A);
    Matr(A,n);
    getchar();
    getchar();
    return 0;
}

@Alexandr1966 · 23.03.2014, 10:15 **[ТС]**

recoder, а можно еще в коде добавить вас попросить cout которые бы поясняли что где вводить например, cout << "Ввод матрицы А" и так далее

@recoder · 23.03.2014, 10:19

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
#include <cmath>
#include <iostream>
 
void Input(int n,int A[][100])
{
    int i,j;
    for (i=0;i<n;i++)
        for (j=0;j<n;j++){
            std::cout << "Enter A[" << i << "][" << j << "]: ";
            std::cin >> A[i][j];
        }
}
 
void Matr(int A[][100],int n)
{
    float w0[100],w[100],summ=0,w0norm[100],e,d,d0;
    int i,j,k;
    for (i=0;i<n;i++)
    w0[i]=0;    
    w0[0]=1;
    do
    {
        for (i=0;i<n;i++)
            summ=summ+w0[i]*w0[i];
        d0=sqrt(summ);
        for (i=0;i<n;i++)
            w0norm[i]=w0[i]/d0;
        for (i=0;i<n;i++)
        {
            w[i]=0;
            for (j=0;j<n;j++)
                w[i]=w[i]+A[i][j]*w0norm[j];
        }
        summ=0;
        for (i=0;i<n;i++)
            summ=summ+w[i]*w[i];
        d=sqrt(summ);
        e=fabs(d-d0);
        for (i=0;i<n;i++)
            w0[i]=w[i];
        summ=0;
    } while(e>0.001);
    std::cout << "\nMatrix:\n";
    for (i=0;i<n;i++){
        for (j=0;j<n;j++)
            std::cout << A[i][j] << " ";
        std::cout << std::endl;
    }
    std::cout << "\nL = " << d << std::endl;  
}
 
int main()
{
    int n,A[100][100];
    std::cout << "Enter N = ";
    std::cin >> n;
    Input(n,A);
    Matr(A,n);
    
    std::cin.get();
    std::cin.get();
    return 0;
}

@Alexandr1966 · 23.03.2014, 14:13 **[ТС]**

recoder, а почему я ввожу только матрицу A?
Где вот это реализовано:
Матрица x1:
0,8077c1
0,7720c1
c1

Добавлено через 3 минуты
И как здесь вообще работает алгоритм подсчета? можете прокомментировать?

Добавлено через 3 часа 43 минуты
recoder, подскажи а

@recoder · 23.03.2014, 14:22

У вас есть задание? Не могли бы вы его полностью написать сюда.

@Alexandr1966 · 23.03.2014, 14:46 **[ТС]**

recoder, Степенной метод нахождения наибольшего по модулю собственного значения и соответствующего ему собственного вектора, написать программу на с++.

Тестовый пример:

Матрица А:
4 2 2
2 5 1
2 1 6

Матрица x1:
0,8077c1
0,7720c1
c1

Ответ:
L1=8,3874

@recoder · 23.03.2014, 15:10

Как-то так!

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
#include <cmath>
#include <iostream>
 
void Input(int n,int A[][100])
{
    int i,j;
    for (i=0;i<n;i++)
        for (j=0;j<n;j++){
            std::cout << "Enter A[" << i << "][" << j << "]: ";
            std::cin >> A[i][j];
        }
}
 
void Matr(int A[][100],int n)
{
    float w0[100],w[100],summ=0,w0norm[100],e,d,d0;
    int i,j,k;
    w0[0]=1;
    for (i=1;i<n;i++)
        w0[i]=0;    
    do
    {
        for (i=0;i<n;i++)
            summ=summ+w0[i]*w0[i];
        d0=sqrt(summ);
        for (i=0;i<n;i++)
            w0norm[i]=w0[i]/d0;
        for (i=0;i<n;i++)
        {
            w[i]=0;
            for (j=0;j<n;j++)
                w[i]=w[i]+A[i][j]*w0norm[j];
        }
        summ=0;
        for (i=0;i<n;i++)
            summ=summ+w[i]*w[i];
        d=sqrt(summ);
        e=fabs(d-d0);
        for (i=0;i<n;i++)
            w0[i]=w[i];
        summ=0;
    } while(e>0.001);
    std::cout << "\nMatrix:\n";
    for (i=0;i<n;i++){
        for (j=0;j<n;j++)
            std::cout << A[i][j] << " ";
        std::cout  << " = " << w0norm[i] << std::endl;
    }
    std::cout << "\nL = " << d << std::endl;  
}
 
int main()
{
    int n,A[100][100];
    std::cout << "Enter N = ";
    std::cin >> n;
    Input(n,A);
    Matr(A,n);
    
    std::cin.get();
    std::cin.get();
    return 0;
}

@Alexandr1966 · 23.03.2014, 16:00 **[ТС]**

recoder,
Матрица x1:
0,8077c1
0,7720c1
c1
вот это где в кодеее???!?!

@recoder · 23.03.2014, 16:19

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

w0norm[100]

Судя по всему это вот это.

@Alexandr1966 · 24.03.2014, 16:32 **[ТС]**

recoder, то есть эта матрица X сама по себе считается?

Добавлено через 22 часа 21 минуту
есть у кого нибудь соображения на этот счет?

Новые блоги и статьи Все статьи Все блоги /
Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек SDL3 и Box2D из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия SDL 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual. . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»

@Alexandr1966 29 / 29 / 2 Регистрация: 17.12.2012 Сообщений: 493
	23.03.2014, 09:28 [ТС]
	ребята отзовитесь) 0

@Alexandr1966 29 / 29 / 2 Регистрация: 17.12.2012 Сообщений: 493
	23.03.2014, 10:15 [ТС]
	recoder, а можно еще в коде добавить вас попросить cout которые бы поясняли что где вводить например, cout << "Ввод матрицы А" и так далее 0

@Alexandr1966 29 / 29 / 2 Регистрация: 17.12.2012 Сообщений: 493
	23.03.2014, 14:13 [ТС]
	recoder, а почему я ввожу только матрицу A? Где вот это реализовано: Матрица x1: 0,8077c1 0,7720c1 c1 Добавлено через 3 минуты И как здесь вообще работает алгоритм подсчета? можете прокомментировать? Добавлено через 3 часа 43 минуты recoder, подскажи а 0

@recoder 135 / 134 / 53 Регистрация: 13.09.2013 Сообщений: 260 Записей в блоге: 2
	23.03.2014, 14:22
	У вас есть задание? Не могли бы вы его полностью написать сюда. 0

@Alexandr1966 29 / 29 / 2 Регистрация: 17.12.2012 Сообщений: 493
	23.03.2014, 14:46 [ТС]
	recoder, Степенной метод нахождения наибольшего по модулю собственного значения и соответствующего ему собственного вектора, написать программу на с++. Тестовый пример: Матрица А: 4 2 2 2 5 1 2 1 6 Матрица x1: 0,8077c1 0,7720c1 c1 Ответ: L1=8,3874 0

@Alexandr1966 29 / 29 / 2 Регистрация: 17.12.2012 Сообщений: 493
	23.03.2014, 16:00 [ТС]
	recoder, Матрица x1: 0,8077c1 0,7720c1 c1 вот это где в кодеее???!?! 0

@Alexandr1966 29 / 29 / 2 Регистрация: 17.12.2012 Сообщений: 493
	24.03.2014, 16:32 [ТС]
	recoder, то есть эта матрица X сама по себе считается? Добавлено через 22 часа 21 минуту есть у кого нибудь соображения на этот счет? 0