水雷——屯

屯【卦三】水雷屯
屯,其卦象为震下坎上
(⊙﹏⊙∥，我不会输入那个卦象，有机会去学习怎么输入，然后给大家展示)
\qquad 屯，元亨，利贞，勿用有攸往，利建侯。
\qquad 《彖》曰：屯，刚柔始交而难生。动乎险中，大亨贞。雷雨之动满盈，天造草昧，宜建侯而不宁。
\qquad 《象》曰：“即鹿无虞”，以从禽也。君子舍之，往吝穷也。
\qquad 屯卦，让我想起来，金庸笔下的一个非常浪漫、非常爱自由的人物——令狐冲。我非常地喜欢他，爱他的浪漫不羁、热爱自由、真性情、不为世俗所困。他在《笑傲江湖》全篇之中，始终是保持着一颗赤子之心，绝对不向所谓的“强者”，如余沧海、任我行之流，屈膝谄媚。但是在小说开始时，他有傲骨、有智慧，但没有能力，勉强成功地英雄救美，但是自己也身负重伤，命悬一线。而后来，他虽然是受到委屈，去思过崖面壁思过，看上去是一段非常煎熬、非常难受的时期，但也是在这段时间中，他发现了华山石洞地秘密，遇到一样潇洒不羁地风清扬，并得其真传——独孤九剑。尽管令狐冲他后来的种种遭遇，确实也一路坎坷，但是相比于他开始的三流水准（武功和见识阅历），他已经能够游刃有余的处理解决各类难题了。（当然，他的情感问题是一件非常棘手的，并且具有扑朔迷离的味道，我只能感慨，金庸真的是大智慧。以后有机会和大家聊一聊他的情感生涯）
\qquad 我想通过令狐冲的遭际说什么呢？当下，可能确实有诸多不顺心之事，抑或是看不惯之人，我们可以选择去“隐居”，不过现在已经2019了，比不了小说中的时代，做一个完全脱节社会的人，不现实也不可能。我想体面地活着，同时也不想处处受他人掣肘，我向往与世隔绝地生活，所以，现在我得苦练七十二技了。学习本身确实蛮痛苦的，不过换个角度来说，比如我喜欢看看《易经》、《庄子》，如果我不学习，我不就看不懂了吗？所以我们把学习当作手段，把兴趣爱好当作学习服务的对象。啊！这样一想，不禁觉得胸中的积郁之气，便大图为快啊！
\qquad 所以，我亲爱的朋友，如果你这会儿依旧在拼命地往上走，拼命地学习，不要“畏途巉岩不可攀”，假使真的遇到了很棘手的问题，不妨说出来，或是找个小说来寄托一下情感，宣泄一下情绪。让我们再想想令狐冲在思过崖上的那段日子吧，“若无漫漫长夜，何来熠熠闪光”，可能我们也在思过崖上，只是风清扬觉得还不是来见我们的时候，所以千万不要放弃哦！d=====(￣▽￣*)b

附一张海中集训班黑板报：

