% Signal parameters:f = [ 440 880 1

% Signal parameters:
f = [ 440 880 1000 2000 ]; % frequencies
M = 256; % signal length
Fs = 5000; % sampling rate

% Generate a signal by adding up sinusoids:
x = zeros(1,M); % pre-allocate 'accumulator'
n = 0:(M-1); % discrete-time grid
for fk = f;
x = x + sin(2*pi*n*fk/Fs);
end
% Filter parameters:
L = 257; % filter length
fc = 600; % cutoff frequency

% Design the filter using the window method:
hsupp = (-(L-1)/2:(L-1)/2);
hideal = (2*fc/Fs)*sinc(2*fc*hsupp/Fs);
h = hamming(L)' .* hideal; % h is our filter
% Choose the next power of 2 greater than L+M-1
Nfft = 2^(ceil(log2(L+M-1))); % or 2^nextpow2(L+M-1)

% Zero pad the signal and impulse response:
xzp = [ x zeros(1,Nfft-M) ];
hzp = [ h zeros(1,Nfft-L) ];

% Transform the signal and the filter:
X = fft(xzp);
H = fft(hzp);
Y = X .* H;
y = ifft(Y);
relrmserr = norm(imag(y))/norm(y) % check... should be zero
y = real(y);

ผลลัพธ์ (ไทย) 2:

% พารามิเตอร์สัญญาณ:
f = [440 880 1000 2000] % ความถี่
M = 256; สัญญาณ% ระยะเวลาใน
Fs = 5000; % อัตราการสุ่มตัวอย่าง% สร้างสัญญาณโดยการเพิ่มขึ้นไซน์: x = ศูนย์ (1, M); % ก่อนการจัดสรร 'สะสม' n = 0 (M-1); % ตารางเวลาที่ไม่ต่อเนื่องสำหรับfk f =; x = x + บาป (2 * * * * * * * * ปี่ n * fk / Fs); ท้ายพารามิเตอร์% กรอง: L = 257; กรอง% ระยะเวลาfc = 600; % ตัดความถี่% การออกแบบตัวกรองโดยใช้วิธีหน้าต่าง: hsupp = (- (L-1) / 2: (L-1) / 2); hideal = (2 * fc / Fs) * sinc (2 * fc * hsupp / Fs); h. = แฮม (L) * hideal; % h คือกรองของเรา% เลือกพลังงานต่อไปของ 2 มากกว่า L + M-1 Nfft = 2 ^ (ceil (log2 (L + M-1))); % หรือ 2 ^ nextpow2 (L + M-1)% ศูนย์แผ่นสัญญาณและกระตุ้นการตอบสนอง: xzp = [ศูนย์ x (1, Nfft-M)]; hzp = [ศูนย์เอช (1, Nfft-L)]% แปลงสัญญาณและตัวกรอง: X = FFT (xzp); H = FFT (hzp);. Y = X * H; y = IFFT (Y); relrmserr บรรทัดฐาน = (imag (y)) / บรรทัดฐาน (y)% การตรวจสอบ ... ควรจะเป็นศูนย์Y = จริง (y);

ผลลัพธ์ (ไทย) 3:

พารามิเตอร์สัญญาณ :
% f = [ 440 880 1000 2000 ] ; % ความถี่
M = 256 ; % สัญญาณความยาว
FS = 5000 ; %

% อัตราการสุ่มสร้างสัญญาณโดยการเพิ่มขึ้น sinusoids :
x = 0 ( , 1 ) ; % ก่อนจัดสรร ' สะสม '
n = 0 : ( m-1 ) ;
% เวลาตารางสำหรับ FK = f ;
x = x sin ( 2 * pi * N * FK / FS ) ;
จบ
% ค่าพารามิเตอร์ตัวกรอง :
L = 257 ; % กรองความยาว
FC = 600 ; %

ตัดความถี่การออกแบบตัวกรองโดยใช้วิธีหน้าต่าง :
hsupp = ( - ( L-1 ) / 2 ( L-1 ) / 2 ) ;
hideal = 2 * FC / FS ( 2 ) * * * * * * hsupp sinc FC / FS ) ;
H = แฮม ( L ) * hideal ; % H คือ %
ของเรากรองเลือกพลังงานต่อไป 2 มากกว่า L M - 1
nfft = 2
( CEIL ( LOG ( L M - 1 ) ) ) ) ) ) ) ; % หรือ 2
nextpow2 ( L M - 1 )
% ศูนย์แผ่นสัญญาณและแรงกระตุ้นการตอบสนอง :
xzp = [ 0 X ( 1 nfft-m ) ] ;
hzp = [ 0 H ( 1 nfft-l ) ] ;

แปลงสัญญาณและตัวกรอง :
x = หน่วย ( xzp ) ;
H = หน่วย ( hzp ) ;
y = x * H ;
y = ifft ( y ) ;
relrmserr = บรรทัดฐาน ( IMAG ( Y ) ( Y ) C ) / เช็ค . . . . . . . ควรจะเป็นศูนย์
Y ;
( Y ) = จริง