陶老师讲解的Simple程序

合集下载

ADAM-6251快速入门手册 - 免费文档

ADAM-6251快速入门手册 - 免费文档

免费文档中心免费文档中心可免积分在线阅读和下载文档∙免费文档∙高等教育∙高中教育∙初中教育∙小学教育∙外语考试∙资格考试∙工作范文∙求职职场∙党团工作∙表格模板∙总结汇报∙经管营销∙工程科技您的位置:免费文档所有分类工程科技信息与通信 ADAM-6251快速入门手册ADAM-6251快速入门手册Adam-6251 快速入门手册 ADAMADAM-6251 快速入门手册第一章产品介绍 .................................................................. ............................................................ 2 1.1 Adam-6200 概述................................................................... ......................................................... 2 1.2 Adam-6251 概述................................................................... ......................................................... 4 1.3 规格说明 .................................................................. ......................................................................4 1.3.1 一般规格 .................................................................. ............................................................. 4 1.3.2环境 .................................................................. .....................................................................5 1.3.3 Adam-6200 系列私有特性 .................................................................. ................................ 5 1.3.4 数字量输入参数 .................................................................. ................................................. 5 第二章 Adam 的软件安装................................................................... ............................................. 6 2.1 初始检查 .................................................................. ....................................................................6 2.2 安装 Advantech Adam/Utility ............................................................. ..................... 6 2.2.1 软件下载路径 .................................................................. ....................................................... 6 2.2.2 软件支持的操作系统 .................................................................. ........................................... 6 2.2.3 安装 Advantech Adam/Utility ............................................................. ................. 6 第三章硬件连接及测试 .................................................................. ................................................ 9 3.1 硬件连接 .................................................................. ......................................................................9 3.1.1 电源连接................................................................... ............................................................ 9 3.1.2 硬件接线................................................................... ............................................................ 9 3.1.3 Adam-6251 数字量输入功能接线 .................................................................. ................... 11 3.2 Adam-6251 的 Modbus TCP 通信协议编程时的地址映射....................................................... 13 3.3 软件测试 .................................................................. .................................................................... 15 3.3.1 Adam 模块通用参数配置 .................................................................. .................................. 15 3.3.2 Adam-6251 模块功能测试................................................................... ................................ 21 (1)数字量输入功能概述 .................................................................. ...................................... 22 (2) DI 通道高级功能设置 .................................................................. ..................................... 24 3.3.3 Adam-6251 GCL 功能 .................................................................. ...................................... 27 第四章例程使用详解 .................................................................. .................................................... 27 4.1 Adam-6251 板卡支持例程列表 .................................................................. ................................ 28 4.2 常用例子使用说明(以 CSharp 例程为例) ................................................................ .............. 28 4.2.1 Adam62XXDIO(数字量输入输出状态瞬时读值) ...................................................... 28 4.2.2 DataStreaming(主动定时上传功能,监测数据流信息) ............................................. 29 4.2.3P2P_UdpEvent(点对点功能和事件触发)..................................................................... 30 第五章遇到问题,如何解决? ................................................................ .................................... 32 1页下一页文档免费下载:ADAM-6251快速入门手册(共32页,当前第1页)你可能喜欢∙快速入门手册∙使用手册∙技术手册∙客户手册∙用户手册∙评价手册∙实施手册∙计价手册∙F2MC-16 SOFTUNE V3 快速入门手册73页∙lingo8.0中文快速入门手册53页∙M218 快速入门手册_V1.2_2010060841页∙TKSA40快速入门手册2页∙通达OA工作流程案例及快速入门手册135页∙会议快速入门手册(浦口区行政服务中心)4页更多与“快速入门手册”相关的内容>>∙DH-DRS录音系统使用手册16页∙OA操作使用手册[1]10页∙电子环使用手册3页∙老板无忧监控管理软件-使用手册26页∙ASA301硬件使用手册10页∙PriciseProject使用手册219页更多与“使用手册”相关的内容>>∙YF17技术手册63页∙R434a技术手册10页∙弯头技术手册24页∙SPE技术手册82页∙项目技术管理手册9页∙A4900A技术手册2页更多与“技术手册”相关的内容>>∙2012年集团客户营销产品手册32页∙装饰管客户服务流程手册18页∙大客户流程手册175页∙客户沟通手册4页∙400终端客户操作手册41页∙金牌客户服务人员手册07072565页更多与“客户手册”相关的内容>>∙PM-300用户手册159页∙TRSWCMV6用户安装手册(标准版)46页∙STM32F103-EVAL+用户手册38页∙CmailPlus用户手册18页∙出纳机用户手册5页∙企业用户操作手册50页更多与“用户手册”相关的内容>>∙英语评价手册(1-6Vocabulaary2011-12)7A30页∙附1嵘昌集团评价手册42页∙扎旗音三中课堂教学评价手册3页∙绩效评价应用手册26页∙评价手册参考11页∙教育教学实践评价手册34页更多与“评价手册”相关的内容>>∙课题管理与实施手册9页∙san实施手册58页∙5S实施手册21页∙erp实施手册73页∙Vmware_vsphere 实施手册48页∙课题实施手册27页更多与“实施手册”相关的内容>>∙计价手册办理程序2页∙计价手册申请范本6页∙《计价手册》登记表3页∙浙江03计价培训手册39页∙《模具计价手册》简介1页∙建设工程计价手册年检认证表6页更多与“计价手册”相关的内容>>ADAM 6251快速入门手册的相关文档搜索∙win10快速入门手册∙起亚k2快速入门手册∙thinkphp快速入门手册∙起亚k5快速入门手册∙起亚k3快速入门手册∙拍立得快速入门手册∙tiny快速入门手册∙adam4000手册ADAM 6251快速入门手册相关文档∙ADAM-6052快速入门手册38 1 Adam-6052 快速入门手册第一章产品介绍 1.1 adamadam-6052概述Adam-6052 是一款基于以太网远程数据采集的智能型 I/O 模块。

scratch怎么编随机计算题

scratch怎么编随机计算题

文章标题:从零开始,学习使用Scratch编写随机计算题一、引言在日常学习和工作中,我们经常需要进行各种计算,而编写随机计算题是一种非常有趣和实用的技能。

Scratch作为一款适合初学者的编程工具,它提供了丰富的功能和图形化的编程界面,非常适合用来编写随机计算题。

本篇文章将带领大家从零开始学习使用Scratch编写随机计算题,让我们一起来深入探讨这个主题。

二、Scratch简介Scratch是麻省理工学院的媒体实验室开发的一款图形化编程工具,它可以让用户轻松地创建自己的互动媒体、动画和游戏等。

Scratch采用积木式编程,用户只需将不同的积木拼接起来,就能完成复杂的编程任务,非常适合初学者使用。

三、什么是随机计算题?在学习和教育中,随机计算题是一种常见的形式,它能够帮助学生巩固知识和提高计算能力。

随机计算题通常包括随机生成的数字和运算符,要求学生进行计算并给出答案,是一种很好的训练方式。

四、使用Scratch编写随机计算题1. 生成随机数在Scratch中,我们可以使用"随机数"积木块来生成随机数。

在舞台中新建一个角色,然后在该角色的脚本中使用"当绿旗被点击"积木块和"重复"积木块来生成随机数,并将结果保存在变量中。

2. 生成随机运算符类似地,我们可以使用"选择"积木块来生成随机的加减乘除运算符。

