Форум программистов, компьютерный форум CyberForum.ru

Решение системы уравнений методом Гаусса. Исправить код - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Найти максимальное значение радиуса шара r и соответствующие значения радиуса основания r3 и высоты h шарового сектора http://www.cyberforum.ru/cpp-beginners/thread1124481.html
Здравствуйте))))) помогите пожалуйста!!!! надежда только на вас! не могу сделать математическую часть курсовой((( Найти максимальное значение радиуса шара r (r изменяется от r1 до r2 c шагом k) и соответствующие значения радиуса основания r3 и высоты h шарового сектора, для которых площадь поверхности шарового сектора не превышает площади поверхности прямоугольного параллелепипеда со сторонами...
C++ Stack around the variable 'stat_c1' was corrupted Добрый вечер форумчане. Написал програмку, все отлично, все работает. Но при завершении программы, т.е. когда выполняеться строчка return 0; Выскакивает ошибка Stack around the variable 'stat_c1' was corrupted и Stack around the variable 'stat_c2' was corrupted Я знаю суть этих ошибок, это обращение к неизвестной итерации масива, но тут то программа закончила уже работу, и только тогда... http://www.cyberforum.ru/cpp-beginners/thread1124479.html
C++ Класс памяти по умолчанию для глобальной переменной
Здравствуйте! Прочитала такую вещь тут Объявление переменных на глобальном уровне - это или определение переменных, или ссылки на определения, сделанные в другом месте программы. Объявление глобальной переменной, которое инициализирует эту переменную (явно или неявно), является определением переменной. Определение на глобальном уровне может задаваться в следующих формах: 1. Переменная...
C++ Заполнить рандомно массив и все отрицательные элементы отправить в конец
В массиве A=(a1, а2, ..., аn) все отрицательные элементы отправить в «хвост» массива. ( Это задание) Вот сам код. // ConsoleApplication4234452.cpp: определяет точку входа для консольного приложения. // #include "stdafx.h" #include <iostream> #include "conio.h" #include <time.h>
C++ Реализовать метод приведения времени в минуты http://www.cyberforum.ru/cpp-beginners/thread1124444.html
Реализовать внешнюю функцию с именем make_mun(), где mun – тип реализуемой структуры. Функция должна получать в качестве аргументов значения для полей структуры и возвращать структуру требуемого типа. При передаче ошибочных параметров следует выводить сообщение и заканчивать работу. Поле first – целое положительное число, часы; поле second – целое положительное число, минуты. Реализовать метод...
C++ Структура-пара. Реализовать метод rangecheck() – проверку заданного числа Реализовать внешнюю функцию с именем make_mun(), где mun – тип реализуемой структуры. Функция должна получать в качестве аргументов значения для полей структуры и возвращать структуру требуемого типа. При передаче ошибочных параметров следует выводить сообщение и заканчивать работу. Поле first – дробное число, левая граница диапазона; поле second – дробное число, правая граница диапазона.... подробнее

Показать сообщение отдельно
Katan0073
0 / 0 / 0
Регистрация: 26.02.2014
Сообщений: 9
19.03.2014, 23:43     Решение системы уравнений методом Гаусса. Исправить код
Добрый день!
Столкнулся с небольшой проблемой, но не удается от неё избавиться. Программа решает систему уравнений методом гаусса.
Суть проблемы: при вводе конкретных коэффициентов уравнений:
C++
1
2
3
0,3  0,1  2 | 2
1   3   0 |1
0.9  0.3  6 |0
при приведении матрицы к стандартному виду последний коэффициент (a33) не обращается в ноль, хотя должен. Причину нашел, но я не знаю, как с ней справиться. Попытался показать проблему на скрине. Как я понял, в этом последнем случае коэффициент k не округляется до 3, а остается 2,99999976. Почему?
Проблемный участок помечен: "приводим к стандартному виду"
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
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#pragma warning(disable: 4996)
using namespace std;
void vvodsohr(); void chtreshsohr(); void nevyazka();
int n, i, j;
 
FILE *ish, *rez;
 
int main()
{
    while (true){
        cout << "MENU:" << endl << "1. Zapis' matricy v file." << endl << "2. Reshenie matricy." << endl << "3. Nevyazka." << endl << "0. Vyhod." << endl;
        cin >> i;
        switch (i){
        case 1: vvodsohr(); break;
        case 2: chtreshsohr(); break;
        case 3: nevyazka(); break;
        case 0: return 0;
        }
    }
}
 
void vvodsohr(){
 
    ish = fopen("ishodnye_dannye.txt", "wb");
    cout << "Vvedite razmer matricy (nxn): ";
    cin >> n;
    fwrite(&n, sizeof(int), 1, ish);
    float **a = new float*[n], *b = new float[n]; //инициализация динамических массивов
    for (i = 0; i < n; i++) a[i] = new float[n];
    for (i = 0; i<n; i++){//заполнение массива и запись в файл
        cout << "Vvodite elementy i svobodny chlen stroki " << i + 1 << ":";
        for (j = 0; j < n; j++){
            cin >> a[i][j];
            fwrite(&a[i][j], sizeof(float), 1, ish);
        }
        cin >> b[i];
        fwrite(&b[i], sizeof(float), 1, ish);
    }
 
    cout << endl << "File sozdan, vvodimye dannye sohraneny." << endl << endl;
    fclose(ish);
    for (i = 0; i < n; i++) delete[]a[i];
    delete[]a;
    delete[]b;
}
 
void chtreshsohr(){
    ish = fopen("ishodnye_dannye.txt", "rb");
    fread(&n, sizeof(float), 1, ish);
    float **a = new float*[n], *b = new float[n];
    for (i = 0; i < n; i++) a[i] = new float[n];
 
    for (i = 0; i < n; i++){//считывание из файла в массив.
        for (j = 0; j < n; j++)
            fread(&a[i][j], sizeof(float), 1, ish);
        fread(&b[i], sizeof(float), 1, ish);
    }
 
 
 
 
    int f, jmax; float k, max, *p = NULL, l; //приводим к стандартному виду
    for (i = 0; i < n; i++){
        max = a[i][i];
        jmax = i;
        for (j = i; j < n; j++)
        if ((a[j][i])>max) { max = a[j][i]; jmax = j; }
 
        p = a[i];
        a[i] = a[jmax];
        a[jmax] = p;
        l = b[i];
        b[i] = b[jmax];
        b[jmax] = l;
 
        for (f = i + 1; f < n; f++){
            k = a[f][i] / a[i][i];
            b[f] = b[f] - k*b[i];
            for (j = i; j < n; j++)
                a[f][j] = a[f][j] - a[i][j] * k;
        }
    }
 
    cout << "Standartnaya:" << endl;
    for (i = 0; i < n; i++){//вывод приведенной матрицы
        for (j = 0; j < n; j++)cout << setw(9) << a[i][j] << " ";
        cout << "|" << b[i] << endl;
    }
 
    for (i = 0; i < n; i++)
    if (a[i][i] == 0){
        cout << "Sistema ne imeet resheniy" << endl;
        cout << "Press any key to exit...";
        getch();
        exit(1);
    }
 
    float s = 0, *x = new float[n];//находим x
    for (i = 0; i<n; i++) x[i] = 0;
    for (i = n - 1; i >= 0; i--){
        s = b[i];
        for (j = n - 1; j>i; j--) s -= x[j] * a[i][j];
        x[i] = s / a[i][i];
    }
 
    for (i = 0; i < n; i++) cout << "x" << i + 1 << "=" << x[i] << " ";
    cout << endl << endl;
    fclose(ish);
    rez = fopen("rezult.txt", "wb");//запись иксов в файл
    for (i = 0; i < n; i++)
        fwrite(&x[i], sizeof(float), 1, rez);
    fclose(rez);
    for (i = 0; i < n; i++) delete[]a[i];
    delete[]a;
    delete[]b;
    delete[]x;
}
Миниатюры
Решение системы уравнений методом Гаусса. Исправить код  
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 20:01. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru