 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1239 | 回復: 11 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
呂京毛衣銀蟲 (小有名氣)
|
[求助]
自己編寫了一段Matlab代碼，運行時出現(xiàn)bug 已有1人參與
|
||
|
一段模擬化學反應的代碼，運行時出現(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];%儲存單體轉化率 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) %轉化率-時間曲線 subplot(224); plot(U,Pn) %平均鏈長-時間曲線 |
回復此樓» 猜你喜歡
哈爾濱理工大學物理系招收物理學考研調劑
已經(jīng)有7人回復
-
0702一志愿吉大B區(qū)求調劑
已經(jīng)有5人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有256人回復
-
求調劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
-
光學工程學碩調劑信息
已經(jīng)有31人回復
-
歡迎加入課題組
已經(jīng)有0人回復
-
散金幣，求好運，祝面上順利！
已經(jīng)有68人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經(jīng)有1人回復
-
調劑
已經(jīng)有0人回復
» 本主題相關商家推薦: (我也要在這里推廣)
1049569501
新蟲 (初入文壇)
- 應助: 3 (幼兒園)
- 金幣: 22.5
- 帖子: 20
- 在線: 24.9小時
- 蟲號: 4306914
- 注冊: 2015-12-23
- 性別: GG
- 專業(yè): 有機高分子功能材料
3樓2016-04-07 18:59:41
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學
- 管轄: 計算模擬區(qū)
2樓2016-04-07 13:47:30
呂京毛衣
銀蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 461.3
- 散金: 7
- 帖子: 121
- 在線: 15.1小時
- 蟲號: 3245265
- 注冊: 2014-05-30
- 專業(yè): 高分子合成化學
4樓2016-04-08 08:36:12
呂京毛衣
銀蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 461.3
- 散金: 7
- 帖子: 121
- 在線: 15.1小時
- 蟲號: 3245265
- 注冊: 2014-05-30
- 專業(yè): 高分子合成化學
|
我運行時第一次會進入到選擇P（i）＝1呀，我試過把這塊程序拿出來寫成一個小程序，就是可以運行的，不需要定義初始值或數(shù)組。這讓我很疑惑是不是問題出在我整個程序里的其他地方？ 發(fā)自小木蟲IOS客戶端 |
5樓2016-04-08 08:46:20
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 311求調劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 277材料科學與工程080500求調劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 298求調劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
|
[考研] 306求0703調劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
|
[考研] 299求調劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[考研] 332求調劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 085700資源與環(huán)境308求調劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 296求調劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 294求調劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 319求調劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調劑 +3 | 領帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 0856調劑，是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 0703化學調劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 328求調劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調劑 2026-03-16 | 12/600 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 085601專碩，總分342求調劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 070300化學學碩求調劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|