将不同的运算符放入"选择"积木块中,并通过随机数来选择其中一个运算符。

3. 显示计算题在Scratch中,我们可以使用角色的气泡对话框来显示生成的随机计算题。

将随机生成的数字和运算符组合起来,并显示在气泡对话框中,让学生进行计算。

4. 检查答案我们可以添加判断逻辑,让学生输入答案并进行比较,从而检查计算的准确性。

如果答案正确,则给出相应的奖励,否则提示错误并给出正确答案。

五、总结与展望通过本篇文章的学习,我们了解了如何使用Scratch来编写随机计算题。

SIMPLE算法及计算例子

SIMPLE算法及计算例子
The effect of shallow water depth on a big FPSO(Floating Production,Storage and Offioading System)is numerically simulated in case three.It has been shown that the viscous plays a very'small role and so it carl be neglected in this condition.The varying law of pressure produced by shallow water to FPSO is given.
第二个问题:压力的一阶导数是以源项的形式出现在动簧方程中。采用分离式求解各 变量的离散方程时,由于压力没有独立的方程,压力与速度的关系隐含于连续性方程中,如 果压力场是正确的,则据此压力场求得的速度场一定满足连续性方程。如何构造求解压力场 的方程,或者说在假定初始乐力分布_|亓如何构造计算压力改进值得方程,就成了分离式求解 法中的一个关键问题。
在交错网格系统中,关『.“、v的离散方程可通过对U、v番白的控制容积作积分而 得出。这时压力梯度的离散形式对虬为(PF—P,,)/dx,对k为(P.。一P,,)/ay,亦即相邻 两点间的压力羞构成了印/彘、印/砂,这就从根本上解决了采川一般网格系统时所遇到 的困难,这也是交错网格的成功之处。
Keywords:SIMPLE algorithm;turbulence modeling;gravity current; numerical simulation;computational fluid mechanics
上海交通大学博十后出站报告
第一章SIMPLE算法概述

2020高考数学第五章不等式、推理与证明、算法初步与复数考点测试38算法初步文(含解析)

2020高考数学第五章不等式、推理与证明、算法初步与复数考点测试38算法初步文(含解析)

