关于本站
人大经济论坛-经管之家：分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
经管之家新媒体交易平台
提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易，真正实现买卖双方的共赢。【请点击这里访问】
TOP热门关键词
求解方程的数值解，已经有了第一部计算，现在的问题是迭代计算，因为迭代计算我不太懂，所以程序也不会编f(x)=0x(k+1)=x(k)-{J(x(k))}^(-1)*f(x(k)),k=0,1,2….其中J和f已经计算出来了，这个Newtons迭代函数编写如下 ...
坛友互助群 |
扫码加入各岗位、行业、专业交流群 |
f(x)=0
x(k+1)=x(k)-{J(x(k))}^(-1)*f(x(k)), k=0,1,2….
其中J和f已经计算出来了，这个Newtons迭代函数编写如下？
Newtons<-function (funs, y, ep=1e-5, it_max=100){
index<-0; k<-1
while (k<=it_max){
y1 <- y; obj <- funs(y);
y<- y-solve(obj$J, obj$f);
........
其中要迭代的函数：
funs<-function(x){
。。。
f<- c(。。,。。。) #这是f(.)
J<- matrix(,..,….,…,
nrow=2, byrow=T)
# Jacobi矩阵 j(.)
return(
list(f=f, J=J))
}
要迭代部分函数是正确的，运行问题在于：
错误于drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) :
Lapack例行程序dgesv: 系统正好是奇异的
是不是迭代函数存在问题啊，谢谢了
扫码或添加微信号：坛友素质互助
「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
人气文章
本文标题：一个方程的求解，数值解 （R语言编写的）
本文链接网址：https://bbs.pinggu.org/jg/ruanjianpeixun_ryuyanpeixun_840134_1.html
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
