当前位置：文档之家 › FISH库函数

FISH库函数

  • 格式：pdf
  • 大小：172.04 KB
  • 文档页数：14

下载文档原格式

  / 14

合集下载

语言(FISH)解释

语言(FISH)解释
msg=' tunnel "height" from base to center (height) '
height = in_def
;交互式输入模型底部到隧道中心的的高度，默认为0.5
;
default = 2.0
msg=' extent of inner primitive in x-direction (dx: extent = x0-dx to x0+dx) '
;计算放射性网格点—调整nx,nu 单元数
x5=(x0*nu+x7*nx)/(nx+nu)
x5c=(nu*x0+nx*x7)/(nx+nu)
x5b=(nu*x0+2*nx*x7)/(2*nx+nu)
x6=(x0*nu+x2*nx)/(nx+nu)
;Compute z coordinates
end_if
;隧道底部的单元数量，这里用if 语句控制nz1 的输入
;
if xr-x0<1e-3 then
nx=0
else
default = 5
msg=' number of zones on far right region (nxr) '
nxr = in_def
end_if
;此if 语句主要是控制隧道上部右侧的单元数量，和网格的合理性
ny2 = in_def
;用来修正y 轴线的单元数量
;
default = 1.2
msg=' radial zoning ratio for ny2 zones (raty) '

Fish总结

Fish总结

Fish是嵌入在FLAC3D中一种编程语言，它允许使用者定义新的变量和函数。

函数被用来扩展FLAC3D的功能和增加用户自定义功能。

例如，输出或者打印新的变量，实行特定的网格生成器，数值测试中应用伺服系统，定义特定的性质以及参数研究的自动化。

伺服系统(servomechanism)又称随动系统，是用来精确地跟随或复现某个过程的反馈控制系统。

伺服系统使物体的位置、方位、状态等输出被控量能够跟随输入目标(或给定值)的任意变化的自动控制系统。

FISH就是简单的嵌入在FLAC3D的数据文件中—DEFINE后面的命令行按函数处理，函数遇到END结束。

函数可以调用其它函数，其它函数也可以调用其它函数。

函数定义的顺序没有影响，只要它们在被使用之前被完全定义就行。

因为函数的编辑形式储存在FLAC3D的存储空间，SAVE命令会保存函数和相关变量的当前值。

2.2 FISH语言规则，变量及函数2.2.1命令行FISH程序包含在普通的FLAC3D数据文件中或者直接从键盘键入。

DIFINE后边的命令行按FISH函数处理；函数以遇到END为结尾。

FISH的有效命令行必须是下述格式中的一种。

1.命令用IF,LOOP等开始。

（见2.3）2.如果一行包含一个或更多用户定义FISH函数的名字，以空格号隔开。

例如fun_1 fun_2 fun_33.命令行包含指派声明（例如等号右边的数学式被运算且其值赋予给等号左边的函数名称或是变量）。

4.命令行由FLAC3D（在FIAH中调用FLAC3D命令）命令组成，通过COMMAND—ENDCOMMAND 分隔FISH指令。

命令行包含在FISH指令的一部分里面。

5.命令行以分号开始或者结束。

函数的变量，函数名称或者声明必须完整拼写，在FLAC3D命令中不能缩写。

不允许有连续的命令行；可以通过引入新变量分解复杂的表示式。

FISH在默认状态下为“对情况不敏感”—即大写与小写字母之间没有区别，所有的名字都转换为小写字母（注意可以通过（SET case sensitivity on）命令改变）。

FISH参考

FISH参考

