Язык С++ как инструмент моделирования на основе решения дифференциальных уравнений в частных производных - page 38

36
в) для трехточечной аппроксимации со вторым порядком
разностная
схема крест примет вид:
 
.
2
5.2
1
,
2
5.2
1
,1 ,1 ,
2
5.2
1
2
, ,0 , 2 exp
,
, ,1 ,
5
4 2
2
2
,0
, ,1 ,
2
5
2 4
2
,1 ,1 ,
5
2
2
0
2
2
2
2
0
1
2
2
1
0
1 2
2
0
0
2
2
2
2
1
0
0
1 2
2
0 2
2
2
0
1 2
2
1
0
2
2
2
2
1
2
2
1
1
0
1 2
2
0
2
2
2
2
1
0
1
0
1 2
2
2
2
2
1
1 2
2
1
N
N
N
j
j
j
j
j
N
k
N
k
N
k
N
k
N
k
k
k
k
k
j
k
j
k
j
k
j
k
j
u
h
h
u
h
u
u
h
u
h
h
u
N j
u
h
u
h
u
h
u
N j hj
u
x при K k u
h h
u
h
u
u
x при K ku
h
u
h h
u
u
N j
u
h
u
h
u u
h
u


   
 

   
 
 

  
 

   
 
 

    
 
 

   
(2.11)
2.1.1. Численное решение задачи гиперболического типа с заданными
начальными и граничными условиями
Программная реализация на С++.
Программная реализация алгоритма решения рассматриваемой
задачи (2.1)-(2.5) на основе ее конечно-разностного представления
начального условия, явной схемы крест для дифференциального
уравнения гиперболического типа с первым начальным условием (2.4), с
аппроксимацией второго начального условия первым порядком точности
по времени и двухточечной аппроксимации граничных условий с первым
порядком точности по координате (2.6).
#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <fstream>
#include <stdio.h>
#include <string>
#define PI 3.1415926535897932384626433832795
using namespace std;
void main()
1...,28,29,30,31,32,33,34,35,36,37 39,40,41,42,43,44,45,46,47,48,...52
Powered by FlippingBook