| 查看: 1385 | 回复: 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为止 ##请问该如何设置循环迭代语句 |
回复此楼» 猜你喜欢
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人