v
向量的1-范数

∥x∥1=∑k=1n∣ξk∣\|x\|_1=\sum_{k=1}^n|\xi_k|x1=k=1nξk

向量的2-范数

∥x∥2=(∑k=1n∣ξk∣2)12\|x\|_2=(\sum_{k=1}^n|\xi_k|^2)^\frac12x2=(k=1nξk2)21

向量的∞\infty-范数

∥x∥∞=max⁡1≤k≤n∣ξk∣\|x\|_\infty=\max_{1\leq k \leq n}|\xi_k|x=max1knξk

列范数=每列模的和取最大

∥A∥1=max⁡1≤j≤n∑i=1n∣aij∣\|A\|_1=\max_{1\leq j \leq n}\sum_{i=1}^n|a_{ij}|A1=max1jni=1naij

行范数=每行模的和取最大

∥A∥∞=max⁡1≤i≤n∑j=1n∣aij∣\|A\|_\infty=\max_{1\leq i \leq n}\sum_{j=1}^n|a_{ij}|A=max1inj=1naij

2-范数导出的方阵A的算子范数

∥A∥2=max⁡∥x∥2=1∥Ax∥2=ρ(AHA)\|A\|_2=\max_{\|x\|_2=1}\|Ax\|_2=\sqrt{\rho(A^HA)}A2=maxx2=1Ax2=ρ(AHA)

frobenius Form

∥A∥F=(∑i=1n∑j=1n∣aij∣2)12\|A\|_F=(\sum_{i=1}^n\sum_{j=1}^n|a_{ij}|^2)^\frac{1}{2}AF=(i=1nj=1naij2)21

谱半径=A中最大的特征值的模

ρ(A)=max{∣λ1∣,∣λ2∣,...,∣λn∣}det(λE−A)=∣λE−A∣\rho(A)=max\{|\lambda_1|,|\lambda_2|,...,|\lambda_n|\} det(\lambda E-A)=|\lambda E-A|ρ(A)=max{λ1,λ2,...,λn}det(λEA)=λEA

正规矩阵 酉对角化

正规矩阵

AHA=AAHA^HA=AA^HAHA=AAH

酉矩阵

AHA=AAH=EA^HA=AA^H=EAHA=AAH=E

正交矩阵,当酉矩阵中

A∈Rn×nA \in \mathbb{R}^{n \times n}ARn×n

Hermite矩阵

AH=AA^H=AAH=A

实对称矩阵,当Hermite矩阵中

A∈Rn×nA \in \mathbb{R}^{n \times n}ARn×n

条件数

cond1A=∥A−1∥1∥A∥1cond_1A=\|A^{-1}\|_1\|A\|_1cond1A=A11A1
cond2A=∥A−1∥2∥A∥2=max⁡{σ(ATA)}min⁡{σ(ATA)}cond_2A=\|A^{-1}\|_2\|A\|_2=\sqrt{\frac{\max\{\sigma(A^TA)\}}{\min\{\sigma(A^TA)\}}}cond2A=A12A2=min{σ(ATA)}max{σ(ATA)}
cond∞A=∥A−1∥∞∥A∥∞cond_\infty A=\|A^{-1}\|_\infty\|A\|_\inftycondA=A1A

严格对角占优

∣aii∣>∑j=1,j≠in∣aij∣|a_{ii}|>\sum_{j=1,j\neq i}^n|a_{ij}|aii>j=1,j=inaij

Jacobi迭代,D为A对角线上元素构成的矩阵

M1=D−1(D−A)M_1=D^{-1}(D-A)M1=D1(DA)

Gauss-Seidel迭代,D为A对角线上元素构成的矩阵

M2=(D−L)−1UM_2=(D-L)^{-1}UM2=(DL)1U

填空题

Lagrange插值多项式

Ln(x)=∑k=0nyklk(x)=∑k=0n(∏i=0,i≠knx−xixk−xi)ykL_n(x)=\sum_{k=0}^ny_kl_k(x)=\sum_{k=0}^n(\prod_{i=0,i\neq k}^n\frac{x-x_i}{x_k-x_i})y_kLn(x)=k=0nyklk(x)=k=0n(i=0,i=knxkxixxi)yk

Newton插值多项式

Nn(x)=f[x0]+∑k=1nf[x0,x1,...,xk]ωk(x)N_n(x)=f[x_0]+\sum_{k=1}^nf[x_0,x_1,...,x_k]\omega_k(x)Nn(x)=f[x0]+k=1nf[x0,x1,...,xk]ωk(x)

