matlab时域路面不平度,程序运行不了（路面不平度）

用三角级数叠加法模拟了4中路面状态下的不平度样本

L=64;         %桥跨

v1=50;      %车辆速度上限

v2=10;      %车辆速度下限

f1=0.2;                %对轮胎的激励时间频率下限

fu=100;                %对轮胎的激励时间频率下限

n1=f1/v1;      %空间频率下限n1=f1/v1,v1=f1/n1;

nu=(fu)/v2;     %空间频率上限

n0=1/2/pi;       %参考空间频率

l=1/2/nu;      %距离采样间隔ΔL= 0.05

N=L/ l;      %采样点数=段

n=1/L;      %空间频率间隔Δn

%路面不平整度系数

i=1;   %位移点数

for i=1:N

i=i+1;

end

for i=1:N

ni=i*n;

g1(i)=4*(abs((ni)/ n0) )^(-2) ;

end

for i=1:N

ni=i*n;

g2(i)=16*(abs((ni)/ n0) )^(-2) ;

end

for i=1:N

ni=i*n;

g3(i)=64*(abs((ni)/ n0) )^(-2) ;

end

for i=1:N

ni=i*n;

g4(i)=256*(abs(