tic;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n=1e5;
m=3e5;
p=m/20;
l=m/p;
A=13.5083*2*12.762*2*10^-20;
T=300;kb=1.38*10^-23;
factor=((4.18*10^13)/(6.022*10^23))^2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear forcex forcey;
for i=n:n+m
forcex(i-(n-1),1)=fx(i,1);
forcey(i-(n-1),1)=fy(i,1);
end
%%%%%%%%
clear myacfx myacfy;
myacfx=autocorr(forcex,p-1);
myacfy=autocorr(forcey,p-1);
%%%%%%%%
clear interx intery inter;
pp=p/200;
for i=pp:pp:p
clear t acx acy;
for j=1:i
t(j,1)=(j-1)*1e-15;
acx(j,1)=myacfx(j,1);
acy(j,1)=myacfy(j,1);
end
interx(i,1)=t(j,1);
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余2页未读,立即下载