现代信号处理方法1-1.pptVIP

现代信号处理方法 主讲： 王忠仁 2012年2月 授课内容 〇、MATLAB语法要点 一、信号的时频分析 二、Radon-Wigner变换 三、分数阶Fourier变换 四、高阶谱估计 参考书 [1] 张贤达. 现代信号处理. 清华大学出版社. 1995 [2] 张贤达,保铮. 非平稳信号分析与处理. 国防工业出版社. 1998 [3] 张贤达. 时间序列分析-高阶统计量方法. 清华大学出版社. 1996 [4] 王忠仁,张静. 复变函数与积分变换. 高等教育出版社. 2006 0．１　MATLAB语言特点 向量与矩阵的运算非常方便 库函数极其丰富 语言简洁，程序简短 图形功能强大 比C、FORTRAN语言更高级的第四代计算机语言： 机器语言 汇编语言 C-语言 MATLAB 0．２　MATLAB语法要点 算术运算符： 标量运算: 加减乘除幂 + - * / ^ 矩阵运算: 加 减 乘 右除 左除 幂 转置 复共轭转置 + - * / \ ^ .’ ’ 若Ax=b 则x=A\b (常用); 若xA=b 则x=A/b (不常用) 标量情况有: 2/5=0.4; 2\5=2.5 矩阵对应元素运算: 乘 右除 左除 幂 .* ./ .\ .^ 数乘矩阵: c*A 冒号运算符: a:h:b 是以a 为初值，h为步长,终值不超 过b（但再增加一个步长就超过b）的行向量 关系运算符: 小于 小于等于 大于 大于等于 等于 不等于 = = == ~= 逻辑运算符: 与 或 非 | ~ 注释与语句分隔: % 开始注释,直到行末; ; 分隔语句,阻止前面变量的结果显示在命令窗口； , 分隔语句,前面变量的结果可显示在命令窗口。 帮助：在命令窗口键入 help 函数名 循环结构: (1) for 语句 格式: for t=表达式1:表达式2:表达式3 语句体 end 其中:表达式1为循环初值; 表达式2为步长,省略时默认步长为1; 表达式3为循环终值. 例 for n=1:5 for m=1:n r(n,m)=m*n; end end (2) while 语句 格式: while 表达式 语句体 end 其中:当表达式值为真时执行语句体; 当表达式值为假时终止循环. 例:找出阶乘超过10100的最小整数. n=1; while prod(1:n)1e100 n=n+1; end 选择(分支)结构: if 语句 格式: if 逻辑表达式1 语句体1 elseif 逻辑表达式2 语句体2 elseif 逻辑表达式3 语句体3 … else 语句体else end 0．３　几个MATLAB常用函数 绘图函数 plot(t,x) 绘制函数x(t)的曲线图; plot(t,x,’.’) 绘制函数x(t)的散点图(点图还有‘o’,‘*’等); stem(t,x) 绘制函数x(t)的火柴杆图(空心); stem(t,x,’.’) 绘制函数x(t)的火柴杆图(实心). FFT fft(x(t),n) 傅立叶变换 Ifft(X(ω),n) 傅立叶反变换 第一章 信号的时频分析 1.1 引言 1.2 Hilbert变换与解析信号 1.3 时频分布及其性质 1.4 二次型时频分布的交叉项 1.5 Wigner-Ville分布及其应用 1.1 　引言 例1.1 用MATLAB编程形成雷克子波,并进行频谱分析. %峰值频率为fp的雷克子波的表达式: % Rc(t)=[1-2(πfp t)2]exp[- (πfp t)2] clc clear tm=6.4; dt=0.1; t=-tm+dt:dt:tm; n=length(t); fp=0.5;