Report 1 and Assignment 1 - 0611866
Report 1 and Assignment 1 - 0611866
KULLIYYAH OF ENGINEERING
DEPARTMENT OF BIOCHEMICAL-BIOTECHNOLOGY ENGINEERING
BTE 4225
COMPUTER SIMULATION IN BIOLOGICAL SYSTEM
MATLAB REPORT
AND
ASSIGNMENT ONE
TUTORIAL 1
Section 1.1: The command Window
>>lookfor sv
This command
is used if we do
not remember
the name of the
function.
The helpwin
command,
invoked
without
arguments,
opens a new
window on the
screen.
Other
example of
using helpwin
command
Here, we can
learn more
about
MATLAB by
opening the
demo
window.
The m-file can also be used as function files. We can directly run the function from the editor window or call
the function from command window.
But then, I found difficulties to calculate the value of y-axis and call other function in other m-file.
Therefore, I stop using the coding as my assignment and I try another way for the assignment. I change the
value of the volume for easier calculation.
ASSIGNMENT 1
Subject: Separation Process for Biochemical Product
Batch Filtration with Change in Pressure
In a filtration of suspended cells under full vacuum using a 0.2 m Buchner Funnel, you collect data for the
volume of filtrate as a function of time. Based on data and equations given, make a graph of batch filtration
operation.
Time (s)
Volume filtered(m3)
0
0
1200
2
1800
2.5
2400
3.25
3000
3.5
3600
3.75
4200
3.9
T = oc (V/A) + oRm
(V/A)
2p
p
The coding:
%ASSIGNMENT 1
%BATCH FILTRATION WITH CHANGE IN PRESSURE
function filtration
t = [0 1200 1800 2400 3000 3600 4200]; % time in minutes
v = [0 2 2.5 3.25 3.5 3.75 3.9 ]; % volume filtered in liters
d = 0.2; % diameter of filter
Area = pi*((d/2)^2)% formula of calculating area of a circle
x=v/Area; % calculating the value of x-axis
y=t./x; % calculating the value of y-axis
Data_x_y_axis = [x',y']
% plotting graph
f=plot (x,y)
grid
% manipulating the graph
set(f,'LineWidth',3, 'Color',[ 1 0.498 1])
xlabel('V/A')
h = get(gca,'xlabel');
set(h,'FontSize',14, 'Color',[ 0.4 0.8 1])
ylabel('t/(V/A)')
h = get(gca,'ylabel');
set(h,'FontSize',14, 'Color',[ 1 0.498 0])
title('Graph of t/(V/A) versus (V/A) of batch filtration')
h = get(gca,'Title');
set(h,'FontSize',14,'Color',[ 0 0.498 0])