0%

作为计算数学专业,在对微分方程算法求解时,经常会用到矩阵的运算,特别是稀疏矩阵、对称矩阵和对角矩阵的运算,本篇文章旨在介绍用spdiags函数构建稀疏的三对角、五对角矩阵.

阅读全文 »

  • 在FreeFEM++中,分段函数是由不同的表达式组成的函数,每个表达式定义了函数在定义域上的一个子区间内的行为。这些子区间可以重叠,并且可以使用不同的函数来定义每个子区间。下面将介绍如何在FreeFEM++中编程实现分段函数。

    A. 假设我们要定义一个分段函数f(x),其定义域为R,并且在x0.5x\leq 0.5时,f(x)=x2f(x)=x^2;在x>0.5x>0.5时,f(x)=xf(x)=x。以下是一种实现方法:

    1
    2
    3
    //分段函数的生成
    //f=(条件 ? 满足输出 : 否则输出)
    f=(x<=0.5 ? x^2 : x);

    B. 多分段实例: 并且在x0.5x\leq 0.5时,f(x)=x2f(x)=x^2;在0.5<x10.5<x\leq 1时,f(x)=x1f(x)=x-1; 在x>1x>1时,f(x)=xf(x)=x

    1
    2
    3
    //多分段函数的生成
    //f=(条件1 ? 满足输出1 : (条件2 ? 满足输出2 : 否则输出)),依次类推
    f=(x<=0.5 ? x^2 : (0.5<x<=1 ? x-1 : x);

Cauchy 中值定理的应用

f(x)f(x)[1,+)[1,+\infty)上连续,在(1,+)(1,+\infty)上可导,已知函数ex2f(x)e^{-x^2}f'(x)(1,+)(1,+\infty)上有界,证明函数xex2f(x)xe^{-x^2}f(x)(1,+)(1,+\infty)上有界.

阅读全文 »

Github有什么用?

Github里边有很多的开源代码,特别是对研究生等研究群体比较使用。你可以用来搜索你需要的代码研究,你也可以建立仓库,分享自己的代码。

Github无法使用?

在国内访问Github,有时候根本无法访问,总是出现网页错误。或者是多次访问进去网页了,想下载一个代码,等ing,太慢了,还时不时下载失败。那么,如何解决这个问题呢?

Github访问工具

下面给大家带来一个简单的方法,亲测有效!

阅读全文 »

  • 关于矩阵的几个常用的结论

证明以下结论:

(1) 秩(A+B)(A+B)\leq(A)+(A)+(B)(B);

(2) A,BA,B均为n×nn\times n矩阵,如果AB=OAB=O, 那么

(A)+(A)+(B)n(B)\leq n;

(3) AAn×nn\times n矩阵, 证明:A=An1|A^*|=|A|^{n-1},这里n2n\geq 2;

(4) AAn×nn\times n矩阵, 这里n2n\geq 2, 证明:

 秩 (A)={n, 秩 (A)=n,1, 秩 (A)=n1,0, 秩 (A)<n1.\text { 秩 }\left(A^{*}\right)= \begin{cases}n, & \text { 秩 }(\boldsymbol{A})=n, \\ 1, & \text { 秩 }(\boldsymbol{A})=n-1, \\ 0, & \text { 秩 }(\boldsymbol{A})<n-1 .\end{cases}

(5) A,BA,B分别为n×mn\times mm×nm\times n矩阵, 证明:

(i)EmBAEn=EnAB=EmBA(i)\left| \begin{matrix} E_m& B\\ A & E_n\\ \end{matrix} \right|=|E_n-AB|=|E_m-BA|

(ii) λEnAB=λnmλEmBA|\lambda E_n-AB|=\lambda^{n-m}|\lambda E_m-BA|.

阅读全文 »

  • 多项式在实数域/复数域上的因式分解及应用

(1)将多项式xn1x^n-1在实数域上因式分解;

(2)如果(x2+x+1)f1(x3)+xf2(x3)(x^2+x+1)|f_1(x^3)+xf_2(x^3), 证明(x1)f1(x),(x1)f2(x)(x-1)|f_1(x), (x-1)|f_2(x).

阅读全文 »