Презентация Метод сеток для решения ДУ в частных производных онлайн
На нашем сайте вы можете скачать и просмотреть онлайн доклад-презентацию на тему Метод сеток для решения ДУ в частных производных абсолютно бесплатно. Урок-презентация на эту тему содержит всего 23 слайда. Все материалы созданы в программе PowerPoint и имеют формат ppt или же pptx. Материалы и темы для презентаций взяты из открытых источников и загружены их авторами, за качество и достоверность информации в них администрация сайта не отвечает, все права принадлежат их создателям. Если вы нашли то, что искали, отблагодарите авторов - поделитесь ссылкой в социальных сетях, а наш сайт добавьте в закладки.
Оцените презентацию от 1 до 5 баллов!
- Тип файла:ppt / pptx (powerpoint)
- Всего слайдов:23 слайда
- Для класса:1,2,3,4,5,6,7,8,9,10,11
- Размер файла:460.50 kB
- Просмотров:68
- Скачиваний:1
- Автор:неизвестен
Слайды и текст к этой презентации:
№1 слайд
Содержание слайда: Тема 6. Метод сеток для решения
ДУ в частных производных
Одномерное нестационарное уравнение теплопроводности:
Явная и неявная схемы
Задача Дирихле для двумерного уравнения Пуассона
Метод простой итерации с релаксацией
Метод Зейделя
Метод продольно-поперечной прогонки
№2 слайд
Содержание слайда: Задача о нагреве стержня, по которому пропускается ток
№3 слайд
Содержание слайда: Одномерное нестационарное уравнение теплопроводности
Область интегрирования
Сетка
Таблица искомого решения
№4 слайд
Содержание слайда: Получение конечноразностной схемы
№5 слайд
Содержание слайда: Явная схема
№6 слайд
Содержание слайда: Реализация явной схемы
u(i)=u0; x(i)=(i-1)*h; gi(i)=G(x(i)+h/2); 1≤i≤N+1
th2= tau/h^2;
t=0;
Plot(x,u);
for k=1:K
for i=2:N
u1(i)=u(i)+th2*(gi(i-1)*u(i-1)-(gi(i-1)+gi(i))*u(i)+
gi(i)*u(i))+tau*f(i);
end
u1(1)=be0;
u1(N1)=be1;
Plot(x,u1);
u=u1;
t=t+tau;
end;
…
№7 слайд
Содержание слайда:
№8 слайд
Содержание слайда:
№9 слайд
Содержание слайда: Реализация метода прогонки
u(i)=u0; x(i)=(i-1)*h; gi(i)=G(x(i)+h/2); 1≤i≤N+1
for k=1:K
c(1)=…; b(1)=…; d(1)-…;
for i=2:N
a(i)=
b(i)=
c(i)=
d(i)= …
end
ks(1)=-c(1)/b(1); et(1)=d(1)/b(1);
for i=2:N1
z=b(i)+a(i)*ks(i-1);
ks(i)=-c(i)/z; et(i)=(d(i)-a(i)*et(i-1))/z;
end;
u1(N1)=be1;
For i=N:-1:1
u1(i)=ks(i)*u1(i+1)+et(i);
End;
Plot(x,u1);
u=u1;
t=t+tau;
end;
№10 слайд
Содержание слайда: Неявная схема второго порядка
Кранка-Николсона
№11 слайд
Содержание слайда:
№12 слайд
Содержание слайда:
№13 слайд
Содержание слайда: Задача Дирихле для двумерного уравнения Пуассона
№14 слайд
Содержание слайда:
№15 слайд
Содержание слайда: Конечно-разностная схема
№16 слайд
Содержание слайда: Метод простой итерации с релаксацией
№17 слайд
Содержание слайда: Реализация метода простой итерации с релаксацией
For k=1:Kit
for i=2:N
for j=2:M
a=… b=… c=… d=… e=…
up= a*u(i-1,j)+b*u(i+1,j)+c*u(i,j-1)+d*u(i,j+1)+e;
u1(i,j)=wr*up+(1-wr)*u(i.j);
end;end; //ij
for i=2:N
u1(i,1)=be0; u1(i,M1)=be1;
for j=2:M
u1(1,j)=al0; u1(N1,j)=al1;
u=u1;
End; //k
surf (x,y,u1);
№18 слайд
Содержание слайда: Метод Зейделя
For k=1:Kit
d=0;
for i=2:N
for j=2:M
a=… b=… c=… d=… e=…
up= a*u(i-1,j)+b*u(i+1,j)+c*u(i,j-1)+d*u(i,j+1)+e;
If abs(u(i,j)-up)>d then d= abs(u(i,j)-up);
u(i,j)=wr*up+(1-wr)*u(i.j);
end;end; //ij
for i=2:N
u(i,1)=be0; u(i,M1)=be1;
for j=2:M
u(1,j)=al0; u(N1,j)=al1;
If d<eps then continue;
End; //k
surf (x,y,u);
№19 слайд
Содержание слайда: Метод продольно-поперечной прогонки
№20 слайд
Содержание слайда: Метод продольно-поперечной прогонки (продолжение)
№21 слайд
Содержание слайда: Метод продольно-поперечной прогонки (продолжение)
№22 слайд
Содержание слайда: Программная реализация
for k=1:Kit
ks(1)=…; et(1)=…;
for j=2:M1 for i=2:N1
a=…b=…c=…d=…;
z=b+a*ks(i-1);
ks(i)=-c/z; et(i)=(d-a*et(i-1))/z;
end;
u(N1,j)=et(N1);
For i=N:-1:1
u(i,j)=ks(i)*u1(i+1,j)+et(i);
End; end; end; (ij)
ks(1)=…; et(1)=…;
for i=2:N1 for j=2:M1
a=…b=…c=…d=…;
z=b+a*ks(j-1);
ks(j)=-c/z; et(j)=(d-a*et(j-1))/z;
end;
u(I,M1)=et(M1);
For j=M:-1:1
u(i,j)=ks(j)*u1(i,j+1)+et(j+1);
End; end; end; (ij)
Plot(x,u1);
u=u1;
t=t+tau;
end;
№23 слайд
Содержание слайда: Конец темы 6
Ваши вопросы
Скачать все slide презентации Метод сеток для решения ДУ в частных производных одним архивом:
