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

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

====================== ======================== ======================

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')

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