猫腻哪本书好看中有dfsane的用法

1797人阅读
R语言(21)
其他(10)
机器学习(11)
首先,dfsane属于BB包,使用时应首先载入library(&BB&)
nleqslv 属于nleqslve包,同样使用时应载入 library(&nleqslve&)
Description
Derivative-Free Spectral Approach for solving nonlinear systems of equations
求解非线性方程组法
dfsane(par, fn, method=2, control=list(),
quiet=FALSE, alertConvergence=TRUE, ...)
看这个公式很多人都头疼,不知道你们是不是,反正我是,那么就那代码样例说话:
library(BB)
fun &- function(x) {
f &- numeric(length(x))
x[1]/(1+exp(x[2]*0.215+x[3]))-134.03
x[1]/(1+exp(x[2]*5.06+x[3]))-80.2
x[1]/(1+exp(x[2]*9.85+x[3]))-13.4
startx &- c(132,0.5,0.5)
result = dfsane(startx,fun,control=list(maxit=2500,trace = FALSE))
theta = result$par
fun函数就是 fn,公式中第二个,你的非线性方程组的公式比如,在本例中 y = a/(1+exp(bx+c)
在本方程中x,y为已知坐标点,求a,b,c的值
在fun中,参数只为x,此x并非自变量x,而是一个list,也就是向量即你要求的参数值(a,b,c的值)的表示,也即用dfsane的求得的解,x[1]表示a,x[2]表示b,x[3]表示c。
由于0 = &a/(1+exp(bx+c) - y, 所以即求 f&= &a/(1+exp(bx+c) - y由于要求a,b,c,在原方程中x,y已知,带入。
求解三个参数最少需要三个方程组联立,f[1],f[2],f[3]即为三个数组联立,用dfsane求解,求得的解即为要求a,b,c的值
startx即为参数par,公式中第一个,即为你要求解公式使用方法的初始值,因为求解公式本质是用搜索方法求最优解,比如高斯-牛顿法,所以得有起始点(即初始值),可以设置为要求解的估计值或猜测一个
control=list(maxit=2500,trace = FALSE),maxit搜索的迭代次数,一般100,200,trace默认值为true,把每一步迭代的结果都给列出事,事实上,我们用不到知道每一步的解,只要一个最终结果就行了,所以设置为false, control 里一般使用这两个参数
method说的使用那种搜索方法,除非题目要求,一般没差,可以不写,其他的就更可以不写了 ,啊哈哈哈
运行结果:
par就是要求的a,b,c的解,即为向量(list)x的值,可以用如下获取:
result = dfsane(startx,fun,control=list(maxit=2500,trace = FALSE))
theta = result$par
theta即为向量x的值
在此例中a = theta[1],b= theta[2],c=theta[3]
http://www.inside-r.org/packages/cran/bb/docs/dfsane&
再来一个例子
#y ~ a +(b*exp(c*x^d))
tempFunction1 &- function(modelParams){
F = numeric(length(modelParams))
F[1] = modelParams[1]+(modelParams[2]*exp(modelParams[3]*(exp(2.751514))^modelParams[4])) - 107.20561
F[2] = modelParams[1]+(modelParams[2]*exp(modelParams[3]*(exp(5.509109))^modelParams[4])) - 70.53738
F[3] = modelParams[1]+(modelParams[2]*exp(modelParams[3]*(exp(6.779219))^modelParams[4])) - 42.44393
F[4] = modelParams[1]+(modelParams[2]*exp(modelParams[3]*(exp(8.542089))^modelParams[4])) - 17.10280
theta = dfsane(par= c(100 + 50*runif(1),runif(1),-runif(1),runif(1)), fn=tempFunction1,control = list(maxit = 2500,trace = FALSE))
#那么a = theta$par[1] b = theta$par[2] c = theta$par[3] d= theta$par[4]
博主微博:文二十口
参考知识库
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
访问:22787次
排名:千里之外
原创:46篇
(2)(2)(2)(3)(1)(2)(5)(2)(1)(5)(5)(1)(1)(9)(4)(2)memberane的用法和样例:
The top of bucket is sealed with rubber membrane, and the rubber memberane inflated at a pressure.
水桶顶端用橡皮膜密封,橡皮膜在压力作用下澎胀。
The concentration of cellular memberane protein was measured by Bradford method and the final concentration was modulated above Img/ml.
用Bradford法测定细胞膜蛋白质浓度,膜蛋白最终浓度控制在1mg/ml以上。
Methods Balb /c mice were immunized with out memberane protein,specificity det ermined by ELISA and protertion of mAb against V.alginolyticus was evalunder in BALB /c mice.
方法用溶藻弧菌外膜蛋白免疫BALB/c小鼠,ELISA鉴定mAb特异性,用BALB/c小鼠进行感染保护实验。
Conclusion Optimizing parameters of the techniques in membrane process, using suitable pretreatment method and cleaning method can reduce membrance fouling and make the memberane property recover.
结论优化膜过程工艺参数,采用适当预处理方法和清洗方法可以减少膜的污染,使膜的性能有较大恢复。
Research Progress of Porous Metal Memberane
多孔金属膜研究进展
memberane的海词问答与网友补充:
memberane的相关资料:
相关词典网站:-ane是什么意思,词典释义与在线翻译:
形容词 ane:
not two or more
-ane的用法和样例:
On ne saurait faire boire un ane qui n'a pas soif.
驴不喝水不能强按头。
You may pick them up ane time you wish.
你们可以从它们中挑选一个你们愿意的时间!
At last ,I lost the chance ane I regre that .
实在想不到他有什么好处,那就想想他的坏处吧。
Fighters ane badly needed in our time.
在这个时代,战士是最需要的。
The Huang Pu Frontier Health ane Quarantine Service of P.R.C.
送达地点中华人民共和国黄埔卫生检疫所。
Ane has many books, most of which he has never read before.
安有许多藏书,大多数他都没读过。
-ane的海词问答与网友补充:
-ane的相关资料:
-ane&:&[词尾]饱和 ...
在&&中查看更多...
相关词典网站:3372人阅读
-----c、c++心得-----(55)
写在最前的三点:
1、所谓图的遍历就是按照某种次序访问图的每一顶点一次仅且一次。
2、实现bfs和dfs都需要解决的一个问题就是如何存储图。一般有两种方法:邻接矩阵和邻接表。这里为简单起
见,均采用邻接矩阵存储,说白了也就是二维数组。
3、本文章的小测试部分的测试实例是下图:
一、深度优先搜索遍历
1、从顶点v出发深度遍历图G的算法
② 依次从顶点v未被访问的邻接点出发深度遍历。
2、一点心得:dfs算法最大特色就在于其递归特性,使得算法代码简洁。但也由于递归使得算法难以理解,原因
在于递归使得初学者难以把握程序运行到何处了!一点建议就是先学好递归,把握函数调用是的种种。
3、算法代码:
二、广度优先搜索遍历
1、从顶点v出发遍历图G的算法买描述如下:
②假设最近一层的访问顶点依次为vi1,vi2,vi3...vik,则依次访问vi1,vi2,vi3...vik的未被访问的邻接点
③重复②知道没有未被访问的邻接点为止
2、一点心得:bfs算法其实就是一种层次遍历算法。从算法描述可以看到该算法要用到队列这一数据结构。我这
里用STL中的&queue&实现。该算法由于不是递归算法,所以程序流程是清晰的。
3、算法代码:
4、小测试:
参考知识库
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
访问:22361次
积分:1856
积分:1856
排名:第17957名
原创:158篇
转载:63篇
(5)(43)(69)(36)(4)(2)(7)(11)(44)

我要回帖

更多关于 高二数学学哪几本书 的文章

 

随机推荐