数字序列的表示方法.docVIP

数字序列的表示方法 两种表示方法 1.数学表达式 2.图形表示 虽然横坐标画成一条连续的直线，但是仅仅对于整数的值才有意义。 采用FFT函数给出两种典型信号的时域/频域图。 clf fs=100; N=128; n=0:N-1; t=n/fs; x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t)+cos(2*pi*30*t); y=fft(x,N); mag=abs(y); f=(0:length(y)-1)*fs/length(y); subplot(211) stem(n,x) xlabel(n); ylabel(x); title(时域图) subplot(212) plot(f,mag); xlabel(f); ylabel(magtude); title(频域图) grid 从上图我们可以看到，信号有三种不同频率频率不同肤质的信号叠加而成 采样定理 clf fs=100; N=16; n=0:N-1; x=[(n-1)==0]; y=fft(x,N); mag=abs(y); f=(0:length(y)-1)*fs/length(y); subplot(211) stem(n,x) xlabel(n); ylabel(x); title(单位脉冲序列) subplot(212) plot(f,mag); xlabel(Frequency (Hz)); ylabel(Magnitude); title(频谱图) grid 绘制数字信号系统H(z)的频率响应、单位冲激响应图形，并加以分析。 频率响应的定义：是一个与系统的特性相关的量，称为单位脉冲相应为的系统的频率响应。利用matlab提供的freqz可以计算序列的傅里叶变换在给定离散频率上的取样值 两边求z变换 求 用freqz函数 对进行逆z变换的 b=1; a=[1,-1,0.75]; N=256; [H,w]=freqz(b,a,N,whole); magH=abs(H(1:N)); w=w(1:N); subplot(2,1,1) plot(w/pi,magH); ylabel(mag); xlabel(frequency unit: pi) title(magnitude response); [h,n]=impz(b,a,N) subplot(2,1,2) plot(n,h) xlabel(n) ylabel(h(n)); title(单位冲击响应) clf fs=100; N=128; n=0:N-1; t=n/fs; x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); b=-2; a=[1,-0.75,-0.25,0.1875]; y=filter(b,a,x); subplot(2,2,1) plot(t,x) xlabel(n) ylabel(x) subplot(2,2,2) plot(t,y) xlabel(n) ylabel(y) x1=fft(x,N); mag=abs(x1); f=(0:length(x1)-1)*fs/length(x1); subplot(223) plot(f,mag); xlabel(Frequency (Hz)); ylabel(Magnitude); title(x-fft) grid y1=fft(y,N); mag=abs(y1); f=(0:length(y1)-1)*fs/length(y1); subplot(224) plot(f,mag); xlabel(Frequency (Hz)); ylabel(Magnitude); title(y-fft) grid 经过对比可以发现，输入信号经过处理幅值和相位都发生了变化，通过频谱图我们可以看到输出信号的幅值要略微大于输入信号。