Lab_data.xlsx contains the experimental drag force data for three different airfoils - flat plate, symmetric, and cambered. The drag coefficient (Cd) can be calculated from the following equation:
Cd = 2Fp / (p * v^2 * S)
Here, Fp (N) is the drag force, p (kg/m^3) is the air density, v (m/s) is the air velocity, and S (m^2) is the reference area. Assume that p = 1.205 kg/m^3, v = 25 m/s, and S = 0.02011 m^2.
(a) Import the data from the "Lab_data.xlsx" in a MATLAB script. Assign different variables for each airfoil.
(b) Create a function named "drag_calc" that accepts Fp, p, v, and S as input arguments and computes the Cd value. An example of the function header can be:
function Cd = drag_calc(Fp, p, v, S)
(c) Call the function in your MATLAB script to calculate the coefficient (Cd) for all airfoils. Plot the Cd values for all airfoils against the angle of attack (a) in the same figure. Use different color symbols for each airfoil to identify the trend of the drag coefficients. Include a figure title, axis labels, and appropriate legends.
(d) Calculate the minimum drag force and the corresponding angle of attack for all airfoils using MATLAB's built-in min() function.
(e) Print the minimum drag coefficient Cd for each airfoil and the corresponding angle of attack (a) values using fprintf().