Chapter 8, Problem 27E

To determine

To write:

A function to read and write the data from the data file, calculate intensities and write another function to print all of the information in a nearly organized table, add a function to calculate the average intensity of the storms, add a function to the program to print all of the information given on the most intense storm.

Answer to Problem 27E

Solution:

The file is,

%MATLAB code to load the data and create a data structure file.

%structure file.

stormdata = [321 2.4 1.5]

save stormdata.dat stormdata -ascii

stormdata = [111 3.3 12.1]

save stormdata.dat stormdata -ascii

clear

load stormdata.dat

%end of function file

%The function file should be placed in the same folder

The function file is,

%MATLAB code to print a table between stormdatavector and intensity.

%function file.

function tableprint(stormdatavector, intensity)

for i = 1:length(stormdatavector)

fprintf('%20.1f\t%20.1f\t%20.1f|t%20f\n', stormdatavector(i).locationcode, stormdatavector(i).rainfallamount, stormdatavector(i).stormduration, intensity(i))

end

end

%end of function file

%The function file should be placed in the same folder.

The function file is,

%MATLAB code to calculate the average intensity.

%function file.

function avgintensity(intensity)

intensitysum = 0;

for i = 1:length(intensity)

intensitysum = intensitysum + intensity(i);

end

avgintensity = intensitysum/length(intensity);

fprintf('Average intensity of storm is %f\n', avgintensity)

end

%end of function file

%The function file should be placed in the same folder.

The function file is,

%MATLAB code to find the information about mostintensestorm.

%function file.

function mostintensestorm(stormvector, intensity)

for i = 1:length(intensity)

for j = i + 1:length(intensity)

if intensity(i)>intensity(j)

temp = intensity(i);

intensity(i) = intensity (j);

intensity(j) = temp;

end

end

end

fprintf('the information of most intense storm is:%f\n', intensity(length(intensity)))

tempindex = index(stormvector, intensity);

fprintf('the index of most intense storm is:%d\n', tempindex);

fprintf('%20s|t%20s\t%20s\t%20s\n','location code','rainfall amount','storm duration','intensities')

fprintf('%20.1f\t%20.1f\t%20.1f\t%20f\n', stormvector(tempindex).locationcode, stormvector(tempindex).rainfallamount, stormvector(tempindex).stormduration, intensity(length(intensity)))

end

function index = index(stormvector, intensity)

for i = 1:length(stormvector)

if intensity(length(intensity)) ==stormvector(i).rainfallamount/stormvector(i).stormduration

index = i;

break;

end

end

end

%end of function file

%The function file should be placed in the same folder.

The function file is,

%MATLAB code to find the output.

%function file.

fprintf('%20s\t%20s\t%20s\t%20s\n','location code','rainfall amount','storm duration','intensities')

for i = 1:size(stormdata, 1)

for j = 1:length(stormdata(i))

stormvector(i) = struct('locationcode', stormdata(i, j),'rainfallamount', stormdata(i, j+1),'stormduration', stormdata(i, j+2));

intensity(i) = stormvector(i).rainfallamount/stormvector(i).stormduration;

end

end

tableprint(stormvector, intensity)

avgintensity(intensity)

mostintensestorm(stormvector, intensity)

%end of function file

%The function file should be placed in the same folder.

Explanation of Solution

The output is,

The average intensity is,

$\text{intensities}=0.2727$

The location code is,

$\text{location}\text{code}=111.0$

The rainfall amount is,

$\text{rainfall}\text{amount}=3.3$

MATLAB Code:

clc

clear all

close all

%MATLAB code to load the data and create a data structure file.

%structure file.

stormdata = [321 2.4 1.5]

save stormdata.dat stormdata -ascii

stormdata = [111 3.3 12.1]

save stormdata.dat stormdata -ascii

clear

load stormdata.dat

%end of structure file

%The file should be placed in the same folder

output;

%The function file should be placed in the same folder

%MATLAB code to print a table between stormdatavector and intensity.

%function file.

function tableprint(stormdatavector, intensity)

for i = 1:length(stormdatavector)

fprintf('%20.1f\t%20.1f\t%20.1f|t%20f\n', stormdatavector(i).locationcode, stormdatavector(i).rainfallamount, stormdatavector(i).stormduration, intensity(i))

end

end

%end of function file

%The function file should be placed in the same folder.

%MATLAB code to calculate the average intensity.

%function file.

function avgintensity(intensity)

intensitysum = 0;

for i = 1:length(intensity)

intensitysum = intensitysum + intensity(i);

end

avgintensity = intensitysum/length(intensity);

fprintf('Average intensity of storm is %f\n', avgintensity)

end

%end of function file

%The function file should be placed in the same folder.

%MATLAB code to find the information about mostintensestorm.

%function file.

function mostintensestorm(stormvector, intensity)

for i = 1:length(intensity)

for j = i + 1:length(intensity)

if intensity(i)>intensity(j)

temp = intensity(i);

intensity(i) = intensity (j);

intensity(j) = temp;

end

end

end

fprintf('the information of most intense storm is:%f\n', intensity(length(intensity)))

tempindex = index(stormvector, intensity);

fprintf('the index of most intense storm is:%d\n', tempindex);

fprintf('%20s|t%20s\t%20s\t%20s\n','location code','rainfall amount','storm duration','intensities')

fprintf('%20.1f\t%20.1f\t%20.1f\t%20f\n', stormvector(tempindex).locationcode, stormvector(tempindex).rainfallamount, stormvector(tempindex).stormduration, intensity(length(intensity)))

end

function index = index(stormvector, intensity)

for i = 1:length(stormvector)

if intensity(length(intensity)) ==stormvector(i).rainfallamount/stormvector(i).stormduration

index = i;

break;

end

end

end

%end of function file

%The function file should be placed in the same folder.

%MATLAB code to find the output.

%function file.

fprintf('%20s\t%20s\t%20s\t%20s\n','location code','rainfall amount','storm duration','intensities')

for i = 1:size(stormdata, 1)

for j = 1:length(stormdata(i))

stormvector(i) = struct('locationcode', stormdata(i, j),'rainfallamount', stormdata(i, j+1),'stormduration', stormdata(i, j+2));

intensity(i) = stormvector(i).rainfallamount/stormvector(i).stormduration;

end

end

tableprint(stormvector, intensity)

avgintensity(intensity)

mostintensestorm(stormvector, intensity)

%end of function file

%The function file should be placed in the same folder.

Save the MATLAB script with name, chapter8_54793_8_27E.m and function files with names tableprint.m, avgintensity.m, mostintensestorm.m and output.m in the current folder. Execute the functions by typing the functions name at the command window to generate output.

Result:

The result is,

Therefore, the result is stated above.