R语言求拟合函数图形的交点坐标

最新推荐文章于 2023-08-28 00:48:11 发布

SCIEISSCI

最新推荐文章于 2023-08-28 00:48:11 发布

阅读量1.9k

点赞数

文章标签： r语言

本文链接：https://blog.youkuaiyun.com/SCIEISSCI/article/details/117082789

版权

该博客介绍了一种使用R语言通过二分类无限逼近方法求解多个拟合函数图形交点坐标的技巧。代码中首先创建了四个曲线的数据结构，然后应用`approxfun`函数进行拟合，并利用`optim`函数结合`combn`进行交点的逼近求解。最终找到交点坐标，确保结果精度达到1e-10。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

R语言求拟合函数图形的交点坐标

图中的交点坐标

通过二分类无限逼近的方法获得的

代码如下：
1）四个曲线的数据结构
v_ac=log(vir_1)
op_log=log(opp)
test1=data.frame(v_ac1,op_log)
v_ac=log(vir_2)
op_log=log(opp)
test2=data.frame(v_ac2,op_log)
v_ac=log(vir_3)
op_log=log(opp)
test3=data.frame(v_ac3,op_log)
v_ac=log(vir_4)
op_log=log(opp)
test4=data.frame(v_ac4,op_log)
2）开始求解：

options(digits=20)
funcLs <- lapply(list(test1, test2, test3,test4), function(DF) approxfun(DF $op_log, DF$ v_ac))
X <- test1$op_log

Filter(Negate(is.null), unlist(combn(funcLs, 2L, FUN=function(funcs) {
mapply(function(lower, upper) {
res <- optim((lower+upper)/2, function(x) abs(funcs[1L] - funcs[2L]),
method=“Brent”, lower=lower, upper=upper, control=list(reltol=1e-10))
if (res $KaTeX parse error: Expected 'EOF', got '&' at position 14: value < 1e-4 &̲& res$ convergence==0L) return(res$par)
}, X[-length(X)], X[-1L])
}, simplify=FALSE), recursive=FALSE))