 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1234 | 回復(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)有6人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有284人回復(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)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
lifang2009
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 321.4
- 散金: 35
- 帖子: 204
- 在線: 33.2小時
- 蟲號: 878282
- 注冊: 2009-10-20
- 專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)
9樓2016-04-08 14:49:35
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強(qiáng)帖: 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è): 有機(jī)高分子功能材料
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ā)表 | |
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +6 | 墨墨漠 2026-03-20 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|