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

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

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

返回列表

下一頁(yè)

呂京毛衣

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 461.3
散金: 7
帖子: 121
在線: 15.1小時(shí)
蟲號(hào): 3245265
注冊(cè): 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ù)組儲(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í)間曲線
  subplot(224);
  plot(U,Pn) %平均鏈長(zhǎng)-時(shí)間曲線

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有231人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(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)商家推薦: (我也要在這里推廣)

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

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

回帖支持 ( 顯示支持度最高的前 50 名 )

lifang2009

銀蟲 (小有名氣)

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

P(i)被賦值時(shí)候沒有問題，Matlab會(huì)自動(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的回帖

普通回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(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í)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

P的初始值沒有定義。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(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í)
蟲號(hào): 4306914
注冊(cè): 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í)
蟲號(hào): 3245265
注冊(cè): 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的回帖

呂京毛衣

銀蟲 (小有名氣)

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

引用回帖:

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

我運(yùn)行時(shí)第一次會(huì)進(jìn)入到選擇P（i）＝1呀，我試過把這塊程序拿出來寫成一個(gè)小程序，就是可以運(yùn)行的，不需要定義初始值或數(shù)組。這讓我很疑惑是不是問題出在我整個(gè)程序里的其他地方？

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

贊一下

回復(fù)此樓

5樓2016-04-08 08:46:20

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

呂京毛衣

銀蟲 (小有名氣)

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

引用回帖:

3樓: Originally posted by 1049569501 at 2016-04-07 18:59:41
沒有定義函數(shù)或數(shù)組P，在循環(huán)前定義一下就行了。

里面的R（i）數(shù)組也沒有定義什么，為什么就沒有問題呢？

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

贊一下

回復(fù)此樓

6樓2016-04-08 09:17:12

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

鑫鑫0001

祝福

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

7樓2016-04-08 10:11:59

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

lifang2009

銀蟲 (小有名氣)

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

可以在z=find(l==i)加斷點(diǎn)調(diào)試，肯定先執(zhí)行第一個(gè)分支了

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

贊一下

回復(fù)此樓

9樓2016-04-08 14:49:35

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

呂京毛衣

銀蟲 (小有名氣)

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

引用回帖:

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

這個(gè)代碼就是可以運(yùn)行的

I=[];t=0
while (t<=3)
i=ceil(rand(1));
z=find(I==i);
if z>=1
P(i)=P(i)+1;
else
P(i)=1;
end
I=[I;i],t=t+1
end
P

贊一下

回復(fù)此樓

10樓2016-04-08 16:04:57

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主呂京毛衣的主題更新

返回列表

下一頁(yè)

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

[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 279分求調(diào)劑一志愿211 +11	chaojifeixia 2026-03-19	12/600	2026-03-21 01:49 by 星空星月
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 319求調(diào)劑 +3	小力氣珂珂 2026-03-20	3/150	2026-03-20 19:47 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 266求調(diào)劑 +5	陽(yáng)陽(yáng)哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽(yáng)陽(yáng)哇塞
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu

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

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

» 猜你喜歡

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

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

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