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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

revisionscc

銅蟲 (初入文壇)

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

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

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

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

%以完全彈性碰撞模型為基礎(chǔ)模擬N個粒子在初始狀態(tài)為同速率的情況下，經(jīng)過一段時間的速度以及速率分布
%
%設(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);                                        %求分子非對心碰撞模式
         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ù)此樓

» 猜你喜歡

生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 已經(jīng)有6人回復(fù)
321求調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有3人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑已經(jīng)有6人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)
學(xué)碩274求調(diào)劑已經(jīng)有8人回復(fù)
壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！已經(jīng)有7人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
22408 359分調(diào)劑已經(jīng)有4人回復(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小時
蟲號: 7733204
注冊: 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 的主題更新