$clc$;

clear all;

close all;

vid$90=$ VideoReader$("$SAE$\_90\_$teja.MOV"$)$;

no$\_$of$\_$frames $=$ vid$90.$NumFrames;

time $=$ vid$90.$Duration;

frame$\_$rate $=$ vid$90.$FrameRate;

for $i=1$:$1$ :no$\_$of$\_$frames

filename $=$ strcat$("$frame$90-",$ num$2$str$($i$),\text{'}.$jpg$\text{'})$;

$a=$read$(vid90,i)$;

imwrite$($a$,$filename$)$;

$a=$ imrotate $(a,1\mathrm{.}13)$;

$g=$imcrop$(rgb2$gray$(a),[580,0,105,1600])$;

imwrite$(g,$filename$)$:

end

base$\_$img $=$ imread$(\text{'}$frame$90-13.$jpg$\text{'})$;

base$\_$img $=$ histeq$($base$\_$img$)$;

threshold $=16$;

binaryMask $=$ base$\_$img $>=$ threshold;

base$\_$img $=$ base$\_$img $.*$ uint$8($binaryMask$)$;

bw$\_$img $=$ imbinarize$($base$\_$img, "global"$)$;

in$\_$frame $=320$;

fin$\_$frame $=390$;

$k=$ in$\_$frame$-1$;

for $j=$ in$\_$frame:$1$:fin$\_$frame

img$\_$read $=$ imread$([\text{'}$frame$90-\text{'},$ num$2$str$($j$),\text{'}.$jpg $\text{'}])$;

img$\_$read $=$ histeq $($img$\_$read$)$;

threshold $=16$;

binaryMask $=$ img$\_$read $$ threshold;

img$\_$read $=$ img$\_$read.$**$ uint$8($binaryMask$)$;

bw$\_$img$1=$ imbinarize$($img$\_$read, "global"$)$;

diff$\_$img $=$ bw$\_$img$1-$ bw$\_$img;

diff$\_$img $=$ diff$\_$img;

$SE=$ strel$(\text{'}$disk$\text{'}\{$:${?}^{-},5)$;

dil$\_$img $=$ imdilate$($diff$\_$img,SE$)$;

subplot $(1,2,1)$;

imshow$($dil$\_$img$)$;

Image $=$ dil$\_$img;

hold on

axis on

B$1=$ sum$($dil$\_$img, $2)$;

$[$minValue$,$ minIndex$]=$min$(B1)$;

$y_{-}$loc $(j-k,1)=j$;

$y_{-}$loc $(j-k,2)=$ minIndex $**(\frac{0\mathrm{.}612}{1000})$;

$t(j-k,1)=\frac{j}{30}$;

pause$(0\mathrm{.}05)$

end

pos$\_$time $=[t,y_{l}oc($:$,2)]$;

degree $=1$;

coefficients $=$ polyfit$($pos$\_$time$((1$:$31),1),$ pos$\_$time$((1$:$31),2),$

xFit $=$ linspace $(10\mathrm{.}6667,11\mathrm{.}667,100)$;

yFit $=$ polyval$($coefficients$,$ xFit;

figure$(2)$

subplot $(2,1,1)$

plot$($pos$\_$time $($ : $,1),$ pos$\_$time$($ : $,2))$;

hold on;

plot$($xFit$,$ yFit, $\text{'}$r$-\text{'},$ 'LineWidth', $2)$;

hold off;

xlabel$(\text{'}$Time$($s$)\text{'})$;

ylabel$(\text{'}$Location of ball'$)$;

title$(\text{'}$Curve fitting for Location$-$Time'$)$;

legend$(\text{'}$Original Points', 'Fitted Curve'$)$;

slopeCoefficients $=$ polyder$($coefficients$)$;

slopeCurve $=$ polyval$($slopeCoefficients$,$ xFit$)$;

subplot $(2,1,2)$;

plot$($xFit$,$ slopeCurve, $\text{'}$g$-\text{'},$ 'LineWidth', $2)$;

xlabel$(\text{'}$Time$\text{'})$;

ylabel$(\text{'}$Velocity$\text{'})$;

title$(\text{'}$Velocity Time Curve'$)$;

velocity $=$ polyval$($slopeCoefficients$,21)$;

$\%$ Viscocity Calculation

dia $=6\mathrm{.}3e-3$;

rho$\_1=870$;

rho$\_$b $=8000$;

$g=9\mathrm{.}81$;

viscocity rho$\_$b $-$ rho$\_1\{$:$)**(g)**di{a}^{{?}^{{?}^{?}}}2$ velocity

I want this programme in matlab with some changes