2018年7月

4阶泊松-费米方程数值求解

常见于离子液体等带电软物质体系的4阶泊松-费米方程:

\begin{equation} \frac{d^2\phi}{dx^2}-\delta_c^2\frac{d^4\phi}{dx^4}=\frac{\sinh \phi}{1+2\gamma \sinh^2 \phi/2}=\rho(\phi) \label{Poisson-Fermi2} \end{equation}

边界条件:

\begin{equation} \begin{split} \phi(0)=&V_0 \\ \phi'''(0)=&0 \\ \phi(\infty)=&0\\ \phi'(\infty)=&0 \end{split} \label{BC2} \end{equation}

下面我们用bvp4c 解方程,重复出文献 Double Layer in Ionic Liquids: Overscreening versus Crowding中 FIG. 2(a) 中的虚线。

依次取 $V_0=1, 10, 100$,解方程,程序如下:

- 阅读剩余部分 -

泊松-费米方程数值求解

常见于离子液体等带电软物质体系的泊松-费米方程

\begin{equation} \frac{d^2\phi}{dx^2}=\frac{\sinh \phi}{1+2\gamma \sinh^2 \phi/2}=\rho(\phi) \label{Poisson-Fermi2} \end{equation}

边界条件:

\begin{equation} \begin{split} \phi(0)=&V_0 \\ \phi(\infty)=&0 \end{split} \label{BC2} \end{equation}

下面我们用bvp4c 解方程,重复出文献 Double Layer in Ionic Liquids: Overscreening versus Crowding中 FIG. 2(a) 中的虚线。

依次取 $V_0=1, 10, 100$,解方程,程序如下:

- 阅读剩余部分 -

matlab 合并绘图,如何加图例?

用语句 legend('-DynamicLegend') 和 hold all ,请看下面的例子:

1
2
3
4
5
x=0:.01:10;
plot(x, sin(x), 'DisplayName','sin');
legend('-DynamicLegend');
hold all; % add new plot lines on top of previous ones
plot(x, cos(x), 'DisplayName','cos');

世界上那些最美的方程



数学方程不仅实用——很多还非常优美。许多科学家承认,他们经常喜欢特定的公式,不仅仅是因为它们的功能强大,还因为它们的形式优雅、简洁及其中所蕴涵的诗一般的真理。

当某些特别著名的方程,比如爱因斯坦的质能方程 $E=mc^2$, 在公众面前享誉极盛时,许多公众不那么熟悉的方程在科学家那儿却拥者甚众。 LiveScience 咨询了许多物理学家、天文学家和数学家,将他们喜爱的数学公式罗列如后:

- 阅读剩余部分 -