《DSP using MATLAB》示例Example5.8

代码：

n = [0:1:99]; x = cos(0.48*pi*n) + cos(0.52*pi*n);
n1 = [0:1:9]; y1 = x(1:1:10);                          % N = 10

figure(‘NumberTitle‘, ‘off‘, ‘Name‘, ‘Exameple5.8 x sequence‘)
set(gcf,‘Color‘,‘white‘);
subplot(3,1,1); stem(n1,y1); title(‘signal x(n), 0<= n <=9, N = 10‘); axis([0,10,-2.5,2.5]);
xlabel(‘n‘); ylabel(‘x(n)‘); grid on;

Y1_DFT = dft(y1,10);                                   % DFT of y1

magY1_DFT = abs(Y1_DFT(1:1:6));
%magY1_DFT = abs(Y1_DFT);
%phaY1_DFT = angle(Y1_DFT)*180/pi                        % degrees
phaY1_DFT = angle(Y1_DFT(1:1:6))*180/pi                  % degrees
%realX_DFT = real(X_DFT); imagX_DFT = imag(X_DFT);
%angX_DFT = angle(X_DFT);                               % radias

k1 = 0:1:5; w1 = 2*pi/10*k1;                            % [0,2pi] axis divided into 501 points.
%k = 0:500; w = (pi/500)*k;                             % [0,pi] axis divided into 501 points.
%X_DTFT = x * (exp(-j*pi/500)) .^ (n‘*k);               % DTFT of x(n)

%magX_DTFT = abs(X_DTFT); angX_DTFT = angle(X_DTFT); realX_DTFT = real(X_DTFT); imagX_DTFT = imag(X_DTFT);
subplot(3,1,2); stem(w1/pi,magY1_DFT); title(‘Samples of DTFT Magnitude, Not Enough‘); %axis([0,N,-0.5,1.5]);
xlabel(‘frequency in \pi units‘);
%ylabel(‘x(n)‘);
grid on;
subplot(3,1,3); stem(w1/pi,phaY1_DFT); title(‘Samples of DTFT Phase, Not Enough‘); %axis([0,N,-0.5,1.5]);
xlabel(‘frequency in \pi units‘);
%ylabel(‘x(n)‘);
grid on;

%% ------------------------------------------------------
%%         zero-padding coperation,  Append  90 zeros
%%           To obtain a dense spectrum
%% ------------------------------------------------------
n2 = [0:1:99]; y2 = [x(1:1:10) zeros(1,90)];           % zero-padding, N = 100

figure(‘NumberTitle‘, ‘off‘, ‘Name‘, ‘Exameple5.8 x sequence‘)
set(gcf,‘Color‘,‘white‘);
subplot(2,1,1); stem(n2,y2); title(‘signal x(n), 0<= n <=9 + 90 zeros, N = 100‘); axis([0,100,-2.5,2.5]);
xlabel(‘n‘); ylabel(‘x(n)‘); grid on;

Y2_DFT = dft(y2,100);                                   % DFT of y2

magY2_DFT = abs(Y2_DFT(1:1:51));
phaY2_DFT = angle(Y2_DFT)*180/pi                        % degrees

k2 = 0:1:50; w2 = (2*pi/100)*k2;                        % [0,2pi] axis divided into 501 points.
subplot(2,1,2); stem(w2/pi,magY2_DFT); hold on; plot(w2/pi,magY2_DFT,‘--r‘);
title(‘DTFT Magnitude‘); axis([0,1,0,10.5]); hold off;
xlabel(‘frequency in \pi units‘);
%ylabel(‘x(n)‘);
grid on;

%% -----------------------------------------------------------
%%          Exameple5.8b   take first 100 samples of x(n)
%%         determine the DTFT
%% -----------------------------------------------------------
figure(‘NumberTitle‘, ‘off‘, ‘Name‘, ‘Exameple5.8b ‘)
set(gcf,‘Color‘,‘white‘);
subplot(2,1,1); stem(n,x); axis([0,100,-2.5,2.5]);
title(‘signal x(n), 0 <= n <= 99, N = 100‘); xlabel(‘n‘); grid on;

