热传导方程(共5页).doc

《热传导方程(共5页).doc》由会员分享，可在线阅读，更多相关《热传导方程(共5页).doc（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上*三点对称计算格式*#include #include#include using namespace std;void main()ofstream outFile;outFile.open(3PT_SYM.txt); int m=5,n=10;/时间坐标分数n，和空间坐标分数m *m，n可以修改double T1111,err1111,x1111;/存放差分法近似结果T【n】【m】,误差err【】【】及真实解x【】double a=1,b=1,c=1,c1=100,c2=100; /空间坐标界a，和时间坐标界b,导热系数c,边界条件c1,c2double delt

2、aX=a/n; /空间步长double deltaT=b/m; /时间步长 double s=c*deltaT/(deltaX*deltaX);double r=1-2*s;/boundary conditionsfor(int i=0;i=n;i+)Ti0=exp(i*deltaT);Tim=exp(i*deltaT+m*deltaX);/Generate first rowfor (int i=1;im;i+)T0i=exp(i*deltaX);/generate remaining rows of Tfor (int i=1;i=n;i+)for (int j=1;jm;j+)Tij=s

3、*(Ti-1j+1+Ti-1j-1)+r*Ti-1j;cout i*deltaT j*deltaX 精确值xik 近似值Tik 误差errik endl;outFile i*deltaT j*deltaX 精确值xik 近似值Tik 误差errik endl;double maxerr=0;for(int i=0;i=n;i+)for(int k=0;kTik?xik-Tik:Tik-xik;couti*deltaX k*deltaT xik Tik errikendl;outFilefixed;outFile.precision(6);outFilei*deltaX k*deltaT xik

4、 Tik errikmaxerr)maxerr=errik;cout最大误差为maxerrendl;coutendl;outFileendlendl;/output the outcomefor (int i=0;i=n;i+)cout按时间推进第i层温度近似值T endl;for (int j=0;j=m;j+)coutTij ;coutendl;/Out put f_filefor (int i=0;i=n;i+)/for (int j=0;j=m;j+) /此时m=5/outFilefixed;outFile.precision(6);/outFile按时间推进第i层近似值温度T end

5、l;outFileTi0 Ti1 Ti2 Ti3 Ti4 Ti5endl;/outFile.close();*五点对称计算格式*#include #include#include using namespace std;void main()ofstream outFile;outFile.open(5PT_SYM.txt);int m=10,n=10;/时间坐标分数n，和空间坐标分数m *m，n可以修改double T1111,err1111,x1111;/存放差分法近似结果T【n】【m】,误差err【】【】及真实解x【】【】double a=1,b=1,c=1,c1=100,c2=100;

6、 /空间坐标界a，和时间坐标界b,导热系数c,边界条件c1,c2double deltaX=a/n; /空间步长double deltaT=b/m; /时间步长double s=c*deltaT/(deltaX*deltaX);/double r=1-2*s;/boundary conditionsfor(int i=0;i=n;i+)Ti0=c1;Tim=c2;/ Ti0=exp(i*deltaT);/ Tim=exp(i*deltaT+m*deltaX);/Generate first rowfor (int i=1;im;i+)T0i=0;/T0i=exp(i*deltaX);/gene

7、rate remaining rows of Tfor (int i=1;i=n;i+)Ti1=11/12*s*Ti-10+(1-5/3*s)*Ti-11+0.5*s*Ti-12+1/3*s*Ti-13-1/12*s*Ti-14;for (int i=1;i=n;i+) Tim-1=11/12*s*Ti-1m+(1-5/3*s)*Ti-1m-1+0.5*s*Ti-1m-2+1/3*s*Ti-1m-3-1/12*s*Ti-1m-4;for (int i=1;i=n;i+)for (int j=2;jm-1;j+)Tij=s*(-1/12*Ti-1j-2+4/3*Ti-1j-1+4/3*Ti-1j

8、+1-1/12*Ti-1j+2)+(1-5/2*s)*Ti-1jcout i*deltaT j*deltaX 精确值xik 近似值Tik 误差errik endl;outFile i*deltaT j*deltaX 精确值xik 近似值Tik 误差errik endl;double maxerr=0;for(int i=0;i=n;i+)for(int k=0;kTik?xik-Tik:Tik-xik;couti*deltaX k*deltaT xik Tik errikendl;outFilefixed;outFile.precision(6);outFilei*deltaX k*delta

9、T xik Tik errikmaxerr)maxerr=errik;cout最大误差为maxerrendl;coutendl;outFileendlendl;/output the outcomefor (int i=0;i=n;i+)cout按时间推进第i层温度近似值T endl;for (int j=0;j=m;j+)coutTij ;coutendl;/Out put f_filefor (int i=0;i=n;i+)/for (int j=0;j=m;j+) /此时m=10/outFilefixed;outFile.precision(6);/outFile按时间推进第i层近似值温度T endl;outFileTi0 Ti1 Ti2 Ti3 Ti4 Ti5 Ti6 Ti7 Ti8 Ti9 Ti10endl;/outFile.close();专心-专注-专业