![](https://img-blog.csdnimg.cn/20201014180756926.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数值分析
文章平均质量分 73
yujunbeta
这个作者很懒,什么都没留下…
展开
-
多变量非线性方程求解问题(牛顿迭代法)
关于多变量的多元函数,我们求解考虑解的问题与单变量是极其类似的。考虑在x处的一阶泰勒展开,我们有0=f(x0)=f(x)+j(x)(x-x0)(后面的等号近似成立,j表示jacobian矩阵. 我们以一个例题来说明牛顿迭代法:考虑这个问题:求解方程组x^2+y^2-5=0 (x+1)*y-3*y-1=0假设给定初值(0,1),那么牛顿迭代法的R程序如下:f1原创 2013-03-09 15:37:45 · 14218 阅读 · 0 评论 -
优化与求解非线性方程组(单变量问题)
求函数极值的问题通常被化简为求解导数为0的点的问题。所以优化问题通常与解非线性方程组联系起来。在前面写点估计中的mle时,我们介绍了R中求解方程极值的函数nlm(),optim(). 我们以一元函数f(x)=ln(x)/(1+x)为例求解函数的极值。f<-function(x) -log(x)/(1+x) #(1)optimize(f,c(0,10))原创 2013-03-06 22:59:46 · 5252 阅读 · 0 评论 -
R语言里的矩阵处理学习笔记
关于矩阵,通常都会使用matlab来做处理。其实使用R也可以对矩阵做出一些简单的处理。而R语言中提供的matrix,matlab包也提供了不少关于矩阵处理的东西(可以通过??matlab来查看具体函数)。一、矩阵的输入通常我们使用函数matrix来创建矩阵,函数的介绍如下:matrix(data = NA,nrow = 1, ncol = 1, byrow = FALSE, dimna原创 2012-10-07 16:09:56 · 12926 阅读 · 2 评论