魔术公式轮胎模型

该文件中包含了魔术公式轮胎模型，代码全部用matlab编写，具有很高的参考价值和实际应用，非常方便！

function [F_xF_yM_z]=combined（alphagamma）
% gamma=0;
% Magic Formula
F_z=3000;
F_z0=3000;
R_0=0.3;
R_e=0.3;
p_Cx1=1.65;
p_Kx1=12;
p_Cy1=1.3;
p_Ey1=-1;
p_Ky1=10;
p_Vy1=0;
p_Dx1=1;
p_Kx2=10;
p_Dy1=1;
p_Ey2=0;
p_Ky2=1.5;
p_Vy2=0;
p_Dx2=0;
p_Kx3=-0.6;
p_Dy2=0;
p_Ey3=0;
p_Ky3=0;
p_Vy3=0.15;
p_Ex1=-0.5;
p_Hx1=0;
p_Dy3=0;
p_Ey4=0;
p_Hy1=0;
p_Vy4=0;
p_Ex2=0;
p_Hx2=0;
p_Hy2=0;
p_Ex3=0;
p_Vx1=0;
p_Hy3=0.25;
p_Ex4=0;
p_Vx2=0;
q_Bz1=6;
q_Bz2=-4;
q_Bz3=0.6;
q_Bz4=0;
q_Bz5=0;
q_Bz9=0;
q_Bz10=0.7;
q_Cz1=1.05;
q_Dz1=0.12;
q_Dz2=-0.03;
q_Dz3=0;
q_Dz4=0;
q_Dz6=0;
q_Dz7=0;
q_Dz8=0.6;
q_Dz9=0.2;
q_Ez1=-10;
q_Ez2=0;
q_Ez3=0;
q_Ez4=0;
q_Ez5=0;
q_Hz1=0;
q_Hz2=0;
q_Hz3=0;
q_Hz4=0;
r_Bx1=5;
r_Bx2=8;
r_Cx1=1;
r_Hx1=0;
r_By1=7;
r_By2=2.5;
r_By3=0;
r_Cy1=1;
r_Hy1=0.02;
r_Vy1=0;
r_Vy2=0;
r_Vy3=-0.2;
r_Vy4=14;
r_Vy5=1.9;
r_Vy6=10;
s_sz1=0;
s_sz2=-0.1;
s_sz3=-1.0;
s_sz4=0;
p_Dxphi1=0.4;
p_Dxphi2=0;
p_Dxphi3=0;
p_Kyphi=1;
p_Dyphi=0.4;
p_Dyph2=0;
p_Dyph3=0;
p_Dyph4=0;
p_Hyphi1=1;
p_Hyphi2=0.15;
p_Hyphi3=0;
p_Hyphi4=-4;
p_epsilongammaphi1=0;
p_epsilongammaphi2=0;
q_Dtphi1=10;
q_Ctphi1=0.2;
q_Crphi2=0.1;
q_Brphi1=0.1;
q_Drphi1=1;
q_Drphi2=-1.5;
A_mu=10;
%---------------------
lambda_Fz0=1;
lambda_mux=1;
lambda_muy=1;
lambda_muV=0;
lambda_Kxkappa=1;
lambda_Kya=1;
lambda_Cx=1;
lambda_Cy=1;
lambda_Ex=1;
lambda_Ey=1;
lambda_Hx=1;
lambda_Hy=1;
lambda_Vx=1;
lambda_Vy=1;
lambda_Kygamma=1;
lambda_Kzgamma=1;
lambda_t=1;
lambda_Mr=1;
lambda_xa=1;
lambda_yk=1;
lambda_Vyk=1;
lambda_s=1;
lambda_Cz=1;
lambda_Mx=1;
lambda_My=1;
r_Ex1=0.3;
r_Ex2=0.3;
% alpha=20*pi/180;

V_cx=30;
epsilon_x=0.001;
V_0=sqrt（9.8*0.3）;
kappa=-1:0.01:1;
zeta_i=1;
F_z0_apostrophe=3000;
zeta_5=1;
q_Bz6=1;
cos_alpha_apostrophe=cos（alpha）;
zeta_2=1;
zeta_0=1;
zeta_8=1;
zeta_7=1;
lambda_Ky=1;
zeta_6=1;

%---------------------
%longgitudinal force
F_z0_star=lambda_Fz0*F_z0;%1
df_z=（F_z-F_z0_star）/F_z0_star;%2
alpha_star=tan（alpha）;%3                             unknown
gamma_star=sin（gamma）;%4                                  unknown
% kappa=-V_sx/abs（V_cx）;%5                                  unknown
% cos（alpha_star）=V_cx/（V_c+epsilon_V）;%66a                unknown
zeta_i=1;
lambda_mux_star=lambda_mux;%7      unknown
lambda_muy_star=lambda_muy;%       unknown
% lambda_mux_star=lambda_mux/（1+lambda_muV*V_s/V_0）;%7      unknown
% lambda_muy_star=lambda_muy/（1+lambda_muV*V_s/V_0）;%       unknown
lambda_mux_apostrophe=A_mu*lambda_mux_star/（1+（A_mu-1）*lambda_mux_star）;%8
lambda_muy_apostrophe=A_mu*lambda_muy_star/（1+（A_mu-1）*lambda_muy_star）;
C_x=p_Cx1*lambda_Cx;%11   >0
mu_x=（p_Dx1+p_Dx2*df_z）*lambda_mux_star;%13   >0
D_x=mu_x*F_z*zeta_i;%12   >0
S_Hx=（p_Hx1+p_Hx2*df_z）*lambda_Hx;%17
K_xk=F_z*（p_Kx1+p_Kx2*df_z）*exp（p_Kx3*df_z）*lambda_Kxkappa;%15       =C_Fk
kappa_x=kappa+S_Hx;%10
E_x=（p_Ex1+p_Ex2*df_z+p_Ex3*df_z^2）*（1-p_Ex4*sign（kappa_x））*lambda_Ex;%14    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      10452  2008-03-30 21:26  extension.m

     文件       6517  2008-03-30 22:13  combined.m

     文件       6768  2008-03-30 22:43  factor.m

     文件       7511  2008-03-31 01:55  drawing.m

----------- ---------  ---------- -----  ----

                31248                    4