考点测试 38算法初步高考概览本考点是高考必考知识点,常考题型为选择题、填空题,分值5分,中、低等难度考纲研读1. 认识算法的含义,认识算法的思想2.理解程序框图的三种基本逻辑结构:次序、条件、循环3.认识几种基本算法语句——输入语句、输出语句、赋值语句、条件语句、循环语句的含义一、基础小题1.给出如图程序框图,其功能是()A .求 a - b 的值B .求 b - a 的值C .求 | a - b | 的值D .以上都不对答案C分析求 |a -b | 的值.2.已知一个算法:① m = a ;②假如③假如b <m ,则c <m ,则m = b ,输出m = c ,输出m ,结束算法;不然履行第 m .3 步;假如 a = 3, b =6, c = 2,那么履行这个算法的结果是( )A 3B 6C 2D m答案C分析 当 a = 3, b = 6, c =2=3,∴ m =c = 2,即输出 m 的值为时,依照算法设计,履行后,2. 应选 C.m = a = 3<b = 6, c =2<m = a3.阅读下边的程序:INPUT xIF x<0 THENx=- xEND IFPRINT xEND则程序履行的目的是()A.务实数x 的绝对值 B .务实数x 的相反数C.求一个负数的绝对值 D .求一个负数的相反数答案A分析由程序可知,当输入的x<0时,取其相反数再赋值给x,其余状况x 不变,而后输出x,则程序履行的目的是务实数x 的绝对值,应选 A.4.阅读程序框图,该算法的功能是输出()A.数列 {2 n- 1} 的第 4 项B.数列 {2 n- 1} 的第 5 项C.数列 {2 n- 1} 的前 4 项和D.数列 {2 n- 1} 的前 5 项和答案B分析依程序框图,有下表:A1371531i23456因为5.当6>5,跳出循环,故输出A=31,而31=25-1,选m=5, n=2时,履行图中所示的程序框图,输出的B.S值为()A.20 B .42 C .60 D .180答案C分析当 m=5, n=2时,程序框图的运算过程以下表所示:k5432S152060故输出 S=60,应选 C.6.以下图程序框图的功能是:给出以下十个数:5, 9, 80, 43, 95, 73, 28, 17,60, 36,把大于60 的数找出来,则框图中的①②应分别填入的是()A.x>60?,i=i- 1 B .x<60?,i=i+1C.x>60?,i=i+ 1 D .x<60?,i=i-1答案C分析关于 A, D,因为i=i-1,则会进入死循环,而关于B,选出的数小于60. 应选C.7.在十进制中,2004=4×10 0+0×10 1+0×10 2+2×10 3,那么在五进制中数码2004折合成十进制为()A. 29 B .254 C .602 C . 2004答案B分析0123,故 B. 2004=4×5+0×5+0×5+2×5= 2548.当x= 0.2 ,用秦九韶算法算多式 f ( x)=3x6+4x5+5x4+6x3+7x2+8x+1的,需要做乘法和加法的次数分是()A.6,6 B . 5,6 C .5,5 D .6, 5答案A分析由f (x) =(((6+ 5)x+4)x+⋯+1)x+0,所以共需要 6 次加法和 6 次乘法,a x a a a a故 A.9.已知一个算法的程序框如所示,当出的果0 ,入的数x 的()A.-3 B .-3或 9C.3或-9 D .-9或-3答案B1x- 8,x≤0,分析本算法框的本求函数y=22- log 3x,x>0的零点,分状况求此分段函数的零点,易解得x=-3或 x=9,故 B.10.如所示的程序框的算法思路源于我国古代有名的“ 子节余定理”,此中“Mod(N,m) =n”表示正整数N除以正整数 m后的余数 n,比如:Mod(10,3)=1.行程序框,出的 i =()A.23 B .38 C .44 D .58答案A分析查验选项A:i=23, Mod(23, 3) =2, Mod(23,5) = 3,Mod(23, 7) = 2,知足题意,应选 A.11.如图是“二分法”解方程的流程图,在①~④处应填写的内容分别是()A.f ( a) f ( m)<0 ;a=m;是;否B.f ( b) f ( m)<0 ;b=m;是;否C.f ( b) f ( m)<0 ;m=b;是;否D.f ( b) f ( m)<0 ;b=m;否;是答案B分析因为题图是“二分法”解方程的流程图,所以判断框的内容是根的存在性定理的应用,所以填 f ( b) f ( m)<0;是,则直接考证精度,否,则先在赋值框中实现b= m 的互换,再考证精度,知足精度则输出结果,结束程序,所以③处填“是”,④处填“否”,在①~④处应填写的内容分别是f ( )( )<0;=;是;否.b f m b m12.下列图是用模拟方法预计圆周率π值的程序框图, P 表示预计结果,则图中空白框内填入 ()A.P=NB .P=4N 10001000 M4MC.P=1000 D .P=1000答案D分析利用几何概型,结构一个 1 的正方形及其内一个半径1、心角 90°M 12M4M的扇形,易知扇形的面S≈1000,又由面公式得S=4π ×1≈1000,解得π ≈1000,故D.二、高考小13.(2018 ·全国卷Ⅱ ) 算11111=1-+-+⋯+-,了下边的程序框,S2*******在空白框中填入()A.i=i+ 1 B .i=i+ 2C.i=i+ 3 D .i=i+ 4答案B1 1 111分析 由 S = 1-2+ 3- 4+⋯+ 99- 100,知程序框 先 奇数 累加,偶数 累加,最后再相减.所以在空白框中 填入i =i + 2,B.14.(2018 ·北京高考) 行如 所示的程序框 , 出的s ()1 5 7 7 A.2 B. 6 C. 6 D.12答案B11 1 1 1 21分析 k = 1,s = 1;s = 1+ ( - 1) ×1+ 1= 1- 2= 2,k = 2,2<3;s =2+ ( - 1)×1+ 2=1 1 55+ = , k =3,此 跳出循 ,所以 出 .故 B.2 36615.(2018 ·天津高考 ) 下 的程序框 ,运转相 的程序,若 入N 的 20,出 T 的 ()A .1B .2C .3D .4答案 B分析第一次循环T=1,i= 3;第二次循环T=1,i= 4;第三次循环T=2,i= 5,满足条件i≥5,结束循环.应选 B.16. (2017 ·全国卷Ⅰ) 右边程序框图是为了求出知足3n- 2n>1000的最小偶数n,那么在和两个空白框中,能够分别填入()A.A>1000?和n=n+ 1B.A>1000?和n=n+ 2C.A≤1000?和n=n+ 1D.A≤1000?和n=n+ 2答案D分析此题求解的是知足3n- 2n>1000 的最小偶数n,可判断出循环结构为当型循环结构,即知足条件要履行循环体,不知足条件要输出结果,所以判断语句应为A≤1000?,另外,所求为知足不等式的偶数解,所以中语句应为n= n+2.应选 D.17.(2017 ·全国卷Ⅲ) 履行下边的程序框图,为使输出S 的值小于91,则输入的正整()数 N的最小值为A.5 B.4 C.3 D.2答案D分析要求的是最小值,察看选项,发现选项中最小的为2,不如将 2 代入查验.当输入的N 为 2 时,第一次循环,= 100,=- 10,= 2;第二次循环,= 90,= 1,= 3,S M t S Mt此时退出循环,输出S=90,切合题意.应选 D.18.(2017 ·天津高考 ) 阅读下边的程序框图,运转相应的程序,若输入N 的值为24,则输出 N的值为()A.0 B.1 C.2 D.3答案C分析履行程序框图,输入N的值为24 时, 24 能被 3 整除,履行是,N=8,8≤3不可立,持续履行循环体;8 不可以被3 整除,履行否,N=7,7≤3不建立,持续履行循环体;7不可以被3 整除,履行否,N=6,6≤3不建立,持续履行循环体; 6 能被3 整除,履行是,N=2,2≤3建立,退出循环,输出N的值为2. 应选C.19.(2017 ·山东高考) 履行两次以下图的程序框图,若第一次输入的x 的值为7,第二次输入的x 的值为9,则第一次、第二次输出的 a 的值分别为()A.0,0 B . 1,1 C .0,1 D .1, 0答案D分析第一次输入 x=7,判断条件,4>7不建立,履行否,判断条件,7÷2=7, 7 不2能被 2 整除,履行否,b= 3,判断条件, 9>7 建立,履行是,输出a=1.第二次输入 x=9,判断条件,4>9不建立,履行否,判断条件,99÷2=2, 9不可以被2整除,履行否, b=3,判断条件,9>9不建立,履行否,判断条件,9÷3=3,9能被3整除,履行是,输出a=0.应选D.三、模拟小题20.(2018 ·衡阳二模)1927年德国汉堡大学的学生考拉兹提出一个猜想:关于每一个正整数,假如它是奇数,对它乘 3 再加1,假如它是偶数,对它除以2,这样循环,最后结果都能获得 1. 固然该猜想看上去很简单,但有的数学家以为“该猜想任何程度的解决都是现代数学的一大进步”.如图是依据考拉兹猜想设计的一个程序框图,则①处应填写的条件及输出的结果分别为()A.a是偶数?6B .a是偶数?8C.a是奇数?5D .a是奇数?7答案D分析阅读考拉兹提出的猜想,联合程序框图可得①处应填写的条件是“ a 是奇数?”,运转状况为a105168421i1234567所以输出的结果为i =7.应选D.21.(2018 ·郑州质检一) 我国古代数学文籍《九章算术》“盈不足”中有一道两鼠穿墙问题:“今有垣厚十尺,两鼠对穿,初日各一尺,大鼠日自倍,小鼠日自半,问几何日相遇?”现用程序框图描绘,以下图,则输出结果n=()A.5 B.4 C.3 D.2答案B分析初始 a=1, A=1,S=0,n=1,第一次循环: S=0+1+1=2, S 小于10,进入下一次循环;第二次循环:119n= n+1=2,a=, A=2, S=2++ 2=,S小于 10,进入下22219135一次循环;第三次循环:n= n+1=3,a=4,A=4,S=2+4+4=4,S 小于10,进入下一次循环;第四次循环:n =n1351n= 4,+1=4,=,=8,=4++8≥10,循环结束,此时a 8AS8应选 B.22.(2018 ·合肥质检一 ) 履行以下图程序框图,若输入的n 等于10,则输出的结果是()1 1A.2 B .-3 C.-2 D. 3答案C1+2分析a=2,i =1,知足 i ≤n=10,进入循环体,第一次循环:a==-3,i=2;1-2足 i ≤ n = 10,第二次循 : a =1+- 3=-1,i = 3; 足 i ≤ n = 10,第三次循 : a =1- -32111+-11+2 = 4; 足i≤ =10,第四次循 : =3= 5;⋯可看出a 的取 周= , = 2,13ina1i1--1-23期性 化, 且周期 4. 可知当 i = 11 与 i = 3 a 的取 同样,即1a =- ,此 ,不 足21i ≤ n = 10,跳出循 体, 出a =- 2,故 C.23.(2018 · 阳模 ) 我国明朝数学家程大位著的 《算法 宗》 里有一道 名世界的目:“一百 一百僧大僧三个更无争,小僧三人分一个,大、小和尚各几丁?”如 所示的程序框 反应了此 的一个求解算法, 出n 的 ( )A .20B .25C .30D .35 答案 B分析开始: = 20;第一步:= 80, = 60+ 80n =21;第二步:=79, =≠100,nm S3 mS797863+ 3 ≠100, n = 22;第三步: m = 78, S = 66+ 3 =92≠100, n = 23;第四步: m = 77, S=69 + 77 ≠100, = 24;第五步:= 76, = 72+76≠100, n = 25;第六步:= 75, =753 n mS3mS75= 100,此 S = 100 退出循 , 出n = 25. 故 B.+ 324.(2018 ·南昌摸底 ) 行如 所示的程序框 , 出 n 的 ()A.1 B.2 C.3 D.4答案C分析依照框图,可知n=1时, f ( x)=( x)′=1,它是偶函数,知足 f ( x)= f (-x),又方程 f ( x)=0无解,则 n=1+1=2;此时, f ( x)=( x2)′=2x,不知足 f ( x)= f (-x),则 n=2+1=3;再次循环, f ( x)=( x3)′=3x2,知足 f ( x)= f (-x),且方程 f ( x)=0有解x=0,跳出循环体,则输出n 的值为3,应选 C.25.(2018 ·深圳调研) 九连环是我国一种传统的智力玩具,其结构如图 1 所示,要将9个圆环所有从框架上解下( 或套上) ,不论是哪一种情况,都需要按照必定的规则.解下( 或套上) 所有9 个圆环所需的最少挪动次数可由如图 2 所示的程序框图获得,履行该程序框图,则输出的结果为()A.170 B .256 C . 341 D .682答案C分析由算法框图,可知i , S 的变化状况以下:i23456789S2510214285170341应选 C.26.(2018 ·邯郸摸底) 我国古代名著《庄子·天下篇》中有一句名言“一尺之棰,日取其半,万世不断”,其意思为:一尺的木棍,每日截取一半,永久都截不完.现将该木棍依此规律截取,以下图的程序框图的功能就是计算该木棍被截取7 天后所剩的长度( 单位:尺) ,则①②③处可分别填入的是()①②③1A i ≤7?s= s- i i =i +1B≤128?1= 2i s = s - ii i1Ci ≤7? s = s - 2ii =i + 1 Di ≤128?1 i = 2is = s -2i答案B分析程序框 的功能是 算木棍被截取7 天后节余部分的 度,在程序运转 程111中, 有:第1 次循 ,s = 1- 2,i = 4;第2 次循 ,s = 1-2- 4,i = 8;第3 次循 ,1 1 1s = 1- 2-4- 8, i= 16;⋯;第7 次循 ,1 1 1s = 1- 2- 4-⋯- 128, i= 256,此 跳出循1体,据此判断可知在判断框① 填入“ i≤128?”, 行框② 填入“s = s - i ”,③填入“ i = 2i ”,故 B.本考点在近三年高考取未波及此 型.。

