ingenieria electrica

1- Definición de variables:

%numero de muestras

n=0:20000;

%frecuencias de las señales de 90 y 150 Hz, frecuencias de muestreo y de portadora

fm1=150;
fm2=90;
fp=fm1*100;
fs=fp*12;wfm1=(2*pi*fm1)/fs;
wfm2=(2*pi*fm2)/fs;
x150=sin(wfm1.*n);
x90=sin(wfm2.*n);
xs=x150+x90;
xr= x90-x150
2.1- Representación señal PBL

%representación señal PBL

Ac=10;
m=0.2;wp=(2*pi*fp)/fs;
sPBL= Ac.*(1+m.*xs).*cos(wp.*n);
plot(n,sPBL), title(‘señal PBL’)


2.2- Representación señal DBL

%representación señal DBL
%siendo r; r=Ess/Em y suponiendo r=0.5
r=0.5;Ess=linspace(-2,2,20001);
Ecs=Ess/r;
sDBL=Ess.*xr.*cos(wp.*n);
plot(n,sDBL), title(‘señal DBL’)

sumpbldbl=sPBL+sDBL;
plot(n,sumpbldbl), title('suma en el espacio de PBL y DBL');



3- Detector deenvolvente:
%suma de señales PBL+DBL---> verde(y la parte azul rectificada)

%detección de la envolvente de la señal: rectificado+FPB

%rectificación de la señal--->azul

% plot(n,sumpbldbl);

holdon

a=find(sumpbldblroja

[b a]=butter(1,wfm1);
sdet=filter(b,a,sumpbldbl);

plot(n,sdet,'r');
hold on


4- Espectro de la señal antes y despues del diezmado:

%diezmado de la señal,siendo M=factor de diezmado
%y d=intervalo de muertras para el diezmado

M=10;
d=0:2000;
sdiezm=decimate(sdet,M);
plot(d,sdiezm), title('señal diezmada')




%espectro de señales antes ydespués del diezmado

%antes del siezmado
wr=0:n:2*pi;
espdetectada=fft(sdet,wr);
subplot(2,1,1)
plot(abs(espdetectada)), axis ([1600 1700 0 1500]) 


%después del diezmado
wd=0:d:2*pi;espdiezmada=fft(sdiezm,wd);
subplot(2,1,2)
plot(abs(espdiezmada)), axis([0 30 0 1500])



%uencias de corte de los filtros
wc150i=(2*pi*120)/fs;
wc150s=(2*pi*180)/fs;
wc90i=(2*pi*65)/fs;wc90s=(2*pi*120)/fs;

wc90=M.*[wc90i/pi wc90s/pi];
[g f]=ellip(4,0.2,40,wc90);
y1=filter(g,f,sdiezm);
[H90,wh90]=freqz(g,f,2000);
plot(wh90,H90), axis ([0 0.2 -1 1]), title('FT del FPB de la señal de...