本节的核心是将常系数微分方程转化为线性代数问题. \[\frac{du}{dt}=\lambda u \quad 的解为 \quad u(t) = Ce^{\lambda t}\] 代入 \(t=0\),可得 \(u(0) = C\),因此有 \(u(t) = u(0)e^{\lambda t}\).这是只有一个变量的情况,在线性代数里,我们扩展到 \(n\) 个方程的情况. \[\frac{d\boldsymbol u}{dt}=A \boldsymbol u \quad 初始条件为向量 \q…
原文:https://mp.weixin.qq.com/s/COpYKxQDMhqJRuMK2raMKQ 微分方程指含有未知函数及其导数的关系式,解微分方程就是找出未知函数.未知函数是一元函数的,叫常微分方程:未知函数是多元函数的,叫做偏微分方程.常微分方程有时也简称方程.微分方程是一门复杂的学科,对于常微分方程来说,可以使用特征值和特征向量的知识求解. 相关前置知识: 微分方程:单变量微积分11——常微分方程和分离变量 泰勒公式:单变量微积分30——幂级数和泰勒级数 泰勒公式在0点展开的原因:…
title: [线性代数]6-3:微分方程的应用(Applications to Differential Equations) categories: Mathematic Linear Algebra keywords: Eigenvalues Eigenvectors Differential Equations toc: true date: 2017-11-22 15:09:04 Abstract: 本文主要介绍线性代数在微分方程中的应用 Keywords: Eigenvalues,E…
基本运算: 一,矩阵的生成 clc ; clear all; close all; 1.直接输入 A = [ 1 ,2 ,3,4;2,3,4,5;3,4,5,6] A = 1 2 3 4 2 3 4 5 3 4 5 6 1)矩阵元素可以是表达式 B = [sin(pi/2),4*8] B = 1 32 2)矩阵元素可以是复数 C = [3+2i,2+6i,4-3i] C = 3.0000 + 2.0000i 2.0000 + 6.0000i 4.0000 - 3.0000i 2.使用函数建立矩阵…
背景: Blei的hlda的C语言实现需要使用C语言的科学计算包GSL,因此决定安装.由于在windows下安装极其繁琐,先在Linux上安装之. 系统环境: Linux version 2.6.35-22-generic (buildd@allspice)  (gcc version 4.4.5 (Ubuntu/Linaro 4.4.4-14ubuntu4) )  #33-Ubuntu SMP Sun Sep 19 20:32:27 UTC 2010 gsl安装过程 1, 下载 gsl-2.0…
GSL(GNU Scientific Library)作为三大科学计算库之一,除了涵盖基本的线性代数,微分方程,积分,随机数,组合数,方程求根,多项式求根,排序等,还有模拟退火,快速傅里叶变换,小波,插值,基本样条,最小二乘拟合,特殊函数等.下面介绍一下GSL的安装和使用. 方法一: 首先从官网下载到源代码(我用的版本是 gsl-1.9)压缩包,解压后进入目录,执行 ./configure make make install 这个过程需要几分钟.这里还有一点需要注意的是,执行 make inst…
注: 此系列为自己之前所搭建网站内容. 由于论文数据处理的需要,需要使用libeemd这个包,需要安装gsl科学库,windows下没有办法,只能转战ubuntu进行科学计算. GSL(GNU Scientific Library)作为三大科学计算库之一,除了涵盖基本的线性代数,微分方程,积分,随机数,组合数,方程求根,多项式求根,排序等,还有模拟退火,快速傅里叶变换,小波,插值,基本样条,最小二乘拟合,特殊函数等 当然最希望通过命令来安装GSL科学库. 方法一: GSL源码包提供了以下二进制包…
Part A: The Issues 议题 第一章 软件品质 第二章 面向对象的标准 Part B: The Road To Object Orientation 通向面向对象之路 第三章 模块性 第四章 复用性方法 第五章 走进对象技术 第六章 抽象数据类型 Part C: Object-Oriented Techniques 面向对象技术 第七章 静态结构: 类 第八章 运行时结构: 对象 第九章 内存管理 第十章 泛型 第一章 软件品质 1.1 外部和内在的因素 1.2 外部因素介绍 1.…
MATLAB版本:R2015b 1.求解符号矩阵的行列式.逆.特征值.特征向量 A = sym('[a11, a12; a21, a22]');deltaA = det(A)invA = inv(A)[V, D] = eig(A) %V的列向量为特征向量,D的主对角线元素为相应的特征值 2.求解代数方程的解析解 syms a b cx = solve('a * x^2 + b * x + c = 0', 'x') 3.求解微分方程(组)的解析解 syms x yY1 = dsolve('D2y…
举例:分别用欧拉法和龙哥库塔法求解下面的微分方程 我们知道的欧拉法(Euler)"思想是用先前的差商近似代替倒数",直白一些的编程说法即:f(i+1)=f(i)+h*f(x,y)其中h是设定的迭代步长,若精度要求不高,一般可取0.01.在定义区间内迭代求解即可.龙哥库塔法一般用于高精度的求解,即高阶精度的改进欧拉法,常用的是四阶龙哥库塔,编程语言如下:y(i+1)=y(i)+h*(k1+2*K2+2*k3+k4)/6;k1=f(xi,yi)k2=f(xi+h/2,yi+h*k1/2);…