Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
0% found this document useful (0 votes)
9 views

Matlab

This document provides information about an assignment for a process stimulation lab conducted by Rushikesh R Mankar as part of a Materials and Systems Engineering Design course at the National Institute of Technology Warangal Department of Mechanical Engineering. The assignment was guided by Prof. Dr. Veera Venkata Nagaraju and involves solving problems using symbolic mathematics software.
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
9 views

Matlab

This document provides information about an assignment for a process stimulation lab conducted by Rushikesh R Mankar as part of a Materials and Systems Engineering Design course at the National Institute of Technology Warangal Department of Mechanical Engineering. The assignment was guided by Prof. Dr. Veera Venkata Nagaraju and involves solving problems using symbolic mathematics software.
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 32

National Institute of Technology Warangal

Department of Mechanical

Materials And Systems Engineering Design

PROCESS STIMULATION LAB


Assignment-5

Guided By
Prof. Dr. Veera Venkata Nagaraju

Rushikesh R Mankar(23MEM6R05)
Slide 12
Problem 6
syms x y z a
f1=2*x+3*y-4*z;
f2=x-2*y+5*z;
a=f2-f1
Solution-
a = 9*z - 5*y – x
%problem 7
syms a b c
f1=2*a+3*b;
f2=a+b;
c=f1*f2
Solution-
c = (a + b)*(2*a + 3*b)

%problem 8
syms x y
f1=3*x-2*y^2+4*x*y;
f2=2*x-5*y;
c=f1*f2
Solution-
c = (2*x - 5*y)*(- 2*y^2 + 4*x*y + 3*x)

%problem 1
syms x y z b a
b=(x^2*y^3*z)*(x^3*y*z^2)
a=simplify(b)
subs(a,[x y z],[0.5 2 3])
a = x^5*y^4*z^3
ans = 27/2

%problem 29
syms x y b a
factor(a*x-a*y+b*x-b*y)
ans =[x - y, a + b]

%problem 30
syms x y b a
factor(2*a*x-3*a*y+2*b*x-3*b*y)
ans = [2*x - 3*y, a + b]

%problem 31
syms x
factor(x^3+3*x^2-x-3)
ans = [x + 3, x - 1, x + 1]

%problem 3
syms x y
simplify((x^3+y^3)/(x+y))
ans = x^2 - x*y + y^2

Slide 13
%problem 1
syms x y a b
y=7*x^4;
a=diff(y,x,2);
b=diff(y,x,3);
a = 84*x^2
b = 168*x

%problem 2
syms x y a b
y=sqrt(x);
a=diff(y,x,2)
b=diff(y,x,3)

a = -1/(4*x^(3/2))

b = 3/(8*x^(5/2))

%problem 3
syms t y a b
y=sqrt(t^3);
a=diff(y,t,2)
b=diff(y,t,3)
a = (3*t)/(t^3)^(1/2) - (9*t^4)/(4*(t^3)^(3/2))
b = 3/(t^3)^(1/2) - (27*t^3)/(2*(t^3)^(3/2)) + (81*t^6)/(8*(t^3)^(5/2))

%problem 4
syms x y a b
y=6+(1/x^3)
a=diff(y,x,2)
b=diff(y,x,3)
%problem 5
syms x y a b
y=3*x-(1/sqrt(x))+(1/x)
a=diff(y,x,2)
b=diff(y,x,3)

%problem 6
syms x y a b
y=(5/x^2)-(1/sqrt(x^7))+2
a=diff(y,x,2)
b=diff(y,x,3)

%problem 7
syms t y a b
y=3*(t-2)^2
a=diff(y,t,2)
a=6

%problem 9
syms x y a b
y=(2/5)*x^3-(4/x^3)+4*sqrt(x^5)+7
a=diff(y,x)
a=

(10*x^4)/(x^5)^(1/2) + (6*x^2)/5 + 12/x^4

%problem 10
syms t y a b
y=5*t+(1/sqrt(t^3))
a=diff(y,t)

%problem 11
syms x y
y=((x+2)^2)/x
diff(y,x)

%problem 13
syms x y
y=7*sin(2*x)-3*cos(4*x)
diff(y,x)

Slide 14
%problem 1
syms x y a b
y=3*x^4+2*x^3-3*x+2
a=diff(y,x,2)
b=diff(y,x,3)

%problem 2
syms t y a
y=((2*t^2)/5)-(1/t^3)+(3/t)-sqrt(t)+1
a=diff(y,t,2)
subs(a,t,1)

%problem 3
syms t y1 y2 a b
y1=3*sin(2*t)+cos(t)
y2=2*log(4*t)
a=diff(y1,t,2)
b=diff(y2,t,2)

