许鑫华老师课堂作业答案.doc

解：布拉格方程：2dsinθ=λ 点阵参数a与d的关系为：d=1/ 有关参数列表如下 1 2 3 4 5 6 7 8 θi 0.3083π 0.3094π 0.3206π 0.3219π 0.3765π 0.3783π 0.4387π 0.4429π 0.3650 0.3557 0.3143 0.3055 0.1379 0.1335 0.03197 0.02762 1.60699 Yi 0.4075 0.4074 0.4171 0.4074 0.4077 0.4077 0.4078 0.4078 3.2704 Xi2 0.1230 0.1265 0.09878 0.09333 0.01902 0.01782 1.022×10-3 7.629×10-4 0.4802349 Yi2 0.1661 0.1660 0.1740 0.1660 0.1662 0.1662 0.1663 0.1663 1.3371 XiYi 0.1469 0.1449 0.1311 0.1246 0.05622 0.05443 0.01304 0.01126 0.68245 取Y=a+bX为一条直线，对于每个数据点（Xi,Yi）则每个点的偶然误差为ei，有ei=a+b Xi,-Yi 所有实验数据点误差平方和为： =（a+b X1,-Y1）2+（a+b X2,-Y2）2+（a+b X3,-Y3）2+（a+b X4,-Y4）2+（a+b X5,-Y5）2+（a+b X6,-Y6）2+（a+b X7,-Y7）2+（a+b X8,-Y8）2 根据最小二乘法原理，误差平方和最小的直线为最佳直线。 求最小值的条件是： =0及=0 可得 故可得所求方程为： a=0.409039-5.746×10-4× 第一题：写出边值问题 解：划分网格，对节点及界点进行编号，如图所示  SHAPE \* MERGEFORMAT  节点 1 2 3 4 5 6 7 8 9 U1 U2 U3 U4 U5 U6 U7 U8 U9 U11 U21 U31 U12 U22 U32 U13 U23 U33 其中i=0,1,2,3,4 j=0,1,2,3,4 对于每一个（xi，yi）∈D，利用数值微分公式可得 =[u(xi+1,yi)-2u(xi,yi)+ u(xi-1,yi)]+Rx =[u(xi,yi+1)-2u(xi,yi)+ u(xi,yi-1)]+Ry 由已知条件可知：+=0 且h1=和h2=1 ∴+= u(xi+1,yi) + u(xi-1,yi)+ u(xi,yi+1) + u(xi,yi-1) -4u(xi,yi)=0 再由边界条件可得： -16 -20+6 -22+8 -12 -34+8 -20-6 -40 -36-6 -34-8 -44 -22-8 -36+6 对于1,3,7,9四个非正则内点，近似为边界上的点： = = = = 对于2,4,5,6,8,五个正则内点，利用差分方程可求得： u1+u3+u4+u5-4u2=0 +u1+u5+u7-4u4=0 u2+u4+u6+u8-4u5=0 u5+u3++u9-4u6=0 u5+u7+u9+-4u8=0 代入边值以及非正则内点近似值可得 4u2 -u4-u5=-56+16 4u4 -u5 =-60 4u5-u2-u4-u6-u8=0 4u6-u5= -108 4u8-u5= -100-16 = 由以上三个方程组联立可求得边值问题。 第二题：写出边值问题 取步长h1=h2=1 解：划分网格，节点与编号如下 节点 5 6 7 10 11 12 U5 U6 U7 U10 U11 U12 U11 U21 U31 U12 U22 U32 由边界条件 可得十个边界点的值为： u1= u2=1 u3= u4=2 u9=2 u8=u13=u14=u15=u16=0 其中u1、u2、u3均指x轴上的点。  对于每一个（xi，yi）∈D，利用数值微分公式可得 =[u(xi+1,yi)-2u(xi,yi)+ u(xi-1,yi)]+Rx =[u(xi,yi+1)-2u(xi,yi)+ u(xi,yi-1)]+Ry 由已知条件可知：+=0 且h1=和h2=1 ∴+= u(xi+1,yi) + u(xi-1,yi)+ u(xi,yi+1) + u(xi,yi-1) -4u(xi,yi)=0 节点5、6、7、10、11、12均为正则内点，利用差分方程可求得： u1+u4+u6+u10-4u5=0 u2+u5+u7+u11-4u6=0 U3+u6+u8+u12-4u7=0 u5+u9+u11+u14-4u10=0 u6+u10+u12+u15-4