myconvpp.m

来自「matlab实习卷积运算。很高兴与大家分享此代码。」· M 代码 · 共 28 行

function y = myconvpp( n )

%maxvalue = 100;
minvalue = 0;
%epsilon = 0.0001;
%MYCONV Summary of this function goes here
%   Detailed explanation goes here
y = inf;
maxvalue = minvalue;
while (maxvalue < 10000000000)
    a = quadv(@(x) funcA(x) .* funcB(n-x), -maxvalue, maxvalue);
    b = quadv(@(x) funcA(x) .* funcB(n-x), -1.2.*maxvalue, 1.2.*maxvalue);
    c = quadv(@(x) funcA(x) .* funcB(n-x), -0.8.*maxvalue, 0.8.*maxvalue);
    epsilon = 0.01 .* a;
    
    if ((a - b < epsilon) && (b - a < epsilon) && (a - c < epsilon) && (c - a < epsilon) && a > 0)
        y = a;
        break;
    end;
        
    if (maxvalue == 0) maxvalue = 1;
    else maxvalue = 2 .* maxvalue;
    end;
end;
maxvalue

%else y = inf;
%end;