| 查看: 1221 | 回复: 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 ] |
回复此楼» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有177人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有118人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】 气质工作站安装以后不能用，请高手帮忙解决一下！急！
已经有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
