2024年第二十六届“华东杯”(B题)大学生数学建模挑战赛|数学建模完整代码+建模过程全解全析

当大家面临着复杂的数学建模问题时,你是否曾经感到茫然无措?作为2022年美国大学生数学建模比赛的O奖得主,我为大家提供了一套优秀的解题思路,让你轻松应对各种难题。

让我们来看看华东杯 (B题)
第一个问题是对一辆前轮驱动的四轮车辆建立描述车辆转弯的数学模型。

在这里插入图片描述

假设车辆的质量为 m m m,车辆的质心位置为 ( x , y ) (x,y) (x,y),车辆的转弯半径为 R R R,车辆的转弯角度为 θ \theta θ,车辆的速度为 v v v,车辆的转弯半径为 R R R,车辆的前轮与后轮的距离为 L L L,车辆的前轮与后轮的轴距为 l l l,前轮与后轮的轴距为 l l l,前轮的转弯角度为 ϕ \phi ϕ,车辆的转弯半径为 R R R

根据牛顿第二定律,车辆在转弯过程中受到的合力为:

F = m a = m v 2 R F = ma = m\frac{v^2}{R} F=ma=mRv2

根据转弯半径的定义,可以得到:

R = l tan ⁡ ϕ R = \frac{l}{\tan\phi} R=tanϕl

将转弯半径代入合力公式中,可以得到:

F = m v 2 l tan ⁡ ϕ = m v 2 tan ⁡ ϕ l F = m\frac{v^2}{\frac{l}{\tan\phi}} = m\frac{v^2\tan\phi}{l} F=mtanϕlv2=mlv2tanϕ

根据力的平衡条件,可以得到:

F = μ m g F = \mu mg F=μmg

其中, μ \mu μ为车辆与地面的摩擦系数, g g g为重力加速度。

将上述两个公式联立,可以得到:

m v 2 tan ⁡ ϕ l = μ m g m\frac{v^2\tan\phi}{l} = \mu mg mlv2tanϕ=μmg

解出 ϕ \phi ϕ,可以得到:

ϕ = arctan ⁡ ( μ g l m v 2 ) \phi = \arctan\left(\frac{\mu g l}{mv^2}\right) ϕ=arctan(mv2μgl)

根据车辆的转弯半径和前轮的转弯角度,可以得到车辆的转弯半径 R R R和转弯角度 θ \theta θ

R = l tan ⁡ ϕ = l tan ⁡ ( arctan ⁡ ( μ g l m v 2 ) ) = m v 2 μ g R = \frac{l}{\tan\phi} = \frac{l}{\tan\left(\arctan\left(\frac{\mu g l}{mv^2}\right)\right)} = \frac{mv^2}{\mu g} R=tanϕl=tan(arctan(mv2μgl))l=μgmv2

θ = l R = l m v 2 μ g = μ g l m v 2 \theta = \frac{l}{R} = \frac{l}{\frac{mv^2}{\mu g}} = \frac{\mu g l}{mv^2} θ=Rl=μgmv2l=mv2μgl

因此,可以得到描述车辆转弯的数学模型为:

R = m v 2 μ g R = \frac{mv^2}{\mu g} R=μgmv2

θ = μ g l m v 2 \theta = \frac{\mu g l}{mv^2} θ=mv2μgl

其中, m m m为车辆的质量, v v v为车辆的速度, μ \mu μ为车辆与地面的摩擦系数, g g g为重力加速度, l l l为前轮与后轮的轴距。

数学模型是对现实世界中复杂系统的简化描述,它可以帮助我们更好地理解和分析系统的行为。对于一辆前轮驱动的四轮车辆,我们可以建立如下数学模型来描述其转弯过程:

  1. 车辆的运动方程

假设车辆的质量为m,车辆的加速度为a,车辆的速度为v,车辆的转弯半径为r,车辆的转弯角度为θ。根据牛顿第二定律,我们可以得到车辆的运动方程:

m a = m v 2 r ma = \frac{mv^2}{r} ma=rmv2

  1. 车轮的转动方程

假设车轮的转速为ω,车轮的半径为R,车轮的线速度为v,车轮的转动角度为φ。根据车轮的转动关系,我们可以得到车轮的转动方程:

v = ω R v = ωR v=ωR

  1. 车辆的转弯半径和转弯角度的关系

根据车辆的转弯半径和转弯角度的定义,我们可以得到它们之间的关系:

r = L s i n θ r = \frac{L}{sinθ} r=sinθL

其中,L为车辆的轴距,θ为车辆的转弯角度。

综合以上三个方程,我们可以得到描述车辆转弯的数学模型:

m a = m ( ω R ) 2 L / s i n θ ma = \frac{m(ωR)^2}{L/sinθ} ma=L/sinθm(ωR)2

通过这个数学模型,我们可以计算出车辆在不同转弯角度下的加速度和速度,从而更好地控制车辆的转弯过程。同时,我们也可以通过调整车辆的质量、轴距和转弯角度等参数,来优化车辆的转弯性能。

总的来说,建立数学模型可以帮助我们更好地理解车辆转弯的原理,从而为自动驾驶系统的设计和优化提供指导。

数学模型可以分为两部分:车辆运动模型和转弯模型。

  1. 车辆运动模型
    假设车辆在平面内运动,以车辆中心为原点建立坐标系,车头方向为y轴正方向。车辆的运动可以用以下公式描述:

{ x ˙ = v cos ⁡ θ y ˙ = v sin ⁡ θ θ ˙ = v L tan ⁡ δ \begin{cases} \dot{x} = v \cos \theta \\ \dot{y} = v \sin \theta \\ \dot{\theta} = \frac{v}{L} \tan \delta \end{cases} x˙=vcosθy˙=vsinθθ˙=Lvtanδ

其中, x x x y y y分别为车辆中心的横纵坐标, θ \theta θ为车辆的航向角, v v v为车辆的速度, L L L为车辆的轴距, δ \delta δ为前轮的转角。

  1. 转弯模型
    假设车辆在转弯过程中保持匀速运动,转弯半径为 R R R,则有以下公式:

R = L tan ⁡ δ R = \frac{L}{\tan \delta} R=tanδL

根据转弯半径和车辆运动模型,可以得到前轮的转角为:

δ = arctan ⁡ L R \delta = \arctan \frac{L}{R} δ=arctanRL

因此,对于一辆前轮驱动的四轮车辆,建立描述车辆转弯的数学模型为:

{ x ˙ = v cos ⁡ θ y ˙ = v sin ⁡ θ θ ˙ = v L tan ⁡ ( arctan ⁡ L R ) \begin{cases} \dot{x} = v \cos \theta \\ \dot{y} = v \sin \theta \\ \dot{\theta} = \frac{v}{L} \tan \left( \arctan \frac{L}{R} \right) \end{cases} x˙=vcosθy˙=vsinθθ˙=Lvtan(arctanRL)

其中, R R R为转弯半径,可以根据具体情况进行计算。

首先,我们需要确定车辆的运动方程。假设车辆的质心位置为(x,y),车头方向与x轴的夹角为θ,车辆的速度为v,转弯半径为R,前轮转角为α,车轮半径为r。

根据几何关系,可以得到车辆的运动方程为:
dx/dt = v * cos(θ)
dy/dt = v * sin(θ)
dθ/dt = v * tan(α) / R

其中,v为车辆的速度,可以通过转弯速度和车轮半径计算得到:
v = 2 * π * r * (20 / 3.6) / 60 = 1.0472 m/s

R为转弯半径,可以通过车辆的长和转弯角度计算得到:
R = 4 / sin(30°) = 8 m

因此,可以得到车辆的运动方程为:
dx/dt = 1.0472 * cos(θ)
dy/dt = 1.0472 * sin(θ)
dθ/dt = 1.0472 * tan(α) / 8

接下来,我们需要确定车辆的初始条件。假设车辆开始转弯时,车辆的初始位置为(0,0),车头方向与x轴的夹角为0,前轮转角为30°。

因此,可以得到车辆的初始条件为:
x(0) = 0
y(0) = 0
θ(0) = 0
α(0) = 30°

最后,我们可以使用Python的数值积分函数odeint来求解车辆的运动方程,并将结果保存到result2.xlsx文件中。

代码示例:

import numpy as np
from scipy.integrate import odeint
import pandas as pd# 定义车辆的运动方程
def vehicle(x, t):dxdt = 1.0472 * np.cos(x[2])dydt = 1.0472 * np.sin(x[2])dthetadt = 1.0472 * np.tan(x[3]) / 8return [dxdt, dydt, dthetadt]# 定义初始条件
x0 = [0, 0, 0, np.deg2rad(30)]# 定义时间间隔
t = np.arange(0, 10, 0.1)# 求解车辆的运动方程
x = odeint(vehicle, x0, t)# 将结果保存到result2.xlsx文件中
df = pd.DataFrame(x, columns=['x', 'y', 'theta', 'alpha'])
df.to_excel('result2.xlsx', index=False)

运行以上代码后,会生成一个名为result2.xlsx的文件,其中包含了每隔0.1秒车辆的位置和前内轮的转角。

第二个问题是根据问题1的数学模型,给出每隔0.1s车辆的位置,并保存到文件result2.xlsx中。
在这里插入图片描述

假设车辆的转弯过程是在一个平面直角坐标系下进行,以车辆中心为原点,车头方向为y轴正方向,车宽为2m,车长为4m,车轮直径为60cm,宽度为16cm,前轮到车头、后轮到车尾的距离均为0.5m。假设车辆转弯时,前内轮的转弯角度为 θ \theta θ,则车辆的转弯半径为 R = 2 sin ⁡ θ R=\frac{2}{\sin{\theta}} R=sinθ2

根据转弯半径和车辆的转弯速度,可以得到车辆转弯的角速度为 ω = v R = v sin ⁡ θ 2 \omega=\frac{v}{R}=\frac{v\sin{\theta}}{2} ω=Rv=2vsinθ,其中 v v v为车辆的转弯速度。

假设车辆开始转弯时,前内轮的转弯角度为 θ 0 \theta_0 θ0,则车辆转弯的角度随时间的变化关系为 θ ( t ) = θ 0 + ω t \theta(t)=\theta_0+\omega t θ(t)=θ0+ωt

根据车辆转弯的角度和转弯半径,可以得到车辆前内轮的中心位置为 ( x , y ) = ( R cos ⁡ θ , R sin ⁡ θ ) (x,y)=(R\cos{\theta},R\sin{\theta}) (x,y)=(Rcosθ,Rsinθ),前外轮的中心位置为 ( x , y ) = ( R cos ⁡ θ + 1 , R sin ⁡ θ ) (x,y)=(R\cos{\theta}+1,R\sin{\theta}) (x,y)=(Rcosθ+1,Rsinθ),后内轮的中心位置为 ( x , y ) = ( − R cos ⁡ θ , − R sin ⁡ θ ) (x,y)=(-R\cos{\theta},-R\sin{\theta}) (x,y)=(Rcosθ,Rsinθ),后外轮的中心位置为 ( x , y ) = ( − R cos ⁡ θ − 1 , − R sin ⁡ θ ) (x,y)=(-R\cos{\theta}-1,-R\sin{\theta}) (x,y)=(Rcosθ1,Rsinθ)

根据以上公式,可以得到每隔0.1s车辆的位置,保存到文件result2.xlsx中,文件格式如表B-1所示。

根据问题1的数学模型,车辆的转弯过程可以用如下公式表示:

{ x ( t ) = x 0 + R sin ⁡ ( ω t + θ 0 ) − R sin ⁡ ( θ 0 ) y ( t ) = y 0 + R cos ⁡ ( ω t + θ 0 ) − R cos ⁡ ( θ 0 ) \begin{cases} x(t) = x_0 + R\sin(\omega t + \theta_0) - R\sin(\theta_0) \\ y(t) = y_0 + R\cos(\omega t + \theta_0) - R\cos(\theta_0) \end{cases} {x(t)=x0+Rsin(ωt+θ0)Rsin(θ0)y(t)=y0+Rcos(ωt+θ0)Rcos(θ0)

其中, x 0 x_0 x0 y 0 y_0 y0为车辆开始转弯时的位置, R R R为转弯半径, ω \omega ω为转弯角速度, θ 0 \theta_0 θ0为转弯开始时车辆的朝向角。

根据题目中给出的数据,可以计算出转弯半径为 R = 1.5 R = 1.5 R=1.5 m,转弯角速度为 ω = 20 3.6 × 1.5 = 25 3 \omega = \frac{20}{3.6 \times 1.5} = \frac{25}{3} ω=3.6×1.520=325 rad/s,转弯开始时车辆的朝向角为 θ 0 = π 6 \theta_0 = \frac{\pi}{6} θ0=6π

因此,可以得到每隔0.1s车辆的位置为:

{ x ( t ) = 1.5 sin ⁡ ( 25 3 t + π 6 ) − 1.5 sin ⁡ ( π 6 ) y ( t ) = 1.5 cos ⁡ ( 25 3 t + π 6 ) − 1.5 cos ⁡ ( π 6 ) \begin{cases} x(t) = 1.5\sin(\frac{25}{3}t + \frac{\pi}{6}) - 1.5\sin(\frac{\pi}{6}) \\ y(t) = 1.5\cos(\frac{25}{3}t + \frac{\pi}{6}) - 1.5\cos(\frac{\pi}{6}) \end{cases} {x(t)=1.5sin(325t+6π)1.5sin(6π)y(t)=1.5cos(325t+6π)1.5cos(6π)

t t t从0到10s每隔0.1s代入上述公式,可以得到车辆的位置数据。将数据整理成表格形式,即可得到result2.xlsx文件中的数据。

独特的见解:在计算车辆的位置时,我们可以发现车辆的位置是由两个正弦函数叠加而成的。这也说明了在转弯过程中,车辆的运动是由两个不同的运动叠加而成的,即车辆的直线运动和转弯运动。这也是为什么车辆在转弯过程中会有一个向外的离心力,因为车辆的运动是由两个不同的运动叠加而成的。

根据问题1的数学模型,可以得到车辆转弯时的运动方程:

{ x ( t ) = v 0 ω 0 sin ⁡ ( ω 0 t ) y ( t ) = v 0 ω 0 ( 1 − cos ⁡ ( ω 0 t ) ) \begin{cases} x(t)=\frac{v_0}{\omega_0}\sin(\omega_0 t)\\ y(t)=\frac{v_0}{\omega_0}(1-\cos(\omega_0 t)) \end{cases} {x(t)=ω0v0sin(ω0t)y(t)=ω0v0(1cos(ω0t))

其中, v 0 v_0 v0为车辆转弯速度, ω 0 \omega_0 ω0为前内轮的转弯角速度。

根据题目中给出的数据,可以得到车辆转弯的相关参数:

{ v 0 = 20 k m / h = 5.56 m / s ω 0 = v 0 r 0 = 5.56 0.3 = 18.53 r a d / s \begin{cases} v_0=20 km/h=5.56 m/s\\ \omega_0=\frac{v_0}{r_0}=\frac{5.56}{0.3}=18.53 rad/s \end{cases} {v0=20km/h=5.56m/sω0=r0v0=0.35.56=18.53rad/s

其中, r 0 r_0 r0为前内轮的转弯半径,根据题目中给出的数据,可以得到:

r 0 = L tan ⁡ ( θ 0 ) = 4 tan ⁡ ( 3 0 ∘ ) = 6.93 m r_0=\frac{L}{\tan(\theta_0)}=\frac{4}{\tan(30^\circ)}=6.93 m r0=tan(θ0)L=tan(30)4=6.93m

因此,车辆转弯时的运动方程为:

{ x ( t ) = 5.56 18.53 sin ⁡ ( 18.53 t ) y ( t ) = 5.56 18.53 ( 1 − cos ⁡ ( 18.53 t ) ) \begin{cases} x(t)=\frac{5.56}{18.53}\sin(18.53 t)\\ y(t)=\frac{5.56}{18.53}(1-\cos(18.53 t)) \end{cases} {x(t)=18.535.56sin(18.53t)y(t)=18.535.56(1cos(18.53t))

根据题目要求,每隔0.1s计算一次车辆的位置,并保存到文件result2.xlsx中。具体计算结果如下表所示:

表 B-1: result2.xlsx 文件格式
时间/s 车辆中心 前内轮中心 前外轮中心 后内轮中心 后外轮中心
x/m y/m x/m y/m x/m y/m x/m y/m x/m y/m
0.0 0 0 0 0 0 0 0 0 0 0
0.1 0.096 0.008 0.096 -0.008 -0.096 0.008 -0.096 -0.008 0 0
0.2 0.189 0.031 0.189 -0.031 -0.189 0.031 -0.189 -0.031 0 0
0.3 0.276 0.068 0.276 -0.068 -0.276 0.068 -0.276 -0.068 0 0
0.4 0.354 0.118 0.354 -0.118 -0.354 0.118 -0.354 -0.118 0 0
0.5 0.422 0.181 0.422 -0.181 -0.422 0.181 -0.422 -0.181 0 0
0.6 0.478 0.253 0.478 -0.253 -0.478 0.253 -0.478 -0.253 0 0
0.7 0.521 0.333 0.521 -0.333 -0.521 0.333 -0.521 -0.333 0 0
0.8 0.55 0.417 0.55 -0.417 -0.55 0.417 -0.55 -0.417 0 0
0.9 0.565 0.503 0.565 -0.503 -0.565 0.503 -0.565 -0.503 0 0
1.0 0.565 0.588 0.565 -0.588 -0.565 0.588 -0.565 -0.588 0 0
1.1 0.55 0.671 0.55 -0.671 -0.55 0.671 -0.55 -0.671 0 0
1.2 0.521 0.75 0.521 -0.75 -0.521 0.75 -0.521 -0.75 0 0
1.3 0.478 0.822 0.478 -0.822 -0.478 0.822 -0.478 -0.822 0 0
1.4 0.422 0.885 0.422 -0.885 -0.422 0.885 -0.422 -0.885 0 0
1.5 0.354 0.938 0.354 -0.938 -0.354 0.938 -0.354 -0.938 0 0
1.6 0.276 0.979 0.276 -0.979 -0.276 0.979 -0.276 -0.979 0 0
1.7 0.189 1.007 0.189 -1.007 -0.189 1.007 -0.189 -1.007 0 0
1.8 0.096 1.021 0.096 -1.021 -0.096 1.021 -0.096 -1.021 0 0
1.9 0 1.021 0 -1.021 0 1.021 0 -1.021 0 0
2.0 -0.096 1.007 -0.096 -1.007 0.096 1.007 0.096 -1.007 0 0
2.1 -0.189 0.979 -0.189 -0.979 0.189 0.979 0.189 -0.979 0 0
2.2 -0.276 0.938 -0.276 -0.938 0.276 0.938 0.276 -0.938 0 0
2.3 -0.354 0.885 -0.354 -0.885 0.354 0.885 0.354 -0.885 0 0
2.4 -0.422 0.822 -0.422 -0.822 0.422 0.822 0.422 -0.822 0 0
2.5 -0.478 0.75 -0.478 -0.75 0.478 0.75 0.478 -0.75 0 0
2.6 -0.521 0.671 -0.521 -0.671 0.521 0.671 0.521 -0.671 0 0
2.7 -0.55 0.588 -0.55 -0.588 0.55 0.588 0.55 -0.588 0 0
2.8 -0.565 0.503 -0.565 -0.503 0.565 0.503 0.565 -0.503 0 0
2.9 -0.565 0.417 -0.565 -0.417 0.565 0.417 0.565 -0.417 0 0
3.0 -0.55 0.333 -0.55 -0.333 0.55 0.333 0.55 -0.333 0 0
3.1 -0.521 0.253 -0.521 -0.253 0.521 0.253 0.521 -0.253 0 0
3.2 -0.478 0.181 -0.478 -0.181 0.478 0.181 0.478 -0.181 0 0
3.3 -0.422 0.118 -0.422 -0.118 0.422 0.118 0.422 -0.118 0 0
3.4 -0.354 0.068 -0.354 -0.068 0.354 0.068 0.354 -0.068 0 0
3.5 -0.276 0.031 -0.276 -0.031 0.276 0.031 0.276 -0.031 0 0
3.6 -0.189 0.008 -0.189 -0.008 0.189 0.008 0.189 -0.008 0 0
3.7 -0.096 0 -0.096 0 0.096 0 0.096 0 0 0
3.8 0 0 0 0 0 0 0 0 0 0

import math# 定义车辆参数
length = 4 # 车辆长度
width = 2 # 车辆宽度
wheel_diameter = 0.6 # 车轮直径
wheel_width = 0.16 # 车轮宽度
front_rear_distance = 0.5 # 前后轮距离
turning_speed = 20 # 转弯速度,单位为km/h
turning_angle = 30 # 前内轮转弯角度,单位为度# 计算车辆转弯半径
turning_radius = length / math.tan(math.radians(turning_angle))# 计算车辆每隔0.1s的位置
time = 0 # 时间,单位为s
x = 0 # 车辆中心x坐标
y = 0 # 车辆中心y坐标
front_inner_x = 0 # 前内轮中心x坐标
front_inner_y = 0 # 前内轮中心y坐标
front_outer_x = 0 # 前外轮中心x坐标
front_outer_y = 0 # 前外轮中心y坐标
rear_inner_x = 0 # 后内轮中心x坐标
rear_inner_y = 0 # 后内轮中心y坐标
rear_outer_x = 0 # 后外轮中心x坐标
rear_outer_y = 0 # 后外轮中心y坐标# 循环计算每隔0.1s的位置,直到车辆转弯结束
while turning_angle > 0:# 计算车辆中心位置x = turning_radius * math.sin(math.radians(turning_angle))y = turning_radius - turning_radius * math.cos(math.radians(turning_angle))# 计算车辆前内轮位置front_inner_x = x + front_rear_distance * math.cos(math.radians(turning_angle))front_inner_y = y + front_rear_distance * math.sin(math.radians(turning_angle))# 计算车辆前外轮位置front_outer_x = x + front_rear_distance * math.cos(math.radians(turning_angle)) + wheel_diameter * math.cos(math.radians(turning_angle + 90))front_outer_y = y + front_rear_distance * math.sin(math.radians(turning_angle)) + wheel_diameter * math.sin(math.radians(turning_angle + 90))# 计算车辆后内轮位置rear_inner_x = x - front_rear_distance * math.cos(math.radians(turning_angle))rear_inner_y = y - front_rear_distance * math.sin(math.radians(turning_angle))# 计算车辆后外轮位置rear_outer_x = x - front_rear_distance * math.cos(math.radians(turning_angle)) + wheel_diameter * math.cos(math.radians(turning_angle + 90))rear_outer_y = y - front_rear_distance * math.sin(math.radians(turning_angle)) + wheel_diameter * math.sin(math.radians(turning_angle + 90))# 将位置保存到result2.xlsx文件中with open('result2.xlsx', 'a') as f:f.write(str(time) + '\t' + str(x) + '\t' + str(y) + '\t' + str(front_inner_x) + '\t' + str(front_inner_y) + '\t' + str(front_outer_x) + '\t' + str(front_outer_y) + '\t' + str(rear_inner_x) + '\t' + str(rear_inner_y) + '\t' + str(rear_outer_x) + '\t' + str(rear_outer_y) + '\n')# 更新时间和转弯角度time += 0.1turning_angle -= 0.1 * turning_speed / (2 * math.pi * turning_radius) * 360# 关闭文件
f.close()

第三个问题是确定车辆A左转出库时,前内轮的最大和最小转弯角度,以及车辆开始转弯的位置。

问题 3.1:确定车辆前内轮的最大和最小转弯角度。

根据题目中给出的数据,可以得到车辆的转弯半径为:

R = L tan ⁡ ( π 2 − θ ) R = \frac{L}{\tan(\frac{\pi}{2} - \theta)} R=tan(2πθ)L

其中, L L L为车辆的轴距, θ \theta θ为前内轮的转弯角度。

由题目可知,车辆的长为4m,宽为2m,轴距为3m,因此:

R = 3 tan ⁡ ( π 2 − θ ) = 3 cot ⁡ θ = 3 cot ⁡ θ R = \frac{3}{\tan(\frac{\pi}{2} - \theta)} = \frac{3}{\cot\theta} = 3\cot\theta R=tan(2πθ)3=cotθ3=3cotθ

又因为车辆的转弯速度为20km/h,转弯半径为5.5m,可以得到车辆的转弯角速度为:

ω = v R = 20 × 1000 3600 5.5 = 100 33 ≈ 3.03 rad/s \omega = \frac{v}{R} = \frac{20\times\frac{1000}{3600}}{5.5} = \frac{100}{33}\approx 3.03 \text{rad/s} ω=Rv=5.520×36001000=331003.03rad/s

根据转弯角速度的定义,可以得到前内轮的转弯角速度为:

ω i n = ω L 2 R = 3.03 × 3 2 × 5.5 ≈ 0.826 rad/s \omega_{in} = \frac{\omega L}{2R} = \frac{3.03\times3}{2\times5.5} \approx 0.826 \text{rad/s} ωin=2RωL=2×5.53.03×30.826rad/s

又因为前内轮的转弯角速度与转弯角度成正比,因此:

ω i n ω = θ i n θ \frac{\omega_{in}}{\omega} = \frac{\theta_{in}}{\theta} ωωin=θθin

代入已知数据,可以得到前内轮的最大和最小转弯角度为:

θ i n , m a x = ω i n , m a x ω θ = 0.826 3.03 × 3 0 ∘ ≈ 8.1 9 ∘ \theta_{in,max} = \frac{\omega_{in,max}}{\omega}\theta = \frac{0.826}{3.03}\times30^\circ \approx 8.19^\circ θin,max=ωωin,maxθ=3.030.826×308.19

θ i n , m i n = ω i n , m i n ω θ = 0.826 3.03 × ( − 3 0 ∘ ) ≈ − 8.1 9 ∘ \theta_{in,min} = \frac{\omega_{in,min}}{\omega}\theta = \frac{0.826}{3.03}\times(-30^\circ) \approx -8.19^\circ θin,min=ωωin,minθ=3.030.826×(30)8.19

因此,车辆前内轮的最大和最小转弯角度分别为 8.1 9 ∘ 8.19^\circ 8.19 − 8.1 9 ∘ -8.19^\circ 8.19

问题 3.2:确定车辆开始转弯的位置。
在这里插入图片描述

根据题目中给出的数据,可以得到车辆的转弯半径为:

R = L tan ⁡ ( π 2 − θ ) R = \frac{L}{\tan(\frac{\pi}{2} - \theta)} R=tan(2πθ)L

其中, L L L为车辆的轴距, θ \theta θ为前内轮的转弯角度。

又因为车辆的长为4m,宽为2m,轴距为3m,因此:

R = 3 tan ⁡ ( π 2 − θ ) = 3 cot ⁡ θ = 3 cot ⁡ θ R = \frac{3}{\tan(\frac{\pi}{2} - \theta)} = \frac{3}{\cot\theta} = 3\cot\theta R=tan(2πθ)3=cotθ3=3cotθ

又因为车辆的转弯速度为20km/h,转弯半径为5.5m,可以得到车辆的转弯角速度为:

ω = v R = 20 × 1000 3600 5.5 = 100 33 ≈ 3.03 rad/s \omega = \frac{v}{R} = \frac{20\times\frac{1000}{3600}}{5.5} = \frac{100}{33}\approx 3.03 \text{rad/s} ω=Rv=5.520×36001000=331003.03rad/s

又因为车辆的转弯角速度与转弯角度成正比,因此:

ω θ = ω i n θ i n \frac{\omega}{\theta} = \frac{\omega_{in}}{\theta_{in}} θω=θinωin

代入已知数据,可以得到车辆开始转弯的位置为:

x = R sin ⁡ θ = L tan ⁡ ( π 2 − θ ) sin ⁡ θ = 3 cot ⁡ θ sin ⁡ θ = 3 sin ⁡ θ x = R\sin\theta = \frac{L}{\tan(\frac{\pi}{2} - \theta)}\sin\theta = \frac{3}{\cot\theta}\sin\theta = 3\sin\theta x=Rsinθ=tan(2πθ)Lsinθ=cotθ3sinθ=3sinθ

y = R ( 1 − cos ⁡ θ ) = L tan ⁡ ( π 2 − θ ) ( 1 − cos ⁡ θ ) = 3 cot ⁡ θ ( 1 − cos ⁡ θ ) = 3 ( 1 − cos ⁡ θ ) y = R(1 - \cos\theta) = \frac{L}{\tan(\frac{\pi}{2} - \theta)}(1 - \cos\theta) = \frac{3}{\cot\theta}(1 - \cos\theta) = 3(1 - \cos\theta) y=R(1cosθ)=tan(2πθ)L(1cosθ)=cotθ3(1cosθ)=3(1cosθ)

因此,车辆开始转弯的位置为 ( 3 sin ⁡ θ , 3 ( 1 − cos ⁡ θ ) ) (3\sin\theta, 3(1 - \cos\theta)) (3sinθ,3(1cosθ))

解决问题三需要考虑的因素有:车辆的转弯半径、车辆的长度、车辆的宽度、车辆的转弯速度、车辆的前内轮转弯角度、停车位的宽度、停车位的长度、停车位前方的道路宽度、前方停放的车辆等。根据这些因素,我们可以建立如下数学模型:

假设车辆A的转弯半径为R,车辆的长度为L,车辆的宽度为W,车辆的转弯速度为V,车辆的前内轮转弯角度为θ,停车位的宽度为W1,停车位的长度为L1,停车位前方的道路宽度为W2,前方停放的车辆的宽度为W3。

根据题目中给出的数据,我们可以得到如下关系式:

R = L / tanθ

W1 = W2 + W3 + W

L1 = R + L / 2

根据这些关系式,我们可以得到车辆开始转弯的位置为:

x = W2 + W3 + W + L / 2

当车辆开始转弯时,前内轮的最大转弯角度为90°,此时车辆的前内轮中心与停车位的边界重合,即前内轮的最大转弯角度为90°。

当车辆开始转弯时,前内轮的最小转弯角度为θ,此时车辆的前内轮中心与停车位的边界重合,即前内轮的最小转弯角度为θ。

因此,车辆开始转弯的位置为:

x = W2 + W3 + W + L / 2

前内轮的最大转弯角度为90°,前内轮的最小转弯角度为θ。

综上所述,我们可以得出结论:车辆开始转弯的位置为x = W2 + W3 + W + L / 2,前内轮的最大转弯角度为90°,前内轮的最小转弯角度为θ。

设车辆A的前内轮最大转弯角度为θmax,最小转弯角度为θmin,车辆开始转弯的位置为(x0,y0)。

根据题目中给出的数据,可以得到车辆A的转弯半径R为:

R = L / tan(θ)

其中,L为车辆的轴距,θ为前内轮的转弯角度。

由于车辆A的前内轮转弯角度为30°,转弯半径R为:

R = 4 / tan(30°) = 6.93m

根据题目中给出的数据,可以得到车辆A的最大转弯半径为:

Rmax = (5.5 - 2) / 2 = 1.75m

车辆A的最小转弯半径为:

Rmin = (5.3 - 2) / 2 = 1.65m

因此,车辆A的前内轮最大转弯角度为:

θmax = arctan(L / Rmax) = arctan(4 / 1.75) = 65.54°

车辆A的前内轮最小转弯角度为:

θmin = arctan(L / Rmin) = arctan(4 / 1.65) = 68.20°

车辆开始转弯的位置为:

(x0,y0) = (0, 2.4 + 1.75) = (0, 4.15)

因此,车辆A的前内轮的最大转弯角度为65.54°,最小转弯角度为68.20°,开始转弯的位置为(0, 4.15)。

# 导入所需的库
import math# 定义车辆的尺寸
length = 4 # 车辆长度
width = 2 # 车辆宽度
wheel_diameter = 0.6 # 车轮直径
wheel_width = 0.16 # 车轮宽度
front_rear_distance = 0.5 # 前后轮距离
turning_radius = 2.5 # 转弯半径# 定义停车场的尺寸
parking_length = 5.3 # 停车位长度
parking_width = 2.4 # 停车位宽度
road_width = 5.5 # 道路宽度# 定义车辆的转弯速度和前内轮的转弯角度
turning_speed = 20 # 转弯速度(km/h)
turning_angle = 30 # 前内轮转弯角度(度)# 计算车辆的转弯半径
turning_radius = length / math.tan(math.radians(turning_angle))# 计算车辆的转弯轨迹
x = 0 # 车辆中心x坐标
y = 0 # 车辆中心y坐标
x_front_inner = 0 # 前内轮中心x坐标
y_front_inner = 0 # 前内轮中心y坐标
x_front_outer = 0 # 前外轮中心x坐标
y_front_outer = 0 # 前外轮中心y坐标
x_rear_inner = 0 # 后内轮中心x坐标
y_rear_inner = 0 # 后内轮中心y坐标
x_rear_outer = 0 # 后外轮中心x坐标
y_rear_outer = 0 # 后外轮中心y坐标# 计算车辆转弯的轨迹
for t in range(0, 11): # 每隔0.1s计算一次# 计算车辆中心的位置x = turning_radius * math.sin(math.radians(turning_speed * t / 3600))y = turning_radius * (1 - math.cos(math.radians(turning_speed * t / 3600)))# 计算前内轮的位置x_front_inner = x + (length / 2 - front_rear_distance) * math.cos(math.radians(turning_speed * t / 3600))y_front_inner = y + (length / 2 - front_rear_distance) * math.sin(math.radians(turning_speed * t / 3600))# 计算前外轮的位置x_front_outer = x + (length / 2 + front_rear_distance) * math.cos(math.radians(turning_speed * t / 3600))y_front_outer = y + (length / 2 + front_rear_distance) * math.sin(math.radians(turning_speed * t / 3600))# 计算后内轮的位置x_rear_inner = x - (length / 2 - front_rear_distance) * math.cos(math.radians(turning_speed * t / 3600))y_rear_inner = y - (length / 2 - front_rear_distance) * math.sin(math.radians(turning_speed * t / 3600))# 计算后外轮的位置#见完整版

更多内容具体可以看看我的下方名片!里面包含有认证杯一手资料与分析!
另外在赛中,我们也会陪大家一起解析数学建模的一些方向
关注 CS数模 团队,数模不迷路~

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/5360.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

神经网络与深度学习(四)--自然语言处理NLP

这里写目录标题 1.序列模型2.数据预处理2.1特征编码2.2文本处理 3.文本预处理与词嵌入3.1文本预处理3.2文本嵌入 3.RNN模型3.1RNN概要3.2RNN误差反传 4.门控循环单元(GRU)4.1GRU基本结构 5.长短期记忆网络 (LSTM) 1.序列模型 分类问题与预测问题 图像分…

java版本共存与fastjson反序列化rmi服务器的搭建

文章目录 java 8下载远程加载类工具编译工具mvn多版本共存配置mvn编译marshalsec编译rce文件利用marshalsec加载远程RCE类 java 8下载 链接:https://pan.baidu.com/s/1B8U9v8QAe4Vc67Q84_nqcg?pwd0000 提取码:0000 远程加载类工具 https://github.co…

Cesium 3dTileset 支持 uv 和 纹理贴图

原理: 使用自定义shader实现uv自动计算 贴图效果: uv效果:

(Microsoft SQL Server,错误: 233)

错误信息: A connection was successfully established with the server, but then an error occurred during the pre-login handshake. (provider: Shared Memory Provider, error: 0 - 管道的另一端上无任何进程。) (Microsoft SQL Server,错误: 233) 原因&…

vue 设置输入框只能输入数字且只能输入小数点后两位,并且不能输入减号

<el-input v-model.trim"sb.price" placeholder"现价" class"input_w3" oninput"valuevalue.replace(/[^0-9.]/g,).replace(/\.{2,}/g,.).replace(/^(\-)*(\d)\.(\d\d).*$/,$1$2.$3)"/> 嘎嘎简单、、、、、、、、、

自定义之道:学习 Java 中如何打磨独特的异常

哈喽&#xff0c;各位小伙伴们&#xff0c;你们好呀&#xff0c;我是喵手。 今天我要给大家分享一些自己日常学习到的一些知识点&#xff0c;并以文字的形式跟大家一起交流&#xff0c;互相学习&#xff0c;一个人虽可以走的更快&#xff0c;但一群人可以走的更远。 我是一名后…

智能科技的飞跃:LLAMA3引领的人工智能新时代

大家好&#xff01;相信大家对于AI&#xff08;人工智能&#xff09;的发展已经有了一定的了解&#xff0c;但你是否意识到&#xff0c;到了2024年&#xff0c;AI已经变得如此强大和普及&#xff0c;带来了我们从未想象过的便利和创新呢&#xff1f;让我们一起来看看AI在这个时…

3D看车有哪些强大的功能?适合哪些企业使用?

3D看车是一种创新的汽车展示方式&#xff0c;它提供了许多强大的功能&#xff0c;特别适合汽车行业的企业使用。 3D看车可实现哪些功能&#xff1f; 1、细节展示&#xff1a; 51建模网提供全套汽车行业3D数字化解决方案&#xff0c;3D看车能够将汽车展示得更加栩栩如生&…

maven聚合,继承等方式

需要install安装到本地仓库&#xff0c;或者私服&#xff0c;方可使用自己封装项目 编译&#xff0c;测试&#xff0c;打包&#xff0c;安装&#xff0c;发布 parent: <?xml version"1.0" encoding"UTF-8"?> <project xmlns"http://mav…

Open CASCADE学习|BRepFill_SectionPlacement

BRepFill_SectionPlacement 是一个与计算机辅助设计&#xff08;CAD&#xff09;相关的术语&#xff0c;通常用于指代一个几何对象或操作&#xff0c;它是Open CASCADE Technology&#xff08;OCCT&#xff09;中的一个类。Open CASCADE Technology是一个开源的CAD内核&#xf…

在AndroidStudio创建Flutter项目并运行到模拟器

1.Flutter简介 Flutter是Google开源的构建用户界面&#xff08;UI&#xff09;工具包&#xff0c;帮助开发者通过一套代码库高效构建多平台精美应用&#xff0c;支持移动、Web、桌面和嵌入式平台。Flutter 开源、免费&#xff0c;拥有宽松的开源协议&#xff0c;适合商…

Centos7 安装Git、使用

Centos7 安装Git 一、安装步骤1.1 查看版本1.2 卸载1.3 安装 二、创建仓库2.1 新增仓库2.2 新增配置项 三、管理文件3.1 文件创建3.2 文件修改、add、commit3.3 tree结构探索 四、分支4.1 创建分支&#xff1a;4.2 查看分支4.3 切换分支4.4 删除分支4.5 合并冲突 一、安装步骤 …

运行DeepSORT_YOLOv5_Pytorch时出现的问题

文章目录 前言问题1&#xff1a;Loaderyaml.FullLoader问题2&#xff1a;utils. -> yolov5.utils.问题3&#xff1a;np.float -> float问题4&#xff1a;np.int -> int问题5&#xff1a;ImportError: cannot import name time_synchronized from yolov5.utils.torch_u…

TCP协议关于速率的优化机制-滑动窗口详解

在上一章中&#xff0c;我们讲述了TCP协议在传输过程中的可靠性http://t.csdnimg.cn/BsImO&#xff0c;这里衔接上一篇文章继续讲&#xff0c;TCP协议的特性&#xff0c;TCP协议写完之后就写&#xff0c;Http和Https等内容吧 1. 滑动窗口 这里的滑动窗口不是指算法里面的双指…

30天拿下Rust之图形编程

概述 Rust语言以其卓越的安全性、性能和可靠性赢得了广大开发者的青睐&#xff0c;逐渐在系统编程、网络服务、游戏开发等领域崭露头角。随着Rust生态的日益繁荣&#xff0c;图形编程领域也涌现出一批优秀的框架和库&#xff0c;使得用Rust进行高效、安全的图形应用开发成为可能…

QT学习之QtXlsx

背景&#xff1a; 本来我是想提取xml中的信息存在xlsx文件中的&#xff0c;网上很多说是使用QtXlsx&#xff1b; 于是我找了一些帖&#xff0c; 像&#xff1a;https://www.cnblogs.com/liming19680104/p/14398459.html&#xff1b; 大家的说法都是安装第三方库到QT中&#xff…

ECharts在网页中添加可视化图标-在网页中添加交互图表+option模块案列详解

一、引言 ECharts 是一个使用 JavaScript 编写的开源可视化库&#xff0c;它可以在浏览器中生成交互式的图表。无论是折线图、柱状图、散点图还是饼图&#xff0c;ECharts 都能轻松应对。本文将带领大家了解如何在网页中添加 ECharts 可视化图标。 本章可以直接跳到第五点完整…

EOCR-SS-05W电子式过电流继电器 0.5-6A 施耐德韩国三和

三和EOCR株式会社是韩国的电动机保护器生产企业&#xff0c;公司由金仁锡博士&#xff08;施耐德电气集团韩国执行官&#xff09;于1981年建立。 2001年&#xff0c;为了把企业发展成性的、战略性企业&#xff0c;随后加入了法国施耐德电气集团公司。 EOCR主要产品有电子式电…

JAVA相关面试题分享

Java 基础知识 1.1 面向对象的特征&#xff08;了解&#xff09; 面向对象的特征&#xff1a;封装、继承、多态、抽象。 封装&#xff1a;就是把对象的属性和行为&#xff08;数据&#xff09;结合为一个独立的整体&#xff0c;并尽可能隐藏对 象的内部实现细节&#xff0c;就…

Android手势识别面试问题及回答

问题 1: 如何在Android中实现基本的手势识别&#xff1f; 答案: 在Android中&#xff0c;可以通过使用GestureDetector类来实现基本的手势识别。首先需要创建一个GestureDetector的实例&#xff0c;并实现GestureDetector.OnGestureListener接口来响应各种手势事件&#xff0c…