儿童编程scratch 数学出题——个性化软件定制教案

儿童编程scratch 数学出题——个性化软件定制教案

少儿编程教学资料,仅供参考,需要可下载并修改后使用!第16课口算生成器任务导航有一天你的数学老师找到你,期望你能够帮助她编写一个软件,给低年级的同学出口算题,想必这难不倒你吧。

同学们来帮助增加口算学习的兴趣,提高口算学习的成绩,我们为小学生们编写一个口算练习的小程序吧!如图16.1所示。

图16.1 口算生成器界面Follow Me1. 简单口算生成器任务分析简单口算生成器的功能如表16.1所示。

表16.1 简单口算练习的功能功能按空格键,在第一个链表中出题,在第二个链表中给出答案。

按空格键时,先清空两个链表的所有信息;使用重复n次控制程序出题数目;在“出题”链表中显示100道加法题,加数分别为0到99之间随机数,随机数可使用“数字和逻辑运算”分类中的来产生,加法算式可使用来连接生成;在第二个链表中显示出两个加数的和,如图16.2简单口算生成器界面。

图16.2简单口算生成器界面第1步:新建Scratch程序,保存文件到指定位置,文件名为“简单口算生成器”。

1.第2步:删除“角色1”。

2.第3步:新建两个链表分别命名为“出题”和“答案”,分别调整好大小,放置到舞台中合适位置,如图16.3所示。

图16.3简单口算生成器设计界面3.第4步:新建三个变量分别命名为“加数1”、“加数2”和“结果”,三个变量都不显示在舞台上,如图16.4。

图16.4 隐藏变量4.第五步:添加脚本如图16.5所示。

图16.5简单口算生成器的脚本5.第6步:按下“空格”键,运行并调试程序。

DIY试修改“简单口算生成器”,尝试出60道加数范围是0到50的加法题。

2. 加减法口算生成器任务分析刚才的“简单口算生成器”只能出加法题,不能出减法题,所以现在来改进简单口算生成器,使其能够随机生成加减法算术题,如图16.6所示。

功能按“空格”键,在第一个链表中随机生成100道加法、减法口算题,在第二个链表中产生答案。

可以再使用一个随机数,使其随机产生1或2,当该数为1时,生成加法题,为2时生成减法题。

Simple软件--快速掌握手册

Simple软件--快速掌握手册

Simple企业管理系统快速上手指南广州斯盟派网络科技有限公司版权信息本文件的版权属于广州斯盟派网络科技有限公司任何形式的散发都必须事先得到广州斯盟派网络科技有限公司的书面许可目录第一章系统介绍 (4)1.1SIMPLE 企业管理系统概述 (4)1.2采购管理介绍 (4)1.3库存管理介绍 (5)1.4销售管理介绍 (5)1.5生产管理介绍 (5)1.6应收管理介绍 (6)1.7应付管理介绍 (6)第二章建立资料 (7)2.1货物分类 (7)2.2货物资料 (8)2.3仓库资料 (8)2.4采购员资料 (9)2.5供应商资料 (9)2.6业务员资料 (9)2.7客户资料 (9)2.8其它基础资料 (9)第三章数据初始化 (10)3.1库存数量初始化 (10)3.2应收款初始化 (10)3.3预收款初始化 (11)3.4应付款初始化 (11)3.5预付款初始化 (11)3.6现金银行初始化 (11)3.7科目余额初始化 (11)第四章采购流程 (13)4.1采购订单 (13)4.2采购订单查询 (14)4.3采购收货单 (16)4.4采购收货单查询 (17)4.5库存汇总表 (19)4.6付款单 (20)4.7应付款报表 (21)第五章销售管理 (24)5.1销售订单 (24)5.2销售发货单 (25)5.3销售退货单 (27)5.4收款单 (28)5.5.应收款报表 (29)第六章应收应付管理 (30)6.1先收货后付款业务 (30)6.2先预付款后收货业务 (30)6.3先发货后收款业务 (33)6.4先收款后发货业务 (34)第七章权限管理 (37)7.1操作员资料 (37)7.2新的操作员进入系统 (39)第八章月末结账 (41)8.1月末库存盘点 (41)8.2月末结账 (43)第一章系统介绍1.1simple 企业管理系统概述Simple企业管理系统是通过在多个行业众多客户的成功应用,形成的一套成熟稳定的ERP 产品。

Simple C语言词法分析器的探讨与实践

Simple C语言词法分析器的探讨与实践
号> =。
> < > = =
G T L T G E L E
EQ
NE
[ ] ( )
/卑
年/
[ ] { )
[ zE _ A— r
][ -Za— A zO-9]
id

号 项
}l
向 符
表 表
请 注 意 . 果按 照 从 s r状 态 到 达 接 受状 态 如 tt a 3的边 的顺 序 . 意 味着 > 则 和一 个 与 之无 关 的字 符 已经 被读 过 。 由于这 个无 关字 符 不是关 系操作 符> 的一 部 分 .而是 下 一 个词 素 的一 部 分 .所 以 f - 0 r W r 针 必 须 回滚 一 个 字 符 。状 态 上 的 水表 示 ad指
wh i e 1 v i o d
( )
( )
符 , 果 该 字 符 是 > 则 转 向状 态 1 否 则 便 告 识 如 , , 别 > : 败 。到达 状 态 1时 。 下 一 个 字符 , 或> 失 遇 如
果 它 是=。 转 向状 态 2 否 则 标 有 o e 的边 表 明 则 . tr h 已经 转 向状态 3 。在 状态 2上 有双 圈 . 示它 是接 表 受状 态 。 当进 入这 个 状态 时 . 态转 换 图识别 了记 状
+ 一 木/ < < > >= == != = =
;, { ) 【 ]() 厂 水 水 /
图 31状态 转 换 图 .
( 下转 第 15页 ) 6
21 0 2年第 4期
福 建 电