X = dft(x,100); magX = abs(X(1:1:51));
k = 0:1:50; w = 2*pi/100*k;
subplot(2,1,2); plot(w/pi, magX); title(‘DTFT Magnitude‘);
xlabel(‘frequency in \pi units‘); axis([0,1,0,60]); grid on;
%% -----------------------------------------------------------
%%          END  Exameple5.8b
%% -----------------------------------------------------------

　　结果：

时间： 2024-11-05 03:14:25

《DSP using MATLAB》示例Example5.8的相关文章

DSP using MATLAB 示例Example3.21

代码: % Discrete-time Signal x1(n) % Ts = 0.0002; n = -25:1:25; nTs = n*Ts; Fs = 1/Ts; x = exp(-1000*abs(nTs)); Ts = 0.001; n = -5:1:5; nTs = n*Ts; Fs = 1/Ts; x = exp(-1000*abs(nTs)); % Analog Signal Dt = 0.00005; t = -0.005:Dt:0.005; xa = x * sinc(Fs*

DSP using MATLAB示例Example3.18

代码: % Analog Signal Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t)); % Continuous-time Fourier Transform Wmax = 2*pi*2000; K = 500; k = 0:1:K; % index array k for frequencies W = k*Wmax/K; % freqency between 0 and +pi, [0,pi] axis divided i

DSP using MATLAB 示例 Example3.19

代码: % Analog Signal Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t)); % Discrete-time Signal Ts = 0.0002; n = -25:1:25; x = exp(-1000*abs(n*Ts)); % Discrete-time Fourier Transform %Wmax = 2*pi*2000; K = 500; k = 0:1:K; w = pi*k/K; % index ar

DSP using MATLAB 示例 Example3.15

上代码: subplot(1,1,1); b = 1; a = [1, -0.8]; n = [0:100]; x = cos(0.05*pi*n); y = filter(b,a,x); figure('NumberTitle', 'off', 'Name', 'Input and Output sequence'); set(gcf,'Color','white'); subplot(2,1,1); stem(n,x); title('Input sequence'); xlabel('n'

DSP using MATlAB 示例Example2.10

上代码 % noise sequence 1 x = [3, 11, 7, 0, -1, 4, 2]; nx = [-3:3]; % given signal x(n) [y,ny] = sigshift(x,nx,2); % obtain x(n-2) set(gcf,'Color','white') subplot(2,1,1);stem(nx,x); title('sequence x(n)'); xlabel('n');ylabel('x(n)'); grid on subplot(2,

DSP using MATLAB 示例Example3.17

DSP using MATLAB示例Example3.6

代码: n = [-5:5]; x = (-0.9).^n; % x(n) = k = -200:200; w = (pi/100)*k; % [0,pi] axis divided into 101 points. X = x * (exp(-j*pi/100)) .^ (n'*k); magX = abs(X); angX = angle(X); realX = real(X); imagX = imag(X); set(gcf,'Color','white'); subplot(2,2,1

DSP using MATLAB 示例Example2.12

代码: b = [1]; a = [1, -0.9]; n = [-5:50]; h = impz(b,a,n); set(gcf,'Color','white'); %subplot(2,1,1); stem(n,h); title('Impulse Response'); xlabel('n'); ylabel('h(n)'); grid on; x = stepseq(0, -5, 50) - stepseq(10, -5, 50) y = filter(b, a, x); figure;

DSP using MATLAB 示例Example2.4

n = [0:10]; x = stepseq(0,0,10) - stepseq(10,0,10); [xe,xo,m] = evenodd(x,n); set(gcf,'Color',[1,1,1]) % 改变坐标外围背景颜色 stem(n,x); title('Rectangular Pulse'); xlabel('n'); ylabel('x(n)') ; axis([-10,10,0,1.2]) grid on figure set(gcf,'Color',[1,1,1]) stem

DSP using MATLAB 示例Example3.9

用到的性质 上代码: n = 0:100; x = cos(pi*n/2); k = -100:100; w = (pi/100)*k; % freqency between -pi and +pi , [0,pi] axis divided into 101 points. X = x * (exp(-j*pi/100)) .^ (n'*k); % DTFT of x % signal multiplied y = exp(j*pi*n/4) .* x; % signal multiplied