%problem 4
syms x y1 y2 a b
y1=2*cos(x)^2
y2=(2*x-3)^4
a=diff(y1,x,2)
b=diff(y2,x,2)

Slide 15
%problem 1
clear all
syms x y a b
y=(x^2-4*x)*(x^2-4*x+1)-20
a=expand(y)
b=y/(x^2-4*x+4)
c=factor(b)
d=solve(x^2-4*x+5)
e=subs(b,2-j)

Slide 16
%problem 3
syms x y z
exp1=x+3*y-z-2;
exp2=x-y+z-3;
exp3=3*x-5*y-4;
[x,y,z]=solve(exp1,exp2,exp3)
Solution-
x = 33/16
y = 7/16
z = 11/8

%problem 5(i)
syms x
y=exp(-x)
a=int(y,0,Inf)
%problem 5(ii)
syms x
y=exp(-x^2)
b=int(y,0,Inf)

%problem 5(iii)
syms x
y=exp(-x^2)*sin(x)
b=int(y,0,Inf)

a=1

b = -(pi^(1/2)*erf(1i/2)*exp(-1/4)*1i)/2

Slide 17
syms f(x)
syms x
f(x)=x^3+2*x^2-x+3
dfdt=diff(f,x,1)
Solution-
dfdt(x) = 3*x^2 + 4*x – 1

Slide 18
syms f(x)
syms x
f(x)=0.2+25*x-200*x^2+675*x^3-900*x^4+400*x^5
dfdt=diff(f,x,1)
Sol-
dfdt(x) = 2000*x^4 - 3600*x^3 + 2025*x^2 - 400*x + 25

%problem 1
syms x y t
y1=sin(3*x)+sin(x)
a=diff(y1,x)
%problem 2
y2=0.5*(sin(9*t)-sin(7*t))
b=diff(y2,t)
%problem 3
y3=cos(5*t)+cos(3*t)
c=diff(y3,t)
%problem 4
y4=(cos(5*t)-cos(t))/8
d=diff(y4,t)
Slide 22
%problem 1
clc
clear all
y1=[1 4 -32]
a=roots(y1)
%problem 2
y2=[1 0 -16]
b=roots(y2)

a =-8
4
b= 4
-4
%problem 3
syms x
solve((x+2)^2==16)
ans =-6
2
%problem 4
clc
clear all
y4=[2 -1 -3]
a=roots(y4)
%problem 5
y5=[6 -5 1]
b=roots(y5)
%problem 6
y6=[10 3 -4]
c=roots(y6)
%problem 7
y7=[1 -4 4]
d=roots(y7)
%problem 8
y8=[21 -25 -4]
e=roots(y8)
%problem 9
y9=[6 -5 -4]
f=roots(y9)
%problem 10
y10=[8 2 -15]
h=roots(y10)
Sol-
a = 1.5000
-1.0000
b = 0.5000
0.3333
c = -0.8000
0.5000

d =2
2

e = 1.3333
-0.1429

f = 1.3333
-0.5000

h = -1.5000
1.2500
%problem 2
syms x y
a=factor(x^2-6*x+9)
b=factor(4*x^2-25)
Sol-
a = [x - 3, x - 3]

b = [2*x - 5, 2*x + 5]
% problem 3
syms x
a=factor(4*x^2+8*x+3)
b=factor(15*x^2+2*x-8)
Sol-
a = [2*x + 3, 2*x + 1]

b = [5*x + 4, 3*x - 2]

SLIDE 23
clear all
syms x
a=solve(x^3-1==0)
b=[1 -1 -5 2]
c=roots(b)

syms x
a=[1 -2 -2 2];
b=round(roots(a),1)
c=[2 -1 -9.08 8.28];
d=round(roots(c),1)
e=[8 36 54 27];
f=round(roots(e),1)

Sol-
b = 2.5000
-1.2000
0.7000

d = -2.3000
1.7000
1.0000

f =-1.5000
-1.5000
-1.5000

SLIDE 26
clear all
syms x
p1=[1 8 5 4 3 2 1 1];
a=roots(p1)
p2=[1 -7 7 -10 -8 7 15];
b=roots(p2)
p3=[1 -8 10 12 8 -15];
c=roots(p3)
Sol-

clear all
syms x
p4=[1 7 12 -25 8];
a=roots(p4)
p5=[1 15 -23 105];
b=roots(p5)
p6=[1 -18 23];
c=roots(p6)
p7=[1 0 0 7];
d=roots(p7)
SLIDE 29

