24小时热门版块排行榜

>论坛更新日志 (2751)
>虫友互识 (259)
>文献求助 (180)
>导师招生 (157)
>考博 (111)
>休闲灌水 (71)
>论文投稿 (52)
>论文道贺祈福 (45)
>硕博家园 (37)
>博后之家 (35)
>基金申请 (33)
>教师之家 (30)
>公派出国 (23)
>招聘信息布告栏 (22)
>考研 (21)
>找工作 (20)

返回列表

michunrong

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 3713.5
红花: 4
帖子: 253
在线: 56.7小时
虫号: 2198363
注册: 2012-12-21
性别: GG
专业: 资源化工

[求助] 如何使用R软件编写循环迭代程序求大神帮助

##求最终的Tb和Theta_T##
##数据
Tmax=c(35.0,30.0,28.0,25.0,25.5,26.0)
Tmin=c(15.4,15.8,13.2,9.8,9.5,9.4)
X=c(0.233,0.233,0.220,0.065,0.138,0.083)
D=c(12,25,46,41,45,62)
Ts=40
To=30
Tbo=0

##Theta_g 其中一个过程量
if (Tbo<Tmin){
  Theta_g=(Tmax+Tmin)/2
} else if (Tbo<Tmax & Tbo>Tmin){
  Theta_g=(Tmax+Tmin)/2+(Tbo-Tmin)^2/2*(Tmax-Tmin)
} else
  Theta_g=Tbo

##Theta_a 另一个过程量
if (To>Tmin & Tmax<Ts){
  Theta_a=(Ts-Tbo)*(Tmax-To)^2/(2*(Ts-To)*(Tmax-Tmin))
} else if (To>Tmin & Tmax>=Ts){
  Theta_a=(Tmax*(Tmax-2*Tbo)-To*(Ts-Tbo)+Ts*Tbo)/(2*(Ts-To))
} else if (Tmin>=To & Tmin<Ts & Tmax<Ts){
  Theta_a=(Ts-Tbo)*(Tmax+Tmin-2*To)/(2*(Tmax-Tmin))
} else if (Tmin>=To & Tmin<Ts & Tmax>Ts){
  Theta_a=(Tmax+Tmin)/2-Tbo-(To-Tbo)*(Ts-Tmin)^2/(2*(Ts-To)*(Tmax-Tmin))
} else
  Theta_a=(Tmax+Tmin)/2-Tbo

## Theta_s
Theta_s=(Theta_g-Theta_a)*D

## Linear formular: Theta_s/X=Theta_T+Tb(D/X) 最小二乘法回归
y=Theta_s/X
x=D/X

##Least square linear regress simulation
Tb=sum((x-mean(x))*(y-mean(y)))/sum((x-mean(x))^2)
Theta_T=mean(y)-Tb*mean(x)

###判定条件是 |Tb-Tbo|<0.01则停止，如果|Tb-Tbo|>=0.01,则令Tbo=Tb，然后回到上面三个方程重新计算Tb,再进行判定，直到|Tb-Tbo|<0.01为止
##请问该如何设置循环迭代语句

回复此楼