wt10fig07.m

来自「斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善」· M 代码 · 共 48 行

%CAPTION :
fprintf('\n');
disp('Figure 10.7')
disp('Window 1: Blurred signal plus white noise (SNR = 26.8db)')
disp('Window 2: Deconvolution with a Wiener filter (SNR = 28.4db).')

close all;
n=4096;

sig=MakeSignal('Piece-Polynomial',n);

sigma=0.245;
Noise=GWN(n,sigma);
[g,I]=convgauss(Noise,-1);

convsig=iconv(g,sig);
Noisy=convsig+Noise;

[H,out_wien] = Wienerdec(g,Noisy,sig,sigma);

figure(1);clf
set(gcf, 'NumberTitle','off', 'Name','Window 1')
plot(Noisy);
handle_axis = gca; % recupere le handle de l'axe actif
set(handle_axis,'FontSize', 21);
axis([0 n -100 200]);

figure(2);clf
set(gcf, 'NumberTitle','off', 'Name','Window 2')
plot(out_wien);
handle_axis = gca; % recupere le handle de l'axe actif
set(handle_axis,'FontSize', 21);
axis([0 n -100 200]);

% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]

    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%