clear all
%problem 8(a)
syms x
a=solve(8*x+3==0)
% problem 8(b)
format long
b=round(solve(8*x+3==0),15)
%problem 8(c)
syms x p q
c=solve(x^3+p*x+q==0)
%problem 8(d)
d=solve(exp(x)-8*x+4)

SOL-
a = -3/8

b = -0.375

c = root(z^3 + p*z + q, z, 1)
root(z^3 + p*z + q, z, 2)
root(z^3 + p*z + q, z, 3)

d =1/2 - lambertw(0, -exp(1/2)/8)


%problem 9(a)
x=linspace(-4,4);
y=x.^3-x;
plot(x,y);
xlabel('x')
ylabel('x^3 -x')
grid on

%problem 9(b)
syms y x
y=sin(1/x^2);
ezplot(y,[-2,2])
grid on
xlabel('x')
ylabel('y')
title('i) using ezplot')
x=linspace(-2,2)
y=sin(1./x.^2);
plot(x,y)
grid on;
xlabel('x')
ylabel('y')
title('ii) Using plot')

%problem 9(c)
syms y x
y=tan(x/2);
ezplot(y,[-pi,pi,-10,10])
grid on
xlabel('x')
ylabel('y')

%problem 9(d)
x=linspace(-4,4);
y1=exp(-x.^2);
y2=x.^4-x.^2;
plot(y1,x,y2,x)
xlabel('x')
ylabel('y')
grid on

SLIDE 37

syms x t
f=4+(3/7)*x-6*x^2
I=int(f,x)
%problem 3(a)
f1=(2*x^3-3*x)/(4*x)
I1=int(f1,x)
%problem 3(b)
f2=(1-t)^2
I2=int(f2,t)

%problem 4
syms x t
f=(3/x^2)
I=int(f,x)
%problem 6
f1=(-5/9)*(t^(-3/4))
I1=int(f1,t)
%problem 8(a)
f2=4*cos(3*x)
I2=int(f2,x)
%problem 8(b)
f3=5*sin(2*t)
I2=int(f3,t)
%problem 10(a)
syms x t
f=5*exp(3*x)
I=int(f,x)
%problem 10(b)
f1=(2/3)*exp(-4*t)
I1=int(f1,t)

%problem 3(a)
syms x t
f=((3*x^2-5*x)/x)
I=int(f,x)
%problem 3(b)
f1=(2+t)^2
I1=int(f1,t)
Sol-
%problem 4(a)
syms x t
f=(4/3)*x^(-2)
I=int(f,x)

f = 4/(3*x^2)

I -4/(3*x)

SLIDE 38
clear all
%problem 19.2
y=inline('1-exp(-x)')
z=quad(y,0,4)
%problem 19.3
y1=@(x) 8+4*cos(x)
z1=quad(y1,0,pi/2)
%problem 19.4
y2=@(x) 1-x-4*x.^3+2*x.^5
z2=quad(y2,-2,4)
%problem 19.6
clear all
a=@(x,y) x.^2-3*y.^2+x*y.^3
I=dblquad(a,0,4,-2,2)
SOL-
a = function_handle with value:

@(x,y)x.^2-3*y.^2+x*y.^3

I = 21.3333
%problem 19.7
clear all
syms x y z
I=int(int(int(x^3-2*y*z,x,-1,3),0,6),-4,4)
Sol-
I = 960
%problem
clear all
y=@(x) x.*exp(2*x)
I=quad(y,0,3)
Sol-

y = function_handle with value:

@(x)x.*exp(2*x)

I = 504.5360

SLIDE 39

%problem 2
clear all
syms x y
I1=int((3-x^2),x,-1,2)
I2=int((x^2-4*x+3),x,1,3)

Sol-

I1,I2

I1 = 6

I2 = -4/3
%problem 3
clear all
syms x y1
y1=(3/2)*cos(x)
I1=int(y1,0,pi)
y2=4*cos(x)
I2=int(y2,0,pi/2)
Sol-

I1 = 0

I2 = 4

%problem 4
clear all
syms x y1
y1=2*sin(2*x)
I1=int(y1,pi/6,pi/3)
y2=@(t) 3.*sin(t)
I2=quad(y2,0,2)
Sol-
y1 = 2*sin(2*x)

I1 = 1

y2 =function_handle with value:

@(t)3.*sin(t)

I2 = 4.2484

%problem 5
clear all
syms x y1
y1=@(x) 5.*cos(3*x)
I1=quad(y1,0,1)
y2=@(t) 3*sec(2*t).^2
I2=quad(y2,0,pi/6)
Sol-
y1 = function_handle with value:

