西安电子科技大学平时作业-计算方法

习题11. x1=4.8675 x1有5位有效数字；x2=4.08675 x2有6位有效数字；X 3=0.08675 x3有4位有效数字；x4=96.4730 x4有6位有效数字； X 5=96×105 x5有2位有效数字；x6=0.00096 x6有2位有效数字。

8.解: y n =5y n-1-2 n=1,2, (1)y 0= 3在计算y 0时有舍入误差，设为e 0,并设求得的y 0的近似值y 0,即e 0= y 0 -y 0，所以，yn = 5y n-1-2 n=1,2,… (2) y 0=y 0-e 0由（1）-（2）得：y n - y n =5(y n-1- y n-1)所以y n -y n =5ne 0 n=1,2,…所以e 10=510e 0=510( 3 -1.73)=20027.42 所以这个计算过程不稳定。

10.解：f(x)=8x 5-0.4x 4+4x 3-9x+1=(8x 4-0.4x 3+4x 2-9)x+1=((8x 3-0.4x 2+4x)x-9)x+1=(((8x 2-0.4x+4)x-9)x+1 =((((8x-0.4)x+4)x-9)x+1b 0=8；b 1=8x-0.4=8×3-0.4=23.6； b 2= b 1x+4=23.6×3+4=74.8； b 3= b 2x=74.8×3=224.4;b 4= b 3x-9=224.4×3-9=664.2; b 5= b 4x+1=664.2×3+1=1993.6; 所以f(3)= b 5=1993.6.8 -0.4 4 0 -9 1X=3 24 70.8 224.4 673.2 1992.6 8 23.6 74.8 224.4 664.2 1993.6 所以f(3)=1993.6. 习题21. 证明：令f(x)=1-x-sinx,则f ′(x)=-1-cosx>0,所以f （x ）在区间[0,1]中连续且严格单调递增。

寻找多数元素的算法试验报告一． 问题描述令[1...]A n 是一个整数序列，A 中的整数a 如果在A 中出现的次数多于/2n ⎢⎥⎣⎦，那么a 称为多数元素。

例如，在序列1,3,2,3,3,4,3中，3是多数元素，因为7个元素中它出现4次。

现在我们就要讨论如何利用计算机来找出一个序列中的多数元素，当然这个多数元素要么不存在，要么就只有一个。

二． 算法描述有几种方法可以解决这个问题，蛮力方法就是把序列中的每个元素和其他每个元素比较，并且对每个元素计数，如果某个元素的计数大于/2n ⎢⎥⎣⎦，就可以断言它是多数元素；否则，在序列中就没有多数元素。

但是这样的比较次数是2(1)/2()n n n -=Θ，这种方法的代价太昂贵。

另一种比较有效的算法是先对这些元素排序，并且计算每个元素在序列中出现多少次。

这在最坏情况下的代价是(log )n n Θ。

因为在最坏情况下，排序这一步需要(log )n n Ω次比较（选择合并排序或者快速排序）。

还有一种方法就是先对序列排序，再寻找中间元素，就是第/2n ⎢⎥⎣⎦元素。

因为多数元素，在排序的序列中一定是中间元素，可以扫描这个序列来测试中间元素是否确实是多数元素。

由于多数元素可以在()n Θ的时间内找到，这个方法要花费()n Θ的时间，但是对于大数据量的序列，中项寻找算法的时间花费非常大，并且算法很复杂。

上述的几种方法都存在各种弊端，这里我们讨论一种非常漂亮的算法，它用的比较次数要少得多。

在介绍这种由归纳法导出的递归算法之前，我们首先需要了解一个观察结论：在原序列中去除两个不同的元素后，那么在原序列中的多数元素在新序列中还是多数元素。

这个观察结论支持下述寻找多数元素候选者的过程。

将计数器置1，并令[1]c A =，从[2]c A =开始，逐个地扫描元素，如果被扫描的元素和c 相等，则计数器加1；如果元素不等于c ，则计数器减1；如果所有的元素都已经扫描完毕并且计数器大于0，那么返回c 作为多数元素的候选者（注意：这里得到的仅仅是候选者，是否是真正的多数元素还有待验证）。

西安电子科技大学出版社《计算方法》任传祥等编著第九章计算方法上机参考答案实验一，算法一#include <stdio.h>#include <math.h>double I0=log(6)/log(5),I1;int n=1;main (){while(1){I1=1.0/(n)-I0*5.0;printf("%d %lf\n", n,I1);if(n>=20)break;elseI0=I1;n++;}}实验一，算法二#include <stdio.h>#include <math.h>double I0=(1/105.0+1/126.0)/2,I1;int n=20;main (){printf("%d %lf\n", n,I0);while(1){I1=1.0/(5.0*n)-I0/5.0;printf("%d %lf\n", n-1,I1);if(n<2)break;elseI0=I1;n--;}}实验二，二分法#include <stdio.h>#include <math.h>#define esp 1e-3double f(double x);main (){double a=1,b=2,x;while(fabs(b-a)>esp){x=(a+b)/2;printf("x=%lf\n",x);if(f(x)==0)break;elseif(f(x)*f(a)<0)b=x;elsea=x;}}double f(double x){return pow(x,3)-x-1;}实验二，牛顿迭代法#include<stdio.h>#include<math.h>double f(double x);double f1(double x);#define esp 1e-3void main(){double x0 = 1.5, x1;x1 = x0 - f(x0) / f1(x0);printf("x=%lf\n", x1);x0 = x1;x1 = x0 - f(x0) / f1(x0);printf("x=%lf\n", x1);while (fabs(x1 - x0)>esp){x0 = x1;x1 = x0 - f(x0) / f1(x0);printf("x=%lf\n", x1);} }double f(double x){return pow(x, 3) - x - 1;} double f1(double x){return 3 * x*x - 1;}弦割法#include<stdio.h>#include<math.h>double f(double x);#define esp 1e-3void main(){double x0 = 1.5, x1=2.0,x2;do{ x2=x1 - (x1-x0)*f(x1) /(f(x1)-f(x0));x0=x1;x1=x2;printf("x=%lf\n", x1);}while (fabs(x1 - x0)>esp);{printf("x=%lf\n", x1);}}double f(double x){return pow(x, 3) - x - 1;}实验3#include <stdio.h>/*列主元高斯消去法*/#include <math.h>float x[3],temp,max;float A[3][4]={10,-2,-1,3,-2,10,-1,15,-1,-2,5,10},c[3][4]={10,-2,-1,3,-2,10,-1,15,-1,-2,5,10}; int n=3,i,k,j,m;void main(){for(i=0;i<n;i++){max=A[i][i];k=i;for(j=j+1;j<n;j++){{max=fabs(A[j][i]);k=j;}}if(k!=i){for(j=i+1;j<=n;j++){temp=A[i][j];A[i][j]=A[k][j];A[k][j]=temp;}}for(j=i+1;j<n;j++)for(m=i+1;m<=n;m++){c[j][m]=c[j][m]+(-c[j][i]/c[i][i])*c[i][m];}}for(i=n-1;i>=0;i--){temp=0.0;for(j=n-1;j>=i+1;j--)temp=temp+c[i][j]*x[j];x[i]=(c[i][n]-temp)/c[i][i];}printf("x[1]=%f\nx[2]=%f\nx[3]=%f\n",x[0],x[1],x[2]);实验四，拉格朗日插值#include<stdio.h>int n=5,i,j;double l,L=0,X=0.5;main(){double x[5]={0.4,0.55,0.65,0.8,0.9};doubley[5]={0.41075,0.57815,0.69675,0.88811,1.02652}; for(i=0;i<n;i++){l=y[i];for(j=0;j<n;j++){if(j!=i)l=l*(X-x[j])/(x[i]-x[j]); } L=L+l;}printf("%lf\n",L);return 0;} X=0.5 X=0.7 X=0.85牛顿插值法#include<stdio.h>#include<math.h>main(){double x[5]={0.4,0.55,0.65,0.8,0.9};doubley[5]={0.41075,0.57815,0.69675,0.88811,1.02652};int n=5,i,j;double z;printf("input z\n");scanf("%lf",&z);double a[5][5];for(i=0;i<5;i++)a[i][0]=y[i];for(i=1;i<5;i++)for(j=i;j<5;j++)a[j][i]=(a[j][i-1]-a[j-1][i-1])/(x[j]-x[j-i]);double N=a[0][0],temp=1.0;for(i=1;i<n;i++){temp=temp*(z-x[i-1]);N=N+a[i][i]*temp;}printf("N=%lf\n",N);return 0;}实验五曲线拟合#include <stdio.h>#include <math.h>float x[5]={1,2,3,4,5};float y[5]={7,11,17,27,40};float A[2][3],c[2][3];float z[2],temp,max;int i,j,k,m;int n=2;void main(){for(i=0;i<5;i++){c[0][0]=A[0][0]+=1;c[0][1]=A[0][1]+=x[i];c[0][2]=A[0][2]+=y[i];c[1][0]=A[1][0]+=x[i];c[1][1]=A[1][1]+=x[i]*x[i];c[1][2]=A[1][2]+=x[i]*y[i];}/* for(i=0;i<2;i++){printf(" %lf %lf %lf\n",A[i][0],A[i][1],A[i ][2]);}*/for(i=0;i<n;i++){max=A[i][i];k=i;for(j=j+1;j<n;j++){if(fabs(A[j][i])>max){max=fabs(A[j][i]);k=j;}} if(k!=i){for(j=i+1;j<=n;j++){temp=A[i][j];A[i][j]=A[k][j];A[k][j]=temp;}}for(j=i+1;j<n;j++)for(m=i+1;m<=n;m++){c[j][m]=c[j][m]+(-c[j][i]/c[i][i])*c[i][m];}}for(i=n-1;i>=0;i--){temp=0.0;for(j=n-1;j>=i+1;j--)temp=temp+c[i][j]*z[j];z[i]=(c[i][n]-temp)/c[i][i];}printf("a=%f\nxb=%f\n",z[0],z[1]); }实验六数值积分/*梯形*/#include<stdio.h>#include<math.h> double f(double x); main(){double x[10],y[10];double h,b=1,a=0,I;int n,i;printf("n\n");scanf("%d",&n);h=(b-a)/n;for(i=0;i<=n;i++){x[i]=a+(i*h);y[i]=f(x[i]);}I=f(a)+f(b);for(i=1;i<=n-1;i++){I=I+2*y[i];}I=(h/2)*I;printf("%lf",I);}double f(double x){double f;f=1.0/(1.0+(x*x));return(f);}/*辛普森*/#include<stdio.h>#include<math.h>double f(double x);main(){double x[30],y[30];double h,b=1,a=0,I;int n,i;printf("n\n");scanf("%d",&n);//点乘2扩展h=(b-a)/n;x[10]=1;y[10]=f(x[10]);for(i=0;i<n;i++){x[2*i]=a+(i*h);y[2*i]=f(x[2*i]);x[2*i+1]=a+(i+(1.0/2.0))*h;y[(2*i)+1]=f(x[(2*i)+1]);}I=f(a)+f(b);for(i=0;i<n;i++){I=I+4*y[(2*i)+1];}for(i=1;i<n;i++){I=I+2*y[2*i];}I=(h/6)*I;printf("%lf\n",I);}double f(double x){double f;f=1.0/(1.0+(x*x));return(f);}/*梯形*//*辛普森*/。

2022年西安电子科技大学公共课《大学计算机基础》期末试卷B(有答案）一、单项选择题1、下列数值不可能是八进制数的是（）A.137B.138C.276D.10112、二进制数110110111转换为十六进制数是（）A. (1B7)16B. (8B1)16C. (DB8)16D. (DB1)163、与二进制数01000011等值的十进制数是（）A.19B.35C.67D.1314、二进制数10111101等于十进制数（）A.187B.189C.191D.1935、下列不属于计算机病毒的特征的是（）A.免疫性B.寄生性C.传染性D.破坏性6、汉字国标码在两个字节中各占用（）位二进制编码。

A. 6B. 7C. 8D. 97、目前，微型计算机存储器容量最大的是（）A.内存B.光盘C.硬盘D.优盘8、Windows 7中的“任务栏”（）A.只能改变位置不能改变大小B.只能改变大小不能改变位置C.既能改变位置也能改变大小D.既不能改变位置也不能改变大小9、在Windows 7的对话框中，下拉框是一个方框，并在右边有一个（）A.向下的黑三角标志B.横向的省略号标志C.颜色变灰的标志D.向右的三角形标志10、在微型计算机中，操作系统的作用是（）A.把源程序编译成目标程序B.便于进行文件夹管理C.管理计算机系统的软件和硬件资源D.管理高级语言和机器语言11、以下有关Windows 7用户账户管理的叙述中，正确的是（）A.新建用户账户时，必须设置账户密码B.可创建的账户类型有标准用户和管理员两种C.已创建的用户账户类型不能被更改D.标准用户可以安装新软件12、Windows 7把所有的系统环境设置功能都统一到了（）A.“我的电脑”B.“附件”C.“控制面板”D.“资源管理器”13、在Windows 7默认配置下，可以将硬盘选中文件永久删除的操作是（）A.按"Delete"键B.按"Shift+Delete”组合键C.单击鼠标右键，在快捷菜单中选择“删除”命令D.单击当前资源管理器窗口的“文件”菜单，并选择“删除”命令14、编辑Word文档时,为文档的文本提供解释需要插入脚注,脚注一般出现在（）A.文档中每一页的顶部B.文档中每一页的底端C.整个文档的结尾D.文档中每一节的结尾15、在Word中，能同时显示水平标尽和垂直标尺的“视图”方式是（）A.普通B.大纲C.Web版式D.页面16、在Word中插入图片，其默认的环绕方式是（）A.嵌入型B.四周型C.紧密型D.浮于文字下方17、下列关于Word文档分栏叙述正确的是（）A.各栏的宽度可以不同B.各栏的间距是固定的C.最多可以设4栏D.各栏之间不能添加分隔线18、在Word中.段落通常是（）A.以句号结束B.以输人回车键结束C.以空格结束D.以分节符结束19、在Word 2010编辑文本时，可以在标尺上直接进行（）A.段落首行缩进操作B.建立表格C.嵌入图片D.分栏操作20、在Word 2010中，需将当前文档保存到其他目录，正确的操作是（）A.单击快速访问工具栏上的“保存”按钮B.选择“文件”→“保存”命令C.选择“文件”→“退出”命令D.选择“文件”→“另存为”命令21、在Excel 2010中，设A10单元格的数字格式为整数，若输入：33.51，则显示（）A.33.51B.33C.34D.ERROR22、在Excel 2010中，当向Excel 2010工作表单元格输入公式时，使用单元格地址D$2引用D列2行单元格，该单元格的引用称为（）A.交叉地址引用B.混合地址引用C.相对地址引用D.绝对地址引用23、在Excel 2010中，计算平均值的函数是（）A.SUMB.MAXC.COUNTD.AVERAGE24、在Excel 2010中，不能实现为单元格定义名称的是（）A.单击工作表左上角名称框，快速定义名称B.单击单元格，输入新名称C.使用“公式”→“定义的名称”命令，在“新建名称”对话框中创建新名称D.使用“公式”→“名称管理器”命令，在“名称管理器”对话框中创建名称25、在Excel 2010单元格中，将数值-100通过“单元格”格式设置后，下列显示正确的负数是（）A.<100>B.［100］C.(100)D.{100}26、在Power Point2010中，下列有关幻灯片主题叙述错误的是（）A主题可以应用于所有幻灯片B主题可以应用于指定幻灯片C主题可以在幻灯片母版中更改D可以对已使用的主题进行更换27、在PowerPoint 2010中，用于实现超链接功能的“动作按钮”位于（）A."插入”-"SmartArt"B."插入”→“图片”C.“插入”-“剪贴画”D.“插入”→“形状”28、HTTP 的中文含义是（）A.超文本传输协议B.文件传输协议C.传输控制协议D.用户数据报协议29、在互联网接入技术ADSL中，通信信道一般被分为上行和下行两部分，其特点是（）A.上行带宽高于下行B.下行带宽高于上行C.上行带宽和下行带宽相同D.带宽随着网速动态发生变化30、以下网络协议中，属于应用层协议的是（）A.TCPB.FTPC.IPD.ARP.二、填空题31、在Windows 7磁盘管理中，优化磁盘包括：整理磁盘碎片和磁盘_______32、一个字节最多可以表示__________个不同的编码。

《计算方法》课程平时作业（2010-2011学年第一学期）作业(考试前交, 给出证明或计算过程、计算程序及计算结果)1. 对向量()12Tn x x x x = 定义1211,max ,nk k k nk x x xx x ∞≤≤====∑设A 是n n ⨯矩阵，规定1111max x A Ax ==，1max x A Ax∞∞∞==，2221max x A Ax ==证明111112max (),max (),.nnkj jk j nj nk k T A a Aa A A A λ∞≤≤≤≤=====∑∑列范数行范数是最大特征值2. 用简单迭代法（即不动点迭代法）求方程32210200x x x ++-=在1x =附近的根. 要求给出不动点方程、程序和运行结果. 3. 用Newton 迭代法求方程32210200x x x ++-=的一个正根，计算结果精确到7位有效数字. 要求给出程序和运行结果. 4. 用牛顿迭代法求方程310x x --=在01x =附近的根. 要求给出程序和运行结果.5. 证明迭代格式()21233k k k k x x a x x a++=+， 00,0a x >>.6. 编写用全主元Gauss 消去法解线性方程组的程序，并求解12345123451234512345123450.024*******4233433241634418x x x x x x x x x x x x x x x x x x x x x x x x x -+-+=⎧⎪-++++=⎪⎪+++-=⎨⎪-++++=⎪⎪+-++=⎩ 7. 用追赶法解线性方程组12345210001121000012100001210000120x x x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦ 要求给出程序和运行结果. 8．给定线性方程组12122132x x x x +=-⎧⎨+=⎩ 问用雅可比迭代法和Gauss-Seidel 迭代法求解是否收敛？ 9. 设有线性方程组123521121422023103x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦ （1）考察用Jacobi 迭代法和Gauss-Seidel 迭代法解此方程组的收敛性； （2）分别用Jacobi 迭代法和Gauss-Seidel 迭代法解此方程组，要求当(1)()410k k x x +--<时迭代终止. 给出求解程序和迭代次数及结果.10．编写幂法程序求矩阵422251216A ⎛⎫⎪= ⎪ ⎪⎝⎭按模最大的特征值1λ和对应的特征向量1x 。

西安电子科技大学计算机组织与体系结构课程实验报告实验名称计算机组织与体系结构计算机学院班Array姓名陈宁学号同作者罗超实验日期 2017 年 9 月 24 日实验地点E-II-311 实验批次第二批一．实验目的1．深入理解基本模型计算机的功能、组成知识；2．深入学习计算机各类典型指令的执行流程；3．学习微程序控制器的设计过程和相关技术，掌握LPM_ROM的配置方法。

4．在掌握部件单元电路实验的基础上，进一步将单元电路组成系统，构造一台基本模型计算机。

5．定义五条机器指令，并编写相应的微程序，上机调试，掌握计算机整机概念。

掌握微程序的设计方法，学会编写二进制微指令代码表。

6．通过熟悉较完整的计算机的设计，全面了解并掌握微程序控制方式计算机的设计方法。

二．实验原理1．在部件实验过程中，各部件单元的控制信号是人为模拟产生的，而本实验将能在微过程控制下自动产生各部件单元控制信号，实现特定的功能。

实验中，计算机数据通路的控制将由微过程控制器来完成，CPU 从内存中取出一条机器指令到指令执行结束的一个指令周期，全部由微指令组成的序列来完成，即一条机器指令对应一个微程序。

2．指令格式 （1）指令格式采用寄存器直接寻址方式，其格式如下：其中，OP -CODE 为操作码，rs 为源寄存器，rd 为目的寄存器，并规定：1,存储器读操作（KRD ）：下载实验程序后按总清除按键（CLR ）后，控制台SWA 、SWB 为“0 0”时，可对RAM 连续手动读入操作。

2,存储器写操作（KWE ）：下载实验程序后按总清除按键（CLR ）后，控制台SWA 、SWB 为“0 1”时，可对RAM 连续手动写操作。

3、启动程序（RP ）：下载实验程序后按总清除按键（CLR ）后，控制台SWA 、SWB 为“1 1”时，即可转入到微地址“01”号“取指令”微指令，启动程序运行。

表6-2 A 、B 、C 各字段功能说明：图6-1 数据通路框图24位微代码中各信号的功能(1) uA5—uA0：微程序控制器的微地址输出信号，是下一条要执行的微指令的微地址。

西安电子科技大学课程论文数学软件与实验最速下降法求最优解姓名：方正阳学号：07117020班级：07117107112016、最速下降法求最优解1 2 n，然后MATLAB 结课大作业摘要：最速下降法,又称为梯度法，是一种重要的无约束最优化方法。

它是 1847年由著名数学家 Cauchy 给出的，其他解析方法或是它的变形，或是受它 启发而得到，因此它是最优化方法的基础。

该法将 n 维问题转化为一系列 不断迭代过程中沿负梯度方向用一维搜索方法寻优的问题，本次程序设计 利用最速下降法算法，反复迭代，最终收敛于局部最优点，即为解出的二 元函数的无约束非线性规划问题 minf(x,y)。

引言：最优化理论作为运筹学中的一个重要理论方法，在工业生产,金融经济活 动,工商管理,国防建设,计算机应用中，都有着重要的应用。

最优化理论 通过给出生产活动中的各类实际问题的数学模型，通过最优化方法，寻求 该问题的最优解或满意解。

最速下降算法是最优化理论中常见的一个重要 算法，理论证明：最速下降算法在一定条件下是收敛的，它能够有效地求 解一部分无约束最优化问题。

一、 实验目的熟悉最速下降法算法思想和步骤，用 MATLAB 语言编程最速下降法 求最优值。

二、 实验要求在最优化计算方法中，要求解 y = f (x 1, x 2 , , x n ) 的局部最小值，可以采用如下的方法进行迭代计算：先给出初始点 x 0 = (x 0 , x 0 , , x 0)根据其梯度方向 ∇f (x 0)，计算一元函数 y (λ1 ) = min f (x λ≥0 0-λ⋅∇f (x 0 )) ，并1 0 0得到x = x -λ1 ⋅∇f (x ) 。

如此反复迭代，最终收敛于局部最优点。

实现 该算法，求 的最优值，a,b,c,d 自定（非 0）三、 实验假设考虑到参数的随机性、代表性，验证程序的正确性、典型性，在此 我们从两个角度出发，一是在 abcd 值确定的情况下改变初始搜索位置 x0,看函数最优解是否相同；二是初始搜索位置 x0 相同，abcd 值不同的 情况下，看函数最优解是否相同。

题目要求1. 编制条件如图所示，用差分法求区域内的电压值。

0v10v0v0v0v0v解：由题意，我们将不规则部分补全，并进行等效处理，如下图结果所示，图示给出的是对整体补全后做3*3 的有限差分结果，当然网格化点数可以根据需要做改变，这里只是体现方法，故只取了 9 个点。

876 a11o a12=10v o-inf54 a21o a22=0v o a23=0v32 a31o a32o a331-1 0 1 2 3 4 5 6 7 8 9根据拉普拉斯 5 点差分原理，可知得到关于电压变量 a(i, j 1, 2, 3) 的i , j方程如下：4a 1,1 a 2,110;a 1,1 4a 2,1 a 3,1 0; a 2,1 4a 3,1 a 3,2 0; a 3,1 4a 3,2 a 3,3 0; a 3,2 4a 3,30.4 1 0 0 0 10 1 4 1 00 0 写成矩阵的形式： Ax b ; 其中， A 0 1 4 10 ， b 0 。

0 0 1 41 0编写程序可以求得0 01 4a , a , a , a , a , 2.6790.7180.1920.0513 0.0132. 在区域一边有个源，边界为 PML 边界，用 FDTD 法求所研究区域的场分布。

建模说明：二维 TE 波在空间传播，采用 PML 边界吸收，点辐射源验证。

FDTD 基本差分方程Yee 采用矩形网格进行空间离散，将每个节点进行编号，节点的编号和其空 间坐标位置按照下面的方式对应起来()(),,,,i j k i x j y k z =∆∆∆ (2-1) 而该点的任意函数()x,y,z,F t 在时刻n t ∆的值可以表示为：()(),,,,,n F i j k F i x j y k z n t =∆∆∆∆ (2-2)式中x ∆、y ∆、z ∆分别为沿,,x y z 方向上离散的空间步长，t ∆是时间步长。