p1_3.m

来自「Digital Signal Processing-A Computer-Bas」· M 代码 · 共 22 行

% Program P1_3.m
% Generation of a random sinusoidal signal
clear all;close all
N = 5; 
A = 4*rand(1,N);
phase = 2*pi*rand(1,N);
n = 0:30;   
f = 0.1;

for i = 1:N;
arg(i,:) = 2*pi*f*n + phase(i); 
x(i,:) = A(i)*cos(arg(i,:));

figure(i);
stem(n,x(i,:));      % Plot the generated sequence
axis([0 30 -4 4]);
grid; 
title('Random Sinusoidal Sequence');
xlabel('Time index n');
ylabel('Amplitude');
axis;
end