汽车制动系统的匹配设计，关于matlab！

Question

问题：我编了个程序，本来要求是“QQ截图未命名1”的图像的，但做出来的却只有两条线，颜色也一样的，不知道问题出在那里，求解，非常感谢，十万火急
程序如下：
ma_k=1313; % 空载质量（kg）
ma_m =1925; % 满载质量（kg）
g=9.80665; % 重力加速度（m/s2）
hg_k=640; % 空载质心高（mm）
hg_m =742; % 满载质心高（mm）
L=3050; % 轴距（mm）
load_f_k=653; % 空载前轴负荷
load_f_m=660; % 满载前轴负荷
b_k=load_f_k*L/ma_k; % 空载质心到后轴的距离（mm）
b_m =load_f_m *L/ma_m ;% 满载质心到后轴的距离（mm）
a_k=L-b_k; % 空载质心到前轴的距离（mm）
a_m =L-b_m ; % 满载质心到前轴的距离（mm）
G_k=ma_k*g; % 空载重力（N）
G_m =ma_m *g; % 满载重力（N）
Fu1=(0:G_m);% 前制动器制动力
Fu2_k=0.5*(G_k/hg_k* (b_k^2+4*hg_k*L*Fu1/G_k).^0.5-(G_k*b_k/hg_k+2*Fu1)); % 空载I线
plot(Fu1,Fu2_k,'g')
hold on
Fu2_m=0.5*(G_m /hg_m *(b_m ^2+4*hg_m *L*Fu1/G_m ).^0.5-(G_m *b_m /hg_m +2*Fu1)); % 满载I线
plot(Fu1,Fu2_m,'b')
hold on
B_z=Fu1/(Fu1+Fu2_k); %β线
plot(Fu1,Fu2_k)
hold off

Answer 1

%问题出在倒数第二行
clear;
close all;
clc;

ma_k=1313; % 空载质量（kg）
ma_m =1925; % 满载质量（kg）
g=9.80665; % 重力加速度（m/s2）
hg_k=640; % 空载质心高（mm）
hg_m =742; % 满载质心高（mm）
L=3050; % 轴距（mm）
load_f_k=653; % 空载前轴负荷
load_f_m=660; % 满载前轴负荷
b_k=load_f_k*L/ma_k; % 空载质心到后轴的距离（mm）
b_m =load_f_m *L/ma_m ;% 满载质心到后轴的距离（mm）
a_k=L-b_k; % 空载质心到前轴的距离（mm）
a_m =L-b_m ; % 满载质心到前轴的距离（mm）
G_k=ma_k*g; % 空载重力（N）
G_m =ma_m *g; % 满载重力（N）
Fu1=(0:G_m);% 前制动器制动力
Fu2_k=0.5*(G_k/hg_k* (b_k^2+4*hg_k*L*Fu1/G_k).^0.5-(G_k*b_k/hg_k+2*Fu1)); % 空载I线
plot(Fu1,Fu2_k,'g')
hold on
Fu2_m=0.5*(G_m /hg_m *(b_m ^2+4*hg_m *L*Fu1/G_m ).^0.5-(G_m *b_m /hg_m +2*Fu1)); % 满载I线
plot(Fu1,Fu2_m,'b')
hold on
B_z=Fu1/(Fu1+Fu2_k); %β线，这里应该是根据制动系的设计给出的，不应这样求
plot(Fu1,(1-B_z)*(Fu1+Fu2_k))%问题出在这里
hold off