1. 连续函数的数字处理
    1. x(t)
      1. n=t/T+k (平移到n=1开始)
      2. 离散幅值为连续幅值的1/T倍
      3. x[n]
        1. 0≤n<N
        2. FFT
        3. X[k] 有限离散频谱
        4. 横轴用w=2π(k-1)/N
        5. X(jw)
          1. 此时[0,pi]正常,[pi,2pi]应该移再[-pi,0]上
        6. w2=(k-1-(M-1)/2)/M;
        7. X(jw) [-pi,pi]的DTFT频谱
  2. 滤波器设计
    1. IIR
    2. FIR
  3. 画图
    1. plot
      1. ylabel,xlabel
      2. legend 标签说明
      3. hold on 在原图上画 hold off
      4. axis([xmin xmax ymin ymax]) 放在plot 后面
      5. title 标题
      6. figure 开辟另一个窗口
      7. subplot 在一个窗口里面画几个图
      8. 例子
        1. plot(f,y3db,'--r');
      9. grid on 显示网格
    2. stem
    3. bar
      1. 绘制直方图
  4. 矩阵与向量
    1. 查看
      1. length
      2. x(1:10)
      3. 从矩阵中取向量
      4. 从矩阵中取矩阵
    2. 矩阵的产生
      1. rand
      2. randn
      3. ones(10,1)
      4. linspace(X1,X2,N)产生一个有N个位于 X1,X2之间的数的向量
      5. logspace(a,b,n)
    3. 运算
      1. 叉乘:cross
      2. norm:求范数 可用于求模
      3. 在代数运算中,出现0/0 的情况,会使结果为NaN
    4. 数值
      1. 舍/入
        1. ceil 取上界 floor 取下界
        2. fix去向零的数靠拢
      2. 求余数
        1. rem:rem(x,y)求x/y的余数
        2. mod
  5. DSP积累
    1. 卷积
      1. 3.9.txt
    2. y=fliter(b,a,x)
    3. 滤波器
      1. 巴特沃兹
        1. [z,p,k]=buttap(N)
          1. 用于计算N阶巴特沃斯归一化(3dB截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子。
      2. [num,den]=zp2tf(Z,P,K)
        1. zero plot to transfer function
      3. [N,wc]=buttord(wp,ws,αp,αs,‘s')
        1. 加上's'就是模拟的
        2. 如果是数字的,则wp,ws要除以π
      4. [z,p,k]=butter(N,wc]
        1. 加上's'就不用归一化
      5. h=freqs(num,den,w)(w是弧度/s)
        1. 要与freqz区分
      6. [H w]=freqz(numd,dend) H是幅频特性曲线
      7. 切比雪夫
        1. [z,p,k]=cheblap(N,Rp)
        2. [num,den]=cheby1(N,Rp,Wn,'type,'s')
          1. type='high', 'low……
        3. [n,Wp]=cheb1ord(Wp,Ws,Rp,Rs,'s')
      8. 注意!!
        1. Wn,Wp等这些,都是弧度/s,因此如果给出频率的话,要将其乘以2π再计算!!!
    4. fft变换
      1. 若要对一个窗口信号做fft,如果x=ones(1,500),fft(x)则会产生y的第一个为500,其他为0 的情况,相当于对常数1做fft,解决方法是fft(x,10000)
  6. 文件流
    1. 语音
  7. 编程规范
    1. 函数:function [y1,y2……]=name(x1,x2……)
    2. clc 清空工作区
  8. 图像处理
    1. I/O
      1. I=imread('gray.png');
      2. imshow(I)
    2. 信息
      1. 直方图
        1. 先I=rgb2gray(I);
        2. 再imhist(I)
    3. 图像增强
      1. 均衡
        1. histeq
    4. 程序
      1. 好的程序
        1. 圆形区域提取直方图
      2. 示例程序
        1. 直方图均衡化
        2. 图像局部增强