@(x)5.*cos(3*x)

I1 = 0.2352

y2 = function_handle with value:

@(t)3*sec(2*t).^2

I2 = 2.5981
%problem 5
clear all
syms x y1 t
y1=1/(sin(4*t)^2)
I1=int(y1,1,2)
y2=3*sin(2*x)-2*cos(3*x)
I2=int(y2,pi/4,pi/2)
Sol-
y1 = 1/sin(4*t)^2

I1 = Inf
y2 = 3*sin(2*x) - 2*cos(3*x)
I2 = 2^(1/2)/3 + 13/6

%problem 7
y1=@(t) 3*exp(3*t)
I1=quad(y1,0,1)
y2=@(x) (2/3)*exp(-2*x)
I2=quad(y2,-1,2)
y1 = function_handle with value:

@(t)3*exp(3*t)
I1 = 19.0855

y2 = function_handle with value:

@(x)(2/3)*exp(-2*x)

I2 = 2.4569
%problem 1
syms x y
I1=int(int(exp(x+y),0,log(y)),1,log(8));
round(I1,4)
ans = 3.3538

%problem 2
syms x y
I1=int(1/(1+x^2+y^2),0,sqrt(1+x^2))
int(I1,0,1)
clear all
%problem 3
syms x y a
I1=int(int(cos(x+y),0,pi/2),0,pi/2)
%problem 4
I2=int(int(sqrt(a^2-x^2-y^2),0,sqrt(a^2-y^2)),0,a)
Sol-

I1 = 0

I2 = (pi*a^3)/6

%problem 8
syms x y z
I1=int(int(int((x*y*z),0,y*z),1,z),0,2)
Sol-
I1 = 7/2
%problem 9
syms x y z
I1=int(int(int(cos(z/x),0,x*y),x,pi/2),0,pi/2)

SLIDE 48
clear all
%problem 1
ode1 = @(t, x, y, z) -10*x + 10*y;
ode2 = @(t, x, y, z) -y - x*z;
ode3 = @(t, x, y, z) (-8/3)*z + x*y - (224/3);
initial_conditions = [-8, 8, 27];
[t, solution] = ode45(@(t, u) [ode1(t, u(1), u(2), u(3)); ode2(t, u(1), u(2), u(3)); ode3(t, u(1),
u(2), u(3))], [0, 10], initial_conditions);

x_values = solution(:, 1);


y_values = solution(:, 2);
z_values = solution(:, 3);
% Plot z versus x
plot(x_values, z_values);
xlabel('x');
ylabel('z');
title('Plot of z vs x');
grid on;

SLIDE 49

%problem 2
clear all
syms y(x)
cond=y(3)==5
a=diff(y,x)==4*x-2
ysol(x)=dsolve(a,cond)

%problem 3
clear all
syms y(x)
a=diff(y,x)==2*x+3
y(x)=dsolve(a)

Sol-
a(x) = diff(y(x), x) == 2*x + 3

ans = C1 + x*(x + 3)

%problem 4
clear all
syms y(x) a
c=diff(y,x,2)+3*diff(y,x)+a*y==0
y(x)=dsolve(c)
y(x)=simplify(y(x))

%problem 5
clear all
syms y(x) a
a=diff(y,x)==(2*x+3*cos(y))/(3*x*sin(y)-2*y)
y(x)=dsolve(a,'Implicit',true)
Sol-
ans = 3*x*cos(y(x)) + y(x)^2 == - x^2 + C1

SLIDE 50
%problem 6
clear all
syms y(x)
a=diff(y,x,2)+8*diff(y,x)+2*y==cos(x)
Dy=diff(y,x)
cond=[y(0)==0,Dy(0)==1]
y(x)=dsolve(a,cond)
simplify(y(x))

ans = (65^(1/2)*sin(x + atan(1/8)))/65 - exp(-x*(14^(1/2) + 4))/130 - exp(x*(14^(1/2) -


4))/130 + (53*14^(1/2)*exp(x*(14^(1/2) - 4)))/1820 - (53*14^(1/2)*exp(-x*(14^(1/2) +
4)))/1820

%problem 7
clear all
syms y(x)
a=diff(y,x,2)+y==sin(x)
Dy=diff(y,x)
cond=[y(0)==1,Dy(0)==-2]
y(x)=dsolve(a,cond)
simplify(y(x))

%problem 8 i)
syms y(x)
a=diff(y,x)==7*x
y(x)=dsolve(a)
ans = (7*x^2)/2 + C1

syms y(x)
%problem 8 ii)
b=diff(y,x,2)+5*diff(y,x)+2*y==0
ysol(x)=dsolve(b)

