Utilizo matlab para escribir el programa de filtro de abeto. Complete la parte del dibujo (necesito explicar en detalle por qué lo escribo de esta manera).
función xt=xtg(N)
La señal x(t) genera una función y muestra la curva característica amplitud-frecuencia de la señal
xt=xtg genera una longitud N, señal AM de frecuencia única xt con ruido aditivo de alta frecuencia, frecuencia de muestreo Fs=1000Hz
Frecuencia portadora fc=Fs/10=100Hz, frecuencia de onda sinusoidal de modulación f0=fc/10= 10Hz
p>N=2000; Fs=1000; T=1/Fs; Tp=N*T;
t=0: T: (N-1)*T ;
fc=Fs/10; f0=fc/10; Frecuencia portadora fc=Fs/10, la frecuencia de la señal de modulación de frecuencia única es f0=Fc/10
mt=cos( 2*pi*f0 *t); Genera una señal de modulación de onda sinusoidal de frecuencia única mt, la frecuencia es f0
ct=cos(2*pi*fc*t); ct, la frecuencia es f0
xt=mt.*ct se multiplica para generar una señal de modulación de frecuencia única xt
nt=2*rand(1, N)-1 genera aleatoriamente; ruido nt
-- ---------Diseñe un filtro de paso alto para filtrar los componentes de baja frecuencia en el ruido nt y generar ruido de paso alto------ ---
fp=150; fs=200; Rp=0.1; As=70; Índice de filtro
fb=;
dev=[10^(-As/20) , (10^(Rp/20)-1)/(10^(Rp/20) 1)];
[n, fo, mo, W]=remezord(fb, m, dev, Fs); Determinar los parámetros necesarios para la función remez
hn=remez(n, fo, mo, W); la función remez para diseñar y filtrar los componentes de baja frecuencia en nt
yt =filter(hn, 1, 10*nt; Filtrar el ruido de baja frecuencia y generar ruido de alta frecuencia yt); p>
====================== ======================== ======================
xt=xt yt ruido más señal
fst=fft(xt, N ); k=0: N-1; f=k/Tp;
subplot(2, 1, 1); ylabel('x(t)');
axis([0, Tp/5, min (xt), max(xt)]('(a) señal más forma de onda de ruido')
subplot(2, 1, 2); plot(f, abs(fst)/max( abs(fst))); title('(b) señal más espectro de ruido') p>
axis([0, Fs/2, 0, 1.2]); xlabel('f/Hz' );ylabel('amplitud')
-------- -------------------------- ------------------------ ------------
Filtro FIR
Implementación de forma de onda
borrar todo; cerrar todo;
== Llame a xtg para generar la señal xt, longitud de xt N=1000 y muestre xt y su espectro ======= ===
N=1000; xt=xtg(N);
fp=120; fs=150; el índice
Usa el método de función de ventana
wc=(fp fs)/Fs;
B=2*pi*(fs-fp)/Fs
Nb=ceil(11*pi/B);
hn=fir1(Nb-1,wc,blackman(Nb));
ywt = fftfilt(hn, xt, N);
plot(...Cómo escribir lo siguiente...Por favor avise