1. (1) 写出元素为3, -1.5, 3E-10的向量。
c(3, -1.5, 3e-10)
(2) 写出从3开始每次增加3,长度为100的向量。
seq(3,by=3,length=100)
(3) 写出(0, 2)重复10次的向量。
rep(c(0,2),10)
(4) 对向量x,写出其元素大于等于0小于1的条件。
if( all( x >= 0 && x < 1) )
(5) 对向量x,写出其元素都等于0的条件。
if( all(x == 0 ) )
(6) 写出包含12个月份名称的向量。
month.name[1:12]
2. 设x为一个长100的整数向量。比如,x <- floor(100*runif(100))。
floor(100*runif(100))
(1) 显示x第21到30号元素。
x[c(21:30)]
(2) 把x第31,35,39号元素赋值为0。
x[c(31,35,39)]=c(0,0,0)
(3) 显示x中除了第1号和第50号的元素之外的子集。
x[-c(1,50)]
(4) 列出x中个位数等于3的元素。
i=1
while(i<=100){
if(x[i]%%10!=3) {x[-c(i)] }
i=i+1
}
(5) 列出x中个位数等于3的元素的下标位置。Hint:use which( )
for( i in 1 : length(x) ){ if( x[ i ] %% 10 == 3 ){print(i)} }
3. 绘制y=x3 +sin(x)*cos(x)的函数图
x=1:100
y=x*x*x+sin(x)*cos(x)
plot(x,y,type='s')
#plot(x,y,xlab=’x轴’,ylab=’y轴’,main=’绘制y=x3 +sin(x)*cos(x)的函数图’,type=’s’)
4. 随机产生两个矩阵A,B,其维数都为5*6,矩阵中每一个元素服从标准正态分布。(1) 将A,B按行合并成一个矩阵C;(2)计算A’*B; (3) 计算矩阵C的行和及列和
AA=rnorm(30,0,1)
A=array(AA,c(5,6))
B=array(rnorm(30,0,1),c(5,6))
(1)C=rbind(A,B)
(2)t(A)%*%B
(3)
行:sum(C[1,]) #计算矩阵第一行和,其他行同理
列:sum(C[,1])
博客地址:http://blog.yoqi.me/?p=831
这篇文章还没有评论