hezai.m

来自「用matlab便的支持向量机预测方面的界面」· M 代码 · 共 57 行

%计算zuo节点荷载


Q=10;
D=5;
Lb=1;

yspac = D/20;
lthu=21;

qn=Q/D*yspac;
qs=1/6*qn*yspac;
qx=2*qs;
qk=0.5*qn*yspac;

degao=zeros(3,lthu);
for i=1:lthu
    degao(1,i)=qs;
end
for i=1:lthu
    degao(2,i)=qx;
end
for i=1:lthu
    degao(3,i)=qk;
end


ind=0;
suan=0;
jisuan=zeros(4,lthu);
for i=1:lthu
    ind=ind+1;
    suan=suan+degao(3,i);
    jisuan(1,ind)=suan;
end
for i=3:lthu-1
    ind=ind+3;
    suan=suan+degao(3,i);
   jisuan(2,ind)=suan;
end
for i=1:lthu-1
    jisuan(3,i)=degao(1,i);
end
for i=2:lthu
    jisuan(4,i)=degao(2,i);
end

suanmain=zeros(2,lthu);
for i=1:lthu
    suanmain(1,i)=jisuan(1,i)+jisuan(2,i)+jisuan(3,i)+jisuan(4,i);
end

for i=1:lthu
    f(2*i-1)=suanmain(1,i);
    f(2*i)=suanmain(2,i);
end
f