Scribd
Upload
58 views5 pages

467 Lab 8

This document summarizes the results of Lab #8 exercises on signal processing and eigendecomposition. It contains: 1) Eigendecomposition of three matrices A with eigenvalues 100.9098, -97.3713, 2.4615. 2) Verification that eigenvalues of matrices A and B are 0.5 and 3. 3) Plotting of two signals u1 and u2 with a shortened time vector, showing only 2 periodic peaks. 4) Plotting a signal x corrupted with noise y, and its amplitude spectrum using FFT with frequencies up to 5Hz.

Uploaded by

Jordon Hon
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
58 views5 pages

467 Lab 8

This document summarizes the results of Lab #8 exercises on signal processing and eigendecomposition. It contains: 1) Eigendecomposition of three matrices A with eigenvalues 100.9098, -97.3713, 2.4615. 2) Verification that eigenvalues of matrices A and B are 0.5 and 3. 3) Plotting of two signals u1 and u2 with a shortened time vector, showing only 2 periodic peaks. 4) Plotting a signal x corrupted with noise y, and its amplitude spectrum using FFT with frequencies up to 5Hz.

Uploaded by

Jordon Hon
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

MEC E 467 – Lab #8

Jordon Hon 1394829

Exercise 1
a) A=[-149 -50 -154; 537 180 546; -27 9 -25];
P=poly(A)

eig(A)=
100.9098
-97.3713
2.4615

b) Rewritten in standard form

A=[2 -1; -1 1];


B=[1 0; 0 2/3];

Eig(A,B)

Verified that Eigenvalues= 0.5 and 3

Exercise 2

1) 628 data points instead of 1000

Fs = 100; % Sampling frequency (Hz)


T = 1/Fs; % Sampling time (s)
L = 628; % Length of data sequence
t = (0:L-1)*T; % Initialise a time vector

u1 = 1*sin(2.2361*t);
u2 = 1.5*sin(2.2361*t);
plot(t,u1,t,u2)
The x-axis is shortened and only 2 full periodic peaks can be seen. The number changed is the
data sequence which is the length of the data calculated, shortening it means less of the data is
seen.
2)
3)

4)
x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
y=x+2*randn(size(t)); % Sinusoids plus noise

plot(y,t)
x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
y=x+2*randn(size(t)); % Sinusoids plus noise

plot(Fs*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('Time (milliseconds)')

NFFT = 2^nextpow2(L); %Next power of 2 from length

Y=fft(x,NFFT)/L;

f = Fs/2*linspace(0, 1, NFFT/2);
%Plot single-sided amplitude spectrum.
plot(f, 2*abs(Y(1:NFFT/2)))
axis([0 5 0 0.01])
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('Time (milliseconds)')

You might also like