15 6
界 面 和遥 控逻 辑 .实 现一 个通 用 遥 控器 遥 控 多个 间的服 务发 现 问题 .这 些将 在后 续 的工作 中展 开 不 同类 型 的受 控 应 用 , 图 2 图 3所 示 。在 系统 研 究 。 如 ,

SimpleAPP重要函数解析

SimpleAPP重要函数解析

SimpleAPP重要函数解析Z-Stack1.4.3-1.2.1中SampleApp中重要函数解析在上几篇文章中我们分析了,整个OSAL的流程,还有添加自己的任务,以及如何运行到自己定义的任务。

这一篇文章主要是分析一下,自己定义的任务中要完成的功能,需要的事件函数是怎样的。

这个例子就是一个简单的点对点的数据发送,其中涉及到较少的网络配置,其中最主要的两个函数是SampleApp_ProcessEvent(uint8 task_id,uint16 events),和SampleApp_Init(taskID),一个是任务的处理函数,一个是初始化函数。

SampleApp_Init(taskID)这个函数在前面的文章中也已经分析过了,现在主要关注一下SampleApp_ProcessEvent(uint8 task_id,uint16 events)函数的实现。

每个应用任务都通过SampleApp_ProcessEvent()函数来处理任务中的事件。

一旦SampleApp_TaskID任务的某个OSAL 事件发生,那么就可以通过调用SampleApp_ProcessEvent()函数来处理。

在SampleApp_ProcessEvent()中有一个事件处理循环,循环检测是哪个事件发生。

