GAMS用户指南
GAMS 经典教材内容提要GAMS软件是建立和求解大型数学规划问题的优秀软件包之一,在各领域有着广泛的应用。本书是GAMS软件用户指南,全书分两篇,第一篇为 Windows gAMs2.50用户指南,详细地介绍了GAMS语言的各个组成部分,并对一些高级主题进行了专门讨论;第二篇讨论几个常用的求解器,包括 BARON、 Cplex10、 DICOPT、 MINOS、MOSEK、PATH4.6、SBB、 MPSGE书是建立和求解大型数学规划问题的高级计算机软件参考手册,可供各領域从事规划设计和管理的人员參考使用,也可作为高等院校师生的参考工具书在版编目(CIP)数据GAMS用户指南/魏传江等编译.一北京:中国水利水电出版社,2009.11(通用建模软件技术丛书)I.①G…Ⅱ.①魏…Ⅲ.①数学模型一建立模型应用软件,GAM一指南Ⅳ.①022-39中国版本图书馆CIP数据核字(2009)第208573号书名|用速椰软件技术外书作者魏传江王浩谢新民孙秀芬等编译出版发行中国水利水电出版社北京市海淀区玉渊潭南路1号D座100038)网址:www.waterpub.com.cnE-mailsales@waterpub.com.cn电话:(010)68367658(营销中心)经售北京科水图书销售中心(零售)电话:(010)88383994,6320264全国各地新华书店和相关出版物销售网点中国水利水电出版社做机排版中心田刷北京市兴怀印刷厂规格184mm×260m16开本32.25印张765千字版次2009年11月第1版2009年11月第I次印刷定价9800元凡购买我杜图书,如有缺页、倒页、胱页的,本社营铺中心负责调撓版权所有·侵权必究GAMS用户指南〉编译校人员名单主要编译校人员魏传江王浩谢新民孙秀芬参加编译校人员(按章节顺序)杨舒媛俞烜郝春沣苟思陈宁叶勇李云鶴郭海丹高辉孙笑微王成丽王志璋柴福鑫张海涛李玉科杨丽丽前言数学世界是对物理世界的抽象。物理世界是个大千世界,包括声、光热、电、机槭等。尽管物理规律各异、物理现象千差万别,但是物理世界有三类最基本的关系:一是物理量之间的逻輯关系;二是物理量的量值;三是物理量的量值的动态依存关系。GAMS( General Algebraic Modeling System)即通用数学模型系统,是世界银行与美国GAMS公司在20世纪90年代初开发的一种旨在建立和求解大型复杂数学规划问题的高级计算机软件。作为GAMS语言则对应着这些物理量、这些最基本的关系、这些数学描述以及物理量的量值之间的逻辑关系AMS语言提出了集合的描述工具,用集合来表示物理量之间的逻辑关系。对于物理量的量值,GAMS语言分已知值和未知值进行处理,已知值用参数来表示,未知值用变量来表示。对于物理量之间的动态依存关系,GAMS语言用方程来表示。在这三种基本的表示下,首先定义物理量之间的逻辑关系,在逻辑关系的基础上定义物理量的量值或者变量的上下界,在集合和变量的基础上定义方程,在方程的基础上定义目标函数和约東条件,在目标函数和约束条件的基础上定义模型,在模型的基础上定义模型库,在模型库的基础上定义决策攴持系统,这样就完成了整个数学描述的任务。因此,GAMS软件不仅是一个数学工具,更是一个逻辑分析工具,对于理清思雏、进行系统分析,都是非常得力的工具与其他计算机语言相比,GAMS语言更加面向分析解决综合问题的计算机用户。由于它能够用一种用户和计算机都易读懂的语言来描述现实世界中精确的数学问题,因而可以灵活有效地建立各种类型的模型,进行多种数学优化问题的计算,极大地扩展了数学规划在策略研究和决策分析领域里的应用。用户只需要使用简单的GAMS语句,便可建立各种线性规划、非线性规划、混合整数规划、混合整数非线性规划、二次约束规划等问题的数学模型然后由GAMS系统运行求解。GAMS是专门为线性、非线性、混合整数等优化问题而设计的,在大型复杂问题中表现得更为突出。GAMS语言在形式上和通常使用的程序语言相似,饪何有编程经验的人都会觉得热悉且易上手。它具有很强的灵活性和通用性,可以简便安全地改变模型中的细节,允许模型的描述不依赖解题的算法,因为GAMS的设计把关系数据库的思想与数学规划巧妙地结合起来,并且将这个思想融入模型的总体设计,问題的最优解可以表示为不依赖于其使用的数据形式。GAMS软件也是建立和求解水资源配置模型的优秀软件包,通过二次开发可构建水资源配置模型和基于运筹学原理的其他水利模型,这为开发拥有自主知识产权的模型软件提供了技术平台。GAMS软件是基础软件,使用周期长、应用范围广、灵活方便是其主要优点。我国水资源严重短缺,水资源供需矛盾突出,与水有关的生态环境问题严重,长期以来非常重视水资源的开发利用与保护,尤其是在水资源配置理论研究和应用方面保持了一定的优势,取得了一系列重要成果。开发基于优化技术的水资源配置模型,首先要对水资源系统、经济系统、生态系统进行概化,建立它们之间的逻辑关系,然后将水文长系列数据、各类需求以及各种参数作为模型的輪入项。另外,水资源配置方案的求解是多目标问題,需要一个能够较好地解决上述计算问题的数学规划软件包作为支持。GAMS软件正是解决水资源问题的理想软件工具之一。GAMS语言没有专门的编輯器,用户可以自选编輯器。任何一个GAMS程序都由集合、数据、模型、求解和报告5个部分构成。GAMS不断地在原来基础上进行更新和提高,从早期的DOS版本已经发展至 Windows版本Windows GAMS2.50是为 Windows平台(95/98/ME和NT/2K/XP)设计的基于综合开发环境的界面,它可以在个人计算机、工作站、大型计算机上运行。本书分为两篇。第一篇为 Windows GAMS2.50用户指南,第1章是概述;第2章是一个独立的指南,通过一个小的运输模型例子进行详细的指导;第3~17章构成了第一部分的主干,有序地介绍了GAMS语言的各个组成部分,并伴随着一些详细的例子,即那些经常从模型库调出的例子;附录A~K是对高级主题的专门讨论第二篇讨论几个常用的求解器,即 BARON、 Cplex10、 DICOPT、MINOS、 MOSEK、PATH4.6、SBB、 MPSGE求解器。主要描述的选项为:用于控制求解器;GAMS求解器如何分析这些选项;如何分析求解器返回的模型和求解器状态代码。 BARON(分枝减小最优化导航)是整体求解非线性(NLP)和混合整数非线性规划( MINLP)的GAMS求解程序; Cplex优化求解器的设计理念是在最少的用户介入下快速求解大型的复杂问题,针对线性次约束混合整教规划问题,提供访问(合法授权的) Cplex的求解算法DICOPT是一个求解混合整数非线性规划( MINLP)问题的程序,涉及到线性二进制变量、整数变量和线性、非线性连续变量; MINOS是一个通用的非线性规划求解程序; MOSEK是求解线性、混合整数线性、二次、混合整数二次、二次约柬和凸非线性数学优化问题的软件包;PATH求解非线性问题,通过互补性把组合方法增加到经典的非线性方程组,能扩大建模的适用范围;SBB是求解混合整数非线性规划( MINLP)模型的新的GAMS求解器,它基于著名的混合整数线性规划标准的分枝定界法(B8B)和一些GAMS支持的标准NLP求解器的组合; MPSGE是GAMS数学规划语言的子系统,它是一个函数库和便于AGE模型表达和分析的雅可比计算程序本书由水利部“948”计划项目“通用数学模型系统( WindOws GAMs2.50)”(合同编号200701)、国家科技攴撑计划项目“流域/区域水资源全要素优化配置关键技术研究”(2007BAB28B02)、国家水体污染控制与治理科技重大专项“辽河流域水质水量优化调配技术及示范研究课题”(2008ZX07208010)资助。本书的顺利完成与项目组翻译人员的共同努力是分不开的,在此对参加翻译的所有工作人员表示真诫的感谢由于译者水平有限,难免会有疏忽或错误,在此恳请诸位专家、读者批评指正。译者2009年9月于北京目录第一篇 Windows GAMS2.50用户指南1概述1.1开发研究GAMS的动因1.2GAMS的基本特征1.2.1一般原理1.2.2文件33334412.3可移植性…1,2.4用户界面…1.2.5模型库1.3本书的组成2GAMS指南……2.1简介62.2GAMS模型的结构……集合2.4数据2.4.1列表数据输入890012.4.2表格数据输入2.4,3直接赋值数据输入2.5变量2.6方程2.6.1方程声明2.6.2GAMS求和(和求积)符号2.6.3方程定义7目标函数2.8模型和求解语句…2.9显示语句……………2.10‘.lo,,l,,up,,m’数据库2.10.1变量的边界和初值的赋值2.10.2最优值的转换和显示2.11GAMS输出2.11.1返回输出…8gg22.11.2错误信息2.11.3引用映射2.11.4方程列表2.11.5模型统计2.11.6状态报告…2.11.7求解报告3GAMS程序3.1简介3.2GAMS程序的结构2.1GAMS输入的格式32.2GAMS语句的分类3.2.3GAMS程序的组织3.3数据类型和定义3.4语言条目9993033233.4.1字符3.4.2保留字3.4.3标识符3.4.4标签3.4.5文本3.4.6数字33353.4.8注释集合定义4.1简介4.2简单集合4.2.1语法4.2.2集合名称4.2.3集合元素4.2.4关联的文本4.2.5集合元素的顺序33377883304.2.6多重集合的声明4.3 alias语句:集合的多重命名4.4子集和范围检查4.5多维集合……4.5.1一对一映射4.5.2多对多映射………4.6小结5数据输入:参数、标量和表格5.1简介5.2标量5.22实例5.3参数5.3.1语法5.3.2实例5.3.3更高维数的参数数据…5.4表格5,4.1语法5.4.2实例5.4.3连续表格5.4.4二维以上的表格5.4,5压缩表格54.6处理长行标签……495.5缩写5.1语法5,5.2实例506带参数的数据处理1简介6.2赋值语句6.2.1 scalar赋值6.2.2带索引赋值…6.2.3显式标签赋值…6.2.4子集赋值6.2.5控制索引问题526,2.6赋值中的扩展范围标识符6.2.7赋值中的缩写6.3表达式6,3.1标准的算术操作符
- 2020-12-04下载
- 积分：1
STM32移植到GD32的注意事项
GD32的性价比和稳定性都很高。我更喜欢用GD32.。但不是100%兼容的。需要注意一些细节。2.内部 Flash1)芯片设置读保护用法描述由于GD的Fash是自己的专利技术,STM的Fash是第三方提供的,所以GD的 Flash和STM的Fash有些许差异。GD的擦除时间会长一点解决方法在写完KEY序列以后,需要读该位,确认key已生效。所以,这里应该插入While(! (FLASH->CR &0x200);//Wait OPTWRE或可简单插入两个NOPNOPONOPO在ST库中,只有FLASH Status FLaSH Erase Option Bytes(void)FLaSH Status FlaSH ProgramOption Byte Data(uint32 t Address, uint8 t DataFLASH Status FLASH_ Enable WriteProtectionuint32 t FLASH PagesFLASH Status FLASH ReadOutProtection (Functional State New State)四个函数需要修改。2)IAP在应用中编程描述GD32由于有fash的0访问时序,同SM32在Fash的Ease和 Program上存在差别,GD32的 Erase和 Program时间比STM32的稍微长些建议对 Erase和 Progran时间进行修改解决方法将宏定义井 define erase timeout(luint32 t)OX000B000O)#define ProgramTimeout(uint32_t)ox00002000修改为:#define erasetimeout((uint32_t)OX00OFFFFF#define ProgramTimeout(uint32 t)OxOOOOFFFF备注: Erase和 Program时间宏定义在stm3210 x flash. c源文件中路径: braiesSTM32F10 x Std Periph Driversr)3)用IAR下载配置解决方法在批量牛产的时候首先会烧写一个USB的boot,这个boot自动运行后在由上位机软件进行烧写应用程序。如果boot程序不能自动运行则需要重新插拔次电源。给生产造成一些麻烦。产生不能自动运行程序的原因是如果程序设置读保护的话需要等待 FLASH CR的第9[ OPTWRE]位为1.如果没有置位的话继续执行就会出错。由」ST的执行速度慢,程序执行到读 FLASH CR寄存器的时候该位已经置1,GD的执行速度比较快,程序运行到这的时候该位还没置1,因此需要在 FLASH ReadOut Protection函数里面添加一些轮询该位为1或者加延时3.ISP烧写软件1)ISP烧写,建议使用官方烧写软件性述芯片内部同有区别解决方法建以到下载最新版本的另外也有专门的烧写软件(可以到论坛下毂如果使用自制的软件或脱机编程器,实现和完全兼容,建议修改以下参数。页擦除等待超时时间增加至整片擦除等待超时时间增加至左右字编程等待超时时间增加至,臾编程等笭超时时间增加至I/0日1)I0口外部中断使用方法措述在关闭期间,如果外部引脚有电平的变化,在使用MR打开中断后会马上进入中断服务程序。理论是打开中断前,不管管脚是否有电平的变化,都不会影响到打开后的中断响应。解决方法所以解决方法就是通过禁用上升沿或者下降沿检测寄存器来开关中断,不能使用IMR屏蔽奇存器。程序如下关闭沿检测,以达到关闭中断的目的,下降沿使用寄存器,上升沿使用寄存器2)在待机模式,PA8引脚特殊设置描述在使用低功耗的情况下,PA8会被MCU在内部被设置为地PA8复用为MU内部频率输出,超低功耗设时需要悬空解决方法在待机模式,PA8悬空不用3)低功耗下必须注意描述在使用低功耗情况下,把软件全部端口(AF)时钟关掉,无论是否有该端凵。4)当有脉冲群冲击管脚摧述需要在在进入中断后关闭中断4.定时器1)定时器输入捕获模式需要软件清中断描述sTM定吋时器输入捕获模式默认能硬件清中断,GD为了更加严格要求配置,需要做软件清中断解决方法软件清除标志位2)定时器向上脉冲计数模式设置述定时器的用法差异解决方法脉冲计数模式下,装载值必须设置为比预期值大,否则不计数在ST上如果重载值不设置(初始为0)的时候,CNT可以正常计数。在GD上如果重载值不设置保持初始为0的时候,会因为重载值为零,即便是来一个脉冲也会导致所有的寄存器复位从而不能正常计数。型号GD32F1系列MCUF|ah256B8及以上的型号)3)TIM、ADC模块描述Tmer、ADC模块的触发信号宽度要求解决方法|由于内部有高速和低速两条外围总线,Tmer、ADC模块和其他外设共同使用这两个总线。GD32F103/101系列Fash128KB及以下的型号, Timer、ADC等模块识別触发信号的条件是触发信号宽度大于模块所在总线的时钟宽度5.串口 USART1) USART连续发送数据字节有空闲位带述字节间有空闲位解决方法|对于一般的通讯米说,不会有影响,只对于一般在通讯上有特殊协议的,才会产生数据不准确的情况所以,特定情况,修改程序6.I2C总线1)硬件L2C特殊配置述GD的C相对STM的来说要少一个标志位解决、宏地址定义改交方法2、硬件I2C在会在向从机发送7bits地址完成后,从札还没来得及识别。(看客户应用)我们可以在发送完7bits后加个延时,让从机完全识别I2C Send7bitAddress(I2C1, EEPROM ADDRESS, I2C Direction TransmitterintOfffhile(i --3、检测ADDR不能使用I2 C Checkevent函数,因为他会清除ADR,可以使用I2 C Get FlagStatus函数就是把while(! I2C CheckEvent(I2C1I2C EVENT MASTER TRANSMITTER MODE SELECTED))Ey while(! I2C GetFlagStatus(I2C1, I2C FLAG ADDR))4、还有个关于编程步骤的严谨性,跟STW想比,我们是先 Clcar ack,再 Clear Arrd。7.ADC采集1)ADC采样设置述ADC启动解决方法|分三个方面时写入后,需要等待一段时间,如果用库的话就在 ADC CMD后面加20us左右的延付如果采用中断获得采样数据后,需要软件清除中断。8. SDio1) SDIO DAT3pin的在1 bit bus mode和4 bit bus mode下的配置摧述1、SD|O在1 bit bus mode下,DAT3pin是低电平,这样会导致 SD Card进入SP!模式。原因:初始化失败的原因主要是因为GD32的芯片SDO的DAT3∏存在BUG2、在4位模式下,通过上面的方法,程序能止常初始化,但不能正常读写SD卡原因:因为DAT3∏在前面已经配置成推挽输出,所以在4位模式下,不危正常读下。在调用4位模式前,把DAT3的端凵配置成复用推挽输入即可解决问题解决方法「1、1 bit bus mode的解决方法:建议在使能之前,先把配置成推挽输出,)且要置成高电平,使保持高电平即可2、4 bit bus mode的解决方法:在调用4位模式前,把DAT3的端口配置成复用输出即可解决问题。2)程序在刚烧完后能正常读写SD卡,断电再上电后,SD卡初始化失败,需要手动复位一次后才正常描述在某些SD卡中,GD32断电再上电,会引起SD卡上的时钟信号不正常,导致SD卡发送命令失败。解决方法在程序中,打开时钟后,增加一小段延时,以保证下时钟信号稳定。这个延时添加的地方:在即的配置文件中,然后在这个函数中找到就在这个后面加个延时。10. USBA, USB OTG1)客户使用的原工程时需要注意几点解|1、在中,增加如下图红色字体语句for (1=0; 18; i++) EPli= GetEndPoiNT(i)for(i=0:iregs. HC REGS [num]->HCCHAR, hcchar d3 2)pdev->host hc Status =HC NAK而V2.1.0版本的NAK处理过程如下else if (hcint b nak)if(hcchar b. eptype = EP TYPE_ INTR)UNMASK HOST INT CHH(num)USB OTG HC Halt(pdev, numelse if ((hcchar. b. eptype = EP TYPE CTRL)(hcchar b eptype = EP_ TYPE BULK))A re-activate the channel *hcchar, b chen =1hcchar b chris =0USB OTG WRITE REG32(&pdev->regs. HC REGS [num]->HCCHAR, hcchar d32)pdey->host HC Status [num]=HC_NAKCLEAR HC INT(hcreg, nak)唯一的区别就是 CLEAR HC INT( here,nak)的位置,在Ⅵ1.0.0版本中对于CTRL和BUK端点的NAK中断没有清除NAK,我们的芯片会因此产生多次IN传输的请求,导致数据传输错误。改为V2.1.1的写法后传输正常。(注意 HC Status在V2.1.0是数组,在Ⅵ1.0.0是单个数据,直接拷贝的话要去掉后面的[num])B.USB外设的工作频率有限制摧述有最低工作频率的要求,也就是APB1分频后的时钟必须大于12MHz,比如HCLK为56MHz,APB1的最大分频系数为4,56/4=14MHz,可以正常工作。11 SPI1)输入与输出配置要求(STM32不需要如此要求)解决丨GD32在使用SP时,o的配置必须严格遵守主从模式下的输入与输出配置,而方法STM32无此要求,相关代码如下主机模式下|o配置(主机以SP为例):GPIO InitStructure gPio Mode gPio Mode af plGPIO_ Init Structure GPIo Speed GPlO Speed 50MHzGPIO InitStructure. GPio Pin= GPlO Pin 5 GPIO Pin_ 7;GPIO Init(GPIOA, &GPIO InitstructureGPio Init Structure gPio Mode gPio Mode IN floating:GPio InitStructure gpio Pin gpio pin 6GPIO Init(GPIOA, &GPIO InitStructure)从机模式下o配置(从机以SP2为例)GPIo Init Structure GPio Mode gPlo Mode IN FloatingGPIO InitStructure GPIO Speed= GPl Speed 50MHzGPIO_InitStructure GPIO_ Pin GPIO Pin_13 GPIO_ Pin_15GPIO Init(GPIOB, &GPIO InitStructure)gPio Initstructure gpio mode gpio mode af pp.GPIO InitStructure. GPio Pin= GPIo Pin 14:GPIO_Init(GPIOB, &GPIO_Initstructure);3)在GD32的SP的时钟信号,空闲状态需要配置成高电平,以保证数据的稳定性,具体代码如下:红色字体代码解决SPI_ InitStructure SPl Direction =SPI_ Direction_ 2Lines fullDuplex;方法SPI InitStructure SPl Mode SPi Mode master.SPI Initstructure SPl Data Size= SPl Data Size 8bSPlInitStructure SPl_CPOL= SPI CPOL HighSPl Initstructure SPl CPHA= SPI CPHA 2EdgeSPI InitStructure SPI NSS= SPI NSS SoftSPI InitStructure SPl BaudRate Prescaler =SPI Baud Rate Prescaler 256:SPI Initstructure, Spi FirstBit= SPI FirstBit MSBSPI InitStructure SPl CRCPolynomial =7;SPl Init(sPl1, &SPl Initstructure);4)当作为从机时,在GD32中,时钟信号必须为8的整数倍。例如:红色字体代码解决SPI_InitStructure SPl_ Direction =SPl_ Direction_ 2Lines_ Full Duplex;方法SPI InitStructure. SPl Mode= SPl Mode MasterSPI InitStructure SPSPi Data Size 8SPl InitStructure SPl CPOL= SPI CPOL High;SPI Init Structure. SPl CPHA SPI CPHA_ 2EdgeSPI Initstructure SPl NSS= SPI NSs SoftSPI Initstructure Spl BaudRatePrescaler= SPl BaudRatePrescaler 256SPI InitStructure. SPl FirstBit= SPl First Bit MSBSPl Initstructure SPl CRCPolynomial =7;SPI Init(SPI1, &SPl InitStructure)5)在GD32中,不能使用SPLS_FLAG_BSY该位来判断SP总线数据是否接收或发送完成12.看门狗1)进入SToP模式前打开看门狗,通过RTc的ALR唤醒后,程序会不断被复位的现像摧述WDG内部有个 Reload信号,KEY奇器写AAA会使其拉高,过段时间自动拉低。在拉底之前进入STOP状态会使 Reload信号一直为高,等到退出SToP后也保持为高,之后再写AAAA没有办法让 Reload产生上升沿,也就没办法更新计数器了解决方法「进STOP之前不要 Reload,也可以调整下程序的顺序,把WwDG的配置放到RIC配置之前,效果是一样的。
- 2021-05-06下载
- 积分：1
|
联系客服（QQ：3324461878）
