CV老李的博客

运用递归的方法求解对称三对角矩阵的特征值

Neuman边界条件对这样一个问题我的目标：我们有如下的效果：一阶近似t=0.06t = 0.1t = 0.9二阶近似1.t=0.62.t=0.13.t = 0.9误差的差别改变步长取t=0.91.M=20,dt = 0.001一阶二阶差别2.M=40,dt=0.0003，t=0.6第二道题t = 0.9...

微分方程组求解（Euler法、RK法、线性多步法）

解一道题格式的表达代码示例：%% 用三种方法求解公式h = 0.01; t = 0:h:1;N = length(t);y = zeros(1,N);y(1) = 1;df = @(t,y)(-y+t+1);%改进Euler法(迭代方法)w = y(1);for k = 1:N-1 wk = w + h*df(t(k), w); z1 = wk; z...

日期：3030-3-29二阶泰勒展开法求解微分方程1.

日期：2020-3-16显式欧拉格式格式为：我运用显式欧拉格式计算三道例题，然后展现其效果，最后将我的代码展示出来1.2.3.呈现出这样的曲线的原因是该函数由一开始的极其陡峭变为后来的非常平稳，其中斜率的变化率非常之大。而欧拉格式走的是直线，在一个步长内斜率是不发生变化的，也就是当这个折线走完一个步长，与真实的曲线之间在数值上已经有很大的差别了。后面随着曲线的陡峭程...

作者：老李日期：2020-3-15用几种欧拉格式算一道题在这篇报告中，我先使用了显示欧拉格式来求解微分方程，隐式欧拉格式中，我使用了改进的欧拉格式和修正的欧拉格式来进行求解。这道题是这样的：我们要做的是，通过已有的信息，求出 y(3)。我使用matlab来实现这个事情首先设立参数和函数句柄代码如下%% function and parament% we are going t...