program_09_5.m

来自「自编的一些小波源程序,用MATLAB编的.希望对那些喜欢小波的人有所帮助.去解压」· M 代码 · 共 18 行

% Program 9_5
% Computation of Approximate Output round-off Noise
%
num = input('Numerator coefficients = ');
den = input('Denominator coefficients = ');
x = 1;
order = max(length(num),length(den))-1;
si = [zeros(1,order)];
varnew = 0; k = 1;
while k > 0.000001
   [y,sf] = filter(num,den,x,si);
   si = sf;
   varold = varnew;
   varnew = varnew + abs(y)*abs(y);
   k = varnew - varold;
   x = 0;
end
disp('Output Noise Variance = ');disp(varnew)