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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

revisionscc

銅蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 38
帖子: 8
在線: 35.5小時(shí)
蟲號(hào): 1470687
注冊(cè): 2011-11-01
性別: GG
專業(yè): 原子和分子物理

[求助] 小弟初學(xué)MATLAB，編了一個(gè)仿真麥克斯韋速率分布的程序，求優(yōu)化識(shí)錯(cuò)

大致就是如題，我已經(jīng)分析出了在分子間碰撞那個(gè)循環(huán)會(huì)對(duì)一次碰撞做兩次運(yùn)算，但不知道怎么解決，并且這個(gè)程序運(yùn)行時(shí)間太長，希望優(yōu)化一下程序，求大牛們指導(dǎo)。。。

以下是文件內(nèi)容以及附帶的.m文件。

%以完全彈性碰撞模型為基礎(chǔ)模擬N個(gè)粒子在初始狀態(tài)為同速率的情況下，經(jīng)過一段時(shí)間的速度以及速率分布
%
%設(shè)模擬的分子數(shù)為10000，二維方形盒子邊長為10^3。初始速度為400
%a為X方向速度矩陣，b為y方向速度矩陣，c為x方向位移矩陣，d為y方向位移矩陣
%
clear
k=rand(1,10000);
a=400*cos(2*pi*k);       %Vx
b=400*sin(2*pi*k);       %Vy
c=10^3*k;                %Sx
k=rand(1,10000);          %重置隨機(jī)變量
d=10^3*k;                %Sy
t=0.01;                   %t
for n=1:1000;
c=c+t*a;
d=d+t*b;
e=(c>=1000&c<=0);    %1為碰壁分子下標(biāo)       17---22 解決分子碰壁問題
f=(d>=1000&d<=0);    %0為不碰壁分子
c=c-e.*a*t;          %重置碰壁分子位移
d=d-f.*b*t;          %
a=a-2*a.*e;          %碰壁分子碰撞方向速度反向
b=b-2*b.*f;          %
g=round(c);          %開始處理分子間碰撞問題，約化碰撞半徑為1       23---39 解決分子間碰撞問題
h=round(d);          %四舍五入分子芯的位置
for i=1:10000       %按順序?qū)ふ曳肿有疚恢孟嗤牧Ｗ?br />       if length(find(and(not((g-g(i))),not((h-h(i))))))~=2    %排除不碰撞分子與兩分子以上的碰撞
      else
         j=find(and(not((g-g(i))),not((h-h(i)))));             %找出相互碰撞的兩分子坐標(biāo)
         u=j(1);                                              %
         v=j(2);
         xx=c(u)-c(v);                                        %求分子非對(duì)心碰撞模式
         yy=d(u)-d(v);                                        %
         theta=atan(yy/xx);                                  %
         a(u)=(a(v)*cos(theta)+b(v)*sin(theta))*cos(theta)+(a(u)*sin(theta)+b(u)*cos(theta))*sin(theta);  %碰后分子速度
         b(u)=(a(v)*cos(theta)+b(v)*sin(theta))*sin(theta)+(a(u)*sin(theta)+b(u)*cos(theta))*cos(theta);  %
         a(v)=(a(u)*cos(theta)+b(u)*sin(theta))*cos(theta)+(a(v)*sin(theta)+b(v)*cos(theta))*sin(theta);  %
         b(v)=(a(u)*cos(theta)+b(u)*sin(theta))*sin(theta)+(a(v)*sin(theta)+b(v)*cos(theta))*cos(theta);  %
      end
end
end
hist(sqrt(a.^2+b.^2),100)       %畫出速率分布直方圖

回復(fù)此樓

» 猜你喜歡

壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！已經(jīng)有4人回復(fù)
081200-314 已經(jīng)有4人回復(fù)
289求調(diào)劑已經(jīng)有9人回復(fù)
085405 考的11408求各位老師帶走已經(jīng)有3人回復(fù)
339求調(diào)劑已經(jīng)有6人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有4人回復(fù)
085701環(huán)境工程，267求調(diào)劑已經(jīng)有15人回復(fù)
材料求調(diào)劑一志愿哈工大324 已經(jīng)有9人回復(fù)
292求調(diào)劑已經(jīng)有15人回復(fù)
304求調(diào)劑已經(jīng)有5人回復(fù)

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

我是好人

1樓 2012-02-26 20:35:03

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

bruceall

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5
帖子: 3
在線: 1.4小時(shí)
蟲號(hào): 7733204
注冊(cè): 2018-01-06

能指導(dǎo)一下嗎

回復(fù)此樓

2樓2018-01-06 17:38:35

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

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