assign函数在循环时候，给变量赋值。
举例说明：
1、
for (i in 1:(length(rowSeq)-1)){
assign(paste("nginx_server_fields7_", i, sep = ""), nginx_server_fields7[(rowSeq[(i-1)+1]):(rowSeq[i+1]), ])
}
2、
for (i in 1:3){
assign(paste("a", i, sep = ""), i:10)
}
ls()
[1] "a1" "a2" "a3" "i"
>a1
[1] 1 2 3 4 5 6 7 8 9 10
>a2
[1] 2 3 4 5 6 7 8 9 10
>a3
[1] 3 4 5 6 7 8 9 10
1、paste函数可用于字符串连接
用法：paste (..., sep = " ", collapse = NULL)，分隔符默认为空格
我们简单举个例子
1）paste("a","b") ##能连接a b
[1] "a b"
paste("a","b","c")
[1] "a b c"
2）设置分隔符paste("a","b",sep="=") ##注意到用等号分隔了
[1] "a=b"
3）连接多个元素paste("a",1:5,sep="") ##会自动每个元素与a相连
[1] "a1""a2""a3""a4""a5"
4）paste("a",1:5,".pdf",sep="") ##比如想批量输出文件名
[1] "a1.pdf""a2.pdf""a3.pdf""a4.pdf""a5.pdf"、
2、collapse参数
谢益辉大佬说这个参数引无数英雄竞折腰啊，大家只会sep
1）paste("a",1,collapse="+") ##注意到了吧，collapse没用上
[1] "a 1"
2）paste(c("a","b","c"),collapse="+") ##这样就很容易理解了，collapse折叠起了元素
[1] "a+b+c"
3）paste("a",1:5,sep="") ##再来更明显的举例
[1] "a1""a2""a3""a4""a5"
4）paste("a",1:5,sep="",collapse="+") ##先在元素间连接，然后折叠
[1] "a1+a2+a3+a4+a5"
3、与paste0函数的区别
paste0函数，默认sep=""是两个函数唯一的区别
paste(c("a","b","c"),1:3)##默认空格符
[1] "a 1""b 2""c 3"
paste(c("a","b","c"),1:3,sep=" ")
[1] "a 1""b 2""c 3"
paste(c("a","b","c"),1:3,sep="")
[1] "a1""b2""c3"
如果两个向量长度不同paste(c("a","b","c"),1:5) ##超出范围后继续从前向后连接
[1] "a 1" "b 2" "c 3" "a 4" "b 5"
paste0(c("a","b","c"),1:3)#默认元素连接为sep=""
[1] "a1""b2""c3"
用法：unlist()函数的作用，就是将list结构的数据du，变zhi成非list的数据，即将list数据变成 字符串向量 或者数字向量的形式
如果是向量的话就直接输出向量
例子：
%in%相当于match()函数的一个缩写。用来判断一个数组或矩阵是否包含在另一个数组或矩阵里。举个例子一目了然：
#首先复制两个变量a和b
>a <- 1:5
>b <- 3:7
>a %in% b #看a的元素是否包含在b中输出结果如下：
[1] FALSE FALSE TRUE TRUE TRUE
用法：apply(x, MARGIN, FUN, ...)
作用：对 矩阵 的行或列使用函数，或者对 数组 的各个维度使用函数
x为数据对象，MARGIN是维度的下标，MARGIN=1表示行，MARGIN=2表示列，FUN是自己指定的任意或自定义函数
注：数据框dataframe也可以使用apply函数，该函数会自动将数据框转化为矩阵，但前提是数据框中各列的数据类型 必须是数值型，否则会报错。
用法：lapply(x,FUN,...)
作用：对列表中的各个元素使用函数
x是格式为列表的数据源，FUN是任意函数。
intersect(data1，data2)：交集
union(data1，data2)：并集
输入了函数对象名称，可以直接看到代码的，如要获得函数对象fivenum的代码，就只需要在Console中键入函数对象名称fivenum就可以得到如下结果：function (x, na.rm = TRUE)
{
xna <- is.na(x)
if (na.rm)
x <- x[!xna]
else if (any(xna))
return(rep.int(NA, 5))
x <- sort(x)
n <- length(x)
if (n == 0)
rep.int(NA, 5)
else {
n4 <- floor((n + 3)/2)/2
d <- c(1, n4, (n + 1)/2, n + 1 - n4, n)
0.5 * (x[floor(d)] + x[ceiling(d)])
}
}
题主的做法主要存在两点问题：
1、quad函数用于计算数值积分，函数表达式中不能包含符号量；
2、被积函数的表达式应该写成关于被积变量的向量化的形式（也就是应该用点运算）。
参考代码：
R=1syms L
rr = 0 : 0.1 : 1
for ii = 1 : length(rr)
r = rr(ii)
f = @(l)(acos((1+l*l-r*r)/(2*l))+r*r*acos((r*r+l*l-1)/(2*r*l))-0.5*sqrt(4*r*r-(1+r*r-l*l)^2))*2*l/(pi*r^4)
fun = @(L) arrayfun(f,L)
J(ii) = quadl(fun,0,r)
end
plot(rr, J)
或者也可以借用楼上 枫箫1 的部分代码，写成：
R=1syms L
rr = 0 : 0.1 : 1
for ii = 1 : length(rr)
r = rr(ii)
SOA=R^2*acos((R^2+L^2-r^2)/(2*R*L))+r^2*acos((r^2+L^2-R^2)/(2*r*L))-...
0.5*sqrt(4*R^2*r^2-(R^2+r^2-L^2)^2)
PAB=SOA/(pi*r^2)
p=2*L/r^2
f=PAB*p
fun = eval(['@(L)' vectorize(f)])
fun = @(l) arrayfun(@(L)eval(f),l)
J(ii) = quadl(fun,0,r)
end
plot(rr, J)
以上代码虽可以运行，但被积函数存在问题——SOA的第一项反余弦的值可能出现复数（因为在r稍小的情况下，acos的参数大于1），请题主再仔细检查一下。
