CURVE FIT

clear allclose all clc % preparing data cpdata=load(\'data\'); temp=cpdata(:,1); cp=cpdata(:,2); % curve fit %cp=a*t+c coeffs=polyfit(temp,cp,1); predictedcp=polyval(coeffs,temp); %compare curve fit with original data plot(temp,cp,\'linewidth\',3); hold on plot(temp,predictedcp,\'linewidth\',3,\'color\',\'r\'); xlabel(\'temp…

% calculation of drag force %inputs   % coefficient of drag cd=0.1; %density in kg/m^3 rho=1.2; % frontal area in m^2 a=5; %velocity of cycle in m/s v=linspace(1,10,20); % fd is drag force %formula fd=0.5*rho*a*v^2*cd; plot(fd,v,\'line width\',2,\'color\',\'r\') xlabel(\'drag force\') ylabel(\'velocity\')   %…

clear allclose all clc % preparing data cpdata=load(\'data\'); temp=cpdata(:,1); cp=cpdata(:,2); % curve fit %cp=a*t+c coeffs=polyfit(temp,cp,1); predictedcp=polyval(coeffs,temp); %compare curve fit with original data plot(temp,cp,\'linewidth\',3); hold on plot(temp,predictedcp,\'linewidth\',3,\'color\',\'r\'); xlabel(\'temp…

function [dtheta_dt] = ode_func(t,theta,b,g,l,m) theta1=theta(1)theta2=theta(2)dtheta1_dt=theta2;dtheta2_dt=-(b/m)*theta2-(g/l)*sin(theta1);dtheta_dt=[dtheta1_dt;dtheta2_dt];end clear allclose allclc b=0.5; g=9.81; l=1; m=1; % initial condition theta_0=[0;3]; % time points tspan=linspace(0,20,500); % solve ODE [t,results]=ode45(@(t,theta)…

% calculation of drag force %inputs   % coefficient of drag cd=0.1; %density in kg/m^3 rho=1.2; % frontal area in m^2 a=5; %velocity of cycle in m/s v=linspace(1,10,20); % fd is drag force %formula fd=0.5*rho*a*v^2*cd; plot(fd,v,\'line width\',2,\'color\',\'r\') xlabel(\'drag force\') ylabel(\'velocity\')   %…

% calculation of drag force %inputs   % coefficient of drag cd=0.1; %density in kg/m^3 rho=1.2; % frontal area in m^2 a=5; %velocity of cycle in m/s v=20; % fd is drag force %formula fd=0.5*rho*a*v^2*cd; plot(fd,v,\'line width\',2,\'color\',\'r\') xlabel(\'drag force\') ylabel(\'velocity\')   % calculation of…

clear allclose all clc % preparing data cpdata=load(\'data\'); temp=cpdata(:,1); cp=cpdata(:,2); % curve fit %cp=a*t+c coeffs=polyfit(temp,cp,1); predictedcp=polyval(coeffs,temp); %compare curve fit with original data plot(temp,cp,\'linewidth\',3); hold on plot(temp,predictedcp,\'linewidth\',3,\'color\',\'r\'); xlabel(\'temp…

function [dtheta_dt] = ode_func(t,theta,b,g,l,m) theta1=theta(1)theta2=theta(2)dtheta1_dt=theta2;dtheta2_dt=-(b/m)*theta2-(g/l)*sin(theta1);dtheta_dt=[dtheta1_dt;dtheta2_dt];end clear allclose allclc b=0.5; g=9.81; l=1; m=1; % initial condition theta_0=[0;3]; % time points tspan=linspace(0,20,500); % solve ODE [t,results]=ode45(@(t,theta)…

function [v]= engine_kinematics(bore,stroke,conrod,cr,startcrank,endcrank) a=stroke/2;R=conrod/a; vs=(pi/4)*bore^2*stroke;vc=vs/(cr-1); theta=linspace(startcrank,endcrank,100)term1= 0.5*(cr-1);term2=R+1-cosd(theta);term3=(R^2-(sind(theta).^2)).^0.5; v=(1+term1*(term2-term3)).*vcend   clear allclose allclc % calculation…

% calculation of drag force %inputs   % coefficient of drag cd=0.1; %density in kg/m^3 rho=1.2; % frontal area in m^2 a=5; %velocity of cycle in m/s v=20; % fd is drag force %formula fd=0.5*rho*a*v^2*cd; plot(fd,v,\'line width\',2,\'color\',\'r\') xlabel(\'drag force\') ylabel(\'velocity\')   % calculation of…

clear allclose all clc % preparing data cpdata=load(\'data\'); temp=cpdata(:,1); cp=cpdata(:,2); % curve fit %cp=a*t+c coeffs=polyfit(temp,cp,1); predictedcp=polyval(coeffs,temp); %compare curve fit with original data plot(temp,cp,\'linewidth\',3); hold on plot(temp,predictedcp,\'linewidth\',3,\'color\',\'r\'); xlabel(\'temp…

function [dtheta_dt] = ode_func(t,theta,b,g,l,m) theta1=theta(1)theta2=theta(2)dtheta1_dt=theta2;dtheta2_dt=-(b/m)*theta2-(g/l)*sin(theta1);dtheta_dt=[dtheta1_dt;dtheta2_dt];end clear allclose allclc b=0.5; g=9.81; l=1; m=1; % initial condition theta_0=[0;3]; % time points tspan=linspace(0,20,500); % solve ODE [t,results]=ode45(@(t,theta)…

% calculation of drag force %inputs   % coefficient of drag cd=0.1; %density in kg/m^3 rho=1.2; % frontal area in m^2 a=5; %velocity of cycle in m/s v=20; % fd is drag force %formula fd=0.5*rho*a*v^2*cd; plot(fd,v,\'line width\',2,\'color\',\'r\') xlabel(\'drag force\') ylabel(\'velocity\')   % calculation of…

a=[1;2;3] is row vector%% is column vector b=[1 2 3] is column vector%% is row vector

function [v]= engine_kinematics(bore,stroke,conrod,cr,startcrank,endcrank) a=stroke/2;R=conrod/a; vs=(pi/4)*bore^2*stroke;vc=vs/(cr-1); theta=linspace(startcrank,endcrank,100)term1= 0.5*(cr-1);term2=R+1-cosd(theta);term3=(R^2-(sind(theta).^2)).^0.5; v=(1+term1*(term2-term3)).*vcend   clear allclose allclc % calculation…

% calculation of drag force %inputs   % coefficient of drag cd=0.1; %density in kg/m^3 rho=1.2; % frontal area in m^2 a=5; %velocity of cycle in m/s v=20; % fd is drag force %formula fd=0.5*rho*a*v^2*cd; plot(fd,v,\'line width\',2,\'color\',\'r\') xlabel(\'drag force\') ylabel(\'velocity\')   % calculation of…

a=[1;2;3] is row vector b=[1 2 3] is column vector

a=[1;2;3] is row vector b=[1 2 3] is column vector

ABSTRACT This project deals with the manufacturing of magnesium composites with three different composition of reinforcements such as Al₂O₃,Silica and Carbon nano tubes. Stir casting process has been employed for producing magnesium matrix composites.To understand and analyse the quality of Magnesium composites,experiments…