代码1:TDSE_cn1d.m % Crank-Nicolson 法解一维薛定谔方程% 等间距网格,稀疏矩阵functionTDSE_cn1d% === 参数设置 ===xmin=-80;xmax=80;Nx=1000;% x 网格tmin=0;tmax=20;Nt=400;% 时间网格Nplot=10;% 画图步数ax=[xmin,xmax,-0.5,0.5];% 高斯波包x0=-17;t0=0;m=1;% 高斯波包的初始时间p0...
Matlab:Crank Nicolson方法求解线性抛物方程 1tic;2clear3clc4M=[10,20,40,80,160,320,640];%x的步数5K=M; %时间t的步数6forp=1:length(M)7hx=1/M(p);8ht=1/K(p);9r=ht/hx^2; %网格比10x=0:hx:1;11t=0:ht:1;12numerical=zeros(M(p)+1,K(p)+1);13numerical(:,1)=exp(x); %...
Crank-Nicolson 方法的离散形式如下，需注意实际实现时需解线性方程组，涉及矩阵运算等操作，具体过程不在此赘述。通过 MATLAB 语言实现上述离散形式，以代码示例的形式呈现 Crank-Nicolson 方法在求解一维非稳态导热微分方程中的应用，旨在展示从理论到实践的转换过程。
同时增加N_x与N_t,保持r=0.5,来看一组结果: 这么一看Heun的精度比Euler好一些,而Crank-Nicolson精度最好 Crank-Nicolson格式作为隐格式,有无条件稳定的特性,也就是不再受到\Delta t不能过大的影响,比如以下实验,固定住\Delta t, 减小\Delta x, 结果如下: CN格式还是非常稳的 代码在这里: importnumpyasnpim...
求解一维热传导方程Crank-Nicolson差分法
A class of alternating segment Crank-Nicolson methods for solving convection-diffusion equations In this paper, we give a class of Alternating Segment Crank-Nicolson (ASC-N) method for the convection-diffusion equation. The method is unconditionally st... Wang WQ - 《Computing Archives for Informat...
目录 摘要1 1.前言3 2.Crank-Nicolson差分法4 2.1)差分法定义4 2.2)差分格式的建立4 2.3)Crank-Nicolson差分格式(六点格式)7 2.4)Crank-Nicolson差分格式的向量表示9 2.5)Crank-Nicolson差分格式的稳定性11 2.6)Crank-Nicolson差分格式的收敛性14 3.数值算例17 3.1)利用Crank-Nicolson方法求解数值算例17 4.总结...
Crank-Nicolson有限差分法对欧式障碍期权定价,以欧式看跌期权为例,运用Matlab编程, 将所得两个结果与基于公式解的结果进行比较,结果表明Crank-Nicolson有限差分法比 Monte-Carlo方法更优。 关键词:Monte-Carlo;障碍期权;有限差分方法;欧式看跌期权 中图分类号:F832.51 ...
最后,建立拟均相二维模型方程,采用Crank-Nicholson隐式差分法,利用MATLAB编程对均四甲苯氧化反应固定床反应器进行了模拟,模拟结果与工业数据较为吻合,进一步模拟分析... 曹茂洪 被引量: 0发表: 2014年 光气合成工艺的研究与建模 应用Matlab编程,克兰克—尼科尔森算法(Crank-Nicolson)数值分析,建立了基于反应动力学方程...
有限差分法求解一维扩散方程。使用Crank-Nicolson格式求解PDE方程,并与真实解析解u(tx)=e^(x+0.1t)比较。 一维扩散方程的有限差分法,采用隐式六点差分格式。 上传者:weixin_42691388时间:2021-09-29 LAB12_EDP:使用 Crank-Nicolson 方法求解抛物线方程-matlab开发 ...