| 查看: 1406 | 回复: 0 | ||
michunrong木虫 (小有名气)
|
[求助]
如何使用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为止 ##请问该如何设置循环迭代语句 |
回复此楼» 猜你喜欢
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
-
调剂310
已经有4人回复
-
北科281学硕材料求调剂
已经有16人回复
-
一志愿上海交大生物与医药专硕324分，求调剂
已经有6人回复
-
一志愿哈工大，085400，320，求调剂
已经有4人回复
-
316求调剂
已经有9人回复
-
085600 材料与化工 329分求调剂
已经有8人回复
-
0856求调剂
已经有6人回复
-
299求调剂
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
采用MATLAB进行数值计算（英文第三版）
已经有145人回复
- 自己的写的matlab跑不动，哪位大神帮我用牛顿迭代编一下
已经有4人回复
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！
已经有11人回复
- 请教matlab如何增加迭代次数
已经有11人回复
- 求助matlab中含有三角矩阵的矩阵方程求解
已经有3人回复
- 什么人可以考注化？
已经有37人回复
- 如何用MATLAB求解脉冲微分方程并画图
已经有5人回复
- 实在是看不出这段matlab程序里哪里出错了，请大家帮忙看看！
已经有10人回复
- 全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题
已经有8人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- FORTRAN新手 求助主程序循环问题
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
