| 查看: 1456 | 回复: 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个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
找到一些相关的精华帖子，希望有用哦~
采用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
10
