Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.93/40: Рейтинг темы: голосов - 40, средняя оценка - 4.93
0 / 0 / 0
Регистрация: 29.04.2013
Сообщений: 7
1

Решение уравнения теплопроводности с помощью неявной разностной схемы

23.05.2015, 16:33. Показов 7893. Ответов 5
Метки нет (Все метки)

День добрый, форумчане. Нужен ваш совет.
Имеется у меня код, это решение уравнения теплопроводности с помощью неявной разностной схемы. и всё бы ничего, и вроде всё работает, за исключением одного момента: пока массив В объявлен размерностью конкретного числа double B[11]; всё работает, как только объявляешь double B[n1]; что мне и требуется, работать всё перестает. Я далеко не мастер в С++, но может быть вы меня наведете на мысль как исправить косяк в задаче

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
#include "math.h"
#include "iostream"
#include"fstream"
using namespace std;
 
 
 
double func( double x , double t )
{ return exp ( - 3.14 * 3.14 * t ) * sin ( 3.14 * x ) ; }
 
 
int main()
{
int n1 = 11 ;
int n,k;
double z;
double A[n1][n1];
double B[11];
double eps[n1];
double X[n1];
double et[n1];
double cou, tau=0.005,h=1.0/(n1-1);
 
cou = 2.0 * tau / ( h * h );
 
for ( int l = 0 ; l < n1 ; l++ )
for( int m = 0 ; m < n1 ; m++ )
A[m][l]=0;
 
for ( int j = 0 ; j < n1 + 1 ; j++ )
{
A[j][j-1] = cou / 2 ;
A[j][j] = - 1 - cou ;
A[j][j+1] = cou / 2 ;
}
 
for ( int i = 0; i < n1; i++ )
{
B[i] = - sin ( i * h * 3.14 );
}
 
n = n1 - 1 ;
eps[0] = - A[0][1] / A[0][0] ;
et[0] = B[0] / A[0][0] ;
 
for(int i=1;i<n;i++)
{
z=A[i][i]+A[i][i-1]*eps[i-1];
eps[i]=-A[i][i+1]/z;
et[i]=(B[i]-A[i][i-1]*et[i-1])/z;
}
 
X[n]=(B[n]-A[n][n-1]*et[n-1])/(A[n][n]+A[n][n-1]*eps[n-1]);
 
for(int i=n-1;i>=0;i--)
X[i]=eps[i]*X[i+1]+et[i];
 
 double max = 0;
for(int i=0; i<n1-1;i++)
if(fabs(func(i*h, tau)-X[i])>max)
max=fabs(func(i*h, tau)-X[i]);
printf("%f", max); 
 
return 0;
}
__________________
Помощь в написании контрольных, курсовых и дипломных работ здесь
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
23.05.2015, 16:33
Ответы с готовыми решениями:

Решение одномерного уравнения теплопроводности (методом сеток с помощью явной схемы)
по численным методам задали задачу, методом сеток с помощью явной схемы решить уравнение...

Задача теплопроводности по неявной разностной схеме
Друзья! Построил программу для решения краевой задачи теплопроводности по неявной разностной...

Применение конечно-разностной схемы для решения дифференциального уравнения турбулентной диффузии
Добрый день, форумчане! Помогите пожалуйста разобраться. Я программист. Сейчас занимаюсь...

Решение уравнения теплопроводности
задача:

5
2423 / 1825 / 402
Регистрация: 15.12.2013
Сообщений: 8,002
23.05.2015, 17:01 2
ОксанаКурбатова, чтобы заработало n1 должно быть константой.Но вы также можете использовать контейнер std::vector<> вместо массивов.
0
0 / 0 / 0
Регистрация: 29.04.2013
Сообщений: 7
23.05.2015, 17:41  [ТС] 3
ну так n1 константа, и остальные массивы прекрасно это принимают
0
2423 / 1825 / 402
Регистрация: 15.12.2013
Сообщений: 8,002
23.05.2015, 17:54 4
Цитата Сообщение от ОксанаКурбатова Посмотреть сообщение
ну так n1 константа, и остальные массивы прекрасно это принимают
Под константой я вот что подразумевал:
C++
1
const int n1 = 11 ;
1
0 / 0 / 0
Регистрация: 29.04.2013
Сообщений: 7
23.05.2015, 18:10  [ТС] 5
Работает) спасибо))
0
0 / 0 / 0
Регистрация: 25.02.2013
Сообщений: 46
27.10.2016, 21:30 6
ОксанаКурбатова, если не сложно, можно узнать условие задачи?
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
27.10.2016, 21:30

Решение уравнения теплопроводности
Всем доброго времени суток! У меня такая задача - решить численным методом уравнение...

Решение уравнения теплопроводности
Люди добрые , помогите пожалуйста вот с таким вопросом. Решаю уравнение теплопроводности с...

Решение уравнения теплопроводности неявным методом
Пытаюсь решить уравнение теплопроводности неявным методом, но ничего не выходит и график какой-то...

Нужно численное решение уравнения теплопроводности
Так сложилось, что с математикой и Матлабом я не в ладах, а налаживать с ними отношения - нет...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
6
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2022, CyberForum.ru