Решить систему алгебраических линейных неоднородных уравнени

@Stas0n · Регистрация: 13.07.2011

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

У меня есть система линейных уравнений. В ней 4000 уравнений.
Киньте плиз код для её решения. Желательно, чтобы он был максимально быстрым.

@Stas0n · 08.08.2011, 18:18 **[ТС]**

Кстати, ты время не засекал? просто даже интересно сколько будет работать....

@grizlik78 · 08.08.2011, 18:20

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

Кстати, ты время не засекал? просто даже интересно сколько будет работать....

Это не так-то просто, весь скрипт отрабатывает за доли секунды.

@Stas0n · 08.08.2011, 18:21 **[ТС]**

to Юра: вот на всякий случай решение системы... я тут проверил - оно неплохо подходит.

@Stas0n · 08.08.2011, 18:24 **[ТС]**

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

Это не так-то просто, весь скрипт отрабатывает за доли секунды.

Разработчики позаботились об этом)
Всего то надо сделать так:
если ты хочешь посмотреть, сколько времени работает конкретной кусок кода в scilsb, то нужно просто написать так:
timer();
// кусок кода, время работы которого ты хочешь узнать
timer()

@grizlik78 · 08.08.2011, 18:29

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

Разработчики позаботились об этом)

Ну хорошо, весь скрипт, с загрузкой и выгрузкой 0,2 с.
Само решение без ввода-вывода и проверки 0,03-0,07 с.
scilab 5.1.1, линукс, x86_64, C2D E6600

@Stas0n · 08.08.2011, 18:31 **[ТС]**

ОК. Завтра протестю на винде. У меня scilsb 5.3.3

@-=ЮрА=- · 09.08.2011, 09:48

Я проверил ваши решения, по вот такому алгоритму

C++
1
2
3
4
5
6
7
8
9
printf("\r\n\tCHECK RESULTS A*X - B BEGIN\r\n");
        for(i = 0; i < m; i++)
        {
            E[i] = 0;
            for(r = 0; r < m; r++)
                E[i] += A[i][r]*X[r];
            E[i] -= B[i];
        }
        printf("CHECK RESULTS DONE\r\n");

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

to Юра: вот на всякий случай решение системы... я тут проверил - оно неплохо подходит.

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

Хорошо, сейчас в соответствующем разделе покажу. А пока вот вам правильное решение, полученное в scilab. Когда программа будет давать что-то похожее, можно будет на неё и посмотреть

- привожу для анализа векторы погрешностей scilab.
Как вы видите все алгоритмы дают погрешность, причём число нулей в векторе погрешностей при решении методом Гаусса, не меньше чем при решении scilab. Также выкладываю исходники алгоритма проверки, не забывайте что для скилаба юзаете полную матрицу А в то время как я удаляю строки и столбцы с нулевыми элементы на главной диагонали, поэтому при запросе на сжатие матрицы вы должны нажать N и У для моего вектора Х.
PS:В заключение скажу что не стоит полностью уповать на решения CAD-систем, ещё по универу помню косяки маткада при решении диффуров...

@Stas0n · 09.08.2011, 10:31 **[ТС]**

Сейчас посмотрим что и как...

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

Добавлено через 5 минут
Опа! Юр, нашел..ну не то чтобы косяк... просто ты считываешь вектор В без какой-либо предварительной подготовки. так?
Прости, просто когда кидал тебе, там короче запятые есть.. мой косяк.
Сейчас посмотимр какие результаты выдаст на "нормальном" столбце.

Добавлено через 18 секунд
Ну чтож.. результаты то и не сильно отличаются

@-=ЮрА=- · 09.08.2011, 10:40

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

Что же касается точности матпакета, тот тут полностью с тобой согласен. Сам как-то неплохо так влетел из-за него.

-Полностью согласен, всё же я за метод Гаусса, т.к. этот код уже проверен в течении 7 лет, как его написал, на СЛАУ небольшого порядка ни разу косо не считал

@Stas0n · 09.08.2011, 10:48 **[ТС]**

Сейчас понять где прога косячит и, тогда ей цены не будет.

@grizlik78 · 09.08.2011, 12:06

Сообщение от -=ЮрА=-

PS:В заключение скажу что не стоит полностью уповать на решения CAD-систем, ещё по универу помню косяки маткада при решении диффуров...

Ошибки бывают везде. Естественно, каким бы не был путь получения решения, надо уметь проверить результат и обязательно делать это, даже если программа/метод раньше работали безотказно.

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

@-=ЮрА=- · 09.08.2011, 13:21

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

Сейчас понять где прога косячит и, тогда ей цены не будет.

Работаю над этим, подумываю выгрузить в файл А и В после сжатия и сравнить с исходными, метод Гаусса там реализован верно, по сему заключаю что всё же ошибка возникает до расчёта, на этапе формирования А или её сжатия, матрица нууу очень большая поэтому баги сложно выявить...

@Stas0n · 09.08.2011, 16:53 **[ТС]**

Вот и мне кажется, что где-то в начале косяк происходит... пока не вижу(

Добавлено через 3 часа 25 минут
Юр, слушай, а как задать точность для вывода?
Вот у меня число 14.56, прога выводит 14.6(

Новые блоги и статьи Все статьи Все блоги /
SUNO Ai - Река Без Дна zorxor 31.07.2026 Автор стихотворения - астрофизик Марина Катыс Ссылка на сгенерированную музыкальную композицию: https:/ / suno. com/ song/ 6f6e5464-b290-4650-be6c-44c85f8d8013 Я говорю, что Время- как вода течет. . .	Из невошедшего на форум (диалог с ИИ-гугла) zorxor 29.07.2026 А вот, что интересно, сказал мне ИИ-гугла: Этот текст — эмоциональный пост пользователя под ником zorxor на интернет-форуме (вероятно, посвященном мистике, непознанному или альтернативной науке). . . .	Был праздник вчера, а я и не знал. kumehtar 28.07.2026 27. 07. 2026г. Intel Core 2 Duo исполнилось 20 лет Новости компьютерного мира и их обсуждение (4) Салют, шампанское, овации! :drink:	Нейтральные знания, чистый код - бла-бла-бла-бла, на самом деле кликбейт и самореклама, плагиат, и вот почему Hrethgir 27.07.2026 То-есть отклонение такой публикации говорит само за себя, и пусть только возьмут на вооружение после отклонения публикации - это будет чистейшим актом плагиата. Отклонял Хабр. Дословно, отклонённая. . .
тв 16 бой ии anaschu 27.07.2026 Великий Перелом ИИ: Как уравнения ОДУ Radau дожали цензурные фильтры Алисы Фиксируем в мемофонде Теории Всего беспрецедентный факт в истории ИИ-зондирования. В затяжном многораундовом. . .	мв 15. непроверенное, возможно, глюк anaschu 27.07.2026 НАУЧНО-АНАЛИТИЧЕСКИЙ ОТЧЕТ. РАЗДЕЛ 1. 1: «НАУКА» (РАСШИРЕННАЯ СТЕХИОМЕТРИЧЕСКАЯ И ГЕНЕТИЧЕСКАЯ ВЕРСИЯ)Тема: Теоретическое обоснование инвариантности 19-мерного тензорного ядра непрерывных ОДУ и. . .	Очистка реквизитов и табличных частей документа при копировании (вариант 2) Maks 26.07.2026 Алгоритм из решения ниже разработан на примере нетипового документа "ЗаявкаНаРаботу", разработанного в КА2. Задача: Заменить алгоритм запрета копирования документов для сотрудников с ролью "Стажер",. . .	Доктрина интенционального знания - Доктрина для портала "Срез". Hrethgir 25.07.2026 Может найдётся кто захочет оценить доктрину. . . Написания правил участия для меня роскошь, требующая лимита времени, поэтому все сообщения не прошедшие модерацию будут видны только участникам портала,. . .

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313

	Решить систему алгебраических линейных неоднородных уравнени 22.07.2011, 15:58. Показов 16393. Ответов 152 Метки нет (Все метки) У меня есть система линейных уравнений. В ней 4000 уравнений. Киньте плиз код для её решения. Желательно, чтобы он был максимально быстрым. 0

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313
	08.08.2011, 18:18 [ТС]
	Кстати, ты время не засекал? просто даже интересно сколько будет работать.... 0

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313
	08.08.2011, 18:31 [ТС]
	ОК. Завтра протестю на винде. У меня scilsb 5.3.3 0

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313
	09.08.2011, 10:31 [ТС]
	Сейчас посмотрим что и как... Добавлено через 6 минут Что же касается точности матпакета, тот тут полностью с тобой согласен. Сам как-то неплохо так влетел из-за него. В этом то и преимущество программы, написанной тобой - можно "пощупать" любой кусок кода и настроить его под себя. Добавлено через 5 минут Опа! Юр, нашел..ну не то чтобы косяк... просто ты считываешь вектор В без какой-либо предварительной подготовки. так? Прости, просто когда кидал тебе, там короче запятые есть.. мой косяк. Сейчас посмотимр какие результаты выдаст на "нормальном" столбце. Добавлено через 18 секунд Ну чтож.. результаты то и не сильно отличаются 0

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313
	09.08.2011, 10:48 [ТС]
	Сейчас понять где прога косячит и, тогда ей цены не будет. 0

@Stas0n 3 / 4 / 1 Регистрация: 13.07.2011 Сообщений: 313
	09.08.2011, 16:53 [ТС]
	Вот и мне кажется, что где-то в начале косяк происходит... пока не вижу( Добавлено через 3 часа 25 минут Юр, слушай, а как задать точность для вывода? Вот у меня число 14.56, прога выводит 14.6( 0