材料学院一哥们找我帮忙做计算，给了数据、公式，和验证结果的excel文档，原理我不懂，不过对数据处理还算熟悉。

数据格式如下：

公式如下：

excel如下：

R代码如下：

后根据材料兄弟的要求又修改了输出格式，代码如下：

#R代码：
w<-read.table("standard.txt",header=TRUE)
#matplot(ww[,1],ww[,-1], type = "b", pch = 21:24, col = 2:5,bg = 2:5)
attach(w)
zrl=data.frame(array(1,dim=c(length(w[,1]),6)))
zrl[,1]<-w[,1]
colnames(zrl)<-c("Frequency.GHz","d1","d2","d3","d4","d5")
for(d in 1:5)
{
c<-3*(10^11)
μ<-μ.-μ..1*1i
ε<-ε.-ε..1*1i
z<-sqrt(μ/ε)*tanh(1i*(0.0205979*f*d)*sqrt(μ*ε))
rl<--20*log10(abs((z-1)/(z+1)))
zrl[,d+1]<-rl
}
detach("w")
write.table(zrl, file="zrl.txt", sep="ttt",col.names=T,row.names = F,quote=F,append=FALSE)