in （var） 若 var 是字符串则输出， 若不是则等待键盘输入。 （返回值依赖于字符的类型。 FISH 首先作为整型输入，然后作为浮点型——若是输入的是可编码为整型或 是浮点型的单值，那么返回值就是整型或是浮点型。这个数字独立占用一行。 但是，若其后有空格、逗号或是括号，那么在改行中其他字符将被忽略。若用 户输入的字符不能编译为单值数字， 那么返回值将会是一系列的字符串。 用户 可以通过函数 type（）确定返回的是什么。 ） 2－14 out（s）将 s 信息显示在屏幕上（若记录文件是打开的，则记录到记录文件里）变量 s 的类型必须是字符型。若无错误函数的返回值为零，若有错误则返回 1（如 s 不是字符型的） string（var）将 var 变为字符型 这些函数的用途之一就是控制输入输出。例 2.7 表明用这些函数给杨氏模量和泊松比 输入参数。 例 2.7 交互输入的控制
变量或函数名必须以非数字开头并且不得包含下列符号（运算符） . , */+-^=<>#()[ ]@;‘“ 用户自定义名可以任意长度，但是他们以为行的有限而截取。一般情况下，名字只要不 与 FISH 表达式（见 2.3 节）和预定义的变量和函数名（见 2.5 节）相同，就可任意选取。 在 FLAC 中还有一些其他的单词需要避免。 表 2.1 包含了所有的可能产生冲突的名字的清单、 但是，潜在的冲突取决与所选名字是如何使用的。例如，单词 gravity 可以用做 FISH 变量， 在 FISH 函数中提供出来。 冲突只有在用 SET 命令设定其值是才出现， 因为 gravity 对于 SET 命令是有争议的。类似的，如果它的名字与 PRINT 命令的参数相同，也会产生冲突。若想 毫无顾虑的使用名字， 请参见表 2.1 或者名字的缩写 （因为 FLAC 允许关键字和命令的缩写） 作为一种可选策略， FLAC 命令 SET safe on 可以用来强制在命令行中的 FISH 变量的得 到认可。@符号在任意 FISH 变量之前可以确认物体并强制命令行避免其他任何翻译。例如 假设 FISH 函数 initial 已建立，如例 2.1 所示。在命令行运行该函数将产生错误，因为命令 INITIAL 的优先权超过了 FISH 的 initial 的优先权。 但是在变量前使用 FISH 确认符号@给出 确定，纠正翻译。

fish初学者指南

fish初学者指南

FISH 初学者指南中国矿中国矿大大 力学力学小杨小杨小杨 QQ QQ QQ：：277502960 yhb20002000@ ，欢迎交流欢迎交流1.1.概述概述概述FISH是一门内植于UDEC的程序设计语言，它能让使用者定义新的变量和功能。

这些功能可以拓展UDEC的用途或者增加用户定义的功能。

例如，可以绘出或者打印新的变量，执行特殊的模拟，数值测试的随动控制，指定不寻常的属性分布和参数研究的自动化。

对使用Itasca软件现有程序结构做工程困难或者不能实现的用户来说，FISH语言发展响应了那些用户的需求. 与其在增加许多新的特殊功能到标准代码中，不如提供一种内置语言使用户能编写他们自己的模块。

一些实用FISH模块已经被写出了，UDEC项目提供一个有关此类的库（见FISH卷的第三节）。

对某些人来说，可能没有经验设计去编写简单的FISH模块或者修改现有的一些简单的模块。

4-2节为非程序员提供一个介绍性的指南。

然而，FISH语言和其他任何程序设计语言一样可以变得非常复杂。

要查阅更多细节，参考FISH卷的第二节。

类似所有规划任务，FISH功能可以被构造按增加的形式，在向前使用更复杂的代码前可以检查运行可行性。

相比较大多数编译器来说，FISH实行错误检测比较少，所以在实际应用它们之前，应该用些小的数据设置测试所有的功能。

FISH程序非常简洁的内植于正常UDEC数据文件，由关键词“DEFINE ”引入一个FISH函数，当遇到关键词“END ”时，该函数结束。

函数可以调用其它函数，被调用的函数又能继续调用其它函数，如此循环。

只要是在使用前都被定义了的，函数的调用次序可以任意（例如：由一个UDEC指令调用）。

FISH函数的编译形式储存在UDEC内存空间，指令“SAVE ”保存函数和各个相关变量当前值。

在第二篇FISH卷提供FISH语言规定的完全定义和内部函数。

它包含了有关语法，数据类型算法变量和函数的各种规范。

在第二篇FISH卷描述了所有FISH语言名，并且在“指令和指令和FISH FISH FISH参考概要参考概要”中提供这些名称的概要。

Fish语言学习笔记

