在工程和科学领域,特别是在热传导、扩散、波动等方程的求解中,Crank-Nicolson方法被广泛应用。 在使用Crank-Nicolson方法进行求解时,通常需要编写相应的数值程序。Matlab作为一种功能强大的数值计算工具,被广泛应用于工程和科学领域。编写Crank-Nicolson方法的Matlab程序,可以帮助工程师和科学家更方便地求解偏微分方程,提高...
Use the Crank-Nicolson Method. We need to discretize the space and time domain. x_i = i h = i \Delta x, \quad t_n = nk = n \Delta t, \quad U^n_i \approx u(x_i,t_n) \frac{U^{n+1}_i - U^{n}_i}{k} = \frac{\kappa}{2h^2}(U^n_{i-1} - 2U^n_i + ...
代码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); %...
Robalo. The Crank-Nicolson-Galerkin finite element method for a non- local parabolic equation with moving boundaries. Numerical Methods for Partial Differential Equations, 2014. doi: 10.1002/num.21957.Rui M. P. Almeida, Jose麓 C. M. Duque, Jorge Ferreira, and Rui J. Robalo. The Crank-...
我们有兴趣使用 CN 方法获得一维热传导方程的稳态解。 边界条件是:在 x=0 和 0.3 m 处 T=300 K,在所有其他内部点处 T=100 K。 α = 〖3*10〗^(-6) m-2s-1 . 这里,t=30 分钟,Δx=0.015m 和Δt=20 秒 (0)踩踩(0) 所需:1积分 ...
In this repository, the lid-driven cavity problem is solved using the Crank-Nicolson/Adams-Bashforth scheme for viscous and convection terms, respectively. For spatial discretization, the second-order central difference scheme is used in a uniform staggered grid. The projection method is used to so...