8.累加、累乘与积分
问题:将向量下标为偶数的分量 (x2, x4, …) 累加, 写出相应表达式.
∑i=1⌊n/2⌋x2i\sum_{i=1}^{\lfloor n/2 \rfloor}x_{2i}∑i=1⌊n/2⌋x2i或∑imod 2=0xi\sum_{i\mod2=0}x_i∑imod2=0xi
问题:各出一道累加、累乘、积分表达式的习题, 并给出标准答案.
∑i=110i=55\sum_{i=1}^{10}i=55∑i=110i=55
该累加表达式的python代码及运行结果如下
sum= 0;
for i in range(1,11):
sum=sum+i
print(sum)
55
[Finished in 104ms]
∏i=110i=3628800\prod_{i=1}^{10}i=3628800∏i=110i=3628800
该累乘表达式的python代码及运行结果如下
multiy= 1;
for i in range(1,11):
multiy=multiy*i
print(multiy)
3628800
[Finished in 99ms]
∫12xdx=3/2\int_{1}^{2}x\mathrm{d}x=3/2∫12xdx=3/2
该积分表达式的python代码及运行结果如下
from sympy import *
x=symbols('x')
print(integrate(x,(x,1,2)))
3/2
[Finished in 1.2s]
问题:你使用过三重累加吗? 描述一下其应用.
程序代码中的三重循环?使用excel时的表合并统计。
问题:给一个常用的定积分, 将手算结果与程序结果对比.
∫12xdx=12x2∣12=2−12=32\int_{1}^{2}x\mathrm{d}x=\frac{1}{2}x^{2}|_{1}^{2}=2-\frac{1}{2}=\frac{3}{2}∫12xdx=21x2∣12=2−21=23
from sympy import *
x=symbols('x')
print(integrate(x,(x,1,2)))
3/2
[Finished in 1.2s]
9.线性回归
问题:自己写一个小例子 (n = 3 , m=1) 来验证最小二乘法.
| x | y |
|---|---|
| 1 | 1 |
| 2 | 2 |
| 3 | 2 |
建模y=ax+by=ax+by=ax+b
则X=[xij]=[111213]\mathbf{X}=[x_{ij}]=\left[\begin{matrix}1&1\\1&2\\1&3\\\end{matrix}\right]X=[xij]=⎣⎡111123⎦⎤
Y=[y1,...,yn]T=[1,2,2]T\mathbf{Y}=[y_1,...,y_n]^\mathrm{T}=[1,2,2]^\mathrm{T}Y=[y1,...,yn]T=[1,2,2]T
w∗=arg min∣∣Xw−Y∣∣22\mathbf{w}^*=\argmin||\mathbf{X}\mathbf{w}-\mathbf{Y}||_2^2w∗=argmin∣∣Xw−Y∣∣22
w=(XTX)−1XTY=[36614]−1×[511]=[1223]\mathbf{w}=(\mathbf{X}^\mathrm{T}\mathbf{X})^{-1}\mathbf{X}^\mathrm{T}\mathbf{Y}=\left[\begin{matrix}3&6\\6&14\\\end{matrix}\right]^{-1}\times\left[\begin{matrix}5\\11\\\end{matrix}\right]=\left[\begin{matrix}\frac{1}{2}\\\frac{2}{3}\\\end{matrix}\right]w=(XTX)−1XTY=[36614]−1×[511]=[2132]
所以模型为y=12x+23y=\frac{1}{2}x+\frac{2}{3}y=21x+32.
问题:写出w\mathbf{w}w的推导过程
∥Xw−Y∥22=(Xw−Y)T(Xw−Y)=(wTXT−YT)(Xw−Y)=wTXTXw−wTXTY−YTXw−YTY\|\mathbf{X}\mathbf{w}-\mathbf{Y\|}_2^2=(\mathbf{X}\mathbf{w}-\mathbf{Y})^\mathrm{T}(\mathbf{X}\mathbf{w}-\mathbf{Y})
=(\mathbf{w}^\mathrm{T}\mathbf{X}^\mathrm{T}-\mathbf{Y}^\mathrm{T})(\mathbf{X}\mathbf{w}-\mathbf{Y})=\mathbf{w}^\mathrm{T}\mathbf{X}^\mathrm{T}\mathbf{X}\mathbf{w}-\mathbf{w}^\mathrm{T}\mathbf{X}^\mathrm{T}\mathbf{Y}-\mathbf{Y}^\mathrm{T}\mathbf{X}\mathbf{w}-\mathbf{Y}^\mathrm{T}\mathbf{Y}∥Xw−Y∥22=(Xw−Y)T(Xw−Y)=(wTXT−YT)(Xw−Y)=wTXTXw−wTXTY−YTXw−YTY
令f(w)=wTXTXw−wTXTY−YTXw−YTYf(\mathbf{w})=\mathbf{w}^{\mathrm{T}}\mathbf{X}^{\mathbf{T}}\mathbf{X}\mathbf{w}-\mathbf{w}^{\mathrm{T}}\mathbf{X}^{\mathrm{T}}\mathbf{Y}-\mathbf{Y}^{\mathrm{T}}\mathbf{X}\mathbf{w}-\mathbf{Y}^{\mathrm{T}}\mathbf{Y}f(w)=wTXTXw−wTXTY−YTXw−YTY,将该式关于 w\mathbf{w}w求导 (使用向量求导法则) 并令其为 0。
根据矩阵求导法则可得:
dJ(w)dw=2XTXw−XTY−XTY−0=2XTXw−2XTY=XTXw−XTY\begin{aligned}\frac{\mathrm{d}J(\mathbf{w})}{\mathrm{d}\mathbf{w}}\end{aligned}=2\mathbf{X}^{\mathrm{T}}\mathbf{X}\mathbf{w}-\mathbf{X}^\mathrm{T}\mathbf{Y}-\mathbf{X}^\mathrm{T}\mathbf{Y}-0=2\mathbf{X}^{\mathrm{T}}\mathbf{X}\mathbf{w}-2\mathbf{X}^\mathrm{T}\mathbf{Y}=\mathbf{X}^{\mathrm{T}}\mathbf{X}\mathbf{w}-\mathbf{X}^\mathrm{T}\mathbf{Y}dwdJ(w)=2XTXw−XTY−XTY−0=2XTXw−2XTY=XTXw−XTY
令 dJ(w)dw=0\begin{aligned}\frac{\mathrm{d}J(\mathbf{w})}{\mathrm{d}\mathbf{w}}\end{aligned}=0dwdJ(w)=0得 w=(XTX)−1XTY\mathbf{w}=(\mathbf{X}^\mathrm{T}\mathbf{X})^{-1}\mathbf{X}^\mathrm{T}\mathbf{Y}w=(XTX)−1XTY
10.Logistics回归
问题:自己推导一遍, 并描述这个方法的特点 (不少于 5 条).
1)1)1)用一个表达式表达了两种情况的概率
2)2)2)采用梯度下降求最优解
3)3)3)使用 sigmoid 函数将距离转成概率
4)4)4)表达方式上( ; )’;'右边为参数,左边为数据
5)5)5)点到超平面的距离带符号,需要分正负

被折叠的 条评论
为什么被折叠?



