数字信号处理MATLAB笔记

全文框架

1.函数笔记

fft():快速傅里叶变换
Y = fft(X,n,dim)
计算x的n点DFT，x长度不够n时补0，比n长时截短，dim为对x做傅里叶变换的维数。

ifft():傅里叶反变换
y = ifft(X,n)
计算X的n点DFT反变换

Stem():火柴杆图
stem(X,Y)
显示离散序列，X与Y大小必须相同，一一对应。如Y是一个长度为6的序列，x可以为0:5，若X省略，程序会自动把横坐标变为Y的大小。
若一个连续序列用stem显示，会变成实心的，因为离散序列每个点下都有一条竖线，因此最好用plot显示连续序列。

Abs():绝对值或幅值
Y = abs(X)
求绝对值或复数的幅度，可用于求幅度谱。

Angle():相位
P = angle(Z)
求复数Z的角度和相位，返回值在±π之间。

Conj():求共轭
ZC = conj(Z)

矩形信号：
矩形序列R(N):可以用zeros和ones构造，如：R20=[zeros(1,40),ones(1,20),zeros(1,41)];
连续矩形信号可用rectpuls()函数，y = rectpuls(t,w)，t为时间序列，w为宽度，宽度以t=0为对称
单位脉冲信号：
只需在原点处设一个脉冲即可

Tripuls():三角波信号
y = tripuls(T,width,skew)
T为时间序列，width为三角波宽度，skew为斜度，范围0-1，默认为0，即对称三角波。
注意三角波以T=0对称，若要完整显示，时间序列要从负到正（如：T=-30:1:30）。
若要移动三角波，对T进行加减（如T+50表示向右移50单位）

Impulse():单位脉冲响应
Y=Impulse(sys,t)
Sys为系统函数，t为时间序列（Ts:dt:Te），y为单位阶跃响应，t省略则自动生成时间间隔

Step():单位阶跃响应
Y=step(sys,t)
Sys为系统函数，t为时间序列（Ts:dt:Te），y为单位阶跃响应，t省略则自动生成时间间隔
[y,t] = step(sys)
返回单位阶跃响应y，时间间隔t

Filter():一维滤波
y = filter(b,a,x)
a,b为滤波器的传递函数的分子分母，x为向量或矩阵，y为输出，大小与x相同。
也可以用来求系统响应，比如系统差分方程： ，则b和a为： ，x即为x(n)，输出y为系统响应

Impz():求离散LTI系统的单位冲激响应
h=impz(b,a,k)
b为x的向量，a为y的向量，k为输出序列的取值范围。h为单位脉冲响应。

tf()：系统函数
sys = tf(Numerator,Denominator)
创建系统函数，numerator为分子向量，denominator为分母向量
后面可加参数，如：‘Variable’,‘p’，意思是把p作为变量

Freqz():离散系统频率响应
[h,w] = freqz(b,a,n)
B为分子，a为分母，返回n个频率点的频率响应，放在h向量中，n个频率放在向量w中，这n个频率均匀设置在（0，pi）上，缺省时默认512

Freqs():连续系统频率响应
h = freqs(b,a,w)
B，a为分子分母，w为角频率向量，返回频率响应h

Ellipord():椭圆滤波器参数
[n,Wp] = ellipord(Wp,Ws,Rp,Rs,’s’)
Rp、Rs分别为通带最大波纹和阻带最小衰减；wp