ASSIGNMENT Course Code Course Name Programme Department Faculty
ECC201A Signals and Systems B.TECH COMPUTER SCIENCE AND ENGINEERING FET
Name of the Student
B.PRAVEEN REDDY
Reg. No
17ETCS002043
Semester/Year
4TH / 2ND
Course Leader/s
Dr T. Christy Bobby
i
Declaration Sheet Student Name
B. PRAVEEN REDDY
Reg. No
17ETCS002043
Programme
B.TECH
Course Code
ECC201A
Course Title
Signals and Systems
Course Date
18/02/19
Course Leader
Dr T. Christy Bobby
Semester/Year 4TH / 2ND
to
18/03/19
Declaration The assignment submitted herewith is a result of my own investigations and that I have conformed to the guidelines against plagiarism as laid out in the Student Handbook. All sections of the text and results, which have been obtained from other sources, are fully referenced. I understand that cheating and plagiarism constitute a breach of University regulations and will be dealt with accordingly.
Signature of the Student
Date
18/03/19
Submission date stamp (by Examination & Assessment Section)
Signature of the Course Leader and date
Signature of the Reviewer and date
ii
Faculty of Engineering and Technology Ramaiah University of Applied Sciences Department Semester/Batch
Electronic and Communication Engineering
B. Tech. in CSE
Programme
rd
4 /2017
ECC201A Course Code Course Leader(s) Ms. Prafulla Kumari K. S & Dr T. Christy Bobby
Signals and Systems
Course Title
Assignment β 02 Reg.No.
Se cti on s Pa rt A Pa rt B
17ETCS002043
Name of Student B. PRAVEEN REDDY
Marks Marking Scheme
A1
MM ax ark s
Essay on DCT for Audio Signal Processing
First Examiner Marks
Moderator
5
Part-A Max Marks
5
B 1.1
Formulation and solution of the difference equation
6
B 1.2
Plotting of traffic
2
B 1.3
Comments on the variability of the traffic
2
B.1 Max
10
2
B 2.3
Computation of Laplace Transforms π(π ), π»(π ) and π(π ) Computation of the response when only secret signal is sent Recovery of the message signal
B 2.4
Modification for the case with convolved inputs
3
Marks
B 2.1 B 2.2
2 3
B.2 Max Marks
10
Total Assignment Marks
25
iii
Course Marks Tabulation Component-1 (B) Assignment A
First Examiner
Remarks
Moderator
Remarks
B.1 B.2 Marks (Max 25 )
Signature of First Examiner
Signature of Moderator
Please note: 1. Documental evidence for all the components/parts of the assessment such as the reports, photographs, laboratory exam / tool tests are required to be attached to the assignment report in a proper order. 2. The First Examiner is required to mark the comments in RED ink and the Second Examinerβs comments should be in GREEN ink. 3. The marks for all the questions of the assignment have to be written only in the Component β CET B: Assignment table. 4. If the variation between the marks awarded by the first examiner and the second examiner lies within +/- 3 marks, then the marks allotted by the first examiner is considered to be final. If the variation is more than +/- 3 marks then both the examiners should resolve the issue in consultation with the Chairman BoE.
iv
PART A Discrete Cosine Transform The Discrete Cosine Transform (DCT) is a transform that is very common when encoding video and audio tracks on computers. Many "codecs" for movies rely on DCT concepts for compressing and encoding video files. The DCT can also be used to analyze the spectral components of images as well. The DCT is very similar to the DFT, except the output values are all real numbers, and the output vector is approximately twice as long as the DFT output. It expresses a sequence of finite data points in terms of sum of cosine functions. Audio Signal Processing Audio signal processing is at the heart of recording, enhancing, storing and transmitting audio content. Audio signal processing is used to convert between analog and digital formats, to cut or boost selected frequency ranges, to remove unwanted noise, to add effects and to obtain many other desired results. This process can be done on an ordinary PC or laptop, as well as specialized recording equipment. An audio signal is a representation of sound, typically as an electrical voltage. Audio signals have frequencies in the audio frequency range of roughly 20 to 20,000 Hz (the limits of human hearing). Audio signals may be synthesized directly, or may originate at a transducer such as a microphone, musical instrument pickup, phonograph cartridge, or tape head. Loudspeakers or headphones convert an electrical audio signal into sound. Digital representations of audio signals exist in a variety of formats. Discrete Cosine Transform for Audio Signal Processing The discrete cosine transform(DCT) signal processing technique come with many applications on a great number of engineering fields. Here I propose to apply techniques to the compression of audio signals. Using spectral analysis and the properties of the DCT, we can treat audio signals as sparse signals in the frequency domain. This is especially true for sounds representing tones. Here I propose the use of DCT to obtain an efficient representation of audio signals, especially when they are sparse in the frequency domain. By using the DCT as signal preprocessor in order to obtain a sparse representation in the frequency domain, here I show that the subsequent application of the signals with less information than the wellknown sampling theorem. This means that our results could be the basis for a new compression method for audio and speech signals. Conclusion I have proposed an efficient joint implementation of DCT, as a method to obtain a sparse audio signal representation, and the application of the compressive sampling algorithm to this sparse signal. The DCT speech signal representation has the ability to pack input data into as few coefficients as possible. This allows the quantizer to discard coefficients with relatively small amplitudes without introducing audio distortion in the reconstructed signal. Although the compressive sampling technique is used primarily for compression sample images, we achieve reasonable results due to the preprocessing of the audio signal.
PART B
v
Solution to Question No. B1: B1.1 Formulation and solution of the difference equation: The given data for this question is follow: ο¨π¦(π) = π1π¦(π β 1) + π2 π¦(π β 2) + π₯(π); ο¨π₯(π) = π΄πππ (ππ) ; π¦(0) = π(0), π¦(1) = π(1) The values of π1 , π2 , π3 , π΄ , π(0) & π(1) are as follow: University ID
a1 a2
A
Student No:
5
6
-4
d(0) d(1) 4
5
So, The given differential equation is follow: ο¨ π¦(π) = π1π¦(π β 1) + π2 π¦(π β 2) + π₯(π) Now, Substitute the value of π1&π2 in the given equation ο¨ π¦(π) = 5π¦(π β 1) β 4π¦(π β 2) + π₯(π) ο¨π¦(π) β 5 π¦(π β 1) + 4 π¦(π β 2) = π₯(π) -----------------------(i) Here, the given equation is discrete signals so, the auxiliary equation is of the given form ο¨1 β 5πβ1 + 4πβ2 = 0 ο¨π2 β 5π + 4 = 0 ο¨π2 β 1π β 4π + 4 = 0 ο¨π(π β 1) β 4(π β 1) = 0 (now, take (π β 1) common in next step) ο¨(π β 1) β (π β 4) = 0 ο¨π = 1
&
π=4
Therefore, the natural response of the given differential equation of discrete signals are as follow, ο¨π¦ π (π) = πΆ1 (1)π + πΆ2 (4)π --------------------------(ii) Now, we have to find the particular solution first then we can find the forced response of the given system: For Particular solution: we have the given value of π₯(π) = π΄πππ (ππ) Now, we can write π₯(π) = 6 β (β1)π [Note: here replace A by 8 and we can write πππ (ππ) = (β1)π ] Now, the value of π¦ π (π) = π β (β1)π -----------(iii) Now, put the value of eq (iii) in eq (i) ο¨π(β1)π β 5π(β1)(πβ1) + 4π(β1)(πβ2) = 6(β1)π 5π(β1)π 4π(β1)π ) + ( ) (β1) (β1)2
ο¨π(β1)π β (
= 6(β1)π [Take (β1)π common in next step]
ο¨(π + 5π + 4π) β (β1)π = 6(β1)π ο¨10 π = 6 3
ο¨π = 5
vi
3 5
Hence, the solution of particular solution is π¦ π (π) = ( ) β (β1)π --------------(iv) Now, we have to write the equation of force response which is follow: ο¨π¦ π (π) = π¦ π (π) + π¦ π (π) ο¨Add eq (ii) and eq (iv) 3
ο¨π¦ π (π) = πΆ1 (1)π + πΆ2 (4)π + (5) β (β1)π ------------------(v) Now, for finding the value of πΆ1 & πΆ2 we have to apply the given initial conditions: Initial conditions are as follow: ο¨ π¦(0) = π(0), π¦(1) = π(1) ο¨π¦(0) = 4, π¦(1) = 5 ο¨Now apply the initial conditions in eq (v) we get the value of πΆ1 & πΆ2 3 5
ο¨π¦(0) = πΆ1 (1)0 + πΆ2 (4)0 + ( ) β (β1)0 3
ο¨π¦(0) = πΆ1 + πΆ2 + 5 3 5
ο¨πΆ1 + πΆ2 + = 4 3
ο¨πΆ1 + πΆ2 = 4 β 5 ο¨πΆ1 + πΆ2 =
20β3 5
ο¨πΆ1 + πΆ2 =
17 5
-------------------(vi)
Now, for π¦(1) 6
ο¨ π¦(0) = πΆ1 (1)1 + πΆ2 (4)1 + (15) β (β1)1 ο¨π¦(1) = 1πΆ1 + 4πΆ2 β
3 5
3 5
ο¨1πΆ1 + 4πΆ2 β = 5 ο¨1πΆ1 + 4πΆ2 = 5 + ο¨πΆ1 + 4πΆ2 =
25+3 5
ο¨πΆ1 + 4πΆ2 =
28 5
3 5
-------------------(vii)
Now, I have to find the value πΆ1 & πΆ2 by using the eq (vi) & (vii) Now, subtract the eq (vi) β eq (vii) ο¨πΆ1 + πΆ2 =
17 5 28
ο¨πΆ1 + 4πΆ2 = 5 -------------------------ο¨
β3πΆ2 =
ο¨
βπΆ2 =
17β28 5
β11 15
11
ο¨πΆ2 = 15
vii
Now, put the value of πΆ2 = ο¨πΆ1 + πΆ2 =
11 in 15
eq (vi)
17 5
ο¨πΆ1 +
11 15
17 5
ο¨πΆ1 =
17 11 β 5 15
ο¨πΆ1 =
51β11 15
=
40
ο¨πΆ1 = 15 8
ο¨πΆ1 = 3 Now substitute the value of πΆ1 & πΆ2 in the forced equation we get the total solution of the differential equation: 6 15
ο¨π¦(π) = πΆ1 (1)π + πΆ2 (4)π + ( ) β (β1)π 8
11
6
ο¨ π¦(π) = (3) (1)π + (15) (4)π + (15) β (β1)π B1.2 Plotting of traffic: Here, we have to plot the traffic as a function of π = 0,1, β¦ . .50. For this we have the total solution of the differential equation 8
11
6
ο¨π¦(π) = (3) (1)π + (15) (4)π + (15) β (β1)π now, we have to plot the discrete signal graph using stem plot in MATLAB: MATLAB CODE: n = [0:50]; y = (8/3).*(1).^n+(11/15).*(4).^n+(6/15).*(-1).^n; fprintf('%f\n',y); stem(n,y,'fill','--') hold on grid on xlabel('value of n = 1..50') ylabel('value of discrete signals') title('Plot of the discrete signal') Command window: 3.800000 5.200000 14.800000 49.200000 190.800000 753.200000 3006.800000 12017.200000
viii
48062.800000 192241.200000 768958.800000 3075825.200000 12303294.800000 49213169.200000 196852670.800000 787410673.200000 3149642686.800000 12598570737.199999 50394282942.799995 201577131761.199980 806308527038.799930 3225234108145.199700 12900936432574.799000 51603745730289.195000 206414982921150.780000 825659931684593.120000 3302639726738366.500000 13210558906953456.000000 52842235627813816.000000 211368942511255260.000000 845475770045021060.000000 3381903080180084200.000000 13527612320720337000.000000 54110449282881348000.000000 216441797131525390000.000000 865767188526101560000.000000 3463068754104406200000.000000 13852275016417625000000.000000 55409100065670500000000.000000 221636400262682000000000.000000 886545601050728000000000.000000 3546182404202912000000000.000000 14184729616811648000000000.000000 56738918467246592000000000.000000
ix
226955673868986370000000000.000000 907822695475945470000000000.000000 3631290781903781900000000000.000000 14525163127615128000000000000.000000 58100652510460510000000000000.000000 232402610041842040000000000000.000000 929610440167368160000000000000.000000 GRAPH: 29
10
Plot of the discrete signal
x 10
9
value of discrete signals
8 7 6 5 4 3 2 1 0
0
5
10
15
20 25 30 value of n = 1..50
35
40
45
50
Figure 1: Discrete plot of the signal obtained as a solution in above question. B1.3 Comments on the variability of the traffic: As we can see in the above graph that traffic is varying with respect to time such that for every increase in time there is some increase in traffic. In the Beginning phases we can see that the increase in traffic with respect to time is gradual but as it reaches the 45 unit mark a very high increase is noticed. Thus beyond this point the increase in traffic/variable in traffic is exponential. According to the values provide to me the variation in traffic results in a positive value. Thus Naturally it is a growth Rate.
x
Question No. B2 Solution to Question No. B2: The data is follows:π(π‘) = π βππ‘ π’(βπ‘) ο¨π€(π‘) = π βππ‘ π’(π‘) ο¨β(π‘) = π π’(π‘) & π (π‘) = π(π‘) + π€(π‘) The value of π , π & π πππ ππ ππππππ€: a
b
c
1
4
8
B2.1 Computation of Laplace Transforms π(π ), π»(π ) and π(π ): Here, we have to find the Laplace transforms of π(π ), π»(π )& π(π ) Now, first I have to find π(π )first, For this we have the given eq of π€(π‘) = π βππ‘ π’(π‘) 1 π‘>0 Here, π’(π‘) is unit step function which is defined as π’(π‘) = { } 0 ππ‘βπππ€ππ π Graph of u(t) is follow: ο¨ π€(π‘) = π βππ‘ π’(π‘) Now, apply Laplace transform β
ο¨πΏ[π€(π‘)] = β«ββ π βπ π‘ π βππ‘ π’(π‘)ππ‘ 0
β
ο¨ πΏ[π€(π‘)] = β«ββ π βπ π‘ π βππ‘ β 0ππ‘ + β«0 π βπ π‘ π βππ‘ β 1ππ‘ β
ο¨π(π ) = β«0 π β(π +π)π‘ ππ‘ π β(π +π)π‘
ο¨ π(π ) = [ β(π +π) ] ο¨π(π ) =
β 0
1 [π ββ β(π +π)
β π0]
1
ο¨π(π ) = π +π [Now, replace the b by 3 in next line we get,] 1
ο¨ π(π ) = π +4 ---------------(i) 1 π +4
Hence, Laplace transform of π(π ) = Now, I have to find the value of π»(π )
For, this I have the given value of β(π‘) = π π’(π‘) Now, apply Laplace transform β
ο¨πΏ[β(π‘)] = β«ββ π βπ π‘ β π β π’(π‘)ππ‘ 0
β
ο¨πΏ[β(π‘)] = π β [β«ββ π βπ π‘ β 0ππ‘ + β«0 π βπ π‘ β 1ππ‘ ] β
ο¨π»(π ) = π β β«0 π βπ π‘ ππ‘
xi
π βπ π‘
ο¨ π»(π ) = π β [ βπ ] ο¨π»(π ) =
π [π ββ βπ
β 0
β π 0]
π
ο¨π»(π ) = π [Now, replace the c by 9 in next line we get,] 8
ο¨ π»(π ) = π --------(ii) Hence, the Laplace transform of π»(π ) =
8 π
Now, we have to find the value of π¦(π ) For finding this we have to draw the L.T.I system and Find π¦(π‘) L.T.I system is follow:
Ζ©
From the above L.T.I system we can calculate π¦(π‘) ο¨π¦(π‘) = β(π‘) β π (π‘) ο¨π¦(π‘) = β(π‘) β [π(π‘) + π€(π‘)] Now, apply Laplace transform ο¨πΏ[ π¦(π‘)] = π»(π ) β [π(π ) + π(π )] ο¨ π(π ) = π»(π ) β π(π ) + π»(π ) β π(π ) ----------------(iii) Here, I have to find the value of π(π‘) first, ο¨ π(π‘) = π βππ‘ π’(βπ‘) Here, I have π’(βπ‘) as a function so the value get reversed in π’(βπ‘) 1 π‘<0 Here, π’(βπ‘) is unit step function which is defined as π’(βπ‘) = { } 0 ππ‘βπππ€ππ π Now, take the value of π(π‘) = π βππ‘ π’(βπ‘) Now, apply Laplace transform β
ο¨πΏ[π(π‘)] = β«ββ π βπ π‘ β π βππ‘ β π’(βπ‘)ππ‘ 0
β
ο¨ πΏ[π(π‘)] = β«ββ π βπ π‘ β π βππ‘ β 1ππ‘ + β«0 π βπ π‘ β π βππ‘ β 0ππ‘ 0
ο¨π(π ) = β«ββ π β(π +π)π‘ ππ‘ [Note: here, I Consider β (π + π) = πΆ] 0
ο¨ π(π ) = β«ββ π πΆπ‘ ππ‘ ο¨π(π ) = [
0 π πΆπ‘ ] πΆ ββ
1
ο¨ π(π ) = πΆ [π 0 β π ββ ] 1
ο¨ π(π ) = β(π +π) [1 β 0] [replace πΆ ππ¦ β (π + π)] 1
ο¨ π(π ) = β(π +1) [πππππππ π ππ¦ 1] 1
ο¨π(π ) = β(π +1) ----------(iv) 1
Hence, the value of π(π ) = β(π +1) Now, From eq (iii)
xii
ο¨ π(π ) = π»(π ) β π(π ) + π»(π ) β π(π ) Now, put the value of π»(π ), π(π )&π(π ) we get the value of π(π )as 8
1
8
1
ο¨ π(π ) = π β β(π +1) + π β π +4 8
1
1
ο¨ π(π ) = [ + ] π β(π +1) π +4 8 (π +1)β(π +4)
ο¨ π(π ) = π [ (π +1)(π +4) ] 8
β3
ο¨ π(π ) = π [(π +1)(π +4)] β24
ο¨π(π ) = [π β(π +1)β(π +4)] ----------(v) Hence, we get the value of π(π ) = [
β24 ] π β(π +1)β(π +4)
B2.2 Computation of the response when only secret signal is sent: ο¨π£(π‘) = π€(π‘) β β(π‘) Now, apply Laplace transform ο¨πΏ[π£(π‘)] = π(π ) β π»(π ) Now, put the value of π(π ) πππ π»(π ) in the above equation: From eq (i) & (ii) ο¨π(π ) = π(π ) β π»(π ) -------(vi) ο¨π(π ) =
1 π +4
β
8 π
8
ο¨π(π ) = (π +4)βπ -----------------(vii) B2.3 Recovery of the message signal: From, equation(iii) we have ο¨ π(π ) = π»(π ) β π(π ) + π»(π ) β π(π ) And, from equation (v) we have ο¨ π(π ) = π(π ) β π»(π ) Now, subtract the eq(iii) β eq(vi) we get, ο¨ π(π ) β π(π ) = π»(π ) β π(π ) + (π»(π ) β π(π )) β (π(π ) β π»(π )) ο¨ π(π ) β π(π ) = π»(π ) β π(π ) ο¨π(π ) =
(π(π )βπ(π )) π»(π )
Now, replace the value of π(π ), π(π )&π»(π ) from eq (vii) , (v) & (ii) we get,
xiii
ο¨π(π ) =
(
β24 8 )β( ) π β(π +1)β(π +4) π β(π +4) 8 π
β3
8
8
[take π common then the π get divided and we get the eq as]
1
ο¨ π(π ) = ((π +1)β(π +4)) β ((π +4)) β3βπ β1
ο¨π(π ) = (π +1)β(π +4) (π +4)
ο¨π(π ) = β (π +1)β(π +4) 1
ο¨π(π ) = β(π +1) Now, we have to find the inverse Laplace of π(π ) to get the value of π(π‘) 1
ο¨πΏβ1 [π(π )] = πΏβ1 [β(π +1)] 1
ο¨π(π‘) = βπ β1π‘ [Note: Laplace inverse of πΏβ1 [β(π +1)] = βπ β1π‘ ] Now, we can replace β1 as π’(βπ‘) ο¨ π(π‘) = π β1π‘ β π’(βπ‘) Hence, we get the correct recovery signal. B2.4 Modification for the case with convolved inputs: Here, in this question we have to Analyze the effect of change in operation between the message and the watermark on the response (π€(π‘) and π(π‘) are convolved instead of being added. Here we have to draw the L.T.I system to find convolved inputs: The L.T.I system is as follow:
Now, here we have to take π (π‘) = π€(π‘) β π(π‘) (as per L.T.I system) Then, From above eq of π¦(π‘) = β(π‘) β π (π‘) Now, apply Laplace transform we get, ο¨πΏ[π¦(π‘)] = πΏ[β(π‘) β π€(π‘) β π(π‘)] ο¨π(π ) = π»(π ) β π(π ) β π(π ) Now, put the value of π»(π ), π(π ) & π(π ) we get, 8
1
1
ο¨ π(π ) = π β π +4 β β(π +1) β8
ο¨ π(π ) = π β(π +4)β(π +1)
xiv
15