<< Chapter < Page | Chapter >> Page > |
Homework, tests, and solutions from previous offerings of this course are off limits, under the honor code.
Form a study group of 3-4 members. With your group, discuss and synthesize the major themes of this week of lectures. Turn in a one page summary of yourdiscussion. You need turn in only one summary per group, but include the names of all group members. Please do notwrite up just a "table of contents."
Construct a WWW page (with your picture ) and email Mike Wakin (wakin@rice.edu) your name (as you want it to appear on theclass web page) and the URL. If you need assistance setting up your page or taking/scanning a picture (both are easy!),ask your classmates.
Follow this learning styles link (also found on the Elec 301 web page ) and learn about the basics of learning styles. Write a short summary of what you learned. Also,complete the "Index of learning styles" self-scoring test on the web and bring your results to class.
Make sure you know the material in Lathi , Chapter B, Sections 1-4, 6.1, 6.2, 7. Specifically, be sureto review topics such as:
Reacquaint yourself with complex numbers by going to the course applets web page and clicking on the Complex Numbers applet (may take a few seconds to load).
(a) Change the default add function to exponential (exp). Click on the complex plane to get a blue arrow, which isyour complex number . Click again anywhere on the complex plane to get a yellow arrow,which is equal to . Now drag the tip of the blue arrow along the unit circle on with (smaller circle). For which values of on the unit circle does also lie on the unit circle? Why?
(b) Experiment with the functions absolute (abs), real part (re), and imaginary part (im) and report your findings.
Reduce the following to the Cartesian form, . Do not use your calculator!
(a)
(b)
(c)
(d)
(e)
Find the roots of each of the following polynomials (show
your work). Use MATLAB to check your answer with the
roots
command and to plot the roots in the
complex plane. Mark the root locations with an 'o'. Putall of the roots on the same plot and identify the
corresponding polynomial (
,
,
(a)
(b)
(c)
(d)
(e)
(f)
is called an Nth Root of Unity .
(a) Why?
(b) Let . Draw in the complex plane.
(c) Let . Draw in the complex plane.
A pair of vectors and are called linearly independent if It is a fact that we can write any vector in as a weighted sum (or linear combination ) of any two linearly independent vectors, where the weights and are complex-valued.
(a) Write as a linear combination of and . That is, find and such that
(b) More generally, write as a linear combination of and . We will denote the answer for a given as and .
(c) Write the answer to (a) in matrix form,
(d) Repeat (b) and (c) for a general set of linearly independent vectors and .
A Julia set is obtained by characterizing points in the complex plane. Specifically,let with complex, and define Then for each in the complex plane, we say if the sequence does not tend to infinity. Notice that if , then each element of the sequence also belongs to .
For most values of , the boundary of a Julia set is a fractal curve - it contains"jagged" detail no matter how far you zoom in on it. The well-known Mandelbrot set contains all values of for which the corresponding Julia set is connected.
(a) Let . Is in ?
(b) Let . What conditions on ensure that belongs to ?
(c) Create an approximate picture of a Julia set in MATLAB.
The easiest way is to create a matrix of complexnumbers, decide for each number whether it belongs to
, and plot the results
using the
imagesc
command. To
determine whether a number belongs to
, it is helpful to define
a limit
on the number of
iterations of
. For a
given
, if the magnitude
remains below some threshold
for all
, we say that
belongs to
. The code
below will help you get started:
N = 100; % Max # of iterations
M = 2; % Magnitude threshold
mu = -0.75; % Julia parameter
realVals = [-1.6:0.01:1.6];
imagVals = [-1.2:0.01:1.2];
xVals = ones(length(imagVals),1) * realVals + ...
j*imagVals'*ones(1,length(realVals));
Jmap = ones(size(xVals));
g = xVals; % Start with g0
% Insert code here to fill in elements of Jmap. Leave a '1'
% in locations where x belongs to J, insert '0' in the
% locations otherwise. It is not necessary to store all 100
% iterations of g!
imagesc(realVals, imagVals, Jmap);
colormap gray;
xlabel('Re(x)');
ylabel('Imag(x)');
This creates the following picture for , , and .
Using the same values for , , and , create a picture of the Julia set for . Print out this picture and hand it in with yourMATLAB code.
imagesc(log(Jmap))
and
colormap jet
for a neat picture.Notification Switch
Would you like to follow the 'Señales y sistemas' conversation and receive update notifications?