matlab实验指导报告,Matlab实验指导书详解.doc

Matlab实验指导书详解.doc

(34页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

19.90 积分

34实验一：Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境。2．学习使用图形函数计算器命令funtool及其环境。二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。Ø 求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)^3*(2-x);Ø 求下列函数的符号积分(1) y=cos(x);(2) y=1/(1+x^2);(3) y=1/sqrt(1-x^2);(4) y=(x-1)/(x+1)/(x+2);Ø 求反函数(1) y=(x-1)/(2*x+3);(2) y=exp(x);(3) y=log(x+sqrt(1+x^2));Ø 代数式的化简(1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2) sin(x)^2+cos(x)^2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。Ø 从y=x^2通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2(4) y4=x^2+2(5) y5=x^4(6) y6=x^2/23．两个函数之间的操作Ø 求和(1) sin(x)+cos(x)(2) 1+x+x^2+x^3+x^4+x^5Ø 乘积(1) exp(-x)*sin(x)(2) sin(x)*xØ 商(1) sin(x)/cos(x);(2) x/(1+x^2);(3) 1/(x-1)/(x-2);Ø 求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设计提示1．初次接触Matlab应该注意函数表达式的文本式描述。2．在使用图形函数计算器funtool时，注意观察1号和2号窗口中函数的图形。四、实验报告要求1．针对图形函数计算器funtool，对每一类型计算记录其中一个图形的曲线。2．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。《Matlab语言及其在电子信息科学中的应用》实验指导书实验二：M文件和Mat文件操作一、实验目的1．定制自己的工作环境。2．编写简单的M文件。3．保存内存工作区中的变量到.mat文件。4．学会只用Matlab帮助。二、实验内容1．使用format命令和File|Peferences菜单定制自己的工作环境。2．编写如下M文件，试调整参数a的大小，观察并记录y1、y2的波形特征。%example1.mt=0:pi/100:4*pi;a=3;y2=exp(-t/a);y1=y2.*sin(a*t);plot(t,y1,'-r',t,y2,':b',t,-y2,':b');3．保存内存工作区变量a、t、y1、y2到example1.mat文件；关闭Matlab，再重新启动；观察内存工作区；重新根据.mat文件恢复原来的工作区变量。4．在命令窗口中查看exp函数的帮助；运行helpwin查看超文本格式的帮助文件，试翻译并记录下信号处理工具箱(Signal Processing Toolbox)中的函数分类(Functions -- Categorical List)。三、设计提示1．可以用命令语句、菜单或按钮等多种方式执行命令。2．用于编辑M文件的文本编辑器还可以执行和调试程序。3．不同的工具箱可能包含同名的函数，查看帮助时应注意在左侧栏选择相应的工具箱类别。四、实验报告要求1．对实验内容2，说明参数a的大小对y1、y2波形特征的影响。2．翻译命令窗口中的exp函数的帮助信息。3．运行helpwin，试翻译并记录下信号处理工具箱(Signal Processing Toolbox)中的函数分类(Functions -- Categorical List)。4．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。实验三：矩阵运算与元素群运算一、实验目的1．掌握数组与矩阵的创建。2．掌握矩阵运算与数组运算。3．掌握基本元素群运算。4．掌握向量与矩阵的特殊处理。二、实验内容1．“：”号的用法。用“：”号生成行向量a=[1 2 3 4 5 6 7 8 9 10]、b=[5 3 1 -1 -3 -5]；2．用线性等分命令linspace重新生成上述的a和b向量。3．在100和10000之间用对数等分命令logspace生成10维的向量c。4．生成范围在[0，10]、均值为5的3×5维的均匀分布随机数矩阵D。5．利用magic函数生成5×5维的魔方矩阵，取其对角向量e，并根据向量e生成一个对角矩阵E。(所谓魔方矩阵就是各行、各列、各对角线元素之和相等。)6．另AA是3×3维魔方矩阵，BB是由A旋转180°得到。CC是一个复数矩阵，其实部为AA，虚部为BB。DD是CC的转置，EE是CC的共轭。分别计算CC和EE的模和幅角。7．ｆ是一个首项为20，公比为0.5的10维等比数列；g是一个首项为1，公差为3的10维等差数列。试计算向量f和g的内积s。8．生成一个9×9维的魔方矩阵，提取其中心的3×3维子矩阵M，利用sum函数检验其各行和各列的和是否相等。9．已知，利用函数生成左上三角矩阵。三、设计提示1．等比数列可利用首项和公比的元素群的幂乘积生成。2．提取子矩阵，可灵活应用“：”号或空阵[ ]。3．尽量用Matlab函数生成上述矩阵或向量，不要用手工逐个输入。四、实验报告要求1．编写实现第二节实验内容中所使用的函数命令，并记录相应的生成结果。2．思考题：是否存在2×2维的魔方矩阵？。3．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五，参考答案1.代码：a=1:10 b=5:-2:-5 a = 1 2 3 4 5 6 7 8 9 10b = 5 3 1 -1 -3 -51. 代码:linspace(1,10,10)linspace(5,-5,6)ans = 1 2 3 4 5 6 7 8 9 10。省略部分。 5 2];a=[1 -1.6 1.3 -.9 .5];freqz(b,a,256,8000) % Fs=8000结果：第6题：k=0.5;z=[1,-3]';p=[-1,-2,-4]';sys3=zpk(z,p,k);impulse(sys3)%频率响应结果：b=[3 5 2];a=[1 -1.6 1.3 -.9 .5];impz(b,a,60)实验九：数字滤波器设计与语音信号处理一、实验目的1．掌握Butterworth滤波器、Chebyshev Ⅰ型、Ⅱ型滤波器的设计，并能用于对语音信号子频带的提取。2．了解FIR滤波器的设计。二、实验内容(边做实验，边将生成结果和图形拷贝到Word文档中)1．将Test1.wav和Test2.wav两个语音文件复制到C:\Matlab6p5\work\目录下。2．设计一个4阶的Chebyshev Ⅱ型模拟高通滤波器，要求截止频率，阻带纹波不超过-30dB。取采样频率。要求结果给出该滤波器的传递函数表达式，并画出该滤波器的频率响应曲线图。3．设计一个6阶的Butterworth数字带通滤波器，要求截止频率，，取采样频率。要求结果给出该滤波器的传递函数表达式，并画出该滤波器的频率响应曲线图。4．设计一个6阶的Butterworth数字低通滤波器，从一段含噪声语音信号Test2.wav(0～4kHz，Fs=8kHz)中提取出500Hz以下的子带信号。在一个图形窗口中分上下两个子窗口分别显示原始语音信号波形和滤波后的500Hz以下子带信号波形。5．设计一个65阶的FIR数字滤波器，通带(归一化)：0～0.45、0.65～0.85。画出该滤波器的频率响应曲线图。三、设计提示1．设计滤波器时，注意频率的归一化。2．Chebyshev Ⅱ型模拟高通滤波器设计，用函数cheby2(n,Rs,Wn,'high','s')。其中Rs=30(即表示：阻带纹波不超过-30dB)。3．设计IIR带通或带阻滤波器时，若给滤波器函数的参数为n，则设计出来的是2n阶的滤波器。4．设计高通滤波器，用'high'关键字；带通用'bandpass'关键字。5．查看模拟滤波器的频率响应曲线用freqs函数，查看数字滤波器的频率响应曲线用freqz函数。6．对于FIR滤波器的设计，第一通带如果从0开始，则应加关键字'DC-1'。四、实验报告要求1．编写实现第二节实验内容所使用的M文件，并在Word文档中记录相应的生成结果和图形。2．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五、参考答案第2题：fc=1000;Fs=10000;Wn=fc/(Fs/2); %归一化n=4;Rs=30;[b,a]=cheby2(n,Rs,Wn,'high','s');sys=tf(b,a)freqs(b,a);结果：s^4 - 1.42e-017 s^3 + 0.04 s^2 - 2.938e-018 s + 0.0002------------------------------------------------------ s^4 + 0.6442 s^3 + 0.2475 s^2 + 0.05581 s + 0.第3题：fc1=1000;fc2=3000;Fs=10000;Wn=[fc1,fc2]/(Fs/2); %归一化n=6;[b,a]=butter(n/2,Wn,'bandpass');sys=tf(b,a)freqz(b,a,512);结果：Transfer function: 0.09853 - 0.2956 z^-2 + 0.2956 z^-4 - 0.09853 z^-6---------------------------------------------------------1 - 1.366 z^-1 + 1.245 z^-2 - 0.8777 z^-3 + 0.6537 z^-4 - 0.2256 z^-5 + 0.0563 z^-6第4题：[x,Fs,Bit]=wavread('Test2.wav');Wn=500/(Fs/2); %归一化n=6;[b,a]=butter(n,Wn);y=filter(b,a,x);subplot(211); plot(x)title('原始语音信号');subplot(212); plot(y)title('500Hz以下子带信号');结果：第5题：n=65;Wn=[0.45,0.65,0.85];b=fir1(n,Wn,'DC-1');a=1;sys=tf(b,a)freqz(b,a,512);结果：实验十：数字图像处理初步一、实验目的1．掌握图像文件的读取与写入。2．掌握图像文件的显示、翻转、缩放、叠加3．掌握黑白和彩色图像在Matlab中矩阵存储格式。二、实验内容(边做实验，边将生成结果和图形拷贝到Word文档中)1．将Test3.jpg和Test4.jpg两个图像文件复制到C:\Matlab6p5\work\目录下。2．读取Test3.jpg到矩阵A中，观察A的维数，并显示该图像。3．对A做一定的处理和变换，使之成为下面的图像，并存储到文件PicOut.bmp中。4．读取Test4.jpg到矩阵X中，去掉红色图层，然后显示出来。三、设计提示1．首先计算A的维数，然后列数放大1倍，再考虑将矩阵左半边的图像数据镜像到右半边去。2．对于彩色图像，矩阵的维数是：行数×列数×图层数(1-红，2-绿，3-蓝)。四、实验报告要求1．编写实现第二节实验内容所使用的M文件，并记录相应的生成结果。2．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。五、参考答案第3题：A=imread('Test3.jpg');figure(1);imshow(A);[x,y,z]=size(A);B=A; %复制左边图像B(:,[y+1:2*y],1)=fliplr(A(:,:,1)); %扩维并复制右边图像的红色图层(第一层)B(:,[y+1:2*y],2)=fliplr(A(:,:,2)); %扩维并复制右边图像的红色图层(第二层)B(:,[y+1:2*y],3)=fliplr(A(:,:,3)); %扩维并复制右边图像的红色图层(第三层)figure(2);imshow(B);imwrite(B,'PicOut.bmp');结果：(略)第4题：X=imread('Test4.jpg');figure(1);imshow(X);Y=X;Y(:,:,1)=0;figure(2);imshow(Y);结果：(略) 关 键 词： matlab 实验 指导 详解

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。