隐式差分法[学习资料]

1、隐式差分法程序及结果#include int Fun(float T);void main()int a=1;float t;for (a;a=1;)printf(输入时间步长：);scanf(%f,&t);a=Fun(t);int Fun(float T)int flag;float a,b,c,h1111,f10=0,A10=0,B10=0,W10=0,Z10=0;int i,k;a=c=1000;b=-(2000+1000/T);for(k=0;k11;k+)h0k=10; h10k=0;for(i=1;i11;i+)hi0=0; for(k=0;k10;k+) for(i=1;i9;i

2、+)fi=-(1000*hik/T);f1=f1-a*h0k+1;f9=f9-c*h10k+1; W1=b;A1=c/W1;for(i=2;i10;i+)Bi=a;Wi=b-Ai-1*Bi;Ai=c/Wi;Z1=f1/W1;for(i=2;i10;i+)Zi=(fi-a*Zi-1)/Wi);h9k+1=Z9;for(i=1;i9;i+)h9-ik+1=(Z9-i-A9-i*h10-ik+1); printf(当时间步长%.2f时：n,T);printf( i,k t i=0 t i=1 t i=2 t i=3 t i=4 t i=5 t i=6 t i=7 t i=8 t i=9 t i=10

3、 t);k=0;for(k=0;k11;k+)printf(k=%dt,k);for(i=0;i11;i+)printf(%.2ft,hik);printf(继续？（1/0）);scanf (%d,&flag);return flag;输入时间步长：0.1当时间步长为0.10时： i,k i=0 i=1 i=2 i=3 i=4 i=5 i=6 i=7 i=8 i=9 i=10k=0 10.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00k=1 10.00 0.84 0.07 0.01 0.00 0.00 0.00 0.00 0.00 0.

4、00 0.00k=2 10.00 1.55 0.19 0.02 0.00 0.00 0.00 0.00 0.00 0.00 0.00k=3 10.00 2.15 0.34 0.05 0.01 0.00 0.00 0.00 0.00 0.00 0.00k=4 10.00 2.67 0.51 0.08 0.01 0.00 0.00 0.00 0.00 0.00 0.00k=5 10.00 3.12 0.70 0.13 0.02 0.00 0.00 0.00 0.00 0.00 0.00k=6 10.00 3.50 0.89 0.18 0.03 0.01 0.00 0.00 0.00 0.00 0.

5、00k=7 10.00 3.84 1.08 0.25 0.05 0.01 0.00 0.00 0.00 0.00 0.00k=8 10.00 4.14 1.27 0.32 0.07 0.01 0.00 0.00 0.00 0.00 0.00k=9 10.00 4.41 1.46 0.39 0.09 0.02 0.00 0.00 0.00 0.00 0.00k=10 10.00 4.64 1.64 0.48 0.12 0.03 0.01 0.00 0.00 0.00 0.00继续？（1/0）1输入时间步长：0.2当时间步长为0.20时： i,k i=0 i=1 i=2 i=3 i=4 i=5 i

6、=6 i=7 i=8 i=9 i=10k=0 10.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00k=1 10.00 1.46 0.21 0.03 0.00 0.00 0.00 0.00 0.00 0.00 0.00k=2 10.00 2.55 0.53 0.10 0.02 0.00 0.00 0.00 0.00 0.00 0.00k=3 10.00 3.37 0.89 0.21 0.04 0.01 0.00 0.00 0.00 0.00 0.00k=4 10.00 4.02 1.26 0.34 0.08 0.02 0.00 0.00

7、0.00 0.00 0.00k=5 10.00 4.53 1.62 0.49 0.13 0.03 0.01 0.00 0.00 0.00 0.00k=6 10.00 4.94 1.95 0.66 0.20 0.05 0.01 0.00 0.00 0.00 0.00k=7 10.00 5.28 2.27 0.83 0.27 0.08 0.02 0.01 0.00 0.00 0.00k=8 10.00 5.57 2.56 1.01 0.35 0.11 0.03 0.01 0.00 0.00 0.00k=9 10.00 5.81 2.83 1.19 0.44 0.15 0.05 0.01 0.00

8、0.00 0.00k=10 10.00 6.02 3.08 1.37 0.54 0.19 0.06 0.02 0.01 0.00 0.00继续？（1/0）1输入时间步长：0.5当时间步长为0.50时： i,k i=0 i=1 i=2 i=3 i=4 i=5 i=6 i=7 i=8 i=9 i=10k=0 10.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00k=1 10.00 2.68 0.72 0.19 0.05 0.01 0.00 0.00 0.00 0.00 0.00k=2 10.00 4.23 1.55 0.53 0.17 0.0

9、5 0.02 0.00 0.00 0.00 0.00k=3 10.00 5.19 2.30 0.93 0.35 0.12 0.04 0.01 0.00 0.00 0.00k=4 10.00 5.83 2.94 1.34 0.57 0.22 0.09 0.03 0.01 0.00 0.00k=5 10.00 6.28 3.48 1.74 0.81 0.35 0.14 0.06 0.02 0.01 0.00k=6 10.00 6.62 3.92 2.11 1.05 0.49 0.22 0.09 0.04 0.01 0.00k=7 10.00 6.89 4.30 2.46 1.30 0.65 0.3

10、1 0.14 0.06 0.01 0.00k=8 10.00 7.10 4.62 2.77 1.55 0.81 0.40 0.19 0.08 0.02 0.00k=9 10.00 7.27 4.89 3.05 1.78 0.98 0.51 0.25 0.11 0.03 0.00k=10 10.00 7.42 5.13 3.31 2.00 1.15 0.62 0.31 0.14 0.04 0.00继续？（1/0）1请输入时间步长：0.52当时间步长为0.52时： i,k i=0 i=1 i=2 i=3 i=4 i=5 i=6 i=7 i=8 i=9 i=10k=0 10.00 0.00 0.00

11、 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00k=1 10.00 2.74 0.75 0.21 0.06 0.02 0.00 0.00 0.00 0.00 0.00k=2 10.00 4.30 1.61 0.56 0.18 0.06 0.02 0.01 0.00 0.00 0.00k=3 10.00 5.26 2.38 0.97 0.37 0.14 0.05 0.02 0.01 0.00 0.00k=4 10.00 5.90 3.03 1.40 0.60 0.25 0.10 0.04 0.01 0.00 0.00k=5 10.00 6.35 3.57 1.81 0.86 0.38 0.16 0.06 0.02 0.01 0.00k=6 10.00 6.68 4.01 2.20 1.11 0.53 0.24 0.10 0.04 0.01 0.00k=7 10.00 6.94 4.38 2.54 1.37 0.69 0.33 0.15 0.06 0.02 0.00k=