- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
SAS/IML模块与矩阵运算
一、SAS/IML模块
1.SAS/IML模块简介
IML是Interactive Matrix Language的缩写,即交互式矩阵语言,也是一种程序设计语言。它提供了程序设计语言所需的条件转移和循环等一系列完整的程序流控制语句。SAS/IML模块提供了一个动态的交互式环境,在此环境下,用户按照IML语法输入语句,然后系统交互的输出计算结果。
SAS/IML模块还提供了用户自定义函数和子程序(过程)的功能。IML语句、函数、子程序、及用户自定义函数和子程序按照一定的规则组成IML程序。IML程序可存放在外部文件中,可随时递交给系统连续执行。
2. IML变量
所有IML变量都为矩阵类型变量,最小为1?1阶矩阵。矩阵元素可为数值型,也可为字符型,缺失值用点号“.”表示。多维矩阵每一行的元素个数(列数)必须相同。
1) 数值型矩阵中的元素可以是整数、小数,也可以用科学计数法表示,如1E-5。所有数值均以双精度格式存储。
2) 字符型矩阵中,若元素中含有空格或某些特殊字符,如“?”、“=”、“*”、“:”、“(”、“)”等,则须将字符串用单引号“”或双引号“”括起来。在同一个字符型矩阵中,每个元素的字符串长度不必相等,但存储的长度等于最大字符串的长度,此长度不超过32676。
3) IML变量的取名规则:
(1) 可由下划线、字母和数字组成,但不能以数字开头。
(2) 不可超过8(32)个字符。
(3) 不区分大小写。
3. 启动和退出IML模块
1) 启动IML模块
PROC IML;
2) 退出IML模块
QUIT;
4. 设置和关闭IML交互模式
1) 设置IML交互模式
RESET PRINT;
2) 关闭IML交互模式
RESET NOPRINT;
二、矩阵变量的建立
矩阵变量可通过定义,赋值,或函数返回来建立。
1.通过矩阵定义建立矩阵变量
1) 定义矩阵的语法
矩阵变量名 = 矩阵初值;
示例
(1) 1×1阶矩阵(数)
a=12;
a=. ;
a=hi there;
a=Hello;
(2) 行向量
x={1 2 3 4 5 6};
(3) 字符型行向量
a={abc defg};
print a;
结果:
A
ABC DEFG
(4) 字符型行向量
a={abc DEFG};
print a;
结果:
A
abc DEFG
(5) 列向量
y={1,2,3,4,5};
(6) 3 ×2阶矩阵
z={1 2, 3 4, 5 6};
2) 重复因子的使用
重复因子符号用表示[] 。例如:
answer={[2] Yes, [2] No};
等价于
answer={Yes Yes, No No};
3) 重新定义或赋新值
a={1 2 3, 6 5 4};
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) 矩阵的算术运算运算符
运算符 含义 + 加 - 减 * 乘 / 除 ** 乘方 # 对应元素相乘 ## 对应元素乘方 - 各元素取相反数
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) 对应元素相乘
(1) 语法
矩阵1#矩阵2
矩阵#
矩阵#(2) 功能
(3) 示例
a={1 2,
3
文档评论(0)