在工程和科学领域,特别是在热传导、扩散、波动等方程的求解中,Crank-Nicolson方法被广泛应用。 在使用Crank-Nicolson方法进行求解时,通常需要编写相应的数值程序。Matlab作为一种功能强大的数值计算工具,被广泛应用于工程和科学领域。编写Crank-Nicolson方法的Matlab程序,可以帮助工程师和科学家更方便地求解偏微分方程,提高...
代码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); %...
MATLAB code: function [error,h,k] = heateq_CN(a,b,ti,tf,m,kappa) % Discretization set up h = abs(b-a)/(m+1); % h = delta x k = 4*h; % time step k x = linspace(a,b,m+2)'; % x(1)= a and x(m+2)= b step = round((tf-ti)/k); % approx. number of time...
我们有兴趣使用 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积分 ...
M. Duque, Jorge Ferreira, and Rui J. Robalo. The Crank-Nicolson-Galerkin finite element method for a nonlocal parabolic equation with moving boundaries. Available from: http://ptmat.fc.ul.pt/arquivo/docs/preprints/pdf/2013/Jorge Ferreira Almeida Duque preprint 017 2013.pdf, 2013....
About Numerical analysis of lid-driven cavity flow using the Crank-Nicolson/Adams-Bashforth scheme Resources Readme Activity Stars 0 stars Watchers 1 watching Forks 0 forks Report repository Releases No releases published Packages No packages published Languages MATLAB 100.0% ...