工程数学基础 考试 要点 复习:矩阵论,数值计算
v向量的1-范数∥x∥1=∑k=1n∣ξk∣\|x\|_1=\sum_{k=1}^n|\xi_k|∥x∥1=∑k=1n∣ξk∣向量的2-范数∥x∥2=(∑k=1n∣ξk∣2)12\|x\|_2=(\sum_{k=1}^n|\xi_k|^2)^\frac12∥x∥2=(∑k=1n∣ξk∣2)21向量的\infty-范数∥x∥∞=max1≤k≤n∣ξk∣\|x\|_\infty=\ma
v
向量的1-范数
∥x∥1=∑k=1n∣ξk∣\|x\|_1=\sum_{k=1}^n|\xi_k|∥x∥1=∑k=1n∣ξk∣
向量的2-范数
∥x∥2=(∑k=1n∣ξk∣2)12\|x\|_2=(\sum_{k=1}^n|\xi_k|^2)^\frac12∥x∥2=(∑k=1n∣ξk∣2)21
向量的∞\infty∞-范数
∥x∥∞=max1≤k≤n∣ξk∣\|x\|_\infty=\max_{1\leq k \leq n}|\xi_k|∥x∥∞=max1≤k≤n∣ξk∣
列范数=每列模的和取最大
∥A∥1=max1≤j≤n∑i=1n∣aij∣\|A\|_1=\max_{1\leq j \leq n}\sum_{i=1}^n|a_{ij}|∥A∥1=max1≤j≤n∑i=1n∣aij∣
行范数=每行模的和取最大
∥A∥∞=max1≤i≤n∑j=1n∣aij∣\|A\|_\infty=\max_{1\leq i \leq n}\sum_{j=1}^n|a_{ij}|∥A∥∞=max1≤i≤n∑j=1n∣aij∣
2-范数导出的方阵A的算子范数
∥A∥2=max∥x∥2=1∥Ax∥2=ρ(AHA)\|A\|_2=\max_{\|x\|_2=1}\|Ax\|_2=\sqrt{\rho(A^HA)}∥A∥2=max∥x∥2=1∥Ax∥2=ρ(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}∥A∥F=(∑i=1n∑j=1n∣aij∣2)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(λE−A)=∣λE−A∣
正规矩阵 酉对角化
正规矩阵
AHA=AAHA^HA=AA^HAHA=AAH
酉矩阵
AHA=AAH=EA^HA=AA^H=EAHA=AAH=E
正交矩阵,当酉矩阵中
A∈Rn×nA \in \mathbb{R}^{n \times n}A∈Rn×n
Hermite矩阵
AH=AA^H=AAH=A
实对称矩阵,当Hermite矩阵中
A∈Rn×nA \in \mathbb{R}^{n \times n}A∈Rn×n
条件数
cond1A=∥A−1∥1∥A∥1cond_1A=\|A^{-1}\|_1\|A\|_1cond1A=∥A−1∥1∥A∥1
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=∥A−1∥2∥A∥2=min{σ(ATA)}max{σ(ATA)}
cond∞A=∥A−1∥∞∥A∥∞cond_\infty A=\|A^{-1}\|_\infty\|A\|_\inftycond∞A=∥A−1∥∞∥A∥∞
严格对角占优
∣aii∣>∑j=1,j≠in∣aij∣|a_{ii}|>\sum_{j=1,j\neq i}^n|a_{ij}|∣aii∣>∑j=1,j=in∣aij∣
Jacobi迭代,D为A对角线上元素构成的矩阵
M1=D−1(D−A)M_1=D^{-1}(D-A)M1=D−1(D−A)
Gauss-Seidel迭代,D为A对角线上元素构成的矩阵
M2=(D−L)−1UM_2=(D-L)^{-1}UM2=(D−L)−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=knxk−xix−xi)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](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+f[x0,x1,x2,x3](x−x0)(x−x1)(x−x2)
迭代法收敛
差商
Hermite插值
一次多项式
S1∗=a0∗+a1∗xS_1^*=a_0^*+a_1^*xS1∗=a0∗+a1∗x
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λE−A
Jordan标准型
初等因子组 y-1
有理标准型 自然法式
不变因子 d1(y) 求相伴矩阵
初值问题
求A的最小多项式 m(y)=dn(y)
∣λE−A∣|\lambda E-A|∣λE−A∣
验证相乘是否为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
deteA=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,23x2−1,25x3−3x
δ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=2b−a[∥g∥22−∑k=0n22k+1∣<g,pk>∣2]
∥g∥22=∫−11∣g(t)∣2dt\|g\|_2^2=\int_{-1}^1|g(t)|^2dt∥g∥22=∫−11∣g(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>|^2∑k=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)dx∫abu(x)v′(x)dx=u(x)v(x)∣ab−∫abu′(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=2b−a[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+21b−a[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+22b−a[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+23b−a[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=4−14T21−T20,S21=4−14T22−T21,S22=4−14T23−T22
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=42−142S21−S20,C21=42−142S22−S21
R20=43C21−C2043−1R_{2^0}=\frac{4^3C_{2^1}-C_{2^0}}{4^3-1}R20=43−143C21−C20
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,...,N−1)
证明题
参考文献
工程数学基础教程

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