数字图像上机二save.docVIP

数字信号与图像处理第二次上机作业 班级：071211班 学号姓名：袁玉凡 典型信号序列的频谱 脉冲序列 利用1行32列全零矩阵 (2) 单位阶跃函数 1到32的单位阶跃函数： （3）矩形序列 矩形序列可以用单位阶跃函数表达，单位矩形信号的脉高为1。这里选-10到40区间上的矩阵序列，脉宽区间为0到10，并作出矩形序列的幅度谱和矩形序列的相位谱。 (4) 指数函数 实指数序列: 参数为：n=0:20 ;a=1.2; x=a.^n; 下面进行周期信号的一个周期和多个周期进行频谱分析，参数如下： fs=10;? N=60;? n=1:N;? t=n/fs;? x?=0.25*sin(2*pi*t)+0.5*sin(2*pi*2*t)+0.25*sin(2*pi*4*t);? 根据算法编出程序运行结果如下： FIR设计（窗函数） 各种窗函数的基本参数：? 窗函数 旁瓣峰值幅度/dB 过渡带宽 阻带最小衰减/dB 矩形窗 -13 4π/N -12 三角形窗 -25 8π/N -25 凯塞窗(α=7.865) -57 10π/N -80 ??? 用窗函数法设计FIR数字滤波器?的基本原理为：窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近.由于往往是无限长序列,而且是非因果的,所以用窗函数将截断,并进行加权处理,得到的函数就作为实际设计的FIR数字滤波器的单位脉冲响应序列。 本题选用凯塞窗，用参数Rs=0.01; 抽样频率fs=8000;得到滤波器的不同频率下的分贝，即增益效应。再选择频率f1=500;f2=1500;f3=2000;f4=3000; 做出了信号滤波前时域图和信号滤波后时域图，更加直观的对比出滤波器的滤波效果。 采用 [M,Wc,beta,ftype]= kaiserord(fcuts,a,dev)函数来估计滤波器阶数等，得到凯塞窗滤波器。这里的函数kaiserord(fcuts,a,dev,fs): (1) fcuts为对应的频率，fs为采样频率；当f用数字频率表示时，fs则不需要写。 (2)a=[1 0]为由f指定的各个频带上的幅值向量，一般只有0和1表示；a和f长度关系为（2*a的长度）- 2=（f的长度） (3)devs=[0.05 10^(-2.5)]用于指定各个频带输出滤波器的频率响应与其期望幅值之间的最大输出误差或偏差，长度与a相等，计算公式： 阻带衰减误差=αs，通带衰减误差=αp，可由滤波器指标得到。 (4) fs缺省为2Hz。 音频信号实验 要求将两段音频混合： Step 1 输入音频信号，将WAV文件转换成变量 Step2 双声道音频信号：产生第一，第二信道。 Step3 混音（归一化） Step4 回声 Step 5 频谱：提取音乐做傅里叶变换进行频率合成。 实验总结： MATLAB是十分强大的数字信号和图像处理平台，其中已经定义好的许多函数十分方便我们的应用，如窗函数中就有矩形窗、三角窗，汉宁窗，海宁窗等多种窗口 西安电子科技大学袁玉凡 1 [键入文字] [键入文字] [键入文字]