close?all
m1=1825;%full
m=1450;
hg=0.95;
hg1=0.73;%full
g=9.8;
a=1.451;
a1=1.947;%full
b1=1.253;%full
b=1.749;
G=m*g;
G1=m1*g;%full
fa=0.1:0.1:1;
f_u1a=1:15000;
f_u2a=0.5*（G/hg*sqrt（b^2+4*hg*（a+b）/G*f_u1a）-（G*b/hg+2*f_u1a））
plot（f_u1‘f_u2‘）;
f_u2=zeros（10400）;
figure（1）
for?j=1:10
????f_u1=1:100:40000;?
????for?i=1:100:40000
????????f_u2（j:）=fa（j）*?m*g-f_u1（1:）;?????
????end
????plot（f_u1f_u2（j:）‘k‘）
????hold?on;
????grid?on;
????axis（[0?40000?0?40000]）
end
hold?on;
figure（2）
f_ur=zeros（10400）
for?j=1:10
????f_uf=1:100:40000;?
????for?i=1:100:40000
????????k2=b+fa（j）*hg;
????????k1=a-fa（j）*hg;
????????k（j）=k1./k2;
????????f_ur（j:）=k（j）*f_uf（1:）;
????end
????plot（f_uff_ur（j:）‘k‘）
????hold?on;