lboutput.m

来自「DSP 系统带有信号生成 滤波等相关功能 下载后 请用解压缩软件打开」· M 代码 · 共 33 行

ws=str2num(get(handles.ws,'string'));
wp=str2num(get(handles.wp,'string'));
rs=str2num(get(handles.rs,'string'));
rp=str2num(get(handles.rp,'string'));
global t;
global t0;
global t1;
global t2;
global y;
global b;
global a;
global n;
global N;
Ts=0.01;
fs=1/Ts;
t=-(t2-t1)/2:Ts:(t2-t1)/2;
fvtool(b,a);
y3=filter(b,a,y);
axes(handles.axes3);
plot(t,y3);
title('滤波后时域输出');
grid on;

y5=xcorr(y3);
N=length(y5);
y4=fft(y5);
y4=fftshift(y4);
t0=-(t2-t1):Ts:(t2-t1);
f0=2*pi/N*t0;
axes(handles.axes4);
plot(f0,y4);
title('滤波后频域输出');
grid on;