![](https://img-blog.csdnimg.cn/20201014180756757.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
文章平均质量分 63
矩阵、数值计算、优化理论
blogZT
这个作者很懒,什么都没留下…
展开
-
非线性系列(三)—— 非线性求解器算法分类
CAE中的非线性方程组求解主要依赖牛顿法(及牛顿法的变体),步骤如下。原创 2024-07-10 20:03:31 · 388 阅读 · 0 评论 -
非线性系列(二)—— 共轭梯度法 Conjugate Gradient Method (线性及非线性)
CG方法第一次被用来解非线性优化问题,是由Fletcher和Reeves提出的。(英语:BiConjugate gradient method)提供了一种处理非对称矩阵情况的推广。共轭梯度法(英语:Conjugate gradient method),是求解系数矩阵为。共轭梯度法中,搜索方向p,是关于A共轭的,即。最后的解为 [1, 1],f(x)最小值为0。其中delta为方向,beta为步长。共轭梯度法是一个迭代方法,它适用于。的线性方程组的数值解的方法。1. 求解线性方程组,首先安装自动微分工具。原创 2024-04-14 08:10:46 · 921 阅读 · 0 评论 -
非线性系列(一)—— 非线性方程组求解与优化的关系 Nonlinear & Optimization
解非线性方程(Nonlinear Equations)的技术在其动机、分析和算法实现方面与优化技术(Optimization)有重叠。原创 2024-04-08 21:10:45 · 555 阅读 · 0 评论 -
奇异值分解(SVD)的应用——图像压缩
我们想寻找一个A的逼近:Ak,使得rank(Ak) = k < n,且|A - Ak|最小。下面的定理(也称为Schmidt-Mirsky, Eckart-Young定理。原创 2024-02-10 20:41:04 · 910 阅读 · 1 评论 -
矩阵的正定(positive definite)性质的作用
注意,本文中正定和半正定矩阵不要求是对称或Hermite的。原创 2024-02-06 20:09:05 · 1314 阅读 · 1 评论 -
特征值求解——重启Arnoldi方法和收缩(deflation)技术
Arnoldi向量转换为Schur向量,注意,只作用在新近收敛的向量和下一个初始向量上。注意上面两个算法的区别,Algorithm 2相比1唯一的不同是,循环是从k+1开始的,也就是只改变Vm和Hm的最后m-k列,即active部分。在下一次的arnoldi过程中,只有m-k个arnoldi向量被计算,也就是active vectors。,我们将这个R的子矩阵取出来记为Rk,对应的Q中的前k列取出来记为Qk,那么我们有。注意,S≡Range(Qk)是A的不变子空间,对应着k个特征值,原创 2023-12-17 15:38:45 · 430 阅读 · 0 评论 -
数值计算project——ODE求解器
一、介绍这篇文章介绍一个ODE solver的实现,可用于一阶线性系统的求解。首先,找到一组ODE方程,可以从网上找算例,或者从comsol等软件中导出。其次,用matlab来求解,看计算结果。最后,自己实现一个ODE solver,并将结果与matlab的结果进行对比。...原创 2023-06-15 23:27:00 · 179 阅读 · 0 评论 -
Krylov子空间与Arnoldi过程
一、应用背景Krylov子空间方法是模型降阶(MOR)的一种,Krylov 子空间方法就是在Krylov子空间中寻找近似解。二、Krylov子空间n是原始矩阵的维度,m为降阶后矩阵的维度,通常m<<n。三、Arnoldi 过程: 计算Km 的一组正交对应上述MGS过程的python代码:import numpy as npdef arnoldi_itera...原创 2023-06-12 23:37:00 · 2087 阅读 · 3 评论 -
指数积分方法(Exponential Integration)求解ODE/DAE问题
1. 简介We can classify an integration method by its stability and computational effort. The y-axis represents the stability and x-axis represents the computational effort.The backward euler is locate...原创 2023-06-23 12:21:00 · 271 阅读 · 0 评论 -
A First course in FEM —— matlab代码实现求解传热问题(瞬态)
这一篇Blog是在A First course in FEM —— matlab代码实现求解传热问题(稳态)基础上更进一步,求解瞬态传热问题。两者的区别如下图所示:1. 问题描述求解下图图所示叶片的温度场在[0-1200s]时间段内的变化,初始条件:T(0)=25℃。控制方程为:2. 模型和网格模型和网格设置详见A First course in FEM —— matlab代码...原创 2023-06-22 17:06:00 · 645 阅读 · 2 评论 -
似然与极大似然估计
一、似然在统计学中,似然性(likelihood)”和“概率”有明确的区分:概率,用于在已知一些参数的情况下,预测接下来在观测上所得到的结果;似然性,则是用于在已知某些观测所得到的结果时,对有关事物之性质的参数进行估值。以高斯分布为例,其可以用参数μ和σ来描述。采样和参数估计是互逆的过程,从分布中采样是已知一些参数,得到观测结果,结果出现的可能性就用“概率”来表示。而在已知猜测结果时,对...原创 2023-08-06 20:15:00 · 63 阅读 · 0 评论 -
FEniCSx入门——解Poisson方程
一、 Poisson方程可转换为线性弱形式其中,u: trial function 是近似解v: test function 是Poisson方程两边同乘的函数,用于转换为弱形式。在Fenics中解PDE时,必须进行下面几个步骤:1、将PDE问题转换成离散变分问题:即寻找u∈V,使得2、选择V和V_hat空间,换句话说就是确定网格和有限元单元的类型。3、定义Inp...原创 2023-05-14 19:59:00 · 470 阅读 · 0 评论 -
A First course in FEM —— matlab代码实现求解传热问题(稳态)
这篇文章会将FEM全流程走一遍,包括网格、矩阵组装、求解、后处理。内容是大三时的大作业,今天拿出来回顾下。1. 问题简介涡轮机叶片需要冷却以提高涡轮的性能和涡轮叶片的寿命。我们现在考虑一个如上图所示的叶片,叶片处在一个高温环境中,中间通有四个冷却孔。假设为稳态,那么叶片内导热微分方程为:内部区域: (扩散方程)边界:(外表面)(内部冷却孔)2.模型2.1几何模型...原创 2023-06-18 15:24:00 · 332 阅读 · 2 评论 -
最优化问题简介及优秀教材《凸优化》介绍
前言最优化广泛应用于科学与工程计算、数据科学、机器学习、人工智能、图像和信号处理、金融和经济、管理科学等众多领域。最优化问题可以归纳为如下定义:最优化问题一般很难求解,除了一些特例。目前已经发展成熟的,能够有效求解的最优化问题可以归为以下三类:最小二乘问题 least-squares problems线性规划问题linear programming problems凸优化问题...原创 2023-07-02 23:28:00 · 325 阅读 · 0 评论