资源描述：
《SASIML模块与矩阵运算.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、SAS/IML模块与矩阵运算一、SAS/IML模块1．SAS/IML模块简介IML是InteractiveMatrixLanguage的缩写，即交互式矩阵语言，也是一种程序设计语言。它提供了程序设计语言所需的条件转移和循环等一系列完整的程序流控制语句。SAS/IML模块提供了一个动态的交互式环境，在此环境下，用户按照IML语法输入语句，然后系统交互的输出计算结果。SAS/IML模块还提供了用户自定义函数和子程序（过程）的功能。IML语句、函数、子程序、及用户自定义函数和子程序按照一定的规则组成IML程
2、序。IML程序可存放在外部文件中，可随时递交给系统连续执行。2.IML变量所有IML变量都为矩阵类型变量，最小为1´1阶矩阵。矩阵元素可为数值型，也可为字符型，缺失值用点号“.”表示。多维矩阵每一行的元素个数（列数）必须相同。1)数值型矩阵中的元素可以是整数、小数，也可以用科学计数法表示，如1E-5。所有数值均以双精度格式存储。2)字符型矩阵中，若元素中含有空格或某些特殊字符，如“?”、“=”、“*”、“:”、“(”、“)”等，则须将字符串用单引号“'”或双引号“"”括起来。在同一个字符型矩阵中，每个
3、元素的字符串长度不必相等，但存储的长度等于最大字符串的长度，此长度不超过32676。3)IML变量的取名规则：(1)可由下划线、字母和数字组成，但不能以数字开头。(2)不可超过8（32）个字符。(3)不区分大小写。3.启动和退出IML模块1)启动IML模块PROCIML;2)退出IML模块QUIT;4.设置和关闭IML交互模式1)设置IML交互模式RESETPRINT;2)关闭IML交互模式RESETNOPRINT;二、矩阵变量的建立矩阵变量可通过定义，赋值，或函数返回来建立。1．通过矩阵定义建立矩阵
4、变量1)定义矩阵的语法<矩阵变量名>＝<矩阵初值>；示例(1)1×1阶矩阵（数）a=12;a=.;a='hithere';a="Hello";(2)行向量x={123456};(3)字符型行向量a={abcdefg};printa;结果：AABCDEFG(4)字符型行向量a={'abc''DEFG'};printa;结果：AabcDEFG(5)列向量y={1,2,3,4,5};(6)3×2阶矩阵z={12,34,56};2)重复因子的使用重复因子符号用表示[]。例如：answer={[2]'Yes',
5、[2]'No'};等价于answer={'Yes''Yes','No''No'};3)重新定义或赋新值a={123,654};a={'Sales''Marketing''Administration'};2.通过赋值建立矩阵变量1)语法<矩阵变量名>＝<表达式>;2)示例x={1,2,3};y={3,2,1};z=x+y;3.通过函数返回建立矩阵变量1)语法<矩阵变量名>＝函数名（<参数表>）;2)示例x=SQRT({1,2,3});三、矩阵的运算及其IML运算符1．矩阵的算术运算1)矩阵的算术运算运
6、算符运算符含义+加-减*乘/除**乘方#对应元素相乘##对应元素乘方-各元素取相反数2)加(1)语法矩阵1+矩阵2矩阵+数(2)功能(3)示例3)减(1)语法矩阵1-矩阵2矩阵-数(2)功能(3)示例4)乘(1)语法矩阵1*矩阵2矩阵*数(2)功能(3)示例5)除(1)语法矩阵1/矩阵2矩阵/数(2)功能对应元素相除，或各元素除以给定的数，结果为同阶矩阵。(3)示例6)乘方(1)语法矩阵**数(2)功能矩阵乘方。其中给定的数为大于等于–1的整数，若非整数将自动取整。A**(-1)等于A的逆矩阵。(3)
7、示例7)对应元素相乘(1)语法矩阵1#矩阵2矩阵#数矩阵#向量(2)功能(3)示例a={12,34};b={48,05};c=a#b;矩阵C中的内容为：C2rows2cols(numeric)416020d={10,100};ad=a#d;矩阵AD中的内容为：AD2rows2cols(numeric)1020300400若改为：d={10100};ad=a#d;则矩阵AD中的内容为：AD2rows2cols(numeric)10200304008)对应元素乘方(1)语法矩阵1##矩阵2矩阵##数(2)
8、功能若矩阵1中的元素有负数，则矩阵2中相应的元素必须为整数。(3)示例例1:a={123};b=a##3;则矩阵B中的内容为：B1row2cols(numeric)1827例2:b=a##.5;则矩阵B中的内容为：B1row3cols(numeric)11.41421361.7320508例3:a={-176,20-8};b=-a;则矩阵B中的内容为：B2rows3cols(numeric)1-7-6-2089)元素取相反数(1)语法–矩阵(2)功能(3
