clc;clear;
m=2045;
I=5428;
a=1.488;
b=1.717;
L=a+b;
Cf=-77.85e3;
Cr=-76.51e3;
u=20;
tao=0.2;
h=0.0025;
LS=40;
%--------------------------------------------------------------------------
%
%??????????????????????直線行駛干擾影響
%
%--------------------------------------------------------------------------
Hs.den=[tao?1];?
Hs.num=[0?h];%駕駛員模型
Hs_sys=tf（Hs.numHs.den）;
%--------------------------------------------------------------------------
const1=-（a*Cf-b*Cr）./u;
const2=-（a^2*Cf+b^2*Cr）./u;
const3=-（Cf+Cr）./u;
const4=a*Cf-b*Cr;
Ys.num=conv（[00-Cf][Iconst2const4（1）]）-conv（[00-a*Cf][0const1Cf+Cr]）;
Ys.den=conv（[mconst30][Iconst2（1）const4]）-conv（[0const10][0const1Cf+Cr]）;
Ys_sys=tf（Ys.numYs.den）;
t=0:0.1:30;
y_Ys_sys=impulse（Ys_syst）*pi/180;
%--------------------------------------------------------------------------
scrsz?=?get（0‘ScreenSize‘）;
h?=?findobj（0?‘Name‘‘Line?Driver?vehicle?close‘）;
if?isempty（h）
h=figure（‘Position‘[scrsz（3）/4?scrsz（4）/4?scrsz（3）/2?scrsz（4）/2]‘name‘...
????‘Line?Driver?vehicle?close‘‘numbertitle‘‘off‘‘color‘‘w‘）;
end
uipanel（‘parent‘h‘position‘[0.5?0?0.002?1]‘backgroundcolor‘[0?0?0]...
????‘shadowcolor‘[0?0?0]‘highlightcolor‘[0?0?0]）;
uipanel（‘parent‘h‘position‘[0?0.5?1?0.002]‘backgroundcolor‘[0?0?0]...
????‘shadowcolor‘[0?0?0]‘highlightc