Week 3 Lecture
Week 3 Lecture
a = 7
E=3
E = 3
d = 1×4
5 10 4 9
g = 2×3
7 49 13
21 1 343
whos
E 1x1 8 double
a 1x1 8 double
d 1x4 32 double
g 2x3 48 double
input function
midterm1=75
midterm1 = 75
midterm2=93
midterm2 = 93
final=68
final = 68
avg_points=(midterm1+midterm2+final)/3
avg_points = 78.6667
table(midterm1,midterm2, final)
1 75 93 68
clear;clc;
1
student1=input('what is your name','s') % characters data type
% 's' means string type data
m1=input('please enter your score') % enter numbers by default
Chapter 5 Plotting
clear;clc;clf
x=1:5:10
y=10:10:100
hold on
plot(x,y,'o')
plot(x,y, '--g')
hold off
plot(x,y,'^r')
function plot
f='x^2+4*sin(2*x)-1'
fplot(f)
fplot(f,[-3 3])
2
plot(year, sales, 'o')
createFit(year, sales)
y=sin(x).*cos(x)
y2=tan(x)
y3=csc(x).*tan(x)
%creating 3x1 plots
subplot(3,1,1)
plot(x,y)
subplot(3,1,2)
plot(x,y2)
subplot(3,1,3)
plot(x,y3)
syms x
y=cos(x)*sin(x)
dy=diff(y)
d2y=diff(y,2)
d3y=diff(y,3)
d4y=diff(y,4)
d5y=diff(y,5)
subplot(3,2,1)
fplot(y)
subplot(3,2,2)
fplot(dy)
subplot(3,2,3)
fplot(d2y)
subplot(3,2,4)
fplot(d3y)
subplot(3,2,5)
fplot(d4y)
subplot(3,2,6)
fplot(d5y)
random graphs
3
clear;clc;clf
syms x
y=cos(x)*sin(x)
dy=diff(y)
d2y=diff(y,2)
d3y=diff(y,3)
d4y=diff(y,4)
d5y=diff(y,5)
subplot(3,3,1)
fplot(y)
subplot(3,3,2)
fplot(dy)
subplot(3,3,3)
fplot(d2y)
subplot(3,3,[4 7])
fplot(d3y)
subplot(3,3,5)
fplot(d4y)
subplot(3,3,[5 6 8 9])
fplot(d5y)
Other graphs
clear;clc;clf
yr=[1988:1994]
sle=[8 12 20 22 18 24 27]
ylabel('sle')
stairs(yr,sle)
stem(yr,sle)
Grade=[11 18 26 9 5]
pie3(Grade)
histogram
y=[58 73 73 53 50 48 56 73 73 66 69 63 74 82 84 91 93 89 91 80 59 69 56 64
63 66 64 74 63 69 ]
4
hist(y)
hist(y,3)
xlabel('T')
ylabel('Frequency')
x=45:10:95
hist(y,x)
t=linspace(0,2*pi,200)
r=3*cos(0.5*t).^2+t
polar(t,r)
Chapter 4 Homework
problem 11
clear;clc;
h: -500:500:10000
p=29.921*(1-6.8753e-6*h)
Tb=49.16*log(p)+44.932
tbl=[h' Tb']
disp(' Boling')
disp(' Altitude Temperature')
disp(' (ft) (degF)')
5
disp(tbl)
Problem 23
clear;clc;
m=1:7
m = 1×7
1 2 3 4 5 6 7
a=[cosd(28.5) 1 0 0 0 0 0
sind(28.5) 0 0 0 0 0 0
-cosd(28.5) 0 -cosd(58.4) 0 cosd(58.4) cosd(28.5) 0
-sind(28.5) 0 -sind(58.4) 0 -sind(58.4) -sind(28.5) 0
0 0 0 -1 -cosd(58.4) 0 1
0 0 0 0 0 sind(28.5) 0
0 0 0 0 0 -cosd(28.5) -1]
a = 7×7
0.87882 1 0 0 0 0
0.47716 0 0 0 0 0
-0.87882 0 -0.52399 0 0.52399 0.87882
-0.47716 0 -0.85173 0 -0.85173 -0.47716
0 0 0 -1 -0.52399 0
0 0 0 0 0 0.47716
0 0 0 0 0 -0.87882
b = 7×1
3000
-6521
-3000
0
0
-7479
0
x=a\b
x = 7×1
-13666
15010
9397.6
10086
7039.6
-15674
13775
tbl=[m' x]
tbl = 7×2
1 -13666
2 15010
3 9397.6
4 10086
5 7039.6
6 -15674
7 13775
6
disp(' Member Force (lb)')
disp(tbl)
1 -13666
2 15010
3 9397.6
4 10086
5 7039.6
6 -15674
7 13775
Problem 27
clear;clc;
A=[1 2 1 1
2 3 0 1
1 4 1 0
1 3 2 0]
A = 4×4
1 2 1 1
2 3 0 1
1 4 1 0
1 3 2 0
B = 4×1
5
12
11
8
Points=A\B
Points = 4×1
4
2
-1
-2
Chapter 5 Homework
Problem 15
clear;clc;
x=linspace(-1,1,200)
x = 1×200
-1 -0.98995 -0.9799 -0.96985 -0.9598 -0.94975
y1=x1.^(2/3)+sqrt(1-x1.^2)
7
Unrecognized function or variable 'x1'.
y2=x1.^(2/3)-sqrt(1-x1.^2)
x1=[0:0.01:1]
x2=[0:-0.01:-1]
plot(x1, y1, x1, y2, x2, -y1, x2,-y2)
plot(x1, y1, x1, y2, x2, y1, x2,y2)
Problem 25
clear;clc;
th=70
th =
70
v0=200
v0 =
200
g=9.81
g =
9.81
syms t
solve((g*t^2)/2 == v0*t*sin(th))
ans =
t = 1×200
0 0.15857 0.31714 0.47571 0.63428 0.79284
ht=v0*t*sin(th)-(g*t.^2)/2
ht = 1×200
0 24.42 48.593 72.519 96.199 119.63
vt=sqrt(v0^2-2*v0*g*t*sin(th)+g^2*t.^2)
vt = 1×200
200 198.8 197.6 196.41 195.22 194.04
subplot(2,1,1)
plot(t,ht)
xlabel('Time')
ylabel('Height')
subplot(2,1,2)
8
plot(t,vt)
xlabel('Time')
ylabel('Velocity')
Problem 37