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