⭐ 欢迎来到虫虫下载站！ | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站
📤 上传资源

📄 thtab66.m

📁 ％ Atomizer Main Directory, Version .802 里面信号含有分解去噪合成过程的代码 ％---------------------------------------
💻 M
📖 第 1 页 / 共 3 页
字号:
🔍
123 
% thtab66: BP Thesis Table 6.6 -- MOF, BOB, MP  Complexity
% Signal:
%	TwinSine-2	4-overcomplete DCT
%	Carbon		WP
%	WernerSorrows	CP D-4
%	FM-Cosine	CP D-4
%	Gong		CP D
%	Noisy TwinSine-2
%	Noisy Carbon
%	Noisy WernerSorrows
%	Noisy FM-Cosine
%	Noisy Gong
% Signal Length:
%	64, 128, 256, 512, 1024
% Parameter Setting: 1e-1, 1e-2, 1e-3
% Observation:
%	CPU times in seconds recorded in Table 6.6
% Warning: This takes hours to run
%

help thtab66

%
% (1) MOF
%

%	
% Sample Size: n = 64
% 
n = 64; D= log2(n);
qmf = MakeONFilter('Symmlet', 8);
natom = n^2;

x = InputSignal('TwinSine2', n);
t = cputime;
c1 = MOF(x, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c3 = MP(x, 'DCT', 4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);

x = InputSignal('Carbon', n);
t = cputime;
c1 = MOF(x, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'WP', D, qmf, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('Carbon')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('WernerSorrows', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('WernerSorrows')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('FM', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('FM-Cosine')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('Gong', n);
t = cputime;
c1 = MOF(x, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('TwinSine2', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'DCT', 4, 0, 0);
cpu3 = cputime - t;
disp('Noisy TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('Carbon', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'WP', D, qmf, 0);
cpu3 = cputime - t;
disp('Noisy Carbon');
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);



SNR = 1;
x = InputSignal('WernerSorrows', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D-4, 0, 0);
cpu3 = cputime - t;
disp('Noisy WernerSorrows')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('FM', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D-4, 0, 0);
cpu3 = cputime - t;
disp('Noisy FM-Cosine')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
t = (1:n)' / n;
x = cos(t * pi/4 * n) ./ ((1:n) .^ .5)';
m = floor(3 *n/4);
x = [zeros(n-m,1);x(1:m)];
%Noise level = 1
randn('seed', 2.055615866000000e+09);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D, 0, 0);
cpu3 = cputime - t;
disp('Noisy Gong')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);

%	
% Sample Size: n = 128
% 

n = 128; D= log2(n);
qmf = MakeONFilter('Symmlet', 8);
natom = n^2;

x = InputSignal('TwinSine2', n);
t = cputime;
c1 = MOF(x, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c3 = MP(x, 'DCT', 4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);

x = InputSignal('Carbon', n);
t = cputime;
c1 = MOF(x, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'WP', D, qmf, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('Carbon')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('WernerSorrows', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('WernerSorrows')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('FM', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('FM-Cosine')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('Gong', n);
t = cputime;
c1 = MOF(x, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('TwinSine2', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'DCT', 4, 0, 0);
cpu3 = cputime - t;
disp('Noisy TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('Carbon', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'WP', D, qmf, 0);
cpu3 = cputime - t;
disp('Noisy Carbon');
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);



SNR = 1;
x = InputSignal('WernerSorrows', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D-4, 0, 0);
cpu3 = cputime - t;
disp('Noisy WernerSorrows')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('FM', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D-4, 0, 0);
cpu3 = cputime - t;
disp('Noisy FM-Cosine')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
t = (1:n)' / n;
x = cos(t * pi/4 * n) ./ ((1:n) .^ .5)';
m = floor(3 *n/4);
x = [zeros(n-m,1);x(1:m)];
%Noise level = 1
randn('seed', 2.055615866000000e+09);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D, 0, 0);
cpu3 = cputime - t;
disp('Noisy Gong')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);

%	
% Sample Size: n = 256
% 

n = 256; D= log2(n);
qmf = MakeONFilter('Symmlet', 8);
natom = n^2;

x = InputSignal('TwinSine2', n);
t = cputime;
c1 = MOF(x, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c3 = MP(x, 'DCT', 4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);

x = InputSignal('Carbon', n);
t = cputime;
c1 = MOF(x, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'WP', D, qmf, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('Carbon')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('WernerSorrows', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('WernerSorrows')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('FM', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('FM-Cosine')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('Gong', n);
t = cputime;
c1 = MOF(x, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('TwinSine2', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'DCT', 4, 0, 0);
cpu3 = cputime - t;
disp('Noisy TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('Carbon', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'WP', D, qmf, 0);
cpu3 = cputime - t;
disp('Noisy Carbon');
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);



SNR = 1;
x = InputSignal('WernerSorrows', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D-4, 0, 0);
cpu3 = cputime - t;
disp('Noisy WernerSorrows')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('FM', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D-4, 0, 0);
cpu3 = cputime - t;
disp('Noisy FM-Cosine')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
t = (1:n)' / n;
x = cos(t * pi/4 * n) ./ ((1:n) .^ .5)';
m = floor(3 *n/4);
x = [zeros(n-m,1);x(1:m)];
%Noise level = 1
randn('seed', 2.055615866000000e+09);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'CP', D, 0, 0);
cpu3 = cputime - t;
disp('Noisy Gong')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOBDN = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);

%	
% Sample Size: n = 512
% 

n = 512; D= log2(n);
qmf = MakeONFilter('Symmlet', 8);
natom = n^2;

x = InputSignal('TwinSine2', n);
t = cputime;
c1 = MOF(x, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c3 = MP(x, 'DCT', 4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);

x = InputSignal('Carbon', n);
t = cputime;
c1 = MOF(x, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'WP', D, qmf, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'WP', D, qmf, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('Carbon')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('WernerSorrows', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('WernerSorrows')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('FM', n);
t = cputime;
c1 = MOF(x, 'CP', D-4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D-4, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D-4, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
disp('FM-Cosine')
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


x = InputSignal('Gong', n);
t = cputime;
c1 = MOF(x, 'CP', D, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
c2 = BOB(x, 'CP', D, 0, 0);
cpu2 = cputime - t;
t = cputime;
c3 = MP(x, 'CP', D, 0, 0, natom, 1e-2);
cpu3 = cputime - t;
fprintf('\n')
fprintf('CPU Running Time of MOF = %8.4e\n', cpu1);
fprintf('CPU Running Time of BOB = %8.4e\n', cpu2);
fprintf('CPU Running Time of  MP = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('TwinSine2', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'DCT', 4, 0, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec3 = MPDN(y, 'DCT', 4, 0, 0);
cpu3 = cputime - t;
disp('Noisy TwinSine-2');
fprintf('\n')
fprintf('CPU Running Time of MOFDN = %8.4e\n', cpu1);
fprintf('CPU Running Time of  MPDN = %8.4e\n', cpu3);


SNR = 1;
x = InputSignal('Carbon', n);
randn('seed', 12345678);
[x y]= NoiseMaker(x, SNR);
t = cputime;
xrec1 = MOFDN(y, 'WP', D, qmf, 0, 1e-10);
cpu1 = cputime - t;
t = cputime;
xrec2 = BOBDN(y, 'WP', D, qmf, 0);
123

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -