S-I 模型用于研究不同情景下传染病的传播。在本题中,假设在整个考虑期间,我们研究的是一个规模恒定为 NN 的封闭社区中的均匀人群。人群被分为两组:一组是患有该疾病的已感染个体,另一组是健康且易感的个体,后者在与已感染个体进行某种形式的接触时会被感染。分别用 xx 和 yy 表示在时间 tt 时的感染个体数量和健康且易感的个体数量,我们可以使用以下微分方程对这种情况进行建模
dxdt=kxy,(I)\frac{\text{d}x}{\text{d}t} = kxy, \quad\quad\quad\quad \text{(I)}
其中 kk 是一个正常数。

  1. (i)

    通过用 NN 和 xx 表示 yy,证明方程 (I) 可以改写为关于 xx 和 tt 的一阶微分方程的形式。已知当 t=0t = 0 时 x=αx = \alpha,求解该微分方程,并证明其解可以表示为
    x=A1+(B−1)e−Nkt,x = \frac{A}{1 + (B - 1)\text{e}^{-Nkt}},
    其中 AA 和 BB 是用 NN 和 \alpha$ 表示的常数。

  2. (ii)

    该模型使用 2003 年新加坡 SARS 疫情的部分数据进行了测试,其中 α=1\alpha = 1,N=206N = 206,kk 的估计值为 8.1835×10−48.1835 \times 10^{-4},时间以天为单位。

    写出相应解曲线的方程,并画出与本情境相关的部分曲线。(你的草图应在坐标轴上进行适当的标注。)说明当 tt 很大时 xx 的变化情况。

先不看答案,让 AI 诊断你卡在哪里?

它不会直接告诉你答案,而是问你卡在哪一步看看 AI 是怎么陪你一步步弄懂这道题的

完整解法

7(i)

1. Express $y$ in terms of $N$ and $x$, and rewrite the differential equation.

Since the population size is constant at NN, the total number of individuals is the sum of infected and healthy individuals:
x+y=Nx + y = N
y=N−xy = N - x
Substituting this into equation (I):
dxdt=kx(N−x)\frac{\text{d}x}{\text{d}t} = kx(N - x)

2. Separate the variables to prepare for integration.

1x(N−x)dxdt=k\frac{1}{x(N - x)} \frac{\text{d}x}{\text{d}t} = k
∫1x(N−x) dx=∫k dt\int \frac{1}{x(N - x)} \,\text{d}x = \int k \,\text{d}t

3. Use partial fractions to split the integrand.

Let 1x(N−x)=Px+QN−x\frac{1}{x(N - x)} = \frac{P}{x} + \frac{Q}{N - x}
1=P(N−x)+Qx1 = P(N - x) + Qx
Comparing coefficients or substituting values for xx gives P=1NP = \frac{1}{N} and Q=1NQ = \frac{1}{N}.
1N∫(1x+1N−x) dx=∫k dt\frac{1}{N} \int \left( \frac{1}{x} + \frac{1}{N - x} \right) \,\text{d}x = \int k \,\text{d}t

4. Integrate both sides and combine logarithmic terms.

1N(ln⁡x−ln⁡(N−x))=kt+C\frac{1}{N} \left( \ln x - \ln(N - x) \right) = kt + C
Note: Modulus signs are not strictly necessary as 0<x<N0 < x < N (the number of infected is positive and cannot exceed the total population).
1Nln⁡(xN−x)=kt+C\frac{1}{N} \ln \left( \frac{x}{N - x} \right) = kt + C
ln⁡(xN−x)=Nkt+NC\ln \left( \frac{x}{N - x} \right) = Nkt + NC

5. Convert to exponential form to isolate the fraction.

xN−x=eNkt+NC\frac{x}{N - x} = \text{e}^{Nkt + NC}
xN−x=DeNkt\frac{x}{N - x} = D\text{e}^{Nkt}
where D=eNCD = \text{e}^{NC} is an arbitrary constant.

6. Substitute the initial condition to find the constant $D$.

Given that x=αx = \alpha when t=0t = 0:
αN−α=De0\frac{\alpha}{N - \alpha} = D\text{e}^{0}
D=αN−αD = \frac{\alpha}{N - \alpha}
Substituting DD back into the equation:
xN−x=(αN−α)eNkt\frac{x}{N - x} = \left( \frac{\alpha}{N - \alpha} \right) \text{e}^{Nkt}

7. Rearrange the equation to make $x$ the subject in the required format.