Fish语言学习笔记

Fish语言学习笔记Fish语言学习笔记FISH 语言中的变量名以及函数名在整个程序中都有效，不管是在FISH 代码中还是在FLAC3D 命令行中。

如果变量名没有被赋值，则默认为0。

我们可以删除或重新定义FISH 函数，方法是用同名的新代码取而代之，如果在DEFINE 后直接跟END 行，那么也就删除了原定义的函数。

当函数被删除了，但原有变量依然存在，因为变量是全局性的，如同在其他位置一样。

控制语句EXIT，会使程序无条件跳到当前函数的结束处。

与节点参数有关命令：1.gp_near(x,y,z)：获得距坐标(x,y,z)最近的节点的指针2.ngp：节点总数。

3.gp_nearall(x, y, z) ：获得距坐标(x,y,z)最近的节点的指针，包括空单元节点在内。

4.gp_id(p_gp)：获得指针为p_gp的节点的id5.gp_head：第一个节点的指针6.gp_next(p_gp)：获得指针为p_gp的节点的下一节点的指针。

7.gp_xpos(p_gp)：获得指针为p-gp的节点的x坐标。

8.gp_xdisp(p_gp)：获得指针为p-gp的节点的x方向位移。

9.gp_xvel(p_gp)：获得指针为p-gp的节点的x方向变形速率。

10.find_gp(id)：获得id为id的节点的指针。

11.gp_extra(p_gp,ind)：获得指针为p_gp的节点的索引（编号）为ind的额外参数。

12.gp_yfunbal(p_gp)：获得指针为p-gp的节点的y方向节点不平衡力。

与单元参数有关命令：1.z_head：第一个单元的指针。

/doc/fa4590171.html,one：单元总数。

3.find_zone(id)：获得编号为id的单元的指针。

4.z_near(x, y, z)：得距坐标(x,y,z)最近的单元的指针。

5.z_nearall(x, y, z) ：获得距坐标(x,y,z)最近的单元的指针，包括空单元在内。

fish 语法

fish 语法

fish 语法
Fish是一种流行的命令行Shell，它具有许多有用的功能和语法。

以下是Fish 语法的一些主要特征：
1. 命令自动补全：Fish具有强大的命令自动补全功能，可以根据您键入的内容自动推断您要输入的命令或选项。

2. 命令历史记录：Fish会自动记录您输入的命令历史记录，并且可以使用方向键或Ctrl+R来搜索和重复以前的命令。

3. 管道和重定向：Fish支持管道和重定向操作，可以将命令的输出发送到文件或其他命令中进行处理。

4. 变量：Fish支持变量，可以将值存储在变量中，并在需要时使用它们。

5. 函数：Fish支持函数，可以将一系列命令组合成一个函数，并在需要时调用它们。

6. 条件语句：Fish支持条件语句，可以根据条件执行不同的命令。

7. 循环语句：Fish支持循环语句，可以重复执行一系列命令。

8. 命令别名：Fish支持命令别名，可以将长命令或常用命令的缩写定义为别名，以便更快地输入它们。

9. 脚本文件：Fish支持脚本文件，可以将一系列命令保存在文件中，并在需要时执行它们。

总之，Fish具有许多有用的功能和语法，使得在命令行中工作更加高效和便捷。

FLAC-3D小技巧

1、如果是外观上看模型，要切换视角的话除了上述的快捷键，命令是：plot set rotation 20 0 60 rotation后面的3个数字是绕x、y、z轴转动的角度。

使用快捷键时，每次操作转动的角度，或者move的大小，可以通过plot set rotincrement和plot set moveincrement来设置。

2、从surfer把数据倒入FLAC3D!不知道各位在做象边坡地形表面这种不规则面的模型时有没有什么好招。

我最近发现了一种方法。

在FLAC3D的fish函数库中有一个叫TOPO的函数，可以生成自由表面模型的网格，关键是它要用到的那个数据表做起来非常麻烦，具体就是在地形图上平行y轴画一系列的等距的直线，然后读出每条与每条直线相交的等高线的x值和高程，然后生成一张表，格式是table y x z。

