陶文铨数值传热学第二版第五章(共6页).docx

精选优质文档-倾情为你奉上精确解：p=1,5,10; x=0:1/19:1;for i=1:1:3 for j=1:1:20 y(i,j)=(exp(p(1,i)*19*x(1,j)-1)/(exp(p(1,i)*19)-1); end plot(x,y(i,:); hold on;end由题对中心差分、一阶迎风、混合格式进行模块编程：他们之间可以通用，只需更改ae关于p的函数即可：程序如下:(1)中心差分p=1,5,10;for i=1:1:3 ae=1-0.5*p(1,i); aw=p(1,i)+ae; ap=ae+aw;for i=1:1:18 for j=1:1:20 a(i,j)=0; endendfor i=1:1:18 j=i; a(i,j)=aw; a(i,j+1)=-ap; a(i,j+2)=ae;endfor i=1:1