版塊導(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)劑公告

返回列表

yali1990

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 469.8
帖子: 53
在線: 17.7小時(shí)
蟲號(hào): 2108135
注冊(cè): 2012-11-05
性別: MM
專業(yè): 高分子物理與高分子物理化

[求助] 自規(guī)避行走鏈

#include "stdafx.h"
#include "math.h"
#include"stdlib.h"//調(diào)用rand()及srand()函數(shù)
#include "time.h"
#define PI 3.1415926
void mass_center_Rg2()
{
double X[10000][3];//鏈段坐標(biāo)
double MC_X[3];//質(zhì)心坐標(biāo)
int l=500;//鏈條數(shù)
int N;//聚合度
srand((unsigned)time(NULL));
double a,b,d;
double r=2.0;
double Sum_Rg2;
double Rg2;
double  X_Dis;
double  Y_Dis;
double  Z_Dis;
FILE *write;
write = fopen("Rg2_compute-1.txt","w"

;
for(N=1000;N<=2000;N+=20)
{
Sum_Rg2=0.0;
for(int j=1;j {
d=0.0;
X[0][0]=0.0;
X[0][1]=0.0;
X[0][2]=0.0;
MC_X[0] = 0.0;
MC_X[1] = 0.0;
MC_X[2] = 0.0;
for(int i=0;i {
d=((float)rand())/RAND_MAX;
a=d*PI;
b=d*PI*2.0;
X[i+1][0]=X[0]+r*sin(a)*cos(b);//鏈單元的坐標(biāo)
X[i+1][1]=X[1]+r*sin(a)*sin(b);//鏈單元的坐標(biāo)
X[i+1][2]=X[2]+r*cos(a);//鏈單元的坐標(biāo)
   if((i+1)>=2)
      {
      for(int counter=i-1;counter>=0;counter--)
      {
   double r0,x0,y0,z0;
            x0=X[i+1][0]-X[counter][0];
            y0=X[i+1][1]-X[counter][1];
            z0=X[i+1][2]-X[counter][2];
            r0=pow(x0,2)+pow(y0,2)+pow(z0,2);
            if(r0>=4.0)
      {
               d=((float)rand())/RAND_MAX;
            a=d*PI;
            b=d*PI*2.0;
            X[i+1][0]=X[0]+r*sin(a)*cos(b);//鏈單元的坐標(biāo)
            X[i+1][1]=X[1]+r*sin(a)*sin(b);//鏈單元的坐標(biāo)
            X[i+1][2]=X[2]+r*cos(a);//鏈單元的坐標(biāo)

         }
         else continue;
   }
}
      }
for(int i=1;i<=N;i++){
MC_X[0]=MC_X[0]+ X[0];
MC_X[1]=MC_X[1]+ X[1];
MC_X[2]=MC_X[2]+ X[2];
   }

      MC_X[0]=MC_X[0]/N;
      MC_X[1]=MC_X[1]/N;
      MC_X[2]=MC_X[2]/N;
         //printf("%f, %f ,%f\n",MC_X[0],MC_X[1],MC_X[2]);

Rg2=0.0;

for(int i=1;i<=N;i++){
X_Dis=X[0]-MC_X[0];
Y_Dis=X[1]-MC_X[1];
Z_Dis=X[2]-MC_X[2];
Rg2+=pow(X_Dis,2)+pow(Y_Dis,2)+pow(Z_Dis,2);
}

Rg2= Rg2/N;
Sum_Rg2+=Rg2;
}

Sum_Rg2=Sum_Rg2/l; //平均，求出均方半徑
fprintf(write,"輸出N=%d,Sum_Rg2=%f\n",N,Sum_Rg2);
printf("N= %d \n", N);
}
fclose(write);
}

int _tmain(int argc, _TCHAR* argv[])
{
    mass_center_Rg2();//調(diào)用函數(shù)求出鏈的質(zhì)心和Rg

      getchar();
      return 0;
}

這個(gè)是我編的代碼，目的是每產(chǎn)生一個(gè)坐標(biāo)，都判斷它是否與前面的坐標(biāo)重合，如果重合，則重新產(chǎn)生一個(gè)新的坐標(biāo)。接下來是算質(zhì)心坐標(biāo)，進(jìn)而計(jì)算均方回轉(zhuǎn)半徑。
目前我的問題是：這個(gè)和我之前計(jì)算無規(guī)行走鏈，得到的均方回轉(zhuǎn)半徑值基本一致，我要怎么改變坐標(biāo)的判斷條件？

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有85人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)

1樓 2012-12-06 17:50:28

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

氧化苯環(huán)側(cè)鏈時(shí)如何避免對(duì)萘環(huán)的影響？已經(jīng)有3人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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