MATLAB的答案

解决线性对流方程(波动方程)松懈Wendroff方案

130(30天)
萨阿德·萨利姆
萨阿德·萨利姆 2018年11月9日
编辑: Torsten 2018年11月9日
这就跟你问声好！
我试图解决这个问题在文本。
我复制我的MATLAB代码解决松懈Wendroff方案。我努力把周期性边界条件。也许有人可以指导吗?
谢谢!
%美5093作业# 5 - 11月13日,2018年
% %问题1
% % Saad Saleem编写的代码
% %这MATLAB代码解决了线性对流方程
清晰的所有
关闭所有
N = 100;%。网格点
达峰时间= 1;%时间= 1
α= 1;%给
h = 0.01;% delX的给定值
解决= 0.005;
maxt =最高温度/解决;%的时间步数
c =α*解决/ h;% c = 0.5
%初始条件
为我= 1:N + 1;
x (i) =(张)* h;
u(我,1)=罪(2 *π* x (i));
结束
%临时在边界
为k = 1: maxt + 1;
u (k) = (N + 1, k);
结束
u = 0 (maxt N);
为k = 1: maxt
为i = 1: N%空间循环
u (i (k + 1) = (i (k) - c / 2 * (u (i + 1 k) - u(张,k)) + (c ^ 2) / 2 * (u (i + 1 k) 2 * u (i (k) +(张,k));
结束
结束

在回答这个问题。

接受的答案

Torsten
Torsten 2018年11月9日
编辑:Torsten 2018年11月9日
N = 100;%。网格点
达峰时间= 1;%时间= 1
α= 1;%给
h = 0.01;% delX的给定值
解决= 0.005;
maxt =最高温度/解决;%的时间步数
c =α*解决/ h;% c = 0.5
u = 0 (N + 1, maxt + 1);
x = 0 (N + 1);
%初始条件
为我= 1:N + 1
x (i) =(张)* h;
u(我,1)=罪(2 *π* x (i));
结束
为k = 1: maxt
情况= u (N, k);
u (1, k + 1) = (1 k) - c / 2 * (u (2 k)情况)+ c ^ 2/2 * (u (2 k) 2 * (k) +情况);
为我= 2:N
u (i (k + 1) = (i (k) - c / 2 * (u (i + 1 k) - u(张,k)) + (c ^ 2) / 2 * (u (i + 1 k) 2 * u (i (k) +(张,k));
结束
uNp2 = u (2 k);
u (N + 1, k + 1) = u (N + 1 k) - c / 2 * (uNp2-u (N, k)) + c ^ 2/2 * (uNp2-2 * u (N + 1, k) + (N, k));
结束
情节(x, u (:, 10), x, u (:, 20), x, u (:, 30), x, u (: 40), x, u (: 50), x, u (:, 60), x, u (:, 70))

更多的答案(0)

类别

找到更多的在可视化在帮助中心和文件交换

标签

下载188bet金宝搏

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!