%problem 9
syms y(x)
a=x*diff(y,x)==2-4*x^3
y(x)=dsolve(a)
Sol-
y(x) = C1 + 2*log(x) - (4*x^3)/3

SLIDE 51
%problem 10
syms y(x)
a=diff(y,x)==(x-1)^2
y(x)=dsolve(a)
Sol-

y(x) = C1 + (x*(x^2 - 3*x + 3))/3

%problem 11
syms y(x)
a=a=diff(y,x)==-(1/exp(3*x))
y(x)=dsolve(a)

%problem 12
syms y1(x) y2(x) y(x)
eqn1=diff(y,x)+2*x*y^2==0
y1(x)=dsolve(eqn1)
%problem 13
eqn2=diff(y,x)==exp(3*x+2*y)
y2(x)=dsolve(eqn2)

SLIDE 52
%problem 1
syms y(x)
a=diff(y,x)==cos(4*x)-2*x
y(x)=dsolve(a)
%problem 2
clear all
syms y(x)
b=2*x*diff(y,x)==3-x^3
y1(x)=dsolve(b)
sol-

%problem 3
syms y(x)
a=diff(y,x)+x==3
y(x)=dsolve(a,'y(1)=2')
%problem 4
clear all
syms y(x)
b=3*diff(y,x)+sin(x)==0
y1(x)=dsolve(b,'y(pi/3)=2/3')
sol-
%problem 5
syms y(x)
a=exp(-x)+2==x-3*diff(y,x)
y(x)=dsolve(a,'y(0)=1')
Sol-

y(x) = exp(-x)/3 - (2*x)/3 + x^2/6 + 2/3

%problem 4
clear all
syms y(x)
b=2*y*(1-x)+x*(1+y)*diff(y,x)==0
y1(x)=dsolve(b,'y(1)=1','Implicit',true)
Sol-
y1(x) =log(y(x)) + y(x) == 2*x - 2*log(x) – 1

%problem 5
clear all
syms y(x)
b=(y/x)*diff(y,x)==(y^2+1)/(x^2+1)
y1(x)=dsolve(b,'Implicit',true)
SOL-
log(y(x)^2 + 1) == C1 + log(x^2 + 1)
(y(x)^2 + 1)/y(x) == 0

%problem 6
clear all
syms y(x)
b=x*y==(1-x^2)*diff(y,x)
y1(x)=dsolve(b,'y(0)=1')
Sol-
y1(x) = 1i/(x^2 - 1)^(1/2)

SLIDE 59
% problem 1
clear all
c=0.517;
tau=3;
[t,y]=ode45(@(t,y) (c-y)/tau,[0 25],0.5)
plot(t,y)
grid on
xlabel('t')
ylabel('x')
% problem 1(b)
clear all
[t,y]=ode45(@(t,y) (y+3*t)/t^2,[1 4],-2)
plot(t,y)
grid on
xlabel('t')
ylabel('x')

% problem 1(c)
clear all
[t,y]=ode45(@(t,y) -y*sin(t),[-1 2],1)
plot(t,y)
grid on
axis('equal')
xlabel('t')
ylabel('y')
% problem 2
clear all
[t,y]=ode45(@(t,y) 2*y+t^(-2),[0.1 6],2.4)
plot(t,y)
grid on
xlabel('t')
ylabel('y')

SLIDE 60
% problem 3
clear all
[t,y]=ode45(@(t,y) (-t*y)/(sqrt(2-y^2)),[0 6],1)
plot(t,y)
grid on
xlabel('t')
ylabel('y')
Problem 4
function solve_ode()
tspan = [0, pi/2];
y0 = [1; -1; 0];
[t, y] = ode45(@odefun, tspan, y0);
plot(t, y)
grid on
xlabel('t')
ylabel('y')
legend('y1', 'y2', 'y3')
end

function dydt = odefun(t, y)


dydt = zeros(3, 1);
dydt(1) = 2 * y(1) + y(2) + 5 * y(3) + exp(-2 * t);
dydt(2) = -3 * y(1) - 2 * y(2) - 8 * y(3) + 2 * exp(-2 * t);
dydt(3) = 3 * y(1) + 3 * y(2) + 2 * y(3) + cos(3 * t);
end

SLIDE 61
syms y x
y=3*sin(x)+7*cos(5*x);
lhs=diff(y,x,2)+y
rhs=-5*cos(2*x);
if(isequal(lhs=rhs))
disp('Yes, the equation is true')
else
disp('No,the equation is false')
end
SOL-

You might also like