版塊導(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í)熱門版塊排行榜

返回列表

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

呂京毛衣

銀蟲 (小有名氣)

應(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 %鏈增長
         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)物鏈長
            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;%平均鏈長
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)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有81人回復(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的回帖

呂京毛衣

銀蟲 (小有名氣)

應(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的回帖

查看全部 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í)
蟲號(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的初始值沒有定義。

我運(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的回帖

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

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

[考研] 求助 +4	夢里的無言 2026-03-21	5/250	2026-03-21 15:53 by ColorlessPI
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 323求調(diào)劑 +3	洼小桶 2026-03-18	3/150	2026-03-20 22:54 by JourneyLucky
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海

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

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

呂京毛衣

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

» 猜你喜歡

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

呂京毛衣

月只藍(lán)

1049569501

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

呂京毛衣

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