Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronics Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

Energy in the signal

Status
Not open for further replies.

praveen450

Member level 5
Joined
May 13, 2011
Messages
89
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,286
Activity points
1,957
Hello,

I'm implementing a matlab code to find energy in the gaussian signal. First I generated a gaussian signal and find the energy in it. Again I differentiated the signal 5 times and find the energy in it. Now, the energy in both the gaussian and 5th derivative gaussian signal should be same. But I'm unable to get this answer. Could anyone help me in this regard.
The matlab code is shown below :


Code ActionScript - [expand]
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
clc
clear
b2=0
c=0.1e-9;
a2=1
x1=linspace(-1,1,2000);
x=x1.*10.^-8;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
in1=exp(-((x)./c).^2);%zero-order gaussian signal
out=abs(in1.^2/50);%power of the signal
%energy =sum(out)%energy of the differentiated signal
energy=trapz(x,out)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%
%der2=-2*in1.*(c^2-2*x.^2)/(c^4);
der2=diff(in1,5)./(x(2)-x(1))^5;
der2=horzcat(der2,0,0,0,0,0)
out=abs(der2.^2/50);%power of the signal
%energy=sum(out)%energy of the differentiated signal
energy=trapz(x,out)

 
Last edited by a moderator:

SunnySkyguy

Advanced Member level 5
Joined
Sep 26, 2007
Messages
6,743
Helped
1,675
Reputation
3,348
Reaction score
1,644
Trophy points
1,413
Location
Richmond Hill, ON, Canada
Activity points
50,732
On what basis do you believe it should be the same energy, if the lower portion of the spectrum is attenuated more than the upper.

I think the Gaussian shaped pulse when differentiated once it has a Guassian like spectrum. But 5 times?

Are you thinking of spectral density or pulse energy?
 

praveen450

Member level 5
Joined
May 13, 2011
Messages
89
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,286
Activity points
1,957
I'm considering pulse energy.
If you take a zero ordered gaussian signal, lets say it has energy E1 in time domain. Now if you differentiate n number of times integrate it over time, the differentiated signal should have same energy. Is it correct?.
 

andre_teprom

Super Moderator
Staff member
Joined
Nov 7, 2006
Messages
9,241
Helped
1,151
Reputation
2,321
Reaction score
1,127
Trophy points
1,403
Location
Brazil
Activity points
53,788
if you differentiate n number of times integrate it over time, the differentiated signal should have same energy. Is it correct?.

Mathematically speaking, not: Integration operation adds a constant factor.
Anyway, differentiation at 5th order seems having no relation to any physical fundament.
 

SunnySkyguy

Advanced Member level 5
Joined
Sep 26, 2007
Messages
6,743
Helped
1,675
Reputation
3,348
Reaction score
1,644
Trophy points
1,413
Location
Richmond Hill, ON, Canada
Activity points
50,732
I'm considering pulse energy.
If you take a zero ordered gaussian signal, lets say it has energy E1 in time domain. Now if you differentiate n number of times integrate it over time, the differentiated signal should have same energy. Is it correct?.

Any pulse, Gaussian or not has a finite DC energy unless it is AC coupled or differentiated causing undershoot on the trail edge where the area above and below zero are now equal. Differentiating the signal attenuates the spectrum 1/s or 6dB per octave or -6dB per harmonic, which changes the energy as repeating pulses have many harmonics up F=/PW50 pulse width then repeats with 2nd harmonics of the repetition rate,

Notice how the 2nd derivative looks like a QRS heart beat ... maybe integrated once. after trace show below :) https://en.wikipedia.org/wiki/File:QRS_complex.png

3878540700_1443205603.jpg


File:QRS_complex.png
 

Status
Not open for further replies.

Similar threads

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Top