%% 线法 D = 扩散率:菲克第二扩散定律% by Prof. Roche C. de Guzman %% 自定义外汇函数 Y1 = F(~,Y,D,nx,dx) c = Y; % 指定浓度为 Y Y1 = 零点(nx,1); % 临时衍生品对于 i = 1:nx-2 % 位置计数器Y1(i+1) = D*(c(i+2)-2*c(i+1)+c(i))/dx^2; %求解导数结尾Y1 = [Y1(2); Y1(2:nx-1); Y1(nx-1)]; % 具有零通量边界的导数结尾
2021-08-30 22:56:32
2KB
matlab
1