冲击电流发生器matlab程序非线性元件

%非线性电阻冲击电流发生器计算程序（8/20）
function cjdlfsq_fxx1（RLu）        
i=zeros（999）;
t=zeros（999）;
dt=0.05;
af=0.25;
ce=700;
nm=450;
nn=540;
C=15;
A=af-1;
t（1）=dt;
di=u*dt/L;
i（1）=di;
i（1）
du=i（1）*dt/C;
u=u-du;
for n=1:1:nn
    k=n+1;
    t（k）=t（n）+dt;
    R1（k）=ce*i（n）^A;
    R2（k）=R+R1（k）;
    di=（u-i（n）*R2（k））*dt/L;
    i（k）=i（n）+di;
    du=i（k）*dt/C;
    u=u-du;
end
im=i（1）;
for j=1:nm
    if（i（j）>im）
        im=i（j）;
    end
end
for j=1:1:nm
    if（i（j）==im）
        m=j;
    end
end
for j=1:nm
    jj=j+m;
    i51=i（jj）-0.5*im;
    i52=i（jj+1）-0.5*im;
    if