N3(x)=f(x0)+f[x0,x1](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+f[x0,x1,x2,x3](x−x0)(x−x1)(x−x2)N_3(x)=f(x_0) +f[x_0,x_1](x-x_0) +f[x_0,x_1,x_2](x-x_0)(x-x_1) +f[x_0,x_1,x_2,x_3](x-x_0)(x-x_1)(x-x_2)N3(x)=f(x0)+f[x0,x1](xx0)+f[x0,x1,x2](xx0)(xx1)+f[x0,x1,x2,x3](xx0)(xx1)(xx2)

迭代法收敛

差商

Hermite插值

一次多项式

S1∗=a0∗+a1∗xS_1^*=a_0^*+a_1^*xS1=a0+a1x

M=span{1,x}M=span\{1,x\}M=span{1,x}

φ0(x)=1,φ1(x)=x\varphi_0(x)=1,\varphi_1(x)=xφ0(x)=1,φ1(x)=x

<φ0,φ0>=∑1×1<\varphi_0,\varphi_0>=\sum1 \times 1<φ0,φ0>=1×1

<φ0,φ1>=<φ1,φ0>=∑1×xk<\varphi_0,\varphi_1>=<\varphi_1,\varphi_0>=\sum1 \times x_k<φ0,φ1>=<φ1,φ0>=1×xk

<φ1,φ1>=∑xk2<\varphi_1,\varphi_1>=\sum x_k^2<φ1,φ1>=xk2

<f,φ0>=∑f(xk)<f,\varphi_0>=\sum f(x_k)<f,φ0>=f(xk)

<f,φ1>=∑xkf(xk)<f,\varphi_1>=\sum x_k f(x_k)<f,φ1>=xkf(xk)

[<φ0,φ0><φ0,φ1><φ1,φ0><φ1,φ1>][a0a1]=[<f,φ0><f,φ1>]\begin{bmatrix} <\varphi_0,\varphi_0> &<\varphi_0,\varphi_1>\\ <\varphi_1,\varphi_0> &<\varphi_1,\varphi_1>\\ \end{bmatrix} \begin{bmatrix} a_0\\ a_1\\ \end{bmatrix} =\begin{bmatrix} <f,\varphi_0>\\ <f,\varphi_1>\\ \end{bmatrix}[<φ0,φ0><φ1,φ0><φ0,φ1><φ1,φ1>][a0a1]=[<f,φ0><f,φ1>]

Euler 稳定性

0<h≤−2λ0<h\leq -\frac{2}{\lambda}0<hλ2

Runge-kutta 稳定性

0<h≤−2.78λ0<h\leq -\frac{2.78}{\lambda}0<hλ2.78

三次样条

初等变换

Smith标准型

λE−A\lambda E-AλEA

Jordan标准型

初等因子组 y-1

有理标准型 自然法式

不变因子 d1(y) 求相伴矩阵

初值问题

求A的最小多项式 m(y)=dn(y)

∣λE−A∣|\lambda E-A|λEA

验证相乘是否为O

eAte^{At}eAt

根据特征值的个数与重根,设A的次数

eAt=a0(t)E+a1(t)A+a2(t)A2e^{At}=a_0(t)E+a_1(t)A+a_2(t)A^2eAt=a0(t)E+a1(t)A+a2(t)A2

eλte^{\lambda t}eλt

T(λt)=a0(t)+a1(t)λ+a2(t)λ2T(\lambda t)=a_0(t)+a_1(t)\lambda+a_2(t)\lambda^2T(λt)=a0(t)+a1(t)λ+a2(t)λ2

计算x(t),其中C为初值

x(t)=eAtCx(t)=e^{At}Cx(t)=eAtC

det⁡eA=etrA\det e^A=e^{trA}deteA=etrA

曲线拟合

最佳平方逼近 Legendre多项式

S3∗(x)=12<f,p0>p0(x)+32<f,p1>p1(x)+52<f,p2>p2(x)+72<f,p3>p3(x)S_3^*(x)=\frac12<f,p_0>p_0(x)+\frac32<f,p_1>p_1(x)+\frac52<f,p_2>p_2(x)+\frac72<f,p_3>p_3(x)S3(x)=21<f,p0>p0(x)+23<f,p1>p1(x)+25<f,p2>p2(x)+27<f,p3>p3(x)

<g,pk>=∫−11g(x)⋅pkdx<g,p_k>=\int_{-1}^{1}g(x)\cdot p_kdx<g,pk>=11g(x)pkdx

pk=1,x,3x2−12,5x3−3x2p_k=1,x,\frac{3x^2-1}{2},\frac{5x^3-3x}{2}pk=1,x,23x21,25x33x

δ2=b−a2[∥g∥22−∑k=0n2k+12∣<g,pk>∣2]\delta^2=\frac{b-a}2[\|g\|_2^2-\sum_{k=0}^n\frac{2k+1}2|<g,p_k>|^2]δ2=2ba[g22k=0n22k+1<g,pk>2]

∥g∥22=∫−11∣g(t)∣2dt\|g\|_2^2=\int_{-1}^1|g(t)|^2dtg22=11g(t)2dt

∑k=0n2k+12∣<g,pk>∣2=12∣<g,p0>∣2+32∣<g,p1>∣2+52∣<g,p2>∣2+72∣<g,p3>∣2\sum_{k=0}^n\frac{2k+1}2|<g,p_k>|^2=\frac12|<g,p_0>|^2+\frac32|<g,p_1>|^2+\frac52|<g,p_2>|^2+\frac72|<g,p_3>|^2k=0n22k+1<g,pk>2=21<g,p0>2+23<g,p1>2+25<g,p2>2+27<g,p3>2

分部积分法

∫abu(x)v′(x)dx=u(x)v(x)∣ab−∫abu′(x)v(x)dx\int_a^b u(x)v'(x)dx=u(x)v(x)|_a^b-\int_a^bu'(x)v(x)dxabu(x)v(x)dx=u(x)v(x)ababu(x)v(x)dx

(u(x)v(x))′=u(x)v′(x)+u′(x)v(x)(u(x)v(x))'=u(x)v'(x)+u'(x)v(x)(u(x)v(x))=u(x)v(x)+u(x)v(x)

Romberg算法 初始计算 迭代计算

I=∫abf(x)dxI=\int_a^bf(x)dxI=abf(x)dx

T20=b−a2[f(a)+f(b)]T_{2^0}=\frac{b-a}2[f(a)+f(b)]T20=2ba[f(a)+f(b)]

T21=T202+b−a21[f(a+b2)]T_{2^1}=\frac{T_{2^0}}2+\frac{b-a}{2^1}[f(\frac{a+b}2)]T21=2T20+21ba[f(2a+b)]

T22=T212+b−a22[f(14(a+b))+f(34(a+b))]T_{2^2}=\frac{T_{2^1}}2+\frac{b-a}{2^2}[f(\frac14(a+b))+f(\frac34(a+b))]T22=2T21+22ba[f(41(a+b))+f(43(a+b))]

T23=T222+b−a23[f(18(a+b))+f(38(a+b))+f(58(a+b))+f(78(a+b))]T_{2^3}=\frac{T_{2^2}}2+\frac{b-a}{2^3}[f(\frac18(a+b))+f(\frac38(a+b))+f(\frac58(a+b))+f(\frac78(a+b))]T23=2T22+23ba[f(81(a+b))+f(83(a+b))+f(85(a+b))+f(87(a+b))]

S20=4T21−T204−1,S21=4T22−T214−1,S22=4T23−T224−1S_{2^0}=\frac{4T_{2^1}-T_{2^0}}{4-1}, S_{2^1}=\frac{4T_{2^2}-T_{2^1}}{4-1}, S_{2^2}=\frac{4T_{2^3}-T_{2^2}}{4-1}S20=414T21T20,S21=414T22T21,S22=414T23T22

C20=42S21−S2042−1,C21=42S22−S2142−1C_{2^0}=\frac{4^2S_{2^1}-S_{2^0}}{4^2-1}, C_{2^1}=\frac{4^2S_{2^2}-S_{2^1}}{4^2-1}C20=42142S21S20,C21=42142S22S21

R20=43C21−C2043−1R_{2^0}=\frac{4^3C_{2^1}-C_{2^0}}{4^3-1}R20=43143C21C20

Runge-Kutta法

{y′=zz′=(1+x2)y+1y(0)=1,z(0)=3\begin{cases} y'=z \\ z'=(1+x^2)y+1 \\ y(0)=1,z(0)=3 \end{cases}y=zz=(1+x2)y+1y(0)=1,z(0)=3

{yn+1=yn+h6(k1+2k2+2k3+k4)zn+1=zn+h6(l1+2l2+2l3+l4)k1=zn,l1=(1+xn2)yn+1k2=zn+h2l1,l2=(1+(xn+h2)2)(yn+h2k1)+1k3=zn+h2l2,l3=(1+(xn+h2)2)(yn+h2k2)+1k4=zn+hl3,l4=(1+(xn+h)2)(yn+hk3)+1y0=1,z0=3(n=0,1,2,...,N−1)\begin{cases} y_{n+1}=y_n+\frac{h}6(k_1+2k_2+2k_3+k_4)\\ z_{n+1}=z_n+\frac{h}6(l_1+2l_2+2l_3+l_4)\\ k_1=z_n,l_1=(1+x_n^2)y_n+1 \\ k_2=z_n+\frac{h}2l_1,l_2=(1+(x_n+\frac{h}2)^2)(y_n+\frac{h}2k_1)+1 \\ k_3=z_n+\frac{h}2l_2,l_3=(1+(x_n+\frac{h}2)^2)(y_n+\frac{h}2k_2)+1 \\ k_4=z_n+hl_3,l_4=(1+(x_n+h)^2)(y_n+hk_3)+1 \\ y_0=1,z_0=3 \end{cases} (n=0,1,2,...,N-1)yn+1=yn+6h(k1+2k2+2k3+k4)zn+1=zn+6h(l1+2l2+2l3+l4)k1=zn,l1=(1+xn2)yn+1k2=zn+2hl1,l2=(1+(xn+2h)2)(yn+2hk1)+1k3=zn+2hl2,l3=(1+(xn+2h)2)(yn+2hk2)+1k4=zn+hl3,l4=(1+(xn+h)2)(yn+hk3)+1y0=1,z0=3(n=0,1,2,...,N1)

证明题

参考文献

工程数学基础教程

Logo

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。

更多推荐