迭代循环

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我想做的事情如下：
语句1；
语句2；
语句3；
有三个未知数，a 是一个向量，B, C均为矩阵。
首先通过语句1，给定B,C的初值可以得到a,
通过语句2，由上步得到的a和初值B可以得到C
通过语句3，由上步得到的a和初值C可以得到B
一直这样循环迭代，一直到收敛

问题是，如何写这个过程呢？
用while吗？
那while里面判断应该是前后两个相差比较小吗？我想同时让a, B, C收敛。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：while Hil 如何写 未知数 如何 未知数

有没有人在？求解惑啊！！！

同求问~~~

坐等答疑者啊！！！！

或者用repeat ……break语句？

size.norm2<-function(s,alpha,d,m){
      t0<-qt(alpha/2,m,lower.tail=FALSE)
     n0<-(t0*s/d)^2
      t1<-qt(alpha/2,n0,lower.tail=FALSE)
     n1<-(t1*s/d)^2
     while(abs(n1-n0)>0.5){
     n0<-(qt(alpha/2,n1,lower.tail=FALSE)*s/d)^2
     n1<-(qt(alpha/2,n0,lower.tail=FALSE)*s/d)^2
}
n1
}

给个例子你，这个是t分布求样本空间的问题， m可以先搞个大点的数 比如1000，再通过迭代，控制while(abs(..)>0)等等，repeat  if ...break应该可以的


你把具体的例子 发出 看看撒

jmpamao 发表于 2012-12-19 21:04
size.norm2

我写的循环如下：
repeat{
         fr<-function(lm0){
             A<-matrix(rep(0,dm^2),ncol=dm)
                   for(i in 1:dm)
                     {
                       B<-lm0*Kx0[i,]%*%t(Et0[,i])
                       A<-B+A
                      }
         X<-M-B
         norm(X,"F")^2-ln*sum(w*abs(lm0))}
         ###optim通过Kx0,Et0初值得到最优lm0

         lm1<-optim(lm0,fr)
         ###奇异值分解利用lm1与Kx0得到Et1
         Et1<-svd(M%*%Kx0%*%diag(lm1[[1]]))$u%*%t(svd(M%*%Kx0%*%diag(lm1[[1]]))$v)
         ###奇异值分解利用lm1与Et0得到Kx1
         Kx1<-svd(M%*%Et0%*%diag(lm1[[1]]))$u%*%t(svd(M%*%Et0%*%diag(lm1[[1]]))$v)
         ##反复迭代直至收敛
         if(sum(abs(lm0-lm1[[1]]))<=eps)  break
         lm0<-lm1[[1]]
         Et0<-Et1
         Kx0<-Kx1   }
其中lm0,Et0,Kx0,是需要迭代至收敛的向量，矩阵，矩阵。
在程序中，我先控制lm0试了下，可是还是报错了，我eps取1
报错如下：
Error in optim(lm0, fr) : non-finite value supplied by optim

水平有限 看不懂