热传导方程的差分格式(共15页).doc

上传人:飞****2 文档编号:15004615 上传时间:2022-05-10 格式:DOC 页数:15 大小:664.50KB
返回 下载 相关 举报
热传导方程的差分格式(共15页).doc_第1页
第1页 / 共15页
热传导方程的差分格式(共15页).doc_第2页
第2页 / 共15页
点击查看更多>>
资源描述

《热传导方程的差分格式(共15页).doc》由会员分享,可在线阅读,更多相关《热传导方程的差分格式(共15页).doc(15页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、精选优质文档-倾情为你奉上热传导方程的差分格式 上机实习报告二零一四年五月专心-专注-专业一维抛物方程的初边值问题分别用向前差分格式、向后差分格式、六点对称格式,求解下列问题:在时刻的数值解,并与解析解。1差分格式形式 设空间步长, 时间步长,网比.(1) 向前差分格式向前差分格式,即 (1) ,其中,以表示网比。(1)式可改写成如下:此格式为显格式。其矩阵表达式如下:(2) 向后差分格式向后差分格式,即(2)其中(2)式可改写成此种差分格式被称为隐格式。其矩阵表达式如下:(3) 六点对称格式六点差分格式: (3)将(3)式改写成其矩阵表达式如下:2利用MATLAB求解问题的过程 对每种差分格

2、式依次取,用MATLAB求解并图形比较数值解与精确解,用表格列出不同剖分时的误差。(1)向前差分格式(2)向后差分格式(3)六点对称格式3 方法总结及分析六点对称方法克服了向前差分格式与向后差分格式方法误差偏大的缺点,更接近准确数值本文向前差分格式,向后差分格式及六点差分格式都是使用三对角系数矩阵,计算简单。根据matlab作,特别明显的是,时,图像解析解波动特别大,与真解差距很大。附件程序(1)向前差分格式源程序:functione=forward(h,t,T)N=1/h;x=zeros(1,N+1);for i=1:Nx(i+1)=i*h;endu=sin(pi.*x);r=t/(h*h)

3、;for j=1:T/t;u=for_climb(u,r,t);endplot(x,u,o)holdu_xt=exp(-pi*pi*T)*sin(pi.*x);plot(x,u_xt,r)e=u_xt-u;function b=for_climb(a,r,t)L=length(a);b=zeros(1,L);for i=1:L-2b(i+1)=r*a(i+2)+(1-2*r)*a(i+1)+r*a(i);end(2)向后差分格式源程序:function e=back(dx,dt,T)M=1/dx; N=T/dt; u1=zeros(1,M+1);x=1:M-1*dx; u1(2:M)=sin(

4、pi.*x);r=dt/dx/dx;r2=1+2*r;for i=1:M-1A(i,i)=r2;if i1 A(i-1,i)=-r;A(i,i-1)=-r;end end u=zeros(N+1,M+1);u(N+1,:)=u1;for k=1:N b=u(N+2-k,2:M); u(N+1-k,2:M)=inv(A)*b;end uT=u(1,:); x1=0,x,1; plot(x1,uT,o)hold u_xt=exp(-pi*pi*T)*sin(pi.*x1);plot(x1,u_xt,r)e=u_xt-uT;(3)六点对称格式:源程序:function e=six(dx,dt,T)

5、M=1/dx;N=T/dt; u1=zeros(1,M+1);x=1:M-1*dx; u1(2:M)=sin(pi*x);r=dt/dx/dx;r2=2+2*r;r3=2-2*r;for i=1:M-1 A(i,i)=r2; if i1 A(i-1,i)=-r; A(i,i-1)=-r;end end for i=1:M-1 B(i,i)=r3;if i1 B(i-1,i)=r;B(i,i-1)=r; end end u=zeros(N+1,M+1);u(N+1,:)=u1;for k=1:N b=B*(u(N+2-k,2:M); u(N+1-k,2:M)=inv(A)*b;end uT=u(1,:); x1=0,x,1; plot(x1,uT,o)hold u_xt=exp(-pi*pi*T)*sin(pi.*x1);plot(x1,u_xt,r)e=u_xt-uT;

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 教案示例

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