sine_art.m

matlab中噪声模型

function [y] =sine_art(N,f)
% [y] = sine_art(N, f);
% N  = length of pulse (in samples)
% f = frequency in Hz
%
%   See also: SINC_ART, SINC SQUARE, SIN, COS, CHIRP, DIRIC, GAUSPULS, 
%             PULSTRAN, RECTPULS SIGMOID_ART and TRIPULS  

%   Author: G. Clifford, 8-10-03
%   Copyright 2003 MIT

if nargin < 2
  f = 5;    
end
if nargin <1
    N = 256;    
end

y=ones(1,N);
i=find(y);
% 'normalise' period
ii = i/N;

% calculate one half of sinc function
y(i)=sin(2*pi*f*ii(i));

if nargin <1
    figure;
    plot(y);
end