Taking the reciprocal on both sides:
N−xx=(N−αα)e−Nkt\frac{N - x}{x} = \left( \frac{N - \alpha}{\alpha} \right) \text{e}^{-Nkt}
Nx−1=(Nα−1)e−Nkt\frac{N}{x} - 1 = \left( \frac{N}{\alpha} - 1 \right) \text{e}^{-Nkt}
Nx=1+(Nα−1)e−Nkt\frac{N}{x} = 1 + \left( \frac{N}{\alpha} - 1 \right) \text{e}^{-Nkt}
x=N1+(Nα−1)e−Nktx = \frac{N}{1 + \left( \frac{N}{\alpha} - 1 \right)\text{e}^{-Nkt}}

8. State the values of constants $A$ and $B$.

Comparing with x=A1+(B−1)e−Nktx = \frac{A}{1 + (B - 1)\text{e}^{-Nkt}}:
A=NA = N
B−1=Nα−1  ⟹  B=NαB - 1 = \frac{N}{\alpha} - 1 \implies B = \frac{N}{\alpha}

答案

dxdt=kx(N−x)\frac{\text{d}x}{\text{d}t} = kx(N - x)
x=N1+(Nα−1)e−Nktx = \frac{N}{1 + \left(\frac{N}{\alpha} - 1\right)\text{e}^{-Nkt}}
A=NA = N, B=NαB = \frac{N}{\alpha}

7(ii)

1. Determine the values of constants A, B, and the exponent coefficient from the given data.

From part (i), A=NA = N and B=NαB = \frac{N}{\alpha}.
Given α=1\alpha = 1 and N=206N = 206:
A=206A = 206
B=2061=206B = \frac{206}{1} = 206
The coefficient for tt is Nk=206×(8.1835×10−4)=0.1685801≈0.169Nk = 206 \times (8.1835 \times 10^{-4}) = 0.1685801 \approx 0.169 (to 3 s.f.).
The equation of the solution curve is:
x=2061+205e−0.169tx = \frac{206}{1 + 205\text{e}^{-0.169t}}

2. Identify key features to sketch the relevant part of the solution curve.

The context considers time t≥0t \ge 0.
At t=0t = 0, x=1x = 1. The curve starts at (0,1)(0, 1).
As t→∞t \to \infty, e−0.169t→0\text{e}^{-0.169t} \to 0, so x→206x \to 206. There is a horizontal asymptote at x=206x = 206.
Since xx models cumulative infections, it is strictly increasing.
The sketch will have a horizontal axis labeled tt and a vertical axis labeled xx. It will show a logistic 'S' shape in the first quadrant, starting at (0,1)(0, 1) and flattening out as it approaches the dotted line x=206x = 206.

3. State the long-term behavior of $x$ for large $t$.

As t→∞t \to \infty, the term e−0.169t\text{e}^{-0.169t} tends to 00.
x=2061+0=206x = \frac{206}{1 + 0} = 206
For large values of tt, xx approaches 206206. (The entire population eventually becomes infected.)

答案

Equation: x=2061+205e−0.169tx = \dfrac{206}{1 + 205\text{e}^{-0.169t}}

For large values of tt, xx approaches 206206.

AI 分析,可能出错

用 AI 批改,圈出写错的那一步,并自动将你的薄弱点存入错题本,形成你的错因库。

同一知识点的其他题

  1. 1曲线上的点 (x,y)(x,y) 满足方程 x=12t2−2tx = \tfrac{1}{2}t^2 - 2t 和 dydx=12t−1\dfrac{dy}{dx} = \dfrac{1}{2t - 1} 其中 tt 为参数,t>2t > 2。 已知当 t=4t = 4…
  2. 2在一个大城镇中,首次发现某种病毒 tt 天后感染该病毒的人数为 xx。假设感染率与 xx 成正比。起初有 5 人感染该病毒,而在首次发现病毒 30 天后,有 5120 人感染该病毒。 在发现病毒 30 天后,人们发现了该病毒的治疗方法…
  3. 3考虑二阶微分方程 d2xdt2+adxdt=b,\frac{d^{2}x}{dt^{2}} + a\frac{dx}{dt} = b, where aa and bb are non-zero constants and dxdt≠ba\frac{dx}{dt} \ne \dfrac{b}{a}…
  4. 4一杯刚泡好的茶,初始温度为 85∘C85^\circ\text{C},被放在一间室温为 25∘C25^\circ\text{C} 的房间里冷却。20分钟后,茶的温度为 55∘C55^\circ\text{C}。已知茶温的变化率与茶温与室温之差成正比。…
  5. 5在一个大城镇中,在首次发现某种病毒 tt 天后,感染该病毒的人数为 xx。假设感染率与 xx 成正比。起初有 5 人感染该病毒,在首次发现该病毒 30 天后,有 5120 人感染。 证明 x=5(2)t3x = 5(2)^{\frac{t}{3}}…