 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1243 | 回復(fù): 11 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
呂京毛衣銀蟲 (小有名氣)
|
[求助]
自己編寫了一段Matlab代碼，運行時出現(xiàn)bug 已有1人參與
|
||
|
一段模擬化學(xué)反應(yīng)的代碼，運行時出現(xiàn)“Undefined function 'P' for input arguments of type 'double'，而我只是想用P（i）作為一個數(shù)組儲存分子數(shù)，請大家?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ū)分單體和初級自由基） 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;%儲存分子數(shù) else P(i)=1;%儲存終止產(chǎn)物鏈長 end I=[I;i]; end end end T=T+t;U=[U;T];%儲存時間 C=(2000-Xm)*100/2000; V=[V;C];%儲存單體轉(zhuǎn)化率 L=length(I);S=sum(I); E=S/L;%平均鏈長 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)化率-時間曲線 subplot(224); plot(U,Pn) %平均鏈長-時間曲線 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有73人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
呂京毛衣
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 461.3
- 散金: 7
- 帖子: 121
- 在線: 15.1小時
- 蟲號: 3245265
- 注冊: 2014-05-30
- 專業(yè): 高分子合成化學(xué)
|
我運行時第一次會進入到選擇P（i）＝1呀，我試過把這塊程序拿出來寫成一個小程序，就是可以運行的，不需要定義初始值或數(shù)組。這讓我很疑惑是不是問題出在我整個程序里的其他地方？ 發(fā)自小木蟲IOS客戶端 |
5樓2016-04-08 08:46:20
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
2樓2016-04-07 13:47:30
1049569501
新蟲 (初入文壇)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 22.5
- 帖子: 20
- 在線: 24.9小時
- 蟲號: 4306914
- 注冊: 2015-12-23
- 性別: GG
- 專業(yè): 有機高分子功能材料
3樓2016-04-07 18:59:41
呂京毛衣
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 461.3
- 散金: 7
- 帖子: 121
- 在線: 15.1小時
- 蟲號: 3245265
- 注冊: 2014-05-30
- 專業(yè): 高分子合成化學(xué)
4樓2016-04-08 08:36:12
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 招08考數(shù)學(xué) +4 | laoshidan 2026-03-20 | 8/400 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|