龙贝格算法主要是不断递推和加速,直到满足精度要求 递推: 加速: 得到T表: MATLAB代码: function I = Romberg(f, a, b, epsilon) I = 0; h = b-a; k = 0; m = 0; T = zeros(5); %下标转换:T^(k)_0 => T(k+1,1) T(1,1) = h/2 * (subs(f,a) + subs(f,b));%即T^(0)_0 delta = 2*epsilon; while delta > epsilon k
取半径=3 用matlab代码实现上式公式: length=3;for Ki = 1:length for Kj = 1:length for Kk = 1:length Ksigma(Ki,Kj,Kk)=exp(-(Ki-2)^2/8-(Kj-2)^2/8-(Kk-2)^2/8); 此公式为:K(u), ρ=3 end endend KONE = convn(ones(size(Img3D_mc,1),size(Img3D_mc,