时间步长与空间步长对稳定性影响_第1页
时间步长与空间步长对稳定性影响_第2页
时间步长与空间步长对稳定性影响_第3页
时间步长与空间步长对稳定性影响_第4页
时间步长与空间步长对稳定性影响_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、对下面的一阶对流扩散方程(1.1)对求解区间Matlab程序function u = yingfeng(a,dt,n,minx,maxx,M)%方程中的常数：a%时间步长：dt%空间节点个数：n%求解区间的左端：minx%求解区间的右端：naxx%时间步的个数：M%求解区间上的数值解：uformat long;h = (maxx-minx)/(n-1);if a0 for j=1:(n+M) u0(j) = IniU(minx+(j-M-1)*h); endelse for j=1:(n+M) u0(j) = IniU(minx+(j-1)*h); endend u1 = u0; for k=

2、1:M if a0 for i=(k+1):n+M u1(i) = -dt*a*(u0(i)-u0(i-1)/h+u0(i); end else for i=1:n+M-k u1(i) = -dt*a*(u0(i+1)-u0(i)/h+u0(i); end end u0 = u1;endif a0 u = u1(M+1):M+n);else u = u1(1:n);endformat long;function ux=IniU(x)format long;if x=-0.1 ux=10*(x+0.1); else ux=0; endelse if x plot(u)u=yingfeng(1,0

3、.005,101,0,1,100);plot(u) u=yingfeng(1,0.006,101,0,1,100); plot(u) u=yingfeng(1,0.007,101,0,1,100); plot(u) u=yingfeng(1,0.008,101,0,1,100); plot(u)u=yingfeng(1,0.009,101,0,1,100);plot(u) u=yingfeng(1,0.01,101,0,1,100); plot(u) u=yingfeng(1,0.011,101,0,1,100); plot(u) u=yingfeng(1,0.012,101,0,1,100); plot(u) u=yingfeng(1,0.013,101,0,1,100); plot(u) u=yingfeng(1,0.014,101,0,1,100); plot(u) u=yingfeng(1,0.017,101,0,1,100); plot(u) u=yingfeng(1,0.020,101,0,1,100); plot(u) u=yingfeng(1,1,101,0,1,100); plot(u)由以上得到的图像我们可以

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论