| 查看: 1289 | 回复: 5 | ||
[求助]
求助高手帮忙用matlab接下面这个算式，具体见附件
|
|
本人在科研过程中找到了这个公式，并做试验验证，公式的计算规则如下，不知道要怎么用matlab计算，请高手帮忙算一下： AN=((0.02105*xN^0.9208)/(1- xN)^0.9208-aN)*0.0002778 1式 BN=(1.515*0.2842* xN /(1+0.2842* xN)-bN)*0.0002778 2式 M=11.8554*A+2.0937*B 3式 总式为0.0010652*(3600*0.8/300-3168*cN/298)-M=0 4式 ①N=1时，a1 =0.007，b1=0.048，x1 =0.0903 将a1、b1、x1代入解出4式中的cN值为c1； ②N=2时a2=0.007+ A1 b2=0.048+ B1 x2= c1 代入取出值c2 ③ N=3时a3=0.007+ A1+ A2 b3=0.048+ B1+ B2 x3= c2 代入取出值c3 ④ N=4时a4=0.007+ A1+ A2+ A3 b4=0.048+ B1+ B2+ B3 x4= c3 代入取出值c4 ⑤ N=5时a5=0.007+ A1+ A2+ A3+ A4 b5=0.048+ B1+ B2+ B3+ B4 x5= c4 代入取出值c5 以此类推。 请高手给出m文件并将计算的结果绘制成曲线，cN-N 或者给出表格形式的cN，N 若不理解请看附件[ Last edited by likuihei on 2011-10-3 at 22:05 ] |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有69人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】 气质工作站安装以后不能用，请高手帮忙解决一下！急！
已经有8人回复
- 求高手帮忙 润色一下我的英语翻译
已经有2人回复
- 【求助】matlab 计算一个二重积分，高手们帮帮忙吧
已经有6人回复
- 求高手帮忙分析一下NMR和IR谱图，第一次做，无从下手，哪怕进来给点建议也好?
已经有43人回复
- 求助高人帮分析下这张孔径分布图，感觉矛盾，解释不通呀！
已经有10人回复
- 【求助/交流】怎么做都提不出DNA，请高手帮一下忙喽！！！！
已经有52人回复
- 【求助】为什么我用PCC氧化醇，得不到想要的醛呢~~~高手帮忙解答下啊
已经有5人回复
- 【求助】请高手帮忙解个方程，用Matlab
已经有6人回复
- 【求助】怎么用MATLAB求数值解啊（高手请帮忙啊）
已经有6人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
【答案】应助回帖
likuihei(金币+20): 2011-10-04 18:41:28
|
N = 10; n = 1:N; a = zeros(1,N); b = zeros(1,N); c = zeros(1,N); A = zeros(1,N); B = zeros(1,N); x = zeros(1,N); a(1) = 0.007; b(1) = 0.048; x(1) =0.0903; for i=1:N A(i) = ((0.02105*x(i)^0.9208)/(1- x(i))^0.9208-a(i))*0.0002778; B(i) = (1.515*0.2842* x(i) /(1+0.2842* x(i))-b(i))*0.0002778; M = 11.8554*A(i)+2.0937*B(i); c(i) = (3600*0.8/300-M/0.0010652)*298/3168; a(i+1) = a(i)+A(i); b(i+1) = b(i)+B(i); x(i+1) = c(i); end plot(n,c,'bo-') |
2楼2011-10-04 10:17:45
3楼2011-10-04 18:41:13
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
|
N = 50; a = zeros(1,N); b = zeros(1,N); c = zeros(1,N); A = zeros(1,N); B = zeros(1,N); x = zeros(1,N); M = zeros(1,N); S = zeros(1,N); a(1) = 0.007; b(1) = 0.048; x(1) =0.0903; for i=1:N A(i) = ((0.02105*x(i)^0.9208)/(1- x(i))^0.9208-a(i))*0.0002778; B(i) = (1.515*0.2842* x(i) /(1+0.2842* x(i))-b(i))*0.0002778; M(i) = 11.8554*A(i)+2.0937*B(i); S(i) = sum(M(1:i)); c(i) = (3600*0.8/300-S(i)/0.0010652)*298/3168; a(i+1) = a(i)+A(i); b(i+1) = b(i)+B(i); x(i+1) = c(i); end plot(1:N,c,'ro-') |
4楼2011-10-04 19:34:00
5楼2011-10-04 20:29:14
change0618
铁杆木虫 (著名写手)
 |
6楼2011-10-04 21:17:35