用surfer可以很方便的生成这张表，在DATA菜单下的convert...可以将GRD文件生成以上要求的那种格式。

3、按找手册中的解释，print gp position 是用来输出gp的位置坐标的；print zone gp是用来输出与zone相关的gp的id值的。

这样我就有一个上面的疑问，两个命令中的gp都有id值，是不是不同命令中的gp的id值相同的话，他们的坐标就可以看作是一样的啊？答：如果要画等值线的话，先输出gp的坐标来，将其另存为一个文件，再输出zone的应力值，将与zone对应的gp的id与前面输出的gp的id值对应起来，找到gp的坐标。

用excel打开上面的文件进行编辑，最后用surfer画等值线。

这样的办法太苯，不知到有没有更好的办法啊？4、plot stensor plane local on out off能够显示的主应力矢量图，但仍有两个问题1）显示的矢量图中应力好像没有变化，如何通过设置符号的大小表示应力的大小变化？2）应力符号好像总在图形没有发生变形前的单元中心，如何使应力符号在发生变形了的单元中心？5、个人经验，无论有限元还是flac，如果支护措施仅仅只考虑其刚度的贡献，岩体的变形和塑性区变化往往都不大，可能需要适当提高岩体的参数，才能反映支护的作用6、＜锚杆支护的数值模拟＞．中国矿业大学岩石力学与土木工程学报前段时间有篇文章专门关于flac锚杆单元应用模型修改的文章，7、如何模拟预应力自由式锚索？首先我想说说我在模拟预应力自由式锚索的做法,然后在说出现的问题,希望大家帮忙解决一下!我采用先整根生成的方法来生成锚索,然后给自由段和锚固段分别赋不同的参数,自由段的灌浆参数统统赋为零,锚固段的按试验值确定.在模拟垫墩时,把自由段外端的交接删除,然后重新生成刚性铰接(之所以这样做是因为模型介质较软,不这样做就收敛不了,介质内发生了流变.而全长粘结式锚索不需这样处理也和试验值拟合的很好),使得该结点不会和周围介质发生相对位移;然后在垫墩作用的面积上再施加一个均布力,大小等于预应力除以垫墩面积.问题是位移值拟合锝还可以,但是锚索的预应力增加值比实测的大的多,有没有其他的方法来模拟自由式锚索???8、Fish函数增加了以下几种新的特性：1.增加了fish变量来获取结点、单元和界面变量2.fish提供了获取结构单元变量的途径3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量5.提供了fish绘图子程序函数够生用户定义的图形内容6.fish函数已经增加了从文件读、写数据的能力FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说注意点：1.fish函数可以嵌套使用;2.以save命令保存模型时，fish函数和变量也同时保存；3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；4.变量或函数名不能以一个数字开头也不能是下列字符：. , * / ^ = > < # ( ) [ ] @ ; “ '5.如果用命令set safe on 指定了编译安全模式，则用户调用fish函数时，函数名前必须加@；6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的类型决定；7.fish函数的调用方法：.可以出项在其他fish函数的单独行中；.可以出现在其他fish函数的表达式中；.出现在flac3d的命令行中；.作为命令set,print,hist的参数。

fish语言 规则

fish语言规则Fish语言规则Fish语言是一种简化的编程语言，其规则包含了以下几个方面：1. 命令规则Fish语言中的命令以英文字母开头，并以换行符或分号作为命令的结束符。

命令可以包含多个参数，参数之间需要使用空格进行分隔。

2. 变量规则Fish语言中的变量以$符号开头，并由字母、数字和下划线组成。

变量可以用于存储和操作数据。

使用变量时，需要在$符号后面加上变量名。

3. 注释规则Fish语言中的注释以#符号开头，可以用于对代码进行解释和说明。

注释内容将被解释器忽略，不会执行。

4. 条件语句规则Fish语言中的条件语句使用if关键字进行定义。

if后面跟上条件表达式，如果条件成立，则执行if后面的代码块。

条件表达式可以包括比较运算符、逻辑运算符和括号。

5. 循环语句规则Fish语言中的循环语句包括for循环和while循环。

