last1.m

来自「在信道仿真中」· M 代码 · 共 29 行

clc;clear;

fd=100;
N=128;              %初始值

f=linspace(1,10,N); % 画图需要

randn('state',1);
t=1:0.01:10;
p(t)=sin(2*pi*10*t)+cos(2*100*pi*t)+0.2*randn(size(t));
P=fft(p(t),N);
PP=abs(P);      % 原始函数


H=lvbo2(fd,N);  %调用函数

Q=PP.*H;

h=ifft(Q,N);    % fft反变换

disp(h)        %输出最终结果


subplot 211     %图示结果
plot (t,p(t))
text(5,1,'滤波之前的函数图像')
subplot 212
plot(f,h);
text (5,0.6,'滤波之后的函数图像')