初學R語言，嘗試着寫了一個函數，調用這個函數，生成一些曲線。
###############################################################################
PictureLine <- function(name.picture,width.picture=1400,height.picture=1400,
width.line=1,x.limit,y.limit,lwd.grid=2,
x.legend,y.legend,cex.legend=3,cex.point=2,
number.grid.major.x,number.grid.minor.x,
number.grid.major.y,number.grid.minor.y,colum.x,colum.y,
x.axis,new.x.axis,y.axis,new.y.axis,picture.legend){
jpeg(name.picture,width=width.picture,height=height.picture*(number.grid.major.y/number.grid.major.x))
oldpar<-par(no.readonly=TRUE)
par(pin=c(width.picture/100,
height.picture/100*(number.grid.major.y/number.grid.major.x)),
lwd=width.line,cex.axis=2.5)
plot( mydate_z20[,colum.x],mydate_z20[,colum.y],type<-"p",pch=0,cex=cex.point,
yaxt="no",xaxt="no",xlim=x.limit,ylim=y.limit,xaxs="i",yaxs="i",
ann=FALSE)
xspline(mydate_z20[,colum.x],mydate_z20[,colum.y],0.2)
lines( mydate_z10[,colum.x],mydate_z10[,colum.y],type<-"p",pch=1,cex=cex.point)
xspline(mydate_z10[,colum.x],mydate_z10[,colum.y],0.2)
lines( mydate_z0[, colum.x], mydate_z0[,colum.y], type<-"p",pch=2,cex=cex.point)
xspline(mydate_z0[, colum.x], mydate_z0[,colum.y], 0.2)
lines( mydate_f10[,colum.x],mydate_f10[,colum.y],type<-"p",pch=5,cex=cex.point)
xspline(mydate_f10[,colum.x],mydate_f10[,colum.y],0.2)
lines( mydate_f20[,colum.x],mydate_f20[,colum.y],type<-"p",pch=6,cex=cex.point)
xspline(mydate_f20[,colum.x],mydate_f20[,colum.y],0.2)
lines( mydate_f30[,colum.x],mydate_f30[,colum.y],type<-"p",pch=8,cex=cex.point)
xspline(mydate_f30[,colum.x],mydate_f30[,colum.y],0.2)
axis(1,at=x.axis,labels=new.x.axis,lty=1,las=1,tck=0.0,font.axis=3,lwd=lwd.grid,padj=1)
axis(2,at=y.axis,labels=new.y.axis,lty=1,las=1,tck=0.0,font.axis=3,lwd=lwd.grid,padj=1)
axis(3,at=x.axis,labels=FALSE, tck=0.0, lwd=lwd.grid)
axis(4,at=y.axis,labels=FALSE, tck=0.0, lwd=lwd.grid)
grid(number.grid.major.x,number.grid.major.y,lty=1,lwd=2,col="black",equilogs=FALSE)
grid(number.grid.minor.x,number.grid.minor.y,lty=3,lwd=1,col="black",equilogs=FALSE)
legend(x.legend,y.legend,inset=0.5,picture.legend,pch=c(0,1,2,5,6,8),bg="white",cex=cex.legend)
par(oldpar)
dev.off()
}
floder.dir <- ""
file.f30 <- paste(floder.dir,"dzf30.dat",sep ="")
file.f20 <- paste(floder.dir,"dzf20.dat",sep ="")
file.f10 <- paste(floder.dir,"dzf10.dat",sep ="")
file.z00 <- paste(floder.dir,"dzz00.dat",sep ="")
file.z10 <- paste(floder.dir,"dzz10.dat",sep ="")
file.z20 <- paste(floder.dir,"dzz20.dat",sep ="")
mydate_f30 <- read.table(file.f30)
mydate_f20 <- read.table(file.f20)
mydate_f10 <- read.table(file.f10)
mydate_z0 <- read.table(file.z00)
mydate_z10 <- read.table(file.z10)
mydate_z20 <- read.table(file.z20)
alpha <- c(-10,-5,0,5,10,15,20,25)
nalpha <- expression(-10,-5,0,5,10,15,20,alpha)
cy<- c(-0.5,0,0.5,1.0,1.5,2.0)
ncy <- expression(-0.5,0,0.5,1.0,1.5,C[y])
cleg <-expression(delta[z]==20^degree,delta[z]==10^degree,delta[z]==0^degree,
delta[z]==-10^degree,delta[z]==-20^degree,delta[z]==-30^degree)
PictureLine(name.picture="cy-a.jpg",
x.limit=c(-10,25),
y.limit=c(-0.5,2.0),
x.legend=10,
y.legend=1.0,
cex.legend=2.5 ,
number.grid.major.x=7,
number.grid.minor.x=35,
number.grid.major.y=5,
number.grid.minor.y=25,
colum.x=8,
colum.y=1,
x.axis=alpha,
new.x.axis=nalpha,
y.axis=cy,
new.y.axis=ncy,
picture.legend=cleg)
生成的圖片如下圖所示：
這個代碼主要是用來繪製風洞實驗數據曲線，繪製不同升降舵偏度下的升力係數，力矩係數，升阻比隨迎角的變化等等。