<< Chapter < Page | Chapter >> Page > |
Let the function y defined by y=\cos x . Plot this function over the interval [-pi,pi]. Use numerical integration techniques to estimate the integral of y over [0, pi/2]using increments of pi/10 and pi/1000.
x=-pi:pi/100:pi;
y=cos(x);plot(x,y),title('Graph of y=cos(x)'),xlabel('x'),ylabel('y'),grid
>>x=0:pi/10:pi/2;>>y=cos(x);>>area1=trapz(x,y)
area1 =0.9918
>>x=0:pi/1000:pi/2;>>y=cos(x);>>area2=trapz(x,y)
area2 =1.0000
Let the function y defined by y=0.04x^{2}-2.13x+32.58 . Plot this function over the interval [3,30]. Use numerical integration techniques to estimate the integral of y over [3,30].
>>x=3:.1:30;>>y=0.04*(x.^2)-2.13.*x+32.58;>>plot(x,y), title('Graph of ...
y=.04*(x^2)-2.13*x+32.58'),xlabel('x'),ylabel('y'),grid
>>area=trapz(x,y)
area =290.3868
A 2000-liter tank is full of lube oil. It is known that if lube oil is drained from the tank, the mass flow rate will decrease from the maximum when the tank level is at the highest. The following data were collected when the tank was drained.
Time [min] | Mass Flow [kg/min] |
---|---|
0 | 50.00 |
5 | 48.25 |
10 | 46.00 |
15 | 42.50 |
20 | 37.50 |
25 | 30.50 |
30 | 19.00 |
35 | 9.00 |
clc
t=linspace(0,35,8) % Data entry for time [min]m=[50 48.25 46 42.5 37.5 30.5 19 9] % Data entry for mass flow [kg/min]% Calculate time intervals
dt=[t(2)-t(1),t(3)-t(2),t(4)-t(3),...t(5)-t(4),t(6)-t(5),t(7)-t(6),t(8)-t(7)]
% Calculate mass outdm=[0.5*(m(2)+m(1)),0.5*(m(3)+m(2)),0.5*(m(4)+m(3)),0.5*(m(5)+...
m(4)),0.5*(m(6)+m(5)),0.5*(m(7)+m(6)),0.5*(m(8)+m(7))]% Calculate differential areas
da=dt.*dm;% Tabulate time and mass flow
[t',m']% Tabulate time intervals, mass out and differential areas
[dt',dm',da']% Calculate the amount of oil drained [kg] in 35 minutesOil_Drained=sum(da)
The output is:
ans =
0 50.00005.0000 48.2500
10.0000 46.000015.0000 42.5000
20.0000 37.500025.0000 30.5000
30.0000 19.000035.0000 9.0000
ans =5.0000 49.1250 245.6250
5.0000 47.1250 235.62505.0000 44.2500 221.2500
5.0000 40.0000 200.00005.0000 34.0000 170.0000
5.0000 24.7500 123.75005.0000 14.0000 70.0000
Oil_Drained =1.2663e+003
A gas is expanded in an engine cylinder, following the law PV ^{1.3} =c. The initial pressure is 2550 kPa and the final pressure is 210 kPa. If the volume at the end of expansion is 0.75 m ^{3} , compute the work done by the gas. Applied Heat for Engineers by W. Embleton and L Jackson, Thomas Reed Publications. © 1999, (p. 80)
clc
disp('A gas is expanded in an engine cylinder, following the law PV^1.3=c')disp('The initial pressure is 2550 kPa and the final pressure is 210 kPa.')
disp('If the volume at the end of expansion is 0.75 m3,')
disp('Compute the work done by the gas.')disp(' ') % Display blank line
n=1.3;P_i=2550; % Initial pressure
P_f=210; % Final pressureV_f=.75; % Final volume
V_i=(P_f*(V_f^n)/P_i)^(1/n); % Initial volumec=P_f*V_f^n;
v=V_i:.001:V_f; % Creating a row vector for volume, vp=c./(v.^n); % Computing pressure for volume
WorkDone=trapz(v,p) % Integrating p*dv
The output is:
A gas is expanded in an engine cylinder, following the law PV^1.3=c
The initial pressure is 2550 kPa and the final pressure is 210 kPa.If the volume at the end of expansion is 0.75 m3,
Compute the work done by the gas.WorkDone =409.0666
A force F acting on a body at a distance s from a fixed point is given by $F=3s+\frac{1}{s^{2}}$ . Write a script to compute the work done when the body moves from the position where s=1 to that where s=10. O. N. Mathematics: 2 by J. Dobinson, Penguin Library of Technology. © 1969, (p. 213)
clc
disp('A force F acting on a body at a distance s from a fixed point is given by')disp('F=3*s+(1/(s^2)) where s is the distance in meters')
disp('Compute the total work done in moving')disp('From the position where s=1 to that where s=10.')
disp(' ') % Display blank lines=1:.001:10; % Creating a row vector for distance, s
F=3.*s+(1./(s.^2)); % Computing Force for sWorkDone=trapz(s,F) % Integrating F*ds over 1 to 10 meters.
The output is:
A force F acting on a body at a distance s from a fixed point is given by
F=3*s+(1/(s^2)) where s is the distance in metersCompute the total work done in moving
From the position where s=1 to that where s=10.WorkDone =149.4000
The pressure p and volume v of a given mass of gas are connected by the relation
(p+\frac{a}{v^{2}})(v-b)=k where a, b and k are constants. Express p in terms of v, and write a script to compute the work done by the gas in expanding from an initial volume to a final volume.
O. N. Mathematics: 2 by J. Dobinson, Penguin Library of Technology. © 1969, (p. 212)
Test your solution with the following input:
a: 0.01
b: 0.001
The initial pressure [kPa]: 100
The initial volume [m3]: 1
The final volume [m3]: 2
clc % Clear screen
disp('This script computes the work done by')disp('The gas in expanding from volume v1 to v2')
disp(' ') % Display blank linea=input('Enter the constant a: ');
b=input('Enter the constant b: ');p_i=input('Enter the initial pressure [kPa]: ');v_i=input('Enter the initial volume [m3]: ');v_f=input('Enter the final volume [m3]: ');k=(p_i+(a/(v_i^2))*(v_i-b)); % Calculating constant k
v=v_i:.001:v_f; % Creating a row vector for volumep=(k./(v-b))-(a./(v.^2)); % Computing pressure for volume
WorkDone=trapz(v,p); % Integrating p*dvdisp(' ') % Display blank line
str = ['The work done by the gas in expanding from ', num2str(v_i),...' m3 to ' num2str(v_f), ' m3 is ', num2str(WorkDone), ' kW.'];disp(str);
The output is:
This script computes the work done by
The gas in expanding from volume v1 to v2Enter the constant a: .01Enter the constant b: .001
Enter the initial pressure [kPa]: 100
Enter the initial volume [m3]: 1
Enter the final volume [m3]: 2The work done by the gas in expanding from 1 m3 to 2 m3 is 69.3667 kW.
Notification Switch
Would you like to follow the 'A brief introduction to engineering computation with matlab' conversation and receive update notifications?