Question

MATLAB HOMEWORK QUESTION DUE TOMORROW: PLEASE HELP!!!

The viscosity of blood depends on many factors (temperature, hematocrit, blood vessel diameter, blood velocity). As such, there isn’t a single formula for calculating blood viscosity. Instead you can interpolate based on formulas and data found in tables.

Your job is to design a program to calculate blood viscosity. The key features of the program are:

• It must ask the user for temperature, hematocrit, blood vessel diameter, and blood velocity.
• For each of these property the program must be able to handle responses in different units (temperature in F, C or K, diameter and velocity in imperial or metric units)
• For each property the program must check that they are reasonable (blood freezes at ~-3C and boils at ~100 C, hematocrit is a percentage so it is ranges from 0 to 1, etc.)
• The program will load in 2 tables: viscosity vs blood vessel size, and viscosity vs velocity.
• The program will calculate the viscosity based on a formula that uses temperature and hematocrit (the Einstein equation applied to blood).
• The program will calculate viscosity by interpolating against the 2 tables and the calculated value.
• The program will present the viscosity with appropriate units (imperial or metric) based on what the user selects. Imperial units for viscosity are lb-s/ft². Metric unit are centipoise (10^-3 Pa-s).

The main program for your design should be called HW5_5. It should be modular (i.e. it calls secondary functions to do the work). These secondary functions are included in HW5_5 as local functions (i.e. in the same file, at the end of the main function). Write the main program and function stubs for all of the secondary functions except for the function that presents the final viscosity to the user. Each function stub should have help text explaining what the function is supposed to do.

Hint 1: the version of HW5_5 that you turn in should run, BUT it does not require any input from the user, data from tables, understanding of the Einstein equation, or how to convert units. The function stubs just need to return appropriate responses. See lecture 10 and/or https://en.wikipedia.org/wiki/Method_stub (a method is similar to a function. You don’t need to worry about the difference for now).

Hint 2: This should take you no more than 60 minutes (10-20 minutes thinking, 20-40 minutes writing). If you are going over this, you are probably misunderstanding the problem. Get help!

Answer 1

This part is of a Translational Biomedical Science Research initiative. It was supported materially by the Albury South West Pathology – a unit of Western Pathology Cluster of NSW Health Australia. 10 years of de-identified archived clinical pathology data for the period of January 1999 to December 2008 constitute the research database. In a 10-year period, the Pathology referred only ten cases for WBV determination, of which five cases that were selected for this evaluation. Selection criteria are as recently published.

The pair of haematocrit (HCT) and serum total protein (TP) values was used in various units to derive four units-converted-extrapolations of the same formulae for:

High shear rate: WBV (208 Sec-1) = (0.12 × HCT) + 0.17(TP - 2.07)

Low shear rate: WBV (0.5 Sec-1) = (1.89 × HCT) + 3.76(TP - 78.42)

• Extrapolation 1: According to the units employed in developing the chart, which converted the units for serum total protein only:

  Where HCT = haematocrit (%) and TP = Serum total proteins (g/L)

• Extrapolation 2: based on units (TP g/dL; HCT %) in referenced material [2].

• Extrapolation 3: based on ‘Prot g/dL’ and ‘HCT fraction’. That is, both units of HCT and TP from the referenced formula were converted.

• Extrapolation 4: WBV 208 Sec1 based on ‘Prot g/L’ and ‘HCT fraction’. That is, only the HCT has been modified.

The fractional contribution of haematocrit or serum proteins to WBV were worked out, in percentages, according to the following formulas

• Haematocrit (HCT %) contribution = (0.12 × HCT)/((0.12 × HCT) + 0.17(Prot – 2.07))

• Protein (TP %) contribution = 0.17(Prot – 2.07)/((0.12 × HCT) + 0.17(Prot – 2.07))