pulse_sine.m

来自「matlab中噪声模型」· M 代码 · 共 42 行

function art = pulse_sine(lenArtS, sampF, freqArt, ampArt)
%
% art = pulse_sine(lenArtS, sampF, freqArt, ampArt);
%
% generate a sine pulse of length lenArtS at freqArt with a sampling
% frequency of freqArt and a normalise amplitude of ampArt.
% Requires sine_art.m
%
% G Clifford 2003 [e-mail] [GNU GPL]

if nargin<4
 ampArt   = 0.5; % normalised amplitude of art 
end
if nargin<3
 freqArt = 5; % hz
end
if nargin<2
 sampF = 256;
end
if nargin<1
 lenArtS  = 0.75; % length in seconds
end

lenArt   = round(sampF*lenArtS);% length in samples

y= sine_art(lenArt,freqArt);  % generate sine wave
% generate a half sine wave to modulate it
x = sine_art(lenArt, (sampF*lenArtS)/(2*lenArt));
%x = sine_art(lenArt, lenArtS/2); 
art = (x.*y)*ampArt;          % modulate and scale

if nargin<1
figure;
subplot(2,3,1)
     plot(y);
subplot(2,3,2)
     plot(x);
subplot(2,3,3)
     plot(art);
subplot(2,1,2)
end