Lax-wendroff.rar

  • ccc_xy
    了解作者
  • matlab
    开发工具
  • 105KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 0
    下载次数
  • 2021-04-09 14:02
    上传日期
使用Lax-wendroff格式解方程 Ut+3Ux=0 , u(x,0)=x^2 ,令 h=0.05,T=0.01 ,并与精确解x=1,t=1比较
Lax-wendroff.rar
  • boundaryconditions.m
    41B
  • uexact.m
    92B
  • result.fig
    92.8KB
  • main.m
    788B
  • contrast.fig
    14KB
  • u0.m
    27B
内容介绍
clear; close all a = 0; b=1; tfinal = 1; m = 20; h = (b-a)/m; k = 0.01; t = 0; n = fix(tfinal/k); y1 = zeros(m+1,1); y2=y1; x=y1; for i=1:m+1, x(i) = a + (i-1)*h; y1(i) = uexact(t,x(i)); y2(i) = 0; end figure(1); plot(x,y1); hold t = 0; for j=1:n, y1(1)=bc(t); y2(1)=bc(t+k); y2(2:m) = 0.5*9*k^2/(h^2)*(y1(1:m-1)+y1(3:m+1)-2*y1(2:m)) - 0.5*3*k/h*(y1(3:m+1)-y1(1:m-1))+y1(2:m); i = m+1; y2(i) = 0.5*9*k^2/(h^2)*(y1(i-2)+y1(i)-2*y1(i-1)) - 3*k/h*(y1(i)-y1(i-1))+y1(i); t = t + k; y1 = y2; plot(x,y2); pause(0.005) end plot(x,y2,'o') u_e = zeros(m+1,1); % for i=1:m+1 % u_e(i) = uexact(t,x(i)); % end % for i=1:m+1 u_e(1:m+1) = uexact(t,x(1:m+1)); % end max(abs(u_e-y2)) figure(2); plot(x,y2,'o',x,u_e)
评论
    相关推荐
    • Lax-Webdroff2D.zip
      斜激波在平面刚壁上反射问题是具有解析解的二维可压缩无黏流动问题。对它采用具有二阶精度 两步差分格式进行数值求解
    • Lax-wendroff.zip
      利用Lax-wendroff两步差分格式求解一维激波管问题Fortran语言
    • lax-wendroff.rar
      二维Fortran编写的激波在壁面反射程序。采用lax wendroff二阶精度编写
    • Lax-Wendroff_Method.rar
      偏微分方程的Lax-Wendroff方法求解数值解
    • 一步Lax-Wendroff.rar
      C++ 一步lax-wendroff格式求解激波管流动
    • lax2.rar
      Numerical analysis example
    • Lax-Wendroff2D-cpp.rar
      利用Lax-wendroff差分格式求解二维平面激波反射问题 (C++语言版本)
    • lax_wendroff_dispersion:证明Lax_Wendroff方法过度分散的实践模型
      lax_wendroff_dispersion 证明Lax_Wendroff方法过度分散的实践模型 LeVeque的“双曲问题的有限体积方法”中的示例模型,其中当试图近似线性双曲PDE的不连续解时,Lax_Wendroff方法的过度分散是显而易见的。 要运行...
    • 一维激波管问题——计算流体力学(Lax-Wendroff格式)
      本文用python3.6对一维激波管lax激波管问题进行了计算,使用了Lax-Wendoff格式,激波管问题的初值采用了lax激波管的初值(具体见代码),推进时间为0.14s,网格数为500,相较于Lax-Friedrichs格式在稀疏波与激波处...
    • matlabcnhelp.rar
      matlab中文帮助很难找的,快速下载