《网格序列法实现等值线生成》实验报告讲述.docVIP

《科学计算可视化》实验报告 实验题目： 采用网格序列法实现等值线生成 姓名 学号 专业 班级 计算机科学与技术学院 2015 年 10 月 29 日 实验内容 第一问：首先生成4*4的网格，在网格顶点进行数组赋值， 利用网格序列法实现等值线生成，等值为4.5，结果用红色标记。 第二问：三角剖分法生成上述等值线 主要技术及实现 （1）利用网格序列法生成等值线，算法如下： 逐行处理网格中的每一单元，对每一单元执行下述操作； 求出单元顶点的Fmin 和Fmax， Fmin=min(F00，F01，F10，F11) Fmax=max(F00，F01，F10，F11) 对于在Fmin和Fmax之间的每一等线值Ft，执行下列操作求出等值线段； 将单元每一顶点与Ft比较，生成记录“+”结点和“-”结点的两个集合{F+}和{F-} 如|F+|=1或|F-|=1，求出该点附近的交点P,Q，连接PQ，成为等值线段。 如|F+|=2：如果两个“+”点相邻，求出两个相对的交点，生成等值线PQ；如果两个“+”点相对，求出四个交点P,Q,R,S,求出中点值Fmid，如果F00与Fmid状态相同，则连接PS，QR，否则，连接PQ，RS。 利用三角剖分法生成上述等值线 1.逐行处理网格中的每一个单元，对每个单元cell执行下述操作。 2. 求出中心点函数Fmid的值，用平均值方法。 3.计算cell四个顶点的Fmax和Fmin. 4.For每个在Fmin和Fmax之间的Ft Do ① If F00=FtF01 Then (AXt,AYt)=lint(F00,F01) /lint为线性插值函数 等值线从(AXt,AYt)开始； If FmidFt Then (BXt,BYt)=lint(F01,Fmid) 连接(BXt,BYt)到等值线上； If F11=Ft Then (BXt,BYt)=lint(F01,F11) 连接(BXt,BYt)到等值线上； Else (BXt,BYt)=lint(Fmid,F11) If F01=Ft Then (BXt,BYt)=lint(F10,F11) 连接(BXt,BYt)到等值线上； Else (BXt,BYt)=lint(Fmid,F10) 连接(BXt,BYt)到等值线上； (BXt,BYt)=lint(F00,F10) 连接(BXt,BYt)到等值线上； END If Else (BXt,BYt)=lint(F00,Fmid) 连接(BXt,BYt)到等值线上； ② Else If F01=FtF11 ③ Else If F11=FtF10 ④ Else If F10=FtF00 Else If 实验结果 在4*4的网格中生成了红色（等线值为4.5）的等值线，网格法生成 三角剖析法生成 实验分析和总结 实验过程中，出现了许多问题，有粗心的原因，比如少写了取余符号，导致结果不对，也有知识欠缺导致的没有思路，经过认真地学习，解决了以上问题。 实验过程中，耐心和细心是很重要的。在算法方面，要注意等值点坐标的公式，很容易出错，还要考虑情况全面一些，以及考虑马鞍点的思路，这是我在实验中的体会。 五、核心源代码 首先将各个点以数组的形式存储 value[0][0] = 0; value[0][1] = 1;value[0][2] = 3;value[0][3] = 2;value[0][4] = 1; value[1][0] = 1;value[1][1] = 3;value[1][2] = 7;value[1][3] = 7;value[1][4] = 2; value[2][0] = 1;value[2][1] = 6;value[2][2] = 9;value[2][3] = 8;value[2][4] = 3; value[3][0] = 3;value[3][1] = 6;value[3][2] = 7;value[3][3] = 6;value[3][4] = 4; value[4][0] = 2;value[4][1] = 3;value[4][2] = 3;value[4][3] = 2;value[4][4] = 3; 网格法 for (i = 0; i4; i++) { for (j = 0; j4; j++) { vertex[0]