基于MATLAB的m文件制动力分配曲线的绘制

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;