基本信息

---

下载

   友情提示：（无需注册或充值，赞助后即可获取资源下载链接）

---

     嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300

本次赞助数额为： 3 元　

支付宝支付

微信支付

   源码介绍

---

clear;clc;close all
i=sqrt(-1);
pi=3.1415926;
dt=power(10,-14);
delta_t=dt;
t=dt:dt:4*power(10,-9);

n=length(t);
N0=3*power(10,24);%m3
E0=power(10,0);
Y(:,1)=[0.000000000000000001;0];

T=power(10,-9);%ns
for k=1:n-1
    z1=C(t(k),Y(:,k));
     z2=C(t(k) delta_t/2,Y(:,k) z1*delta_t/2);
      z3=C(t(k) delta_t/2,Y(:,k) z2*delta_t/2);
       z4=C(t(k) delta_t,Y(:,k) z3*delta_t);
       Y(:,k 1)=Y(:,k) delta_t*(z1 2*z2 2*z3 z4)/6;
end


chao=Y(1,:);
 x1=power(chao,2);

lag=2000;

ACF= autocorr(x1,lag);
ACF = ACF(2:end);   
figure(1)
plot(ACF);