Question

O tittps/ualearn blackboard.com/bbcswebdav/pid plot a new function on the same graph overwriting the orizinal function 1 of11ρ Functions to choose from: y A sin(Bx C)+ D y A sin2(Bx+ C) + D y A In(Bx + C)D y A exp(Bx + C)+ D ž. The file "Test Data.xlsx" contains in-cylinder pressure data from a firing engine. The data is resolved in crank angle degrees (CAD) and is taken every 1/2 degree. A full cycle of data consists of 1440 points. The file can be converted into a 2D array in MATLAB by using the xlsread function. With the file in your working directory, the code to read the data should look something like "data - xlsread('Test Data'),". This will give you a 2D array with CAD in the first column and pressure in the second column. For this problem, write a script file that a, Plots all cycles on a single graph {x = CAD, y-Pressure) b. Finds and reports the maximum and minimum pressure for each cycle c. Finds and plots the average of all the cycles d. Plots the cycle with the maximum peak pressure and the cycle with the minimum peak pressure (largest and smallest max pressure) on the graph with the average of all the cycles L. Make the maximum and minimum cycles dotted lines Hints: a. b. Convert the pressure column into a matrix with each cycle being a new column Remember that you can reference all columns lA(1.)-all columns in row 11 and all rows TA(:,1) all rows in column 2] using the array indexes ME349 Page I O Type here to search

problem2, the file "Test_Data.xlsx" contains in cylinder pressure data from a firing engine. the data is resolved in crank angle degree and is taken every 1/2 degree.

Answer 1

clc
clear
close all

%LOADING DATASET FROM MS EXCEL - Please note: Data for only 3 cycles was used

data = xlsread('Test_Data.xlsx');

CAD = data(:,1);
PRESSURE = data(:,2);

numer_of_data_points_per_cycle = 1440;
number_of_cycles = length(CAD)/numer_of_data_points_per_cycle;
cycles = 1:number_of_cycles;
anges_per_cycle = CAD(1:1440);

for i=1:number_of_cycles
p(i,:) = PRESSURE((i-1)*numer_of_data_points_per_cycle+1 :numer_of_data_points_per_cycle*i);
end

for i=1:number_of_cycles
plot(anges_per_cycle,p(i,:),'color',rand(1,3))
legInfo{i} = sprintf('Cycle #%d', i);
hold on
end
legend(legInfo)
xlabel('CAD')
ylabel('Pressure')
xlim([min(anges_per_cycle) max(anges_per_cycle)])
title('Pressure vs CAD for all cycles')
grid on

for i=1:number_of_cycles
minimum_pressure_for_each_cycle(i) = min(p(i,:));
end
minimum_pressure_for_each_cycle

for i=1:number_of_cycles
maximum_pressure_for_each_cycle(i) = max(p(i,:));
end
maximum_pressure_for_each_cycle

for i=1:number_of_cycles
avg_pressure_each_cycle(i) =mean(p(i,:));
end
avg_pressure_each_cycle
avg_pressure_all_cycles = mean(avg_pressure_each_cycle)

[psort,pos] = sort(maximum_pressure_for_each_cycle,'descend');

cycle_with_max_peak_pressure = pos(1);
cycle_with_min_peak_pressure = pos(end);

figure
plot(anges_per_cycle,p(cycle_with_max_peak_pressure,:),'-- m')
hold on
plot(anges_per_cycle,p(cycle_with_min_peak_pressure,:),'-. c')
hold on
plot([anges_per_cycle(1) anges_per_cycle(end)],[avg_pressure_all_cycles avg_pressure_all_cycles],'k','linewidth',2)
xlim([min(anges_per_cycle) max(anges_per_cycle)])
xlabel('CAD')
ylabel('Pressure')
title('Minimum,maximum peak pressure,overall average pressure')
grid on
legend(['cycle #' num2str(cycle_with_max_peak_pressure)],['cycle #' num2str(cycle_with_min_peak_pressure)],['overall average'])

minimum pressure for each cycle 0.0053 0.0167 0.0062 maximum_pressure_for_each_eycle - 1.0e+03 1.2612 2.8776 1.9240 avg pressure each cycle 248.1078 544.9359 367.8497 avg pressure all cycles 386.9645

Pressure vs CAD for all cycles 3000 Cycle #1 Cycle #2 2500 ___ Cycle #3)- 2000 1500 1000 500 100 200300 400 500 600 700 CAD

Minimum,maximum peak pressure,overall average pressure 3000 Cycle #2 Cycle #1 2500 _overall average 2000 1500 1000 500 0 100 200300 400 500 600 700 CAD

Excel Data Arrangement Used: L6 1 angle pressure 0.5 1.280199909 1 8.006156603 1.5 1.205456398 16.03408284 10.08110479 16.86059957 A030AAK, 2.5 cycle #1 starts at A21 438 718.5 1439 440 719.5 84.72942358 1441 69.78738973 33.02179181 719 6.385776 cycle #1 ends at A1441 1442 0.5 8.453493964 1443 1444 1445 1446 1447 1448 1449 1450 19.16673891 18.98593828 33.30155667 38.71665677 14.22613089 70.27288603 24.38158702 117.5949979 31 ORAAQAsa 1.5 2.5 3.5 4.5 cycle #2 starts at 1442

COMMENT DOWN FOR ANY QUERY RELATED TO THIS ANSWER
PLEASE PLEASE GIVE A THUMBS UP