clc
clear all
close all
% approximation to sin function using infinite series formula
x_deg=input('Enter an angle in degree
')
x_rad=x_deg*pi/180;%convert angle in rad
n=50;
% infinite series uses radians
y = 0;
j=1;
while( j<=n)
y = y + (-1)^(j+1)*(x_rad
ENGR 351-Homework Set III
Fall 2014
Due: Oct 1
B.A. DeVantier
General Instructions: You may print neatly or use word processing for the written portion of your
homework. Graphs may be done by hand or computer spreadsheet software. For any work done with
s
Assignment 1
c
2012
Coventry University
225MSE/1
225MSE: Linear Control Systems
Module Leader: J. Mahtani, 1st Term Lecturer: S. Oleksowicz
Assignment 1: State space form, feedback connections, transfer functions, and poles and zeros
This is an individual
PART A:
Free body diagram of both masses m1 and m2.
(
(
)
)
(
(
)
)
Free body diagram for
(
(
)
Free body diagram for
PART B:
Standard form of differential equations for both masses
From above figures we can see the total forces acting on both masses
Question 1:
We are given following specifications
We have
From the table given in lecture notes we can see that the best window which is suited for this is variable
length window that is also known as Kaiser Window .As the signal have corrupted noise of f
Part 1
Running sum filter of length 7 is defined as
[ ]
[
]
So we can write the difference equation of running sum filter by expanding the summation
We get
[ ]
[ ]
[
]
[
]
[
]
[
]
[
]
[
]
Representing the difference equation of 7 point running sum filter
Problem # 01
The relative frequencies in the data are:
r
0
1
2
3
4
f
11/40
14/40
8/40
6/40
1/40
Number of calls per hour
If the process is Poisson
For r = 0, 1, 2, 3, 4 this gives (rounding to 4 dp) the following probabilities:
r
P
0
0.2592
1
0.3500
Probl
Frequency sampled filter in MATLAB:
We are given following specifications
pass band=00.25 Fs
with attenuation=0 db
stop band frequency =0.4 Fs0.5 Fs
with attenuation at least=35 db
In order to calculate the filter order we need normalize transition band
T
Question: Determine maximum bending stress?
Solution:
Maximum bending stress=
Where
max =
Mc
I
is the bending stress, M is the moment about axis. I is the
moment of inertia of plane axis.
Using
M=
I=
w L2
c
1
b h3
12
Max bending stress=
w L2
c
c
1
b h3
12
Design of the type II Cheby shev Highpass filter:
The specification of the filter are given below
Fp=0.7 KHz
Fs=0.5 KHz
Ft=2 KHz
Epsilon=1 dB
AlphaS=32dB
The Matlab code is given here
%Type 2 Chebyshev Highpass Filetr(Inverse)
clc
clear all
Fp= 0.7 ;%KHz