/******************************************************************** ** @fn SampleApp_ProcessEvent** @brief Generic Application Task event processor. This function* is called to process all events for the task. Events* include timers, messages and any other user defined events.* 这个函数被用来调用处理所有的事件,事件有定时器消息用户自己定义的* @param task_id - The OSAL assigned task ID. 任务ID号* @param events - events to process. This is a bit map and can* contain more than one event. 处理的事件,这是一个位图** @return none*/uint16 SampleApp_ProcessEvent( uint8 task_id, uint16 events ){afIncomingMSGPacket_t *MSGpkt;//系统事件号SYS_EVENT_MSG = 0x8000if ( events & SYS_EVENT_MSG ){//检索收到的命令,没有收到返回NULLMSGpkt = (afIncomingMSGPacket_t*)osal_msg_receive( SampleApp_TaskID );while ( MSGpkt ) //如果不为空时,判断消息的类型{switch ( MSGpkt->hdr.event ) //这里是判断SYS_EVENT_MSG事件类型,不同的SYS_EVENT_MSG类型需要不同的处理。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

C$DEBUGCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC This computer program was copied from the graduate student courseC program of the University of Minnesota. Part of it was re-formu- C lated to meet the microcomputer environment. Some inappropriateC expressions were also corrected. The program is used only for theC teaching purpose. No part of it may be published. You may use itC as a frame to re-develop your own code for research purpose.C --------Instructor of Numerical Heat Transfer, XJTU,1998.12-------C --------注释者:吴志根, NHT -CFD/EHT 科研中心,2003.7,版权所有-------CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C*******************************************************************C -------------------------MAIN PROGRAM----------------------------C*******************************************************************LOGICAL LSTOPCOMMON/CNTL/LSTOPC********************************************************************OPEN(8,FILE='RESULT.txt') !输出数据结果文件名CALL GRID !绘制网格CALL SETUP1 !设置与网格系统有关的,在计算过程中保持不变的几何参数及设置CALL START !设置初值10 CALL DENSE !密度的计算(可压缩流体密度计算时需要)CALL BOUND !边界条件设定CALL OUTPUT !输出打印IF(.NOT.LSTOP) GO TO 15 !由iter 控制CLOSE(8) !计算完毕,关闭打开的文件RESULT.TXTSTOP !停止计算15 CALL SETUP2 !主要计算核心模块GO TO 10 !再来一次计算ENDCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCSUBROUTINE DIFLOW !计算乘方格式中的系数A(|p|),扩导D(参考P151),P142(5-10), P146(5-19)C********************************************************************COMMON/COEF/FLOW,DIFF,ACOFC******************************************************************ACOF=DIFF ! A(|p|),DIF(FLOW.EQ.0.) RETURN !流量为零,不考虑对流,只有扩散TEMP=DIFF-ABS(FLOW)*0.1! /Pe ul ρ=Γ,/P u x ρδ∆=Γ,()x x e a D A P ∆=⋅,/x x e F D P ∆=,0.1D F −ACOF=0. !550 10max[0,(10.1)](10.1) 101 0e E e e e e e P a P P P D P ∆∆∆∆∆ > =−−< = IF(TEMP.LE.0.) RETURN !10e P ∆>TEMP=TEMP/DIFF !10.1e P ∆−ACOF=DIFF*TEMP**5 !5(10.1)e D P ∆⋅−RETURNENDCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CSUBROUTINE SOLVE !采用交替方向先迭代法求解代数方程并补以块修正技术以促进收敛,TDMA ,P281,P269,P100PARAMETER(NI=52,NJ=52,NIJ=NI,NFMAX=10,NFX3=NFMAX+3)C******************************************************************CHARACTER TITLE*8LOGICAL LSOLVE,LPRINT,LBLK,LSTOPCOMMON F(NI,NJ,NFMAX),P(NI,NJ),RHO(NI,NJ),GAM(NI,NJ),CON(NI,NJ),1 AIP(NI,NJ),AIM(NI,NJ),AJP(NI,NJ),AJM(NI,NJ),AP(NI,NJ),2 X(NI),XU(NI),XDIF(NI),XCV(NI),XCVS(NI),3 Y(NJ),YV(NJ),YDIF(NJ),YCV(NJ),YCVS(NJ),4 YCVR(NJ),YCVRS(NJ),ARX(NJ),ARXJ(NJ),ARXJP(NJ),5 R(NJ),RMN(NJ),SX(NJ),SXMN(NJ),XCVI(NI),XCVIP(NI)COMMON DU(NI,NJ),DV(NI,NJ), FV(NI),FVP(NI),1 FX(NI),FXM(NI),FY(NJ),FYM(NJ),PT(NIJ),QT(NIJ)COMMON/INDX/NF,NP,NRHO,NGAM,L1,L2,L3,M1,M2,M3,1IST,JST,ITER,LAST,TITLE(NFX3),RELAX(NFX3),TIME,DT,XL,YL,2IPREF,JPREF,LSOLVE(NFX3),LPRINT(NFX3),LBLK(NFX3),MODE3,NTIMES(NFX3),RHOCONC******************************************************************ISTF=IST-1JSTF=JST-1IT1=L2+IST ! 在Solve 中临时值IT2=L3+IST ! 在Solve 中临时值JT1=M2+JST ! 在Solve 中临时值JT2=M3+JST ! 在Solve 中临时值C******************************************************************DO 999 NT=1,NTIMES(NF) !The number of repetitions of the sweeps inN=NF ! NF 表示所求的变量名称,1=U ,2=V ,3=P ,……….C-------------------------------------------------------------------IF(.NOT.LBLK(NF)) GO TO 10 !当LBLK 为真时,对F(I,J,NF)进行块修正,参看TDMA PT(ISTF)=0. ! TDMA 中系数,i P ,P100;QT(ISTF)=0. ! TDMA 中系数,i Q ,P100;DO 11 I=IST,L2BL=0. !块修正中系数i A 以P100(4-38b )为基准;i B :以P281(7-23)为基准BLP=0. !块修正中系数i B 以P100(4-38b )为基准;i P :以P281(7-23)为基准BLM=0. !块修正中系数i C 以P100(4-38b )为基准;i N :以P281(7-23)为基准 BLC=0. !块修正中系数i D 以P100(4-38b )为基准;i C :以P281(7-23)为基准 DO 12 J=JST,M2 !参考P281(7-24),参考P99(4-38a ),BL=BL+AP(I,J) !p a ∑IF(J.NE.M2) BL=BL-AJP(I,J) !2,P N S J M JST BL a a a ≠=−−∑∑∑时, IF(J.NE.JST) BL=BL-AJM(I,J) !2,P S J M J JST BL a a =≠=−∑∑时,;2,P N J M J JST BL a a ≠==−∑∑时,;2,P J M J JST BL a ===∑时,BLP=BLP+AIP(I,J) !E a ∑BLM=BLM+AIM(I,J) !W a ∑ BLC=BLC+CON(I,J)+AIP(I,J)*F(I+1,J,N)+AIM(I,J)*F(I-1,J,N) 1 +AJP(I,J)*F(I,J+1,N)+AJM(I,J)*F(I,J-1,N)-AP(I,J)*F(I,J,N) !1,1,,1,1,E i j W i j N i j S i j P i j b a a a a a φφφφφ+−+−+++++∑∑∑∑∑∑,P282(7-24a )12 CONTINUEDENOM=BL-PT(I-1)*BLM ! i i i A C P −⋅IF(ABS(DENOM/BL).LT.1.E-10) DENOM=1.E25 !0i i i i A A C P >>−⋅,i i i A C P ≈⋅,不修正;0ii i iA A C P >−⋅,修正PT(I)=BLP/DENOM !若DENOM =2510,PT(I)=0;ii i iB AC P −⋅QT(I)=(BLC+BLM*QT(I-1))/DENOM !若DENOM =2510,QT(I)=0;i i ii i iD C Q A C P +⋅−⋅11 CONTINUEBL=0. !参考P100 BL →余量:,i j R P281DO 13 II=IST,L2I=IT1-II !递归求解BL=BL*PT(I)+QT(I) !余量累乘BL=[[……….]*PT(I-1)+QT(I-1)]*PT(I)+QT(I)DO 13 J=JST,M213 F(I,J,N)=F(I,J,N)+BL ! 同一竖快中采用同一个修正值C------------------------------------------------------------------------PT(JSTF)=0. ! TDMA 中系数,i P ,P100;QT(JSTF)=0. ! TDMA 中系数,i Q ,P100;DO 21 J=JST,M2BL=0. !块修正中系数i A 以P100(4-38b )为基准;i B :以P281(7-23)为基准BLP=0. !!块修正中系数i B 以P100(4-38b )为基准;i P :以P281(7-23)为基准 BLM=0. !!块修正中系数i C 以P100(4-38b )为基准;i N :以P281(7-23)为基准 BLC=0. !!块修正中系数i D 以P100(4-38b )为基准;i C :以P281(7-23)为基准 DO 22 I=IST,L2 !参考P281(7-24),参考P99(4-38a ),BL=BL+AP(I,J) !p a ∑IF(I.NE.L2) BL=BL-AIP(I,J) ! 2,P W W J L IST BL a a a ≠=−−∑∑∑时,IF(I.NE.IST) BL=BL-AIM(I,J) ! 2,P W J L J IST BL a a =≠=−∑∑时,;2,P E J L J IST BL a a ≠==−∑∑时,;2,P J M J JST BL a ===∑时,BLP=BLP+AJP(I,J) !N a ∑BLM=BLM+AJM(I,J) !W a ∑ BLC=BLC+CON(I,J)+AIP(I,J)*F(I+1,J,N)+AIM(I,J)*F(I-1,J,N) 1 +AJP(I,J)*F(I,J+1,N)+AJM(I,J)*F(I,J-1,N)-AP(I,J)*F(I,J,N) !1,1,,1,1,E i j W i j N i j S i j P i j b a a a a a φφφφφ+−+−+++++∑∑∑∑∑∑,P282(7-24a )22 CONTINUEDENOM=BL-PT(J-1)*BLM ! i i i A C P −⋅IF(ABS(DENOM/BL).LT.1.E-10) DENOM=1.E25 !0ii i i A A C P >>−⋅,i i i A C P ≈⋅,不修正;0ii i iA A C P >−⋅,修正PT(J)=BLP/DENOM ! 若DENOM =2510,PT(I)=0;ii i iB AC P −⋅QT(J)=(BLC+BLM*QT(J-1))/DENOM ! !若DENOM =2510,QT(I)=0;i i ii i iD C Q A C P +⋅−⋅21 CONTINUEBL=0. !参考P100 BL →余量:,i j R P281DO 23 JJ=JST,M2J=JT1-JJ !递归求解BL=BL*PT(J)+QT(J) !余量累乘BL=[[……….]*PT(I-1)+QT(I-1)]*PT(I)+QT(I)DO 23 I=IST,L2 ! 同一横条中采用同一个修正值23 F(I,J,N)=F(I,J,N)+BL10 CONTINUEC---------------------------------------!以上是块修正,以下是直接TDMA 交替方向迭代DO 90 J=JST,M2 !参考P99(4-38a),P273(7-19a,b);本块为I 方向,从下到上 PT(ISTF)=0. ! TDMA 中系数,i P ,P100;QT(ISTF)=F(ISTF,J,N) ! TDMA 中系数,i Q ,P100;DO 70 I=IST,L2DENOM=AP(I,J)-PT(I-1)*AIM(I,J) ! i i i A C P −⋅;AP(I,J):P a ,i A ;AIM(I,J):W a ,i C PT(I)=AIP(I,J)/DENOM !i i i iB AC P −⋅;AIP(I,J):E a ,i B TEMP=CON(I,J)+AJP(I,J)*F(I,J+1,N)+AJM(I,J)*F(I,J-1,N) !N N S S b a T a T ++:iD ;P273(7-19a,b) QT(I)=(TEMP+AIM(I,J)*QT(I-1))/DENOM !i i i i i i D C Q A C P +⋅−⋅ 70 CONTINUEDO 80 II=IST,L2I=IT1-II !递归求解80 F(I,J,N)=F(I+1,J,N)*PT(I)+QT(I) !P100(a),1i i i i P Q φφ+=⋅+90 CONTINUEC------------------------------------------------------------------DO 190 JJ=JST,M3 !参考P99(4-38a),P273(7-19a,b);本块为I 方向,从上到下J=JT2-JJPT(ISTF)=0. ! TDMA 中系数,i P ,P100;QT(ISTF)=F(ISTF,J,N) ! TDMA 中系数,i Q ,P100;DO 170 I=IST,L2DENOM=AP(I,J)-PT(I-1)*AIM(I,J) ! i i i A C P −⋅;AP(I,J):P a ,i A ;AIM(I,J):W a ,i C PT(I)=AIP(I,J)/DENOM !i i i iB AC P −⋅;AIP(I,J):E a ,i B TEMP=CON(I,J)+AJP(I,J)*F(I,J+1,N)+AJM(I,J)*F(I,J-1,N) !N N S S b a T a T ++:iD ;P273(7-19a,b) QT(I)=(TEMP+AIM(I,J)*QT(I-1))/DENOM !i i i i i iD C Q A C P +⋅−⋅ 170 CONTINUEDO 180 II=IST,L2I=IT1-II !递归求解180 F(I,J,N)=F(I+1,J,N)*PT(I)+QT(I) !P100(a),1i i i i P Q φφ+=⋅+190 CONTINUEC------------------------------------------------------------------DO 290 I=IST,L2 !参考P99(4-38a),P273(7-19a,b);本块为J 方向,从左到右PT(JSTF)=0. ! TDMA 中系数,i P ,P100;QT(JSTF)=F(I,JSTF,N) ! TDMA 中系数,i Q ,P100;DO 270 J=JST,M2DENOM=AP(I,J)-PT(J-1)*AJM(I,J) ! i i i A C P −⋅;AP(I,J):P a ,i A ;AJM(I,J):S a ,i C PT(J)=AJP(I,J)/DENOM !i i i iB AC P −⋅;AJP(I,J):N a ,i B TEMP=CON(I,J)+AIP(I,J)*F(I+1,J,N)+AIM(I,J)*F(I-1,J,N) !E E W W b a T a T ++:iD ;P273(7-19a,b) QT(J)=(TEMP+AJM(I,J)*QT(J-1))/DENOM !i i i i i i D C Q A C P +⋅−⋅ 270 CONTINUEDO 280 JJ=JST,M2J=JT1-JJ !递归求解280 F(I,J,N)=F(I,J+1,N)*PT(J)+QT(J) ! P100(a),1i i i i P Q φφ+=⋅+290 CONTINUEC-------------------------------------------------------------------DO 390 II=IST,L3 !参考P99(4-38a),P273(7-19a,b);本块为J 方向,从右到左I=IT2-IIPT(JSTF)=0. ! TDMA 中系数,i P ,P100;QT(JSTF)=F(I,JSTF,N) ! TDMA 中系数,i Q ,P100; DO 370 J=JST,M2DENOM=AP(I,J)-PT(J-1)*AJM(I,J) ! i i i A C P −⋅;AP(I,J):P a ,i A ;AJM(I,J):S a ,i C PT(J)=AJP(I,J)/DENOM !i i i iB AC P −⋅;AJP(I,J):N a ,i B TEMP=CON(I,J)+AIP(I,J)*F(I+1,J,N)+AIM(I,J)*F(I-1,J,N) !E E W W b a T a T ++:iD ;P273(7-19a,b) QT(J)=(TEMP+AJM(I,J)*QT(J-1))/DENOM !i i i i i iD C Q A C P +⋅−⋅ 370 CONTINUEDO 380 JJ=JST,M2J=JT1-JJ !递归求解380 F(I,J,N)=F(I,J+1,N)*PT(J)+QT(J) ! P100(a),1i i i i P Q φφ+=⋅+390 CONTINUEC*******************************************999 CONTINUEENTRY RESET !重新设置某些参数初始值DO 400 J=2,M2DO 400 I=2,L2CON(I,J)=0. !baAP(I,J)=0. !P400 CONTINUERETURNEND CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE SETUPC******************************************************************PARAMETER(NI=52,NJ=52,NIJ=NI,NFMAX=10,NFX3=NFMAX+3)CHARACTER TITLE*8LOGICAL LSOLVE,LPRINT,LBLK,LSTOPCOMMON F(NI,NJ,NFMAX),P(NI,NJ),RHO(NI,NJ),GAM(NI,NJ),CON(NI,NJ),1 AIP(NI,NJ),AIM(NI,NJ),AJP(NI,NJ),AJM(NI,NJ),AP(NI,NJ),2 X(NI),XU(NI),XDIF(NI),XCV(NI),XCVS(NI),3 Y(NJ),YV(NJ),YDIF(NJ),YCV(NJ),YCVS(NJ),4 YCVR(NJ),YCVRS(NJ),ARX(NJ),ARXJ(NJ),ARXJP(NJ),5 R(NJ),RMN(NJ),SX(NJ),SXMN(NJ),XCVI(NI),XCVIP(NI)COMMON DU(NI,NJ),DV(NI,NJ), FV(NI),FVP(NI),1 FX(NI),FXM(NI),FY(NJ),FYM(NJ),PT(NIJ),QT(NIJ)COMMON/INDX/NF,NP,NRHO,NGAM,L1,L2,L3,M1,M2,M3,1IST,JST,ITER,LAST,TITLE(NFX3),RELAX(NFX3),TIME,DT,XL,YL,2IPREF,JPREF,LSOLVE(NFX3),LPRINT(NFX3),LBLK(NFX3),MODE3,NTIMES(NFX3),RHOCONCOMMON/CNTL/LSTOPCOMMON/SORC/SMAX,SSUMCOMMON/COEF/FLOW,DIFF,ACOFDIMENSION U(NI,NJ),V(NI,NJ),PC(NI,NJ)EQUIVALENCE(F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3),PC(1,1))DIMENSION COF(NI,NJ,6),COFU(NI,NJ,6),COFV(NI,NJ,6),COFP(NI,NJ,6)EQUIVALENCE(COF(1,1,1),CON(1,1))DIMENSION UHAT(NI,NJ),VHAT(NI,NJ)C******************************************************************1 FORMAT(//15X,'COMPUTATION IN CARTESIAN COORDINATES') !打印设置2 FORMAT(//15X,'COMPUTATION FOR AXISYMMETRIC SITUATION') !打印设置3 FORMAT(//15X,'COMPUTATION IN POLAR COORDINATES') !打印设置4 FORMAT(14X,40(1H*),//) !打印设置DATA LSTOP,LSOLVE,LPRINT/.FALSE.,NFX3*.FALSE.1 ,NFX3*.FALSE./ !默认值设置;数组LSOLVE,LPRINT:FALSEDATA LBLK/NFX3*.TRUE./ !默认值设置;数组LBLK:TRUEDATA MODE,LAST,TIME,ITER/1,5,0.,0/ !默认值设置;DATA RELAX,NTIMES/NFX3*1.,NFX3*1/ !默认值设置;数组RELAX=1,NTIMES=1 DATA DT,IPREF,JPREF,RHOCON/1.E+30,1,1,1./ !在非稳定中需要修改C-------------------------------------------------------------------ENTRY SETUP1 !设置与网格系统有关的,在计算过程中保持不变的几何参数及设置 NP=NFMAX+1 !NFMAX+1; P(I,J)与F(I,J,NP)共用存储器,NFMAX=10,NP=11NRHO=NP+1 !NFMAX+2; RHO(I,J) 与 F(I,J,NRHO) 共用存储器,NRHO=12NGAM=NRHO+1 !NFMAX+3; GAM(I,J) 与 F(I,J,NGAM) 共用存储器,NGAM=13 L2=L1-1 !L2,L3,M2,M3相对位置(参照附图1)L3=L2-1M2=M1-1M3=M2-1X(1)=XU(2) !X(1)=XU(2)=0,相对位置(参照附图1)DO 5 I=2,L2 !参照附图15 X(I)=0.5*(XU(I+1)+XU(I))!网格格点的x轴位置,控制面XU(I)已在UGRID中求出X(L1)=XU(L1)Y(1)=YV(2) !Y(1)=YV(2)=0,相对位置(参照附图1)DO 10 J=2,M210 Y(J)=0.5*(YV(J+1)+YV(J)) !网格格点的y轴位置Y(M1)=YV(M1)DO 15 I=2,L1 !以下到22-3都是x方向设置15 XDIF(I)=X(I)-X(I-1) !网格x轴间距DO 18 I=2,L218 XCV(I)=XU(I+1)-XU(I) !各主控制体积x轴宽度DO 20 I=3,L220 XCVS(I)=XDIF(I) ! 交错网格中U (I,J)在x方向上的控制容积宽度XCVS(3)=XCVS(3)+XDIF(2) !边界定义(参照附图1)XCVS(L2)=XCVS(L2)+XDIF(L1) !边界定义(参照附图1)DO 22 I=3,L3XCVI(I)=0.5*XCV(I) !主控制体积x轴宽度XCV(I)与U(I,J)控制容积宽度重叠部分 22 XCVIP(I)=XCVI(I) !主控制体积x轴宽度XCV(I)与U(I+1,J)控制容积宽度重叠部分XCVIP(2)=XCV(2) !边界定义(参照附图1),参看YL,YLM 计算式XCVI(L2)=XCV(L2) !边界定义(参照附图1)DO 35 J=2,M1 !以下到45-3都是y方向设置35 YDIF(J)=Y(J)-Y(J-1) !网格y轴间距DO 40 J=2,M240 YCV(J)=YV(J+1)-YV(J) !各主控制体积y轴宽度DO 45 J=3,M245 YCVS(J)=YDIF(J) ! 交错网格中V (I,J)在y方向上的控制容积宽度YCVS(3)=YCVS(3)+YDIF(2) !边界定义(参照附图1)YCVS(M2)=YCVS(M2)+YDIF(M1) !边界定义(参照附图1)IF(MODE.NE.1) GO TO 55 !坐标系参数,1直角,2轴对称圆柱,3极坐标(默认1)DO 52 J=1,M1 !直角坐标半径(三种坐标表达式统一书写)RMN(J)=1.0 ! V(I,J)所在点对应的半径52 R(J)=1.0 !直角坐标系中名义半径为1(参考书P90)GO TO 5655 DO 50 J=2,M1 !轴对称圆柱,极坐标半径50 R(J)=R(J-1)+YDIF(J)!R(1)已经定义(参考例4),每个网格点(I,J)对应的半径RMN(2)=R(1) !边界定义(RMN(I)与XU(I)类似),(参照附图1)DO 60 J=3,M260 RMN(J)=RMN(J-1)+YCV(J-1) ! V(I,J)所在点对应的半径RMN(M1)=R(M1) !边界定义(参照附图1)56 CONTINUE !直角,圆柱轴对称坐标东西尺度系数设置DO 57 J=1,M1SX(J)=1. !主节点Y(J)所对应x 方向的尺度系数 (参考书P90)SXMN(J)=1. !由YV(J)点所对应x 方向的尺度系数(参考书P90)IF(MODE.NE.3) GO TO 57 !极坐标东西尺度系数设置(参考书P90)SX(J)=R(J) !参考书P90,SX =rIF(J.NE.1) SXMN(J)=RMN(J) !尺度系数边界定义,SXMN(1)=1,RMN(1)无定义57 CONTINUEDO 62 J=2,M2YCVR(J)=R(J)*YCV(J) ! 主控制体积中r y ∆,东西导热面积 1 y r r r r ∆= ∆ ∆直角圆柱极坐标ARX(J)=YCVR(J) ! 主控制体积垂直于x 方向的面积IF(MODE.NE.3) GO TO 62 !参考书P90ARX(J)=YCV(J) ! r ∆,参考上句注释,P89:r ∆;11YCV(J)1r y ∆×∆×=×=62 CONTINUEDO 64 J=4,M364 YCVRS(J)=0.5*(R(J)+R(J-1))*YDIF(J) ! v(I,J)控制容积的面积r y ∆,东西导热面积YCVRS(3)=0.5*(R(3)+R(1))*YCVS(3) !边界定义,参考附图1,讲义P10图d ,书P90表4-2YCVRS(M2)=0.5*(R(M1)+R(M3))*YCVS(M2) !注意:R(1),R(2),R(3);R(M1),R(M2),R(M3)之间位置关系;V 2,U 2均不求解,V 3控制容积为YCVRS(3),V M2控制容积为YCVRS(M2)IF(MODE.NE.2) GO TO 67 !轴对称圆柱坐标下DO 65 J=3,M3 !the part of ARX(J) that overlaps on the C.V. for V(I,J)ARXJ(J)=0.25*(1.+RMN(J)/R(J))*ARX(J) !主控制体积垂直于x 方向的面积与V(I,J)控制容积垂直于x 方向的面积重叠部分;RMN(J)+R(J)R(J)YCV(J)2R(J)2×××65 ARXJP(J)=ARX(J)-ARXJ(J) ! 主控制体积垂直于x 方向的面积与V(I,J+1)控制容积垂直于x 方向的面积重叠部分GO TO 68 67 DO 66 J=3,M3 !直角、极坐标下 ARXJ(J)=0.5*ARX(J)66 ARXJP(J)=ARXJ(J)68 ARXJP(2)=ARX(2)!边界定义 ARXJ(M2)=ARX(M2) !边界定义DO 70 J=3,M3FV(J)=ARXJP(J)/ARX(J) !质量流量界面插值系数70 FVP(J)=1.-FV(J) !质量流量界面插值系数DO 85 I=3,L2FX(I)=0.5*XCV(I-1)/XDIF(I) !一般界面插值系数,x 方向,P200 85 FXM(I)=1.-FX(I) FX(2)=0. !边界定义FXM(2)=1. !边界定义 FX(L1)=1. !边界定义FXM(L1)=0. !边界定义DO 90 J=3,M2FY(J)=0.5*YCV(J-1)/YDIF(J) ! 一般界面插值系数,x 方向,P200,如上90 FYM(J)=1.-FY(J) !密度RHOM(在V(I,J)位置): 插值系数:FY(J) RHO(I,J)+FYM(J) RHO(I,J-1) FY(2)=0. !边界定义 FYM(2)=1.!边界定义 FY(M1)=1.!边界定义 FYM(M1)=0. !边界定义CON,AP,U,V,RHO,PC AND P ARRAYS ARE INITIALIZED HEREDO 95 J=1,M1 !设置速度、压力、密度等初场 DO 95 I=1,L1PC(I,J)=0. !压力场压力初始修正值U(I,J)=0. !x 轴方向初速度场U.V(I,J)=0. !y 轴方向初速度场VCON(I,J)=0. !在离散化方程中常数项 b 或者在GAMSOR 模块中表示源项S CAP(I,J)=0. !系数a P 或者在GAMSOR 模块中S PRHO(I,J)=RHOCON !初始密度场设置P(I,J)=0. !初始压力场设置95 CONTINUE !目录名输出设置IF(MODE.EQ.1) PRINT 1!参看初始块format 设置 IF(MODE.EQ.1) WRITE(8,1)IF(MODE.EQ.2) PRINT 2I,J ARXJP(J)ARXJ(J)R(1)1R(J)RMN(J)FX FXM I,J FXM 1XDIF(2)FXL1XDIF(L1)Y DIF(2)Y DIF(M1)FY M 1FY FY FY M M1I.JIF(MODE.EQ.2) WRITE(8,2) IF(MODE.EQ.3) PRINT 3 IF(MODE.EQ.3) WRITE(8,3) PRINT 4 WRITE(8,4) RETURNC--------------------------------------------------------------- ENTRY SETUP2CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C 1.建立每一轮迭代中变量U 、V 及P*和其他变量离散方程的系数 a E 、a W 、a N 、a S 、a P 及b 项。

相关文档
最新文档