Реализация метода Зейделя

@4unkur · Регистрация: 16.06.2014

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

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

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
??? Error using ==>
mtimes
Inner matrix dimensions
must agree.
 
Error in ==> gseid at 12
            X(j) =
            (B(j)-A(j,1:j-1)
            * X(1:j-1) -
            A(j,j+1:N) *
            P(j+1:N)) /
            A(j,j);

ВОТ сам код :

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
function X = gseid(A, B, P, eps, max1)
 
N = length(B);
 
for k=1:max1
    for j=1:N
        if j==1
            X(1) = (B(1) - A(1,2:N) * P(2:N)) / A(1,1);
        elseif j==N
            X(N) = (B(N) - A(N,1:N-1) * (X(1:N-1))') / A(N,N);
        else
            X(j) = (B(j)-A(j,1:j-1) * X(1:j-1) - A(j,j+1:N) * P(j+1:N)) / A(j,j);
        end
    end
    err = abs(norm(X'-P));
    relerr = err / (norm(X) + eps);
    P=X';
    if (err < eps) || (relerr < eps)
        break;
    end
end
X = X'

@Nick07 · 25.12.2014, 08:30

А это не то:
He пyтaй мaтpичнoe и пoэлeмeнтнoe yмнoжeниe/дeлeниe/cтeпeнь!
СЕКРЕТНО! Зловещая точка претыкания

@4unkur · 25.12.2014, 15:41 **[ТС]**

Nick07, Похоже не то, судя по тому, что все равно эта ошибка вылазит

@Nick07 · 25.12.2014, 16:51

Судя по сообщение об ошибки, или размерность матриц не согласована (см. умножение матриц) или в цикле идет обращение к не существующему элементе матрицы. Пустите в цикле вывод на монитор промежутков и сразу все станет ясно.

@4unkur · 25.12.2014, 20:03 **[ТС]**

проблема в том, что этот код написал сам учитель ))). Походу надо делать все сначала самому.

@VSI · 25.12.2014, 20:59

Сообщение от 4unkur

проблема в том, что этот код написал сам учитель...

А учитель что, не человек? И не может ошибаться? Зря Вы так думаете... Ошибаться могут ВСЕ

@4unkur · 26.12.2014, 15:30 **[ТС]**

VSI, Да вы правы, ему 70 лет. Достойный возраст. И да я не нашел ошибку, а немного модернизировал этот же код. Теперь все работает и я сдал задание. Всем спасибо

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	25.12.2014, 08:30
	А это не то: He пyтaй мaтpичнoe и пoэлeмeнтнoe yмнoжeниe/дeлeниe/cтeпeнь! СЕКРЕТНО! Зловещая точка претыкания 0

@4unkur 16 / 16 / 10 Регистрация: 16.06.2014 Сообщений: 180
	25.12.2014, 15:41 [ТС]
	Nick07, Похоже не то, судя по тому, что все равно эта ошибка вылазит 0

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	25.12.2014, 16:51
	Судя по сообщение об ошибки, или размерность матриц не согласована (см. умножение матриц) или в цикле идет обращение к не существующему элементе матрицы. Пустите в цикле вывод на монитор промежутков и сразу все станет ясно. 0

@4unkur 16 / 16 / 10 Регистрация: 16.06.2014 Сообщений: 180
	25.12.2014, 20:03 [ТС]
	проблема в том, что этот код написал сам учитель ))). Походу надо делать все сначала самому. 0

@4unkur 16 / 16 / 10 Регистрация: 16.06.2014 Сообщений: 180
	26.12.2014, 15:30 [ТС]
	VSI, Да вы правы, ему 70 лет. Достойный возраст. И да я не нашел ошибку, а немного модернизировал этот же код. Теперь все работает и я сдал задание. Всем спасибо 0