for循环用于遍历一个集合或者一段范围内的值，while循环用于在条件成立时重复执行一段代码。

6. 函数规则Fish语言中的函数以function关键字进行定义。

函数可以接受参数，并返回一个值。

函数的定义以function关键字开始，后面跟上函数名和参数列表，然后是函数体。

7. 输入输出规则Fish语言中的输入输出使用echo关键字进行操作。

echo后面跟上要输出的内容，可以是字符串、变量或者表达式。

8. 异常处理规则Fish语言中的异常处理使用try和catch关键字进行定义。

try后面跟上要执行的代码块，catch后面跟上要处理的异常类型和对应的代码块。

9. 数组规则Fish语言中的数组以@符号开头，并由方括号包围。

数组可以存储多个值，并通过索引进行访问和操作。

10. 字符串处理规则Fish语言中的字符串可以使用单引号或双引号进行定义。

字符串可以进行拼接、截取和替换等操作。

11. 模块引入规则Fish语言中的模块使用source关键字进行引入。

source后面跟上模块的路径。

grads fish函数

grads fish函数如何使用GRADs软件包中的fish函数进行统计分析？GRADs（Grid Analysis and Display System）是一个用于进行大气和海洋科学数据分析和可视化的软件包。

它是一个免费开源的软件，拥有丰富的功能和灵活的用户界面。

GRADs中的fish函数是一个用于执行统计分析的重要工具。

本文将介绍如何使用GRADs软件包中的fish函数进行统计分析，并为读者提供一步一步的操作指南。

1. 准备数据在使用fish函数之前，首先需要准备好要进行统计分析的数据。

GRADs 支持多种数据格式，包括常见的二进制和ASCII格式。

如果数据量较大，建议使用GRIB（GRIdded Binary）格式，以便更高效地处理大量数据。

在这里，假设我们有一组气象观测数据，包括温度、湿度和风速等参数。

这些数据可以是同一站点的连续观测值，也可以是多个站点的同一时刻观测值。

确保数据已经按照要求的格式组织，并且可以在GRADs中正确加载和显示。

2. 加载数据在GRADs中，使用“sdfopen”命令加载数据文件。

例如，如果要加载名为“data.sdf”的数据文件，可以执行以下命令：sdfopen data.sdfGRADs将自动识别并加载数据文件中的数据集。

3. 设置分析区域和时间段在进行统计分析之前，需要定义感兴趣的分析区域和时间段。

GRADs提供了一系列命令用于设置分析区域和时间段，例如“gxout”和“gxadd”命令可以设置经纬度范围，而“d”命令可以设置时间范围。

gxout statgxadd -R<lon1>/<lon2>/<lat1>/<lat2>d <start_time>,<end_time>通过设置这些参数，GRADs将仅在指定的区域和时间范围内执行后续的统计分析。

4. 执行统计分析现在可以使用fish函数执行统计分析了。

fish学数学c语言

fish学数学c语言鱼学数学是一门综合学科，涵盖了多个分支，包括数学分析、线性代数、概率统计、微分方程等等。

在数学中，我们通常会使用编程语言来解决一些复杂的数学问题。

C语言是一种功能强大的高级编程语言，它与数学之间有着密切的联系。

在本文中，我们将探讨如何使用C语言解决一些常见的数学问题。

数学计算是C语言中常用的一个功能。

在C语言中，我们可以使用一些基本的运算符（如加减乘除）来进行数学运算，还可以使用一些数学函数库来实现更复杂的计算。

例如，我们可以使用C语言中的<math.h>库来调用sin、cos、tan等三角函数。

下面是一个例子，展示了如何使用C语言来计算三角函数的值。

```c#include <stdio.h>#include <math.h>int main(){double angle = 3.14159 / 4; // 弧度表示角度double sin_value = sin(angle);double cos_value = cos(angle);double tan_value = tan(angle);printf("sin(%.2f) = %.2f\n", angle, sin_value);printf("cos(%.2f) = %.2f\n", angle, cos_value);printf("tan(%.2f) = %.2f\n", angle, tan_value);return 0;}```以上代码将输出角度为45度（即π/4）的正弦、余弦和正切的值。

