版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

呂京毛衣

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 461.3
散金: 7
帖子: 121
在線: 15.1小時(shí)
蟲號: 3245265
注冊: 2014-05-30
專業(yè): 高分子合成化學(xué)

[求助] 自己編寫了一段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ù)組儲存分子數(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];%儲存時(shí)間
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)化率-時(shí)間曲線
  subplot(224);
  plot(U,Pn) %平均鏈長-時(shí)間曲線

回復(fù)此樓

» 猜你喜歡

光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有174人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有79人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2016-04-07 09:58:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lifang2009

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 321.4
散金: 35
帖子: 204
在線: 33.2小時(shí)
蟲號: 878282
注冊: 2009-10-20
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

P(i)被賦值時(shí)候沒有問題，Matlab會自動(dòng)處理，但是在右側(cè)作為一個(gè)操作數(shù)的時(shí)候必須有確定的值

發(fā)自小木蟲IOS客戶端

贊一下(1人)

回復(fù)此樓

8樓2016-04-08 14:46:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 12 個(gè)回答

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

P的初始值沒有定義。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

2樓2016-04-07 13:47:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

1049569501

新蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 22.5
帖子: 20
在線: 24.9小時(shí)
蟲號: 4306914
注冊: 2015-12-23
性別: GG
專業(yè): 有機(jī)高分子功能材料

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

沒有定義函數(shù)或數(shù)組P，在循環(huán)前定義一下就行了。

贊一下

回復(fù)此樓

3樓2016-04-07 18:59:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

呂京毛衣

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 461.3
散金: 7
帖子: 121
在線: 15.1小時(shí)
蟲號: 3245265
注冊: 2014-05-30
專業(yè): 高分子合成化學(xué)

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2016-04-07 13:47:30
P的初始值沒有定義。

這種數(shù)組的初始值怎么定義呢？應(yīng)該是空的吧

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

4樓2016-04-08 08:36:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 12 個(gè)回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 303求調(diào)劑 +4	元夕元 2026-03-20	5/250	2026-03-25 00:20 by 元夕元
[考研] 08工學(xué)調(diào)劑 +12	用戶573181 2026-03-20	17/850	2026-03-24 21:45 by 番茄燉了西紅柿
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6	小羊睡著了? 2026-03-23	9/450	2026-03-24 20:21 by 小羊睡著了?
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4	17501029541 2026-03-23	6/300	2026-03-24 10:21 by 戴圍脖的小蚊子
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +7	vv迷 2026-03-22	7/350	2026-03-23 23:44 by Txy@872106
[考研] 361求調(diào)劑 +3	Glack 2026-03-22	3/150	2026-03-23 22:03 by fuyu_
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 333求調(diào)劑 +3	ALULU4408 2026-03-23	3/150	2026-03-23 19:04 by macy2011
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 考研調(diào)劑 +4	來好運(yùn)來來來 2026-03-21	4/200	2026-03-22 12:15 by 星空星月
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 求助 +5	夢里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

呂京毛衣

[求助] 自己編寫了一段Matlab代碼，運(yùn)行時(shí)出現(xiàn)bug 已有1人參與

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

lifang2009

月只藍(lán)

1049569501

【答案】應(yīng)助回帖

呂京毛衣

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 自己編寫了一段Matlab代碼，運(yùn)行時(shí)出現(xiàn)bug 已有1人參與