 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1235 | 回復(fù): 11 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
呂京毛衣銀蟲(chóng) (小有名氣)
|
[求助]
自己編寫(xiě)了一段Matlab代碼，運(yùn)行時(shí)出現(xiàn)bug 已有1人參與
|
||
|
一段模擬化學(xué)反應(yīng)的代碼，運(yùn)行時(shí)出現(xiàn)“Undefined function 'P' for input arguments of type 'double'，而我只是想用P（i）作為一個(gè)數(shù)組儲(chǔ)存分子數(shù)，請(qǐng)大家?guī)臀铱纯�，怎么修改？謝謝大家! 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 %鏈引發(fā)（不區(qū)分單體和初級(jí)自由基） 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 %鏈增長(zhǎng) 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;%儲(chǔ)存分子數(shù) else P(i)=1;%儲(chǔ)存終止產(chǎn)物鏈長(zhǎng) end I=[I;i]; end end end T=T+t;U=[U;T];%儲(chǔ)存時(shí)間 C=(2000-Xm)*100/2000; V=[V;C];%儲(chǔ)存單體轉(zhuǎn)化率 L=length(I);S=sum(I); E=S/L;%平均鏈長(zhǎng) Pn=[Pn;E]; end N=sum(P(i)); f(i)=P(i)./N;%數(shù)量分布 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) %轉(zhuǎn)化率-時(shí)間曲線(xiàn) subplot(224); plot(U,Pn) %平均鏈長(zhǎng)-時(shí)間曲線(xiàn) |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有58人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
7樓2016-04-08 10:11:59
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專(zhuān)家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線(xiàn): 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專(zhuān)業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
2樓2016-04-07 13:47:30
1049569501
新蟲(chóng) (初入文壇)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 22.5
- 帖子: 20
- 在線(xiàn): 24.9小時(shí)
- 蟲(chóng)號(hào): 4306914
- 注冊(cè): 2015-12-23
- 性別: GG
- 專(zhuān)業(yè): 有機(jī)高分子功能材料
3樓2016-04-07 18:59:41
呂京毛衣
銀蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 461.3
- 散金: 7
- 帖子: 121
- 在線(xiàn): 15.1小時(shí)
- 蟲(chóng)號(hào): 3245265
- 注冊(cè): 2014-05-30
- 專(zhuān)業(yè): 高分子合成化學(xué)
4樓2016-04-08 08:36:12
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +8 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 299求調(diào)劑 +3 | 某某某某位 2026-03-21 | 3/150 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專(zhuān)碩 +19 | 不知道叫什么！ 2026-03-15 | 21/1050 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 8/400 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 材料工程專(zhuān)碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專(zhuān)碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|