1.求解抽样函数广义积分的数学方法
$${\int }_0^{\infty }\frac{sin(ax)}{x}\mathrm{d}x=\{\begin{array}{ll}\frac{\pi }{2}& a>0\\ -\frac{\pi }{2}& a<0\end{array}$$
证明：
$${\int }_0^{\infty }\frac{sin(ax)}{x}\mathrm{d}x={\int }_0^{\infty }\frac{sin(ax)}{ax}\mathrm{d}(ax).$$
令$t=ax$,
（1）当$a>0$时，$$\begin{gathered} 原式={\int }_0^{\infty }\frac{sint}{t}\mathrm{d}t={\int }_0^{\infty }(e^{-at}\frac{sint}{t}){|}_{\infty }^0\mathrm{d}t={\int }_0^{\infty }({\int }_0^{\infty }{e}^{-at}sint\mathrm{d}a)\mathrm{d}t \\ ={\int }_0^{\infty }({\int }_0^{\infty }{e}^{-at}sint\mathrm{d}t)\mathrm{d}a(1.1) \end{gathered}$$
此处我们单独计算${\int }_0^{\infty }{e}^{-at}sint\mathrm{d}t$
$$\begin{gathered} {\int }_0^{\infty }{e}^{-at}sint\mathrm{d}t={\int }_0^{\infty }\frac{sint}{-a}\mathrm{d}{e}^{-at}=\frac{e^{-at}sint}{-a}{|}_0^{\infty }-{\int }_0^{\infty }\frac{e^{-at}}{-a}\mathrm{d}sint \\ =0-0+{\int }_0^{\infty }\frac{e^{-at}cost}{a}\mathrm{d}t={\int }_0^{\infty }\frac{cost}{-a^2}\mathrm{d}{e}^{-at}=\frac{e^{-at}cost}{-a^2}{|}_0^{\infty }-{\int }_0^{\infty }\frac{e^{-at}}{-a^2}\mathrm{d}cost \\ =0-\frac{1}{-a^2}+{\int }_0^{\infty }\frac{e^{-at}sint}{a^2}\mathrm{d}t=\frac{1}{a^2}+\frac{1}{a^2}{\int }_0^{\infty }{e}^{-at}sint\mathrm{d}t \end{gathered}$$
即
$${\int }_0^{\infty }{e}^{-at}sint\mathrm{d}t=\frac{\frac{1}{a^2}}{1+\frac{1}{a^2}}=\frac{1}{1+a^2}a>0$$
则式$（1.1）$的值为$${\int }_0^{\infty }\frac{1}{1+a^2}\mathrm{d}a=arctana{|}_0^{\infty }=\frac{\pi }{2}$$
故当$a>0$时
$${\int }_0^{\infty }\frac{sin(ax)}{x}\mathrm{d}x=\frac{\pi }{2}$$
(2)当$a<0$时，与（1）同样，可以得到
$${\int }_0^{\infty }\frac{sin(ax)}{x}\mathrm{d}x=-\frac{\pi }{2}$$
具体过程读者可以仿照（1），自己尝试证明。

综合（1）和（2），则
$${\int }_0^{\infty }\frac{sin(ax)}{x}\mathrm{d}x=\{\begin{array}{ll}\frac{\pi }{2}& a>0\\ -\frac{\pi }{2}& a<0\end{array}$$
以下给出上面证明部分的LaTeX源代码：

$$
\int_0^\infty \frac{sin(ax)}{x}\mathrm{d}x\,= \begin{cases}  
 \frac{\pi}{2}\,& a > 0 \\
-\frac{\pi}{2}\,& a < 0
\end{cases}
$$
证明：
$$\int_0^\infty \frac{sin(ax)}{x}\mathrm{d}x\,=\int_0^\infty \frac{sin(ax)}{ax}\mathrm{d}(ax)\,\,. $$
令$t=ax$,
（1）当$a>0$时，$$
原式=\int_0^\infty\frac{sint}{t}\mathrm{d}t\,\,=\int_0^\infty (e^{-at}\frac{sint}{t})\big|^{0}_{\infty}\mathrm{d}t \,=\int_0^\infty(\int_0^\infty e^{-at}sint\mathrm{d}a)\mathrm{d}t\\=\int_0^\infty(\int_0^\infty e^{-at}sint\mathrm{d}t)\mathrm{d}a\quad (1.1)$$
此处我们单独计算$\int_0^\infty e^{-at}sint\mathrm{d}t$
$$
\int_0^\infty e^{-at}sint\mathrm{d}t\,=\int_0^\infty\frac{sint}{-a}\mathrm{d}e^{-at}\,=\frac{e^{-at}sint}{-a}\bigg|_0^\infty-\int_0^\infty\frac{e^{-at}}{-a}\mathrm{d}sint\,\\
\qquad =0-0+\int_0^\infty\frac{e^{-at}cost}{a}\mathrm{d}t\,=\int_0^\infty\frac{cost}{-a^2}\mathrm{d}e^{-at}\,= \frac{e^{-at}cost}{-a^2}\bigg|_0^\infty-\int_0^\infty\frac{e^{-at}}{-a^2}\mathrm{d}cost\,\\
=0-\frac{1}{-a^2}+\int_0^\infty\frac{e^{-at}sint}{a^2}\mathrm{d}t\,=\frac{1}{a^2}+\frac{1}{a^2}\int_0^\infty e^{-at}sint\mathrm{d}t\,
$$
即
$$
\int_0^\infty e^{-at}sint\mathrm{d}t\,=\frac{\frac{1}{a^2}}{1+\frac{1}{a^2}}=\frac{1}{1+a^2}\,\qquad a>0
$$
则式$（1.1）$的值为$$\int_0^\infty\frac{1}{1+a^2}\mathrm{d}a=arctana\big|_0^\infty=\frac{\pi}{2}$$
故当$a>0$时
$$\int_0^\infty \frac{sin(ax)}{x}\mathrm{d}x\,= \frac{\pi}{2}$$

2.继昨日求解球面方程交点的方法
$$\{\begin{array}{l}(x-a_1{)}^2+(y-b_1{)}^2+(z-c_1{)}^2=d_1^2(2.1)\\ (x-a_2{)}^2+(y-b_2{)}^2+(z-c_2{)}^2=d_2^2(2.2)\\ (x-a_3{)}^2+(y-b_3{)}^2+(z-c_3{)}^2=d_3^2(2.3)\end{array}$$首先对该方程组进行加减消元(2.1)-(2.2),(2.2)-(2.3)
$$\{\begin{array}{l}(2a_2-2a_1)x+(2b_2-2b_1)y+(2c_2-2c_1)z=d_1^2-d_2^2+a_2^2+b_2^2+c_2^2-a_1^2-b_1^2-c_1^2(2.4)\\ (2a_3-2a_2)x+(2b_3-2b_2)y+(2c_3-2c_2)z=d_2^2-d_3^2+a_3^2+b_3^2+c_3^2-a_2^2-b_2^2-c_2^2(2.5)\end{array}$$
可以解出
$$\{\begin{array}{l}y=\frac{((2c_2-2c_1)(2a_3-2a_2)-(2c_3-2c_2)(2a_2-2a_1))x+(d_1^2-d_2^2+a_2^2+b_2^2+c_2^2-a_1^2-b_1^2-c_1^2)(2c_3-2c_2)-(d_2^2-d_3^2+a_3^2+b_3^2+c_3^2-a_2^2-b_2^2-c_2^2)(2c_2-2c_1)}{(2b_3-2b_2)(2c_2-2c_1)-(2b_2-2b_1)(2c_3-2c_2)}.\\ \\ z=\frac{((2b_2-2b_1)(2a_3-2a_2)-(2b_3-2b_2)(2a_2-2a_1))x+(d_1^2-d_2^2+a_2^2+b_2^2+c_2^2-a_1^2-b_1^2-c_1^2)(2a_3-2a_2)-(d_2^2-d_3^2+a_3^2+b_3^2+c_3^2-a_2^2-b_2^2-c_2^2)(2a_2-2a_1)}{(2b_3-2b_2)(2c_2-2c_1)-(2b_2-2b_1)(2c_3-2c_2)}.\end{array}$$
将$y$和$z$回代到$（2.1）$式中，解下面这个关于$x$的二元一次方程
$$(x-a_1{)}^2+(y-b_1{)}^2+(z-c_1{)}^2=d_1^2$$
进而求得$x,y,z$的值。

另附Matlab源代码如下：

syms x y z a1 a2 a3 b1 b2 b3 c1 c2 c3 d1 d2 d3 m n;
eqns=[(x-a1)^2+(y-b1)^2+(z-c1)^2==d1^2,(x-a2)^2+(y-b2)^2+(z-c2)^2==d2^2,(x-a3)^2+(y-b3)^2+(z-c3)^2==d3^2];
eqn=[(y-b1)^2+(z-c1)^2-(y-b2)^2-(z-c2)^2==d1^2-d2^2+(x-a2)^2-(x-a1)^2,(y-b2)^2+(z-c2)^2-(y-b3)^2-(z-c3)^2==d2^2-d3^2+(x-a3)^2-(x-a2)^2];
var=[y,z];
[E,g]=equationsToMatrix(eqn,var);
Solution=inv(E)*g; 
soly=Solution(1,:);
solz=Solution(2,:);
f(x)=(x-a1)^2+(soly-b1)^2+(solz-c1)^2-d1^2;
t=solve(f(x),x);
Y(x)=soly;
Z(x)=solz;
r=Y(t);
s=Z(t);
Output=[t;s;r];

关于此代码说明，由于是含有9个参数的二元一次方程组，因此该解的输出受到Matlab自身输出长度的限制，存在截断问题，另外Matlab关于非齐次方程似乎没有通式或是函数直接调用，并且根据线性代数的知识，我们要知道m个n次的方程组，在消去最高次项后，真正有效的n-1次方程数目只有m-1个，所以该例子是无法直接用线性方程的知识来求解的！

非常感谢陆同学的帮助，在她的引导下，有一个更一般性的算法。
其源代码如下：

syms x y z x1 x2 x3 y1 y2 y3 z1 z2 z3 a11 a12 a13 a21 a22 a23 a31 a32 a33 d1 d2 d3;
B=[x1,y1,z1;x2,y2,z2;x3,y3,z3];
A=[a11,a12,a13;a21,a22,a23;a31,a32,a33];
C=[x,y,z;x,y,z;x,y,z];
M1=C-B;
M2=M1.*M1;
M3=A*transpose(M2);
b=[d1 ;d2 ;d3];
eqn0=M3-b;
eqn=transpose(eqn0(:,1));
vars=[x,y,z];
[solx ,soly ,solz]=solve(eqn,vars);
sol=[solx(1,:),solx(2,:);soly(1,:),soly(2,:);solz(1,:),solz(2,:)];

3.傅里叶$（Fourier）$变换补充:
昨天我写了一个很简单粗暴的式子：
$$Sa(\frac{t}{2})\leftrightarrow 2\pi g(-\omega ).$$
这个就是$（Fourier）$
$$f(t)\leftrightarrow F(j\omega ).$$
$（Fourier）$变换的定义式如下：
$$F(j\omega )={\int }_{-\infty }^{+\infty }f(t)e^{-j\omega t}\mathrm{d}t.$$
即单纯的做一个定义，类似于$\Gamma (x)$ 、$\mathrm{B}(x)$一样的。
其具体意义到底是什么呢？今天已经写了这么久呢，就明天再说吧！

（(〃￣︶￣)人(￣︶￣〃)其实是猴哥?忘记了傅里叶变换是怎么来，容我去复习一天，明天再来和大家分享?）