| 查看: 1219 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
自己编写了一段Matlab代码，运行时出现bug 已有1人参与
|
|||
|
一段模拟化学反应的代码，运行时出现“Undefined function 'P' for input arguments of type 'double'，而我只是想用P（i）作为一个数组储存分子数，请大家帮我看看，怎么修改？谢谢大家! Xm=2000;Xa=2;Xr=0; % [Xm]=2 mol/L kd=0.00001;kp=1000;kt=10000000; kd1=0.00001;kp1=1;kt1=10000; % VNa=1000 T=0; j=0; I=[];X=[];U=[];V=[];Pn=[]; while(T<=2400) Rd=kd1*Xa; Rp=kp1*Xm*Xr; Rt=0.5*kt1*Xr*Xr; R=Rd+Rp+Rt; Pd=Rd/R;Pp=Rp/R;Pt=Rt/R; r1=rand(1); t=(1/R)*log(1/r1); r2=rand(1); if r2<=Pd %链引发（不区分单体和初级自由基） Xm=Xm-2; Xa=Xa-1; Xr=Xr+2; j=j+1; X=[X;j]; R(j)=1; j=j+1; X=[X;j]; R(j)=1; else if r2>Pd&&r2<=Pd+Pp %链增长 Xm=Xm-1; b=ceil(Xr*rand(1)); r3=X(b); R(r3)=R(r3)+1; else if r2>Pd+Pp&&r2<=1 Xr=Xr-2; c=ceil(Xr*rand(1)); r4=X(c); X(c)=[]; d=ceil((Xr-1)*rand(1)); r5=X(d); X(d)=[]; i=R(r4)+R(r5); z=find(I==i); if z>=1 P(i)=P(i)+1;%储存分子数 else P(i)=1;%储存终止产物链长 end I=[I;i]; end end end T=T+t;U=[U;T];%储存时间 C=(2000-Xm)*100/2000; V=[V;C];%储存单体转化率 L=length(I);S=sum(I); E=S/L;%平均链长 Pn=[Pn;E]; end N=sum(P(i)); f(i)=P(i)./N;%数量分布 m(i)=i.*P(i); M=sum(m(i)); w(i)=m(i)./M;%重量分布 subplot(221); plot(I,f(i)) subplot(222); plot(I,w(i)) subplot(223); plot(U,V) %转化率-时间曲线 subplot(224); plot(U,Pn) %平均链长-时间曲线 |
回复此楼» 猜你喜欢
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
物理学I论文润色/翻译怎么收费?
已经有163人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
lifang2009
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 321.4
- 散金: 35
- 帖子: 204
- 在线: 33.2小时
- 虫号: 878282
- 注册: 2009-10-20
- 专业: 导航、制导与传感技术
9楼2016-04-08 14:49:35
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68126.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2016-04-07 13:47:30
1049569501
新虫 (初入文坛)
- 应助: 3 (幼儿园)
- 金币: 22.5
- 帖子: 20
- 在线: 24.9小时
- 虫号: 4306914
- 注册: 2015-12-23
- 性别: GG
- 专业: 有机高分子功能材料
3楼2016-04-07 18:59:41
4楼2016-04-08 08:36:12