通过使用C语言的数学库函数，我们可以在数学问题中方便地应用三角函数，这对于解决许多实际问题非常有帮助。

在数学中，矩阵是一个重要的概念。

线性代数是研究矩阵和线性方程组的学科，也是许多数学问题的基础。

在C语言中，我们可以使用二维数组来表示和操作矩阵。

  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
value . . .设置成数据文件格式。var1, var2,是变量名，须设置初始值。
3. 在使用DEF命令之后，通过输入地址的函数名来调用FISH函数。若使用命令
fishcall，不得不给出函数名。
FLAC3D中的FISH函数以不同的方式相互作用。你可以看2.4节，其中具体讲述了 FISH 和FLAC3D.之间的连接。
做FLAC3D图表中功率频谱的傅立叶变换
文件“FFT.FIS”执行图表数据的傅立叶变换，使功率频谱由另一个图表输出。定 义输入表的ID 号为fft in，输出为fft out。输出表使用相同的ID，fft out，在调用文件 “FFT.FIS.”之前并不存在。fft in和fft out须在分离的FISH函数中说明，在 访问文件“FFT.FIS”之前执行。 对功率频谱有几点说明，功率频谱被定义为 N / 2 实数，
1*1*Z柱形区域（z轴方向为垂直方向）。在此区域定义8个坐标。假定底部由n_zon row * n_zon_col矩形单元构成，单元体的维数统一为dx和dy。我们还可以 假设底部所以网格点一起向z_base升高。现在问题归纳为如何定义4个基准点。
有一些方法可以使用，对地形数据进行格式化。等值线程序诸如SURFER， 在一个不规则扰动数据表中插入一个规则网格，插入的数据通过图表输入到 FLAC3D。在这个实例中，我们使用了一种不同的方法。从等高线地图开始（如图 1所示）我们采用y轴几段不同区域，每部分的距离为dy。每部分创建一个列表， 每个节点坐标（ x,z）与等高线相交。编写数据转换器可以相对简单的去快速
图1 输入加速度
图2 位移响应频谱 图3 Pseudo 响应频谱
图 频谱
图 4 Pseudo 加速度响应
绘制地形网格 这节主要描述FISH函数“TOPO.FIS,”，讲解了应用FISH函数来创建给定地形
模型表面单元。பைடு நூலகம்
绘制网格主要依靠一系列GENERATE zone brick命令。每个命令建立一个
参考文献 Craig, Jr., R. R. Structural Dynamics—An Introduction to Computer Methods. New York: John Wiley and Sons, 1981.
数据文件“SPEC.DAT”
def cr tab i=0 p2 = 2.*pi loop while i <= num point xx=end time*float(i)/float(num point) i = i+1 yy = sin(xx*p2/per1) table(1,xx) = yy end loop end set num point 250 end time 3.0 set per1 1.0 cr tab ca spec.fis set pmin = 0.5
3 FISH库函数 这节主要讲FISH库函数，它可以用来进行FLAC3D分析。函数用来建模和计
算解析解，包括网状结构的生成、绘图、设置材料的属性和控制解的范围。
FISH库函数的文件在“\FISH\Library”目录下。
执行FISH库函数可以分为以下三步。 1. 用命令call filename.fis访问FLAC3D数据文件。 2. FISH变量，如果在函数文件的顶端给出，用命令set var1 = value var2 =
（i）增大输入表中数组元素的值。
(ii) 权为2。（数组维数不必计算）
fft的运算法则要求输入数据带有一个恒时间步长，计算时间步长，从图表
中插入数据，存储到数组中。
下面的例子用来检验fft FISH函数，输入正炫波（频率1Hz，振幅为1），余
炫波（频率5Hz，振幅为2），正炫波（频率10Hz，振幅为3）。FISH 函数cr_tab
找出加速度响应频谱 用文件“SPEC.FIS”可得到位移响应频谱，虚拟速度响应频谱和虚拟加速度
响应频谱。Acc_in规定输出图表的ID号，这三个输出表分别由sd_out, sv_out 和 sa_out识别。当前存在的任意图表将被删除或复写。
命令dmp规定响应分析持续的衰减。近似计算衰减响应——dmp越高，响应 越不准确。pmin 和 pmax命令计算周期的范围， n_point规定输出表格的数 值点。
图1 三次输入波的总和
图2 功率频谱；功率与频率（Hz）
输出的浮点数带有特定的精度 FISH文件 “NUMBER.FIS”用来输出带有特定精度的浮点数，精确度小于16位数 字。FISH文件正常输出浮点值的精度小于5位。 提取数字符合精度要求的运算法则：
（a） 把值转化成浮点型（打字类型） （b） 测定k，k = [log(number)-1] if log(number) < 0.0
数据文件“TOPOTAB.DAT”
; Tables of interpolated data from contour map ; (Figure 6.10, Statistics and Data Analysis in Geology, John C. Davis, ; John Wiley & Sons, 1973) ; Table number corresponds to constant y = 40 ; with each x,z pair at point where the line intersects a contour level ; ; Table 1 @ y = 40 table 1 0 950 table 1 20 925 table 1 64 900 table 1 130 875 table 1 160 875 table 1 209 900 table 1 315 900 table 1 353 875 table 1 400 870 ; Table 2 @ y = 80 table 2 0 910 table 2 20 925 table 2 50 925 table 2 110 900 table 2 180 900 table 2 215 925 table 2 248 950 table 2 276 950 table 2 306 925 table 2 337 900 table 2 365 875 table 2 376 875 table 2 400 900 ; Table 3 @ y = 120 table 3 0 890 table 3 122 875 table 3 158 875 table 3 185 900 table 3 245 925 table 3 287 925 table 3 308 900 table 3 365 875 table 3 392 900 table 3 400 910 ; Table 4 @ y = 160 table 4 0 910 table 4 8 900
= [log(number) ] 否则 [ ..] 表示括号内的是整数值。 （c）提取的每个值的精度小于规定精度＋1，储存在一个数组 中。 数位小于规定精度＋1，最后一位数若大于等于5就舍入。因此结果位使用者规定 精度，以字符串类型储存的浮点值。位数以k阶精度输出。 注：如果使用者没有规定精度限制，FISH函数自定义精度为7位。 数据文件“NUMTEST.DAT”阐述了如何打印10位精度的数字。在这个例子中，任 意20位的数字以10位精度输出。 数据文件“NUMTEST.DAT”
执行此程序须大量时间，如果Ni是输入值，Np是输出值，则计算结果相应增
加Np × Ni ×log(Ni )。公式表明结果稍加变动，周期接近0；然而，改进小周期的
精度却增加了计算时间。 Craig (1981)运算法则。一个简单的正炫波作为输入加速度被输入。从0.5—2分成 50个点（见“SPEC.DAT”）,图1显示了输入加速度——周期1.0的正炫波。图2到图 4为不同响应频谱，显示了周期1.0处的峰值点。
set pmax = 2. set dmp = 0. set acc in = 1 set sd out = 2 set sv out = 3 set sa out = 4 set n point = 50 spectra ; table 1 name ’input acceleration’ ; table 2 name ’displacement response spectrum’ ; table 3 name ’pseudo-velocity spectrum’ ; table 4 name ’pseudo-acceleration spectrum’ ; ret
数据文件“FFT.DAT”
def cr tab i=1 p2 = 2.*pi loop while i <= num point xx=end time*float(i)/float(num point) i = i+1 yy = sin(xx*p2/per1)+2.*cos(5.*xx*p2/per1)+3.*sin(10.*xx*p2/per1) table(1,xx) = yy end loop end set num point 1024 end time 12.0 set per1 1.0 cr tab def tab ind fft in = 1 fft out = 2 end tab ind ca fft.fis fftransform ret
( ) P0
=
1 N2
⋅
f0
2
(1 )
[ ] ( ) ( ) Pk
=
1 N2
⋅
fk 2 +
f2 N −k
PN
2
=
1 N2
⋅
⎜⎛ ⎜⎝
fN
2
⎟⎞ 2 ⎟⎠
(2) (3)
N——原始数据区值的一半；
P——功率频谱的输出量
f——原始数据傅立叶变换值
k——从0到N/2变化