版塊導(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)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有124人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒有用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í)候沒有問(wèn)題，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小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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呀，我試過(guò)把這塊程序拿出來(lái)寫成一個(gè)小程序，就是可以運(yùn)行的，不需要定義初始值或數(shù)組。這讓我很疑惑是不是問(wèn)題出在我整個(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ù)組也沒有定義什么，為什么就沒有問(wèn)題呢？

發(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ā)表

[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 生物學(xué)調(diào)劑 +3	Surekei 2026-03-21	3/150	2026-03-21 18:31 by 學(xué)員8dgXkO
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +8	楓橋ZL 2026-03-18	10/500	2026-03-21 15:29 by Shawn0911
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客