romberg1.m

基于极大嫡谱估计准则对变形数据进行预测 基于MATLAB平台编制了相应程序

function q = romberg1(a,b,n,randa)
h=(b-a)/n;
q=0;
x=a;
for i=1:n
    x=x+(n-1)*h;
    x0=x;
    x1=x+h/8;
    x2=x+2*h/8;
    x3=x+3*h/8;
    x4=x+4*h/8;
    x5=x+5*h/8;
    x6=x+6*h/8;
    x7=x+7*h/8;
    x8=x+8*h/8;
    zz=h/5670*(441*fx1(x0,randa)+2048*fx1(x1,randa)+736*fx1(x2,randa)+2048*fx1(x3,randa)+2036*fx1(x4,randa)+2048*fx1(x5,randa)+736*fx1(x6,randa)+2048*fx1(x7,randa)+441*fx1(x8,randa));
    q=q+zz;
end