matlab71739.docVIP

西北工业大学 《基于MATLAB的数字信号处理》 学 院： 电子信息学院 学　　号： 2009201219 姓　　名： 李 慧 敏 专 业： 信号与信息处理 实验时间： 2010年6月 指导教师： 冯 燕 西北工业大学 2010年 06 月实验一 MATLAB基本1．熟悉MATLAB；1．2． 1．， ，求使 2 .利用利用公式计算的值。 3. 在同一坐标系下绘制曲线。 四、设计思想 1. ，则； 2. ，使用循环实现 3. 使用hold on，或者使用plot的多坐标 五、程序代码及注释 1.程序： A=[0 2 1;2 -1 3;-3 3 -4]; B=[1 2 3;2 -3 1]; x=B/A 2.程序： pai=1; for i=1:12 pai=pai+((-1)^i)/(1+2*i); end pai=pai*4 3.程序： t=[0:0.01:2*pi]; plot(t,t.^2,r,t,-t.^2,g-.,t,(t.^2).*cos(t),b--) 六、实验结果 1. 结果为。 2. 得 3. 得图 实验二 基于MATLAB信号处理实验 一、实验目的及要求 1. 回顾数字信号处理的主要内容； 2. 掌握利用MATLAB进行信号处理的方法； 3. 了解信号处理工具箱中一些函数的功能； 二、实验设备（环境）及要求 1．2． 1．，其中，，为白噪声。求信号的DFT及其反DFT变换，画出原始信号、频谱特性和反变换后信号的曲线。 2 . 设计FIR滤波器（低通，四种窗函数，fs＝1200Hz, 截止频率300Hz, 阶数自选, 画频率特性并分析比较)。 四、设计思想 1.先生成信号，再按采样，最后调用fft和ifft进行计算。 2.阶数选为50，计算指标，然后选4个窗函数进行设计。 五、程序代码及注释 1. 程序： fs=10;N=200;n=[1:N]; r=randn(1,N);%白噪声 x=5*sin(2*pi*2*n/fs)+3*sin(2*pi*2.05*n/fs)+r;%信号波形 a1=fft(x);%DFT a2=ifft(a1);%反DFT subplot(3,1,1);plot(n/fs,x);title(原始信号); subplot(3,1,2);stem(n,abs(a1));title(信号频谱); subplot(3,1,3);plot(n/fs,a2);title(反变换后信号); 2. 程序： upper_stop=300/600; n=50;%50阶 figure(1); window=boxcar(n+1);%窗函数 b1=fir1(n,upper_stop,window); freqz(b1);title(矩形窗); figure(2); window=triang(n+1);%窗函数 b1=fir1(n,upper_stop,window); freqz(b1);title(三角窗); figure(3); window=hanning(n+1);%窗函数 b1=fir1(n,upper_stop,window); freqz(b1);title(海宁窗); figure(4); window=bartlett(n+1);%窗函数 b1=fir1(n,upper_stop,window); freqz(b1);title(bartlett窗); 六、实验结果 1.结果如下图所示 2.用四种窗函数的设计结果如下 实验三 基于MATLAB的图像处理实验 一、实验目的及要求 1．；1．2． 1． 2．对一幅图像进行DCT变换，分别保留最大的总变换系数个数的1/4、1/8和1/16个变换系数，进行反DCT得到重构图像，显示原图像和重构图像并计算重构图像的峰值信噪比。 四、设计思想 1.选取一张jpg格式的图片，用rgb2gray函数将之转换为灰度图像，然后分别添加椒盐噪声和高斯白噪声，并分别选取3×3、5×5、7×7的模板进行均值和中值滤波。 2. 选取一张jpg格式的图片，用rgb2gray函数将之转换为灰度图像。进行DCT变换，对所得矩阵进行拉直，再升序排序，取出处于总数的1/4、1/8和1/16处的元素值，然后分别将A中小于这三个数值的灰度值置零并进行反DCT变换，得到恢复图像。峰值信噪比（PSNR）只需按公式计算即可。 五、程序代码及注释 1.程序： I=imread(jord