CSE-CE-ISyE524 Homework 5- Least Squares Solution

• Submit a single pdf image of your Jupyter notebook via Canvas.
• Please denote the start of each question in your Julia workbook using a LARGEFONT.
• Please answer “written” questions in Markdown boxes on your Julia workbook, not by using comments in Code boxes.
• A Julia notebook containing “starter” code for the questions below is supplied.
1. Voltage smoothing. We would like to send a sequence of voltage inputs to the manipulator arm of a robot. The desired signal is shown in the plot below (also available in voltages.csv).

Unfortunately, abrupt changes in voltage cause undue wear and tear on the motors over time, so we would like to create a new signal that is similar to the one above but with smoother transitions. If the voltages above are given by v1,v2,...,v200, one way to characterize smoothness is via the sum of squared differences:
R(v) = (v2 − v1)2 + (v3 − v2)2 + ··· + (v200 − v199)2
When R(v) is smaller, the voltage is smoother. Solve a regularized least squares problem that explores the tradeoff between matching the desired signal perfectly and making the signal smooth. Include a plot comparing the original signal to a few different smoothed versions obtained using regularized least squares and with varying degrees of smoothness.
2. Spline fitting. We are running a series of experiments to evaluate the properties of a new fluorescent material. As we vary the intensity of the incident light, the material should fluoresce different amounts. Unfortunately, the material isn’t perfectly uniform and our method for measuring fluorescence is not very accurate. After testing 200 different intensities, we obtained the result below (also available in xy_data.csv). The intensities xi and fluorescences yi are recorded in the first and second columns of the data matrix, respectively.

1 of 2
The material has interesting nonlinear properties, and we would like to characterize the relationship between intensity and fluorescence by using an approximate model that agrees well with the trend of our experimental data. Although there is noise in the data, we know from physics that the fluorescence must be zero when the intensity is zero. This fact must be reflected in all of our models!
a) Polynomial fit. Find the best cubic polynomial fit to the data. In other words, look for a function of the form y = a1x3 + a2x2 + a3x + a4 that has the best possible agreement with the data. Remember that the model should have (exactly) zero fluorescence when the intensity is zero! Include a plot of the data along with your best-fit cubic on the same axes.
b) Spline fit. Instead of using a single cubic polynomial, we will look for a fit to the data using two quadratic polynomials. Specifically, we want to find coefficients pi and qi so that our data is well modeled by the piecewise quadratic function:
if 0 ≤ x < 4 if 4 ≤ x < 10
These quadratic functions must be designed so that: • as in the cubic model, there is zero fluorescence when the intensity is zero.
• both quadratic pieces have the same value at x = 4.
• both quadratic pieces have the same slope at x = 4.
In other words, we are looking for a smooth piecewise quadratic. This is also known as a spline (this is just one type of spline, there are many other types). Include a plot of the data along with your best-fit model.
3. Hovercraft rendezvous. Alice and Bob are cruising on Lake Mendota in their hovercrafts. Each hovercraft has the following dynamics:
Dynamics of each hovercraft:
a) Find the sequence of thruster inputs for Alice (uA) and Bob (uB) that achieves a rendezvous at t = 60 while minimizing the total energy used by both hovercraft:
total energy =
Plot the trajectories of each hovercraft to verify that they do indeed rendezvous.
b) In addition to arriving at the same place at the same time, Alice and Bob should also make sure their velocity vectors match when they rendezvous (otherwise, they might crash!) Solve the rendezvous problem again with the additional velocity matching constraint and plot the resulting trajectories. Is the optimal rendezvous location different from the one found in the first part?
2 of 2