E7 – Prof. Alam - Fall 2019, UC Berkeley E7 Homework Assignment 10: Least Squares and Linear Regression The purpose of this lab is to improve your familiarity with solving least squares problems and...

Please complete HW10_E7_10 using the .m file for data and check with autograder HW10_ag.p file for 100%.


E7 – Prof. Alam - Fall 2019, UC Berkeley E7 Homework Assignment 10: Least Squares and Linear Regression The purpose of this lab is to improve your familiarity with solving least squares problems and performing linear regressions using Matlab. Note 1: A template will not be provided for this assignment. You still need to publish and submit your m-file, so neatly organize it using code cells and comments as was done in the previous templates. The autograder does not rely on the format of this file. You may run the autograder as many times as you like to check your work as you complete the assignment. Note 2: Remember to upload the following to the bCourses website: → All .m files including the functions that you are asked to create. Specifically: Problem1.m, compInput.m, simRobot.m, MyPolyRegressor.m, FindPolyfitOrder.m, PlotPolyregression.m, MMfit.m, and getBasisCombinations.m. Directions to upload your solutions through bCourses can be found here. The homework is due at 11:59am on Monday, November 25, 2019. NO LATE HOMEWORK WILL BE ACCEPTED. 1. One of the first steps to selecting a cell phone signal booster for your home or office, vehicle, or large building is to determine how strong the outside signal is. When most people talk about cell phone signal strength, they talk about “bars” in reference to the signal strength bar indicator on the phone. While bars are an easy way to talk about signal strength, it turns out that it’s not a very accurate way to signal test. While many phones show the signal strength on a 5-bar scale, some phones only have 4 bars while others have 8. Even among phones with the same number of bars, there’s no standardization, so the strength of a 4 bar signal on one phone can be very different than a 4 bar signal on another. Finally, bars are not very granular so for example, saying 3 bars of signal is not being very specific. When experts discuss cell phone signal strength, they measure the signal in decibels. Decibels are a logarithmic unit of measuring signal strength and are very precise making them ideal for performing a signal test of just how strong of a signal is that you’re currently receiving. Most phones have a setting we call Field Test Mode that can show you useful information about your phone, including the signal strength in decibels. You can find the Field Test Mode on your phone following the instructions here. Sarah has conducted a simple experiment to find out about the relation between the Received Signal Strength (RSS) on her phone versus distance into the Northbrae tunnel. The results of her measurements appear below: 1 of 12 https://guides.instructure.com/m/4212/l/54353-how-do-i-upload-a-file-to-my-assignment-submission https://www.ubersignal.com/field-test-mode#ios_always_on https://en.wikipedia.org/wiki/Northbrae_Tunnel E7 – Prof. Alam - Fall 2019, UC Berkeley Distance into the tunnel (steps) Cell phone received signal strength (dB) 0 -92 5 -93 10 -105 15 -121 20 -107 25 -105 30 -101 (a) Store the presented data for “distance into the tunnel” and “cell phone received signal strength” into row arrays called “Dist” and “RSS”, respectively. Now, plot the measured RSS data versus distance. Label the x and y axes with “Distance into the tunnel (steps)” and “Cell phone received signal strength (dB)”, respectively. Save the plot as “CellData.png” using the print command. (b) Use MATLAB’s built-in functions, polyfit() and polyval(), to add linear and quadratic least squares curves that fit to the data. Plot them on top of the real measured data points, and use legends to differentiate the measured data from its linear and quadratic fits. Save this plot as “CellFit.png” using the print command. (c) As it is apparent from her experimental results, Sarah could not walk through the tunnel beyond 30 steps. Using the quadratic curve-fit from part (b), estimate the cell phone’s received signal strength (RSS) at DistNew = [35, 40, 45] steps into the tunnel. Store your results in a new row array called “RSSEst”. Note: Save your code for this problem as “Problem1.m” and submit it along with all other function files that you are asked to create for other problems. 2. Suppose you work at a robotics company and one of the robots that they make moves back and forth on a frictionless track, moved by a motor that applies a force to the robot. The robot is described by its state, which is its instantaneous position p and velocity v. Your task is to determine the force required to move the robot from its current state to some desired state in a specified amount of time. The motor applies a force for a specified number of time steps. The length of each time step is the same and the force is constant during each time step. Let ∆ be the duration of one time step and nSteps be the number of time steps. (For example, if ∆ = 0.1 seconds and nSteps = 10, then a constant force could be applied to the robot from 0 to 0.1 seconds, a different constant force from 0.1 to 0.2 seconds, and so on.) Let p0 and v0 be the initial position and velocity of the robot and pk and vk be the position and velocity of the robot at the end of the time step k. Then, the position and velocity of the robot can be described by the difference equations given below. pk+1 = ∆2uk 2m + vk∆ + pk vk+1 = ∆uk m + vk (1) 2 of 12 E7 – Prof. Alam - Fall 2019, UC Berkeley where m is the mass of the robot and uk is the constant force applied to the robot during time step k. (Think about what these equations are saying in the most general terms: where the robot is next (pk+1, vk+1) is a function of where the robot is now (pk, vk) and the applied force uk.) These equations can be rewritten as a matrix equation as shown below.[ pk+1 vk+1 ] = [ 1 ∆ 0 1 ][ pk vk ] + [ ∆2 2m ∆ m ] uk (2) Make sure you agree that Equations (1) and (2) are equivalent before continuing. Now we will make Equation (2) more concise by defining the vector xk := [ pk vk ] , the ma- trix C := [ 1 ∆ 0 1 ] , and the vector D := [ ∆2 2m ∆ m ] . Then, Equation (2) can be written as xk+1 = Cxk +Duk (3) For k = 0 and k = 1 this would give x1 = Cx0 +Du0 x2 = Cx1 +Du1 Notice that we can plug the equation for x1 into the right hand side of the equation for x2. This gives x2 = C 2x0 + CDu0 +Du1 Notice that C2 = C · C is the result of multiplying the matrix C by itself, and it is not an element-by-element exponentiation. More precisely, C2 = [ 1 ∆ 0 1 ] · [ 1 ∆ 0 1 ] = [ 1 2∆ 0 1 ] This process can be repeated for each subsequent equation until finally you find an equation for xnSteps (the final state of the robot) that depends only on the initial state x0, and the force inputs uk at each step. For example if nSteps equals 10 then we would have x1 = Cx0 +Du0 x2 = Cx1 +Du1 = C 2x0 + CDu0 +Du1 x3 = Cx2 +Du2 = C 3x0 + C 2Du0 + CDu1 +Du2 ... x10 = Cx9 +Du9 = C 10x0 + C 9Du0 + C 8Du1 + . . .+ CDu8 +Du9 3 of 12 E7 – Prof. Alam - Fall 2019, UC Berkeley This final equation can be put into the Au = b form where A is a 2-by-nSteps array, b is a 2-by-1 array, and u is a nSteps-by-1 array of the unknown input forces: [ C9D C8D · · · CD D ]  u0 u1 ... u8 u9  = x10 − C10x0 (4) As seen in the lectures, solutions to Eq. (4) can be found by using the backslash operator or the pinv command in Matlab. Recall that the backslash operator gives a the solution with the minumum number of nonzero elements and the pinv command gives the minimum norm solution. (a) Write a function named compInput with the following function declaration line begin code 1 function [u1, u2] = compInput(m, delta, nSteps, pi, vi, pf, vf) end code where the input arguments are • m: the mass of the robot; • delta: the time stepsize to be used; • nSteps: the number of time steps; • pi: the initial position of the robot; • vi: the initial velocity of the robot; • pf: the final position of the robot; and, • vf: the final velocity of the robot. The output arguments are both nSteps-by-1 arrays of the input force at each time step to be applied to the robot that takes the robot from the initial to the final position in the specified time. More specifically, • u1 is the input force with the minimum number of nonzero elements (determined using the backslash operator); and, • u2 is the input force with the minimum euclidean norm (determined by using the pinv function). (b) Now that the input force to be applied to the robot has been computed, we want to sim- ulate the robot using Equation 2 to determine its position and velocity at each time step. This allows us to confirm that the input force does indeed move the robot to the desired final position and velocity. Write a function named simRobot with the following function declaration line begin code 1 function [p, v] = simRobot(m, delta, nSteps, pi, vi, u) end code 4 of 12 E7 – Prof. Alam - Fall 2019, UC Berkeley where the input arguments are • m: the mass of the robot; • delta: the time stepsize to be used; • nSteps: the number of time steps; • pi: the initial position of the robot; • vi: the initial velocity of the robot; and, • u: the input force to apply to the robot. The output arguments are • p: an (nSteps+1)-by-1 array of the robots position. The first element is the initial position and the remaining elements are the positions after each time step; and, • v: an (nSteps+1)-by-1 array of the robots velocity. The first element is the initial velocity and the remaining elements are the velocities after each time step. An example of these functions being used to determine the input force and then simulate the robot is given below. begin code 1 % Define input parameters 2 m = 2; 3 delta = 0.01; 4 pi = 1; 5 vi = 5; 6 pf = 1.8; 7 vf = 0; 8 nSteps = 10; 9 10 % Determine input force 11 [u1, u2] = compInput(m, delta, nSteps, pi, vi,