实验一 Matlab软件的使用..

合集下载

Matlab实验第一次实验答案

Matlab实验第一次实验答案

实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法;2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容:1、帮助命令使用help命令,查找sqrt(开方)函数的使用方法;解:sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.RemarksSee sqrtm for the matrix square root.Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.41422、矩阵运算(1)矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B解:A=[1 2;3 4 ];B=[5 5;7 8 ];A^2*B(2)矩阵除法已知A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/B解:A=[1 2 3;4 5 6;7 8 9 ];B=[1 0 0;0 2 0;0 0 3 ];A\B,A/B(3)矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'解:A=[5+1i,2-1i,1;6*1i,4,9-1i ];A1=A.',A2=A'(4)使用冒号表达式选出指定元素已知:A=[1 2 3;4 5 6;7 8 9];求A 中第3列前2个元素;A 中所有列第2,3行的元素;方括号[]解:A=[1 2 3;4 5 6;7 8 9 ];B=A([1,2],[3]),C=A(2:end, : )用magic 函数生成一个4阶魔术矩阵,删除该矩阵的第四列3、多项式(1)求多项式42)(3--=x x x p 的根解:A=[1 0 -2 -4];B=roots(A)(2)已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ,求矩阵A的特征多项式;解:A=[1.2 3 5 .9 ; 5 1.7 5 6 ;3 9 0 1 ;1 2 3 4];A=poly(A);A=poly2sym(A)把矩阵A作为未知数代入到多项式中;4、基本绘图命令(1)绘制余弦曲线y=cos(t),t∈[0,2π]解:t=0:.1:2*pi;y=cos(t);plot(t,y),grid(2)在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5),t∈[0,2π]t=0:.1:2*pi;y1=cos(t-0.25);y2=sin(t-0.25);plot(t,y1,t,y2)grid5、基本绘图控制绘制[0,4π]区间上的x1=10sint曲线,并要求:(1)线形为点划线、颜色为红色、数据点标记为加号;(2)坐标轴控制:显示范围、刻度线、比例、网络线(3)标注控制:坐标轴名称、标题、相应文本;t=0:.01:4*pi;y=10*sin(t);plot(t,y,'-.',t,y,'r')grid6、基本程序设计(1)编写命令文件:计算1+2+…+n<2000时的最大n值;解法1:s=0;i=0;while(s<2000) i=i+1;s=s+i; ends=s-i,i=i-1解法2:s=0;for i=1:1000; s=s+i;if(s>2000) ,break;endends=s-i,i=i-1(2)编写函数文件:分别用for和while循环结构编写程序,求2的0到n次幂的和。

实验1Matlab的基本使用方法

实验1Matlab的基本使用方法
最简单的方法是将用户目录设置为当前目录,设置方法是在命令窗口中输入指令 cd d:\mydir 或者 pwd d:\mydir 另一种方法是将用户目录设置为 Matlab 的搜索路径中,方法是在命令窗口中输入指令:
7
path(path,’d:\mydir’)。用这种方法设置的目录信息只在当前环境中有效,一旦 Matlab 重新启 动,以前的设置就无效。
等于)
1
2) 逻辑运算:&(与)、|(或)、~(非)(非零元素都表示为真) 关系函数和逻辑函数:
函数
功能
函数
功能
all(x)
检查 x 是否全为 1
any(x)
检查 x 是否有不为零的元素
find(x)
找出非零元素的位置标识
isempty(x) 检查 x 是否为空阵
isinf(x)
检查 x 是否为无穷大
Your variables are:
RXYZtxyz
>> whos
Name
Size
Bytes Class
R
31x31
7688 double array
X
31x31
7688 double array
Y
31x31
7688 double array
Z
31x31
7688 double array
t
1x501
(11) 帮助.
3、实验指导 Matlab 是美国 Mathworks 公司推出的一套高性能的数值分析和计算软件,它将矩阵运
算、数值分析、图形处理、编程技术结合在一起,为用户提供了一个强有力的科学及工程问 题分析计算和程序设计的工具。
Matlab 语言可以被认为是一种解释性语言,用户可以在 Matlab 的工作空间中输入一个 指令,也可以在编辑器中编写应用程序,应用程序执行时,Matlab 软件对其中的命令和函 数进行翻译,然后在 Matlab 环境中对它进行处理,最后返回结果。 (1)数值、变量和表达式

实验1_Matlab操作入门

实验1_Matlab操作入门

实验一 Matlab操作入门实验目的:1.了解Matlab的工作环境。

掌握help命令的使用。

了解general函数库中有哪些常用的函数,掌握who, whos, clear, save, load 等命令的使用。

2.掌握矩阵的输入方法;掌握利用elmat函数库中一些函数,如linspace, logspace, zeros, ones, eye, rand, randn, magic等创建矩阵的方法。

3.了解elmat函数库中的预定义变量。

4.掌握elmat函数库中diag, blkdiag, triu, tril, fliplr, flipud, rot90等函数的用法。

5.掌握elmat函数库中ndims, size, numel, length,disp等函数的使用。

6.掌握datafun函数库中的一些基本统计函数max, min, mean, median, sum, prod 等的基本用法。

7.掌握matfun函数库中的一些基本矩阵操作函数的使用:inv(A)—求矩阵A的逆,det(A)—求矩阵A的行列式值, rank(A)-求矩阵的秩,trace(A)-求矩阵A的迹(即主对角线元素之和),norm(A)-求向量或矩阵的范数,cond(A)-求矩阵A的条件数,etc。

实验内容:一、Matlab环境1.启动Matlab,了解Matlab的主要窗口:命令窗口,当前目录窗口,工作空间窗口,历史命令窗口。

2.输入help, 了解你的计算机上安装的Matlab中,有哪些帮助主题(函数库、工具箱)。

3.输入 help general(通用命令函数库), 了解Matlab中有哪些通用的命令。

4.输入 help syntax,阅读具体的内容,了解Matlab语句或命令的格式。

5.输入 help format, 了解format命令的使用格式。

6.输入help timefun(时间、日期函数库), 了解有哪些时间和日期函数。

MATLAB实验报告

MATLAB实验报告

实验一、Matlab 软件的安装与基本使用一、实验目的和要求:学会MA TLAB 软件的安装方法,熟悉MA TLAB 的工作环境及其基本功能的设置,以及基本指令的使用,了解MA TLAB 数值数组及其矩阵操作的使用、复数的输入方法。

二、实验内容和步骤:1、完成MA TLAB 软件的安装。

2、安装完成后启动MA TLAB 软件,熟悉MA TLAB 的工作环境。

打开MA TLAB 7,默认打开的窗口包括: (1)命令窗口(Command Window); (2)命令历史窗口(Command History); (3)工作间管理窗口(Workspace);(4)当前路径窗口(Current Directory)此外,还有捷径键、编译窗口、图形窗口和帮助窗口等其他种类的窗口,结合课本实例,熟悉每个窗口的功能和基本使用。

3、(1)求23)]47(212[÷-⨯+的运算结果,注意命令后有无分号的作用。

命令及结果如下:t=[12+2*(7-4)]/3^2t = 2(2)求20012x x v t at =++的值,其中t=5s 、010x m =、015/v m s =、29.81/sec a m =-。

命令及结果如下:syms t x0 v0 a; t=5,x0=10,v0=15,a=-9.18; x=x0+v0*t+1/2*a*t^2 t = 5 x0 = 10 v0 = 15 x = -29.7500(3) 设3x =, 4v =,用MA TLAB 命令求下列各式的值:1)22log()x v + 命令及结果如下:syms x v x=3,v=4; log(x^2+v^2)x = 3 ans =3.21892)2(2)x v - 命令及结果如下:syms x vx=3,v=4;f=(sqrt(x-3))/(x-2*v)^2 x = 3 f = 0 3) |sin 2|vx e命令及结果如下: syms x v x=3,v=4;t=abs(2*x)*exp(v) x = 3t = 327.58894、用三种不同的方法输入矩阵1234246836912a ⎛⎫ ⎪= ⎪ ⎪⎝⎭,运行并观察实验结果,矩阵的三种基本输入方法如下:a=[1,2,3,4;3,4,6,8;3,6,9,12] a =1 2 3 4 3 4 6 8 3 6 9 12>> a=[1,2,3,43,4,6,8 3,6,9,12]a =1 2 3 4 3 4 6 8 3 6 9 12>> >>利用M 文件创建 1 2 3 4 3 4 6 8 3 6 9125、输入复数矩阵⎥⎦⎤⎢⎣⎡-+++=i ii i24356223B ,并求其实部、虚部、模和幅角(分别用弧度和角度表示),了解MA TLAB 向量化运算的特点。

信息论编码实验一MATLAB的基本操作

信息论编码实验一MATLAB的基本操作

实验一 MATLAB 的基本操作一、实验目的1、掌握Matlab 软件使用的基本方法;2、熟悉Matlab 的数据表示、基本运算方法;3、熟悉Matlab 绘图命令及基本绘图控制。

二、实验仪器与软件1、PC 机 1台2、MATLAB7.0环境三、实验原理MATLAB 环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。

MATLAB 有3种窗口,即:命令窗口(The Command Window )、m-文件编辑窗口(The Edit Window )和图形窗口(The Figure Window ),而Simulink 另外又有Simulink 模型编辑窗口。

1、命令窗口(The Command Window )当MATLAB 启动后,出现的最大的窗口就是命令窗口。

用户可以在提示符“>>”后面输入交互的命令,这些命令就立即被执行。

在MATLAB 中,一连串命令可以放置在一个文件中,不必把它们直接在命令窗口内输入。

在命令窗口中输入该文件名,这一连串命令就被执行了。

因为这样的文件都是以“.m ”为后缀,所以称为m-文件。

2、m-文件编辑窗口(The Edit Window )我们可以用m-文件编辑窗口来产生新的m-文件,或者编辑已经存在的m-文件。

在MATLAB 主界面上选择菜单“File/New/M-file ”就打开了一个新的m-文件编辑窗口;选择菜单“File/Open ”就可以打开一个已经存在的m-文件,并且可以在这个窗口中编辑这个m-文件。

如,计算3[(12)34)]2+⨯-÷,只需在提示符“>>”后输入“((1+2)*3-4)/2^3”,然后按Enter 键。

该命令行涉及加、减、乘、除及幂运算符,MATLAB 运算的执行次序遵循的优先规则为:从左到右执行;幂运算具有最高的优先级,乘法和除法具有相同的次优先级,加法和减法有相同的最低优先级;使用括号可以改变前述优先次序,并由最内层括号向外执行。

MATLAB实验代码

MATLAB实验代码

实验一MATLAB软件使用初步一、实验目的∙学习MATLAB软件的安装过程,熟悉MATLAB软件界面的组成及基本使用方法。

∙理解数组(Array)的分类,及标量(scalar)、矢量(vector)和矩阵(Matrix)的区别,熟悉数组与矩阵的构造方法,掌握数组与矩阵的基本运算法则。

二、实验要求1、掌握MATLAB软件的启动与退出的方法。

2、掌握MATLAB软件的桌面环境。

3、掌握MATLAB软件的基本使用方法。

4、掌握M文件编辑器的使用。

5、学会使用MATLAB的帮助系统。

6、掌握标量、矢量和矩阵的构造方法。

7、掌握数组的四则运算。

8、掌握矩阵的四则运算。

9、掌握基本的数据输入和输出函数。

三、实验内容1.通过桌面快捷图标或“开始—程序”菜单运行MATLAB,熟悉MATLAB的桌面环境,如图1所示。

图1 MATLAB桌面环境2.浏览各个菜单项的内容,试着改变各个窗口的字体与字号。

3.熟悉matlab 中的特殊变量:i 、j 、pi 、ans 等等。

4.在实验报告中回答下面哪些赋值语句是正确的,如果不正确给出理由。

a=1正确a20=1.2;正确 年薪=100000错误不能是汉字_fee=20.0错误开头不能下划线coordinate.x=1正确5.已知,,用MATLAB 分别执行下列语句。

并在实验报告中记录语句和1u =3v =结果。

a. b. c. d. 43u v ()222v u v -+333v v u -243v πu=1,v=3 (4*u)/(3*v) (2*v^-2)/(u+v)^2 (v^3)/(u^3-v^3)(4**pi*v^2)/3ans=0.4444 0.2222 -1.0385 37.69916. 运行以下命令并将运行结果与工作空间浏览器中的状态进行比较。

>> a=sin(pi); b=[1, 2]; c=[1; 2];>> who 所有可用变量 >> whos 显示变量所占的矩阵大小和所占的字节在实验报告中回答Size 下的1x1、1x2、2x1的含义>> whos Name Size Bytes Classa 1x1 8 double arrayb 1x2 16 double arrayc 16 double array矩阵大小7.通过工作目录浏览器(Current Directory Browser )或者“File”菜单中相应菜单项改变当前工作目录。

matlab实验

matlab实验

实验一 MATLAB基本操作一、实验目的1、了解MATLAB应用程序环境2、掌握MATLAB语言程序的书写格式和MATLAB语言程序的结构。

3、掌握在MATLAB应用环境下编写程序4、掌握MATALB调试过程,帮助文件5、掌握MATLAB语言上机步骤,了解运行一个MATLAB程序的方法。

6、本实验可在学习完教材第一章后进行。

二、主要仪器及耗材PC电脑,MATLAB6.5软件三、实验内容和步骤1、MATLAB语言上机步骤:(1)、进入系统在C盘或其他盘上找到MATLAB或MATLAB6.5,然后双击其图标打开文件夹。

然后进行编辑源程序->编译->连接->执行程序->显示结果(2)、常用命令编辑切换(F6),编译(F9),运行(CTRL+F9),显示结果(ALT+F5)其它常用命令见“附录一”。

2、有下面的MATLAB程序。

(1)数值计算功能:如,求方程 3x4+7x3 +9x2-23=0的全部根p=[3,7,9,0,-23]; %建立多项式系数向量x=roots(p) %求根(2)绘图功能:如,绘制正弦曲线和余弦曲线x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));(3)仿真功能:如,请调试上述程序。

3、熟悉MATLAB环境下的编辑命令,具体见附录一。

三、实验步骤1、静态地检查上述程序,改正程序中的错误。

2、在编辑状态下照原样键入上述程序。

3、编译并运行上述程序,记下所给出的出错信息。

4、按照事先静态检查后所改正的情况,进行纠错。

5、再编译执行纠错后的程序。

如还有错误,再编辑改正,直到不出现语法错误为止。

四、实验注意事项1、记下在调试过程中所发现的错误、系统给出的出错信息和对策。

分析讨论对策成功或失败的原因。

2、总结MATLAB程序的结构和书写规则。

五、思考题1、matlab到底有多少功能?2、MATLAB的搜索路径3、掌握使用MATLAB帮助文件实验二 MATLAB 矩阵及其运算一、 实验目的1、了解矩阵的操作,包括矩阵的建立、矩阵的拆分、矩阵分析等2、了解MATLAB 运算,包括算术运算、关系运算、逻辑运算等3、掌握字符串的操作,了解结构数据和单元数据。

实验1 Matlab软件环境的基本使用

实验1 Matlab软件环境的基本使用

实验1 Matlab 的基本运算(1)实验目的:1、 了解Matlab 的发展和主要功能;2、 熟悉Matlab 工作环境的各个窗口;3、 掌握建立矩阵的方法;4、 掌握Matlab 各种表达式的书写规则以及常用函数的使用。

实验内容:1. 在命令窗口依次输入下列命令,根据执行结果分析其功能:helplookfor invhelp invwhich inv2. 建立自己的工作目录,再将自己的工作目录设置到Matlab 搜索路径下。

用help 命令查询自己的工作目录;3. 利用Matlab 的帮助功能,分别查询inv 、plot 、max 、round 等函数的功能及用法;4. 先求下列表达式的值,然后显示Matlab 工作空间的使用情况并保存全部变量。

(1)sin z e =+122851 (2)ln(z x =+212,其中.i x +⎡⎤=⎢⎥-⎣⎦2120455 (3)...sin(.)ln ,.,.,.,...,.,.,.a a e e a z a a --+=++=---03033030330292828293022提示:利用冒号表达式生成a 向量,求各点的函数值时用点乘运算。

(4),,,:.:.,t t z t t t t t t ⎧≤<⎪=-≤<=⎨⎪-+≤<⎩224201112005252123其中提示:用逻辑表达式求分段函数值。

5. 已知:A -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦12344347873657,B -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦131203327 求下列表达式的值。

(1)*A B A B I +-+6和,其中I 为单位矩阵;(2)*A B 和.A B *(3)^A 3和.^A 3(4)/A B 和\A B(5)[],A B 和()^[,],:;A B ⎡⎤⎣⎦1326. 将课件中的相关例题进行演示实验。

思考题:1. 与其他高级语言相比,Matlab 有哪些显著的特点?2. Matlab 操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立的窗口?又怎样恢复? 实验报告提交:将所有实验内容进行总结(实验内容6不必写),要求表述清晰,结果明确。

实验一_MATLAB_软件环境及基本操作

实验一_MATLAB_软件环境及基本操作

实验⼀_MATLAB_软件环境及基本操作实验⼀ MATLAB 软件环境及基本操作⼀、实验⽬的熟悉 MATLAB 软件的⼯作环境,练习 MATLAB 命令窗⼝的基本操作。

⼆、实验内容1、命令窗⼝(Command Window )运⾏⼊门(1)最简单的计算器使⽤法【例 1】求[12 + 2×(7- 4)]÷32的算术运算结果。

1)⽤键盘在 MATLAB 命令窗⼝中命令提⽰符>>后输⼊以下内容:(12+2*(7-4))/3^22)在上述表达式输⼊完成后,按【Enter 】键,执⾏命令。

3)MA TLAB 命令窗⼝中将显⽰以下结果:ans =【例 2】简单矩阵的输⼊步骤。

123A=4567891)在键盘上输⼊下列内容A = [1,2,3; 4,5,6; 7,8,9]2)按【Enter 】键执⾏命令。

3)在命令执⾏后,MATLAB 命令窗⼝中将显⽰以下结果:A =1 2 34 5 67 8 9【例 3】矩阵的分⾏输⼊A=[1 2 34 5 67 8 9](显⽰结果同例2)【例 4】命令的续⾏输⼊。

…是续⾏符号,表明下⼀⾏是本⾏的接续。

S = 1 – 1/2 + 1/3 –1/4 + 1/5 – 1/6 ... + 1/7 – 1/8S =0.6345 (2)复数和复数矩阵【例1】复数z1 = 3+ 4i,z2 =1+ 2i,z3 = 2eπi/6表达,及计算123z z zz(1)直⾓坐标表⽰法z1= 3 + 4iz1 =3.0000 +4.0000i(2)采⽤运算符构成的直⾓坐标表⽰法和极坐标表⽰法>> z2 = 1 + 2 * i %运算符构成的直⾓坐标表⽰法。

z2 =1.0000 +2.0000i>> z3=2*exp(i*pi/6) %运算符构成的极坐标表⽰法z3 =1.7321 + 1.0000i>> z=z1*z2/z3z =0.3349 + 5.5801i【例-2】复数矩阵的⽣成及运算>> A=[1,3;2,4]-[5,8;6,9]*iB=[1+5i,2+6i;3+8*i,4+9*i]C=A*BA =1.0000 - 5.0000i 3.0000 - 8.0000i2.0000 - 6.0000i 4.0000 - 9.0000iB =1.0000 + 5.0000i2.0000 + 6.0000i3.0000 + 8.0000i4.0000 + 9.0000iC =1.0e+002 * 0.9900 1.1600 - 0.0900i1.1600 + 0.0900i 1.3700【例3】求上例复数矩阵C 的实部、虚部、模和相⾓。

matlab实验报告

matlab实验报告

MATLAB程序设计软件实验报告专业及班级____通信中兴131_______姓名____魏增_______________学号_____6102213869________日期_____2015.6.15_________南昌大学实验报告学生姓名: 魏增 学 号: 6102213869 班级: 中兴131班 实验类型:□ 验证 □ 综合 ■ 设计 □ 创新 实验日期: 实验成绩:实验一 MA TLAB 的基本使用一、 实验目的1.了解MA TALB 程序设计语言的基本特点,熟悉MA TLAB 软件的运行环境;2.掌握变量、函数等有关概念,掌握M 文件的创建、保存、打开的方法,初步具备将一般数学问题转化为对应计算机模型处理的能力;3.掌握二维图形绘制的方法,并能用这些方法实现计算结果的可视化。

二、 MATLAB 的基础知识通过本课程的学习,应基本掌握以下的基础知识: 一. MA TLAB 简介二. MA TLAB 的启动和退出 三. MA TLAB 使用界面简介 四. 帮助信息的获取五. MA TLAB 的数值计算功能六. 程序流程控制 七. M 文件八. 函数文件九. MATLAB 的可视化 三、上机练习1. 仔细预习第二部分内容,关于MA TLAB 的基础知识。

2. 熟悉MA TLAB 环境,将第二部分所有的例子在计算机上练习一遍3、已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=123456789,987654321B A 。

求A*B ,A .* B ,比较二者结果是否相同。

并利用MA TLAB 的内部函数求矩阵A 的大小、元素和、长度以及最大值。

解:>> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1];>> A*Bans =30 24 18 84 69 54 138 114 90 >> A.*B ans =9 16 2124 25 2421 16 9 两者结果不同 >> [m,n]=size(A) m =3 n =3 >> b=sum(A) b =12 15 18 >> a=length(A) a = 3 >>max(A)ans =7 8 94、Fibonacci 数组的元素满足Fibonacci 规则:),2,1(,12=+=++k a a a k k k ;且121==a a 。

实验一 MATLAB基本操作及运算(含实验报告)

实验一  MATLAB基本操作及运算(含实验报告)

实验一 MATLAB 基本操作及运算一、 实验目的1、 理解Matlab 数据对象的特点;2、 掌握基本Matlab 运算规则;3、 掌握Matlab 帮助的使用方法;二、 实验的设备及条件计算机一台(带有MATLAB7.0以上的软件环境)。

三、 实验内容要求建立一个名为experiment01.m 的,把与实验内容1-7相关的实验命令都放入该文件中,题与题之间用相应注释分割。

注意对实验中出现的相关函数或变量,请使用help 或doc 查询相关帮助文档,学习函数的用法。

1、 建立以下标量:1) a=102) b=2.5×10233) c=2+3i ,(i 为虚数单位)4) d=3/2πj e ,(j 为虚数单位,这里要用到exp ,pi )2、 建立以下向量:1) aVec=[3.14 15 9 26]2) bVec=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡18228871.2 3) cVec=[5 4.8 … -4.8 -5 ] (向量中的数值从5到-5,步长为-0.2)4) dVec=[100 100.01 … 100.99 101] (产生1到10之间的等对数间隔向量,参考logspace ,注意向量的长度)3、 建立以下矩阵:1)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=2222 aMat aMat 一个9×9的矩阵,其元素全为2;(参考ones 或zeros )2)⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=1000005000001 bMat bMat 是一个9×9的矩阵,除主对角上的元素为[1 2 3 4 5 4 3 2 1]外,其余元素均为0。

(参考diag )。

3)10020109212291111=cMatcMat 为一个10×10的矩阵,可有1:100的向量来产生(参考reshape )4)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=NaN NaN NaN NaN NaN NaN NaN NaNNaN NaN NaN NaNdMatdMat 为3×4的NaN 矩阵,(参考nan )5)⎥⎦⎤⎢⎣⎡---=8710225113eMat 6)产生一个5×3随机整数矩阵fMat ,其值的范围在-3到3之间。

实验一 MATLAB的基本操作

实验一 MATLAB的基本操作

实验一 MATLAB 的基本操作一、实验目的:(1)学会MA TLAB 的基本操作。

(2)掌握MA TLAB 的基本命令。

二、实验环境:安装MATLAB7.0以上版本的计算机三、实验内容:1.文件的建立、保存与运行步骤一:打开MATLAB 软件,建立m 文件(File —New —M-file ),并保存命名为:test.m 。

步骤二:在test.m 文件中输入以下命令:0:1:10;%t 2*1;plot(t,y);%figure;%stem(t,y);%t y t ==+对设定取值连续信号作图打开新的图形窗口离散信号作图步骤三:保存(File--save ),运行程序(debug —run 或者直接按F5)。

观察输出图形。

步骤四:将程序第二句改为:2*1y t =+重复步骤三,观察出现的变化。

步骤五:去掉程序命令:figure ,观察出现的变化。

注意:文件的取名规范;分号(;)的作用;百分号(%)的作用。

2.程序的调试步骤一:观察主程序窗口中的三个部分:Workspace ,Command History ,Command Window 。

在Command Window 窗口的》符号后输入命令:x=2;观察各窗口中内容的变化。

步骤二:建立新的m 文件,命名为test2.m ,保存。

输入以下命令:(1,5);(6)clearclcx zeros y x ==步骤三:保存、运行该程序,观察Command Window 和Workspace 。

Command Window 显示出现错误,双击直接跳转到出现差错的程序行,修改程序为:y=x(4)。

保存、运行程序,再次观察Command Window 和Workspace 。

注意:程序的调试,是MATLAB 学习的必需内容。

灵活运用单步调试、断点设置等调试方法,会对以后的编程给予极大的方便。

3.help 命令的应用步骤一:在在Command Window 窗口的》符号后输入命令:help plot 步骤二:理解窗口显示内容。

实验1 MATLAB基本操作

实验1 MATLAB基本操作

一、实验目的1.熟悉MATLAB 的使用方法及特点;学会建立MATLAB 搜索路径;熟悉MATLAB 工作空间、MATLAB 集成环境、命令窗口;掌握MATLAB 的通用命令、管理命令和函数、管理变量和工作空间的使用方法;2.掌握MATLAB 基本操作及矩阵基础知识,包括:输入矩阵、矩阵的转置、矩阵元素求和、矩阵下际、矩阵连接、矩阵行列删除、矩阵产生和操作、逻辑和关系运算、操作符和特殊字符、基本矩阵和矩阵操作、基本矩阵和阵列;掌握特殊变量和常数;掌握基本数学函数。

二、实验内容1.利用基本矩阵产生 3x3 和15x8 的单位阵,全1 阵,全0 阵,均匀分布的随机阵([-1,1]之间),正态分布随机阵(方差4,均值1)eye(3),eye(15,8)ones(3),ones(15,8)zeros(3),zeros(15,8)1-2*rand(3),1-2*rand(15,8)2*randn(3)+1, 2*randn(15,8)+12.利用diag()函数和rot90()产生下列矩阵:0 0 8 2 0 4a= 0 -7 5 b= 0 5 02 3 0 7 0 8然后求解a 阵的逆矩阵aa 及b 阵的特征值和对应特征向量,并利用reshape 将aa 阵变换成行向量。

X=diag([2 -7 8]),Y=diag([3 5],-1),a=rot90(X+Y)M=diag([2 5 8]),N=diag([7 0 4]),b=M+rot90(N)aa=inv(a) reshape(aa,1,9) d=eig(b) [V,D]=eig(b)aa =-0.1339 0.2143 0.50000.0893 -0.1429 00.1250 0 0ans = -0.1339 0.0893 0.1250 0.2143 -0.1429 0 0.5000 0 0d =-1.082811.08285.0000V = -0.7921 -0.4030 00 0 1.00000.6104 -0.9152 0D = -1.0828 0 00 11.0828 00 0 5.00003产生一均匀分布在(-5,5)随机阵(50x2),精确到小数点后一位x=5-round(100*rand(50,2))/104.编程实现当α∈[-π,π],间隔为1o 时,求解正弦和余弦的值,并利用plot()函数绘制正弦,余弦曲线。

matlab 实验报告

matlab 实验报告

matlab 实验报告Matlab实验报告引言:Matlab是一种强大的数值计算和可视化软件,广泛应用于科学、工程和经济等领域。

本实验报告将介绍我在使用Matlab进行实验过程中的一些经验和结果。

实验一:矩阵运算在这个实验中,我使用Matlab进行了矩阵运算。

首先,我创建了一个3x3的矩阵A和一个3x1的矩阵B,并进行了矩阵相乘运算。

通过Matlab的矩阵乘法运算符*,我得到了一个3x1的结果矩阵C。

接着,我对矩阵C进行了转置操作,得到了一个1x3的矩阵D。

最后,我计算了矩阵C和矩阵D的点积,并将结果输出。

实验二:数据可视化在这个实验中,我使用Matlab进行了数据可视化。

我选择了一组实验数据,包括时间和温度两个变量。

首先,我将数据存储在一个矩阵中,并使用Matlab的plot函数将时间和温度之间的关系绘制成曲线图。

接着,我使用Matlab的xlabel、ylabel和title函数添加了横轴、纵轴和标题。

最后,我使用Matlab的legend函数添加了图例,以便更好地理解图表。

实验三:数值积分在这个实验中,我使用Matlab进行了数值积分。

我选择了一个函数f(x)进行积分计算。

首先,我使用Matlab的syms函数定义了符号变量x,并定义了函数f(x)。

接着,我使用Matlab的int函数对函数f(x)进行积分计算,并将结果输出。

为了验证结果的准确性,我还使用了Matlab的diff函数对积分结果进行了求导操作,并与原函数f(x)进行了比较。

实验四:信号处理在这个实验中,我使用Matlab进行了信号处理。

我选择了一个音频文件,并使用Matlab的audioread函数读取了该文件。

接着,我使用Matlab的fft函数对音频信号进行了傅里叶变换,并将结果绘制成频谱图。

为了进一步分析信号的特征,我还使用了Matlab的spectrogram函数绘制了信号的时频图。

通过对信号的频谱和时频图的观察,我可以更好地理解信号的频率和时域特性。

MATLAB 语言及其用 实验(答案)

MATLAB 语言及其用 实验(答案)

《MATLAB 语言及其用》实验指导书目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法;2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容:1、帮助命令使用help命令,查找 sqrt(开方)函数的使用方法;在 CommandWindowL里输入help,接在在search里输入sqr即可。

sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.TipsSee sqrtm for the matrix square root. Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.4142See Alsonthroot | realsqrt | sqrtm2、矩阵运算(1)矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*BA=[1 2;3 4];B=[5 5;7 8];C=A^2*B>> format compactC =105 115229 251(2)矩阵除法已知 A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/BA=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];C=A\B,D=A/BC =1.0e+016 *0.3152 -1.2609 0.9457-0.6304 2.5218 -1.89130.3152 -1.2609 0.9457D =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000(3)矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'A=[5+i,2-i,1;6*i,4,9-i];B=A.', C=A'B =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000iC =5.0000 - 1.0000i 0 -6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i(4)使用冒号表达式选出指定元素已知: A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素;A中所有列第2,3行的元素;A=[1 2 3;4 5 6;7 8 9];B1=A([1,2],[3])B2=A([2,3],:)B1 =36B2 =4 5 67 8 9方括号[]用magic函数生成一个4阶魔术矩阵,删除该矩阵的第四列A=magic(4)B=A(:,[1,2,3])或A=magic(4)A(:,4)=[]A =16 2 3 135 11 10 89 7 6 124 14 15 1B =16 2 35 11 109 7 64 14 153、多项式(1)求多项式4xxp的根=x)2(3--Y=[1 0 -2 -4];S=roots(Y)S =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i(2)已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ,求矩阵A的特征多项式;把矩阵A作为未知数代入到多项式中;A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]P=poly(A)polyval(P,A)A =1.2000 3.0000 5.0000 0.90005.0000 1.7000 5.00006.00003.0000 9.0000 0 1.00001.00002.00003.00004.0000P =1.0000 -6.9000 -77.2600 -86.1300 604.5500ans =1.0e+003 *0.3801 -0.4545 -1.9951 0.4601-1.9951 0.2093 -1.9951 -2.8880-0.4545 -4.8978 0.6046 0.43530.4353 0.0840 -0.4545 -1.16174、基本绘图命令(1)绘制余弦曲线 y=cos(t),t∈[0,2π](2)在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5),t∈[0,2π](1)t=[0:0.05:2*pi];y=cos(t);plot(t,y)-1-0.8-0.6-0.4-0.20.20.40.60.81(2)t=[0:0.05:2*pi];y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1)hold onplot(t,y2)5、基本绘图控制绘制[0,4π]区间上的x1=10sint曲线,并要求:(1)线形为点划线、颜色为红色、数据点标记为加号;(2)坐标轴控制:显示范围、刻度线、比例、网络线(3)标注控制:坐标轴名称、标题、相应文本;t=[0:0.1:4*pi];x1=10*sin(t);plot(t,x1,'r-.+'); %画图,显示红色、点划线、标记加号;axis([0,15,-10,10]); %定义显示范围,横轴为[0,15],纵轴为[-10,10];title('曲线x1=10sint'); %显示标题;xlabel('T轴');ylabel('X1轴'); %显示坐标轴名称;set(gca,'xminortick','on');set(gca,'yminortick','on'); %显示刻度线;grid on %显示网络线T 轴X 1轴6、基本程序设计(1)编写命令文件:计算1+2+…+n<2000 时的最大n 值; (2)编写函数文件:分别用for 和while 循环结构编写程序,求2的0到n 次幂的和。

实验一MATLAB环境及命令窗口的使用

实验一MATLAB环境及命令窗口的使用

实验⼀MATLAB环境及命令窗⼝的使⽤实验⼀MATLAB环境及命令窗⼝的使⽤⼀、实验⽬的1.熟悉MA TLAB的操作环境及基本操作⽅法。

2.熟悉MA TLAB的通⽤参数设置。

3.熟悉MA TLAB的搜索路径及设置⽅法。

4.熟悉MA TLAB帮助信息的查阅⽅法。

⼆、实验内容及实验结果1.熟悉MATLAB的开发环境,了解主菜单和⼯具栏的内容,⼯作空间窗⼝、历史命令窗⼝、当前路径窗⼝的功能。

2.利⽤菜单设置MATLAB的Command Window中字体的⼤⼩,并更改输出格式。

⽂件夹。

4.在E盘根⽬录创建⽂件夹mypath,⽤菜单⽅法和path函数的⽅法将E:\mypath加⼊到搜路径中,指出两种⽅法的区别。

区别:⽤菜单⽅法,如果单击了“save”按钮则添加的搜索⽬录不会因为MATLAB的关闭⽽消失;利⽤path函数⽅法,path命令的扩展搜索路径仅在当前MATLAB环境下有效。

5.完成下列操作:(1)在MATLAB命令窗⼝下输⼊以下命令:x=0:pi/10:2*pi;y=sin(x);(2)在⼯作区间窗⼝选择变量y,再在⼯作空间窗⼝选择绘图菜单命令或在⼯具栏中点击绘图命令按钮,绘制变量y的图形。

(3)利⽤帮助学习save、load命令的⽤法,将⼯作区中变量全部保存在mydata.mat中,清空⼯作区间,重新载⼊变量x,y的信息,并把它们保存在mydata1.mat中。

6.计算y=1.3^3*sin(pi/3)*sqrt(26),实现(1)结果⽤format命令按不同的格式输出。

(2)观察在进⾏上述计算后命令历史窗⼝的变化,⽤功能键实现回调刚才计算的语句。

(3)回调计算语句,把sin改为sn运⾏,观察反馈信息。

若回调语句在语句后⾯加“;”号,看输出有何不同。

三、实验⼼得1.基本了解MATLAB得操作环境和操作⽅法。

2.学会利⽤菜单设置MATLAB的Command Window中字体⼤⼩及更改输出格式。

matlab实验一实验报告

matlab实验一实验报告

matlab实验一实验报告实验一:Matlab实验报告引言:Matlab是一种强大的数学软件工具,广泛应用于科学计算、数据分析和工程设计等领域。

本实验旨在通过使用Matlab解决实际问题,探索其功能和应用。

一、实验目的本次实验的主要目的是熟悉Matlab的基本操作和常用函数,了解其在科学计算中的应用。

二、实验内容1. 数值计算在Matlab中,我们可以进行各种数值计算,包括基本的加减乘除运算,以及更复杂的矩阵运算和方程求解。

通过编写相应的代码,我们可以实现这些功能。

例如,我们可以使用Matlab计算两个矩阵的乘积,并输出结果。

代码如下:```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A * B;disp(C);```2. 数据可视化Matlab还提供了强大的数据可视化功能,可以将数据以图表的形式展示出来,更直观地观察数据的规律和趋势。

例如,我们可以使用Matlab绘制一个简单的折线图,来展示某个物体在不同时间下的位置变化。

代码如下:```matlabt = 0:0.1:10;x = sin(t);plot(t, x);xlabel('Time');ylabel('Position');title('Position vs. Time');```3. 图像处理Matlab还可以进行图像处理,包括图像的读取、处理和保存等操作。

我们可以通过Matlab对图像进行增强、滤波、分割等处理,以及进行图像的压缩和重建。

例如,我们可以使用Matlab读取一张图片,并对其进行灰度化处理。

代码如下:```matlabimg = imread('image.jpg');gray_img = rgb2gray(img);imshow(gray_img);```三、实验结果与分析在本次实验中,我们成功完成了数值计算、数据可视化和图像处理等任务。

MATLAB实验

MATLAB实验

实验一 MATLAB 集成环境使用与基本操作命令练习一 实验目的熟悉MATLAB 语言编程环境;熟悉MATLAB 语言命令 二 实验仪器和设备装有MATLAB6.0以上计算机一台 三 实验原理MATLAB 是以复杂矩阵作为基本编程单元的一种程序设计语言。

它提供了各种矩阵的运算与操作,并有较强的绘图功能。

1.1 基本规则1.1.1 一般MATLAB 命令格式为[输出参数1,输出参数2,……]=(命令名)(输入参数1,输入参数2,……) 输出参数用方括号,输入参数用圆括号如果输出参数只有一个可不使用 括号。

1.1.2 %后面的任意内容都将被忽略,而不作为命令执行,一般用于为代码加注释。

1.1.3 可用↑、↓键来重现已输入的数据或命令。

用←、→键来移动光标进行修改。

1.1.4 所有MATLAB 命令都用小写字母。

大写字母和小写字母分别表示不同的变量。

1.1.5 常用预定义变量,如pi 、Inf 、NaN 、ans1.1.6 矩阵的输入要一行一行的进行,每行各元素用空格或“,”分开,每行用“;”分开。

如⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321AMATLAB 书写格式为A=[1 2 3 ;4 5 6 ;7 8 9] 在MATLABZ 中运行如下程序可得到A 矩阵 a=[1 2 3;4 5 6;7 8 9] a =1 2 3 4 5 6 7 8 91.1.7 需要显示命令的计算结果时,则语句后面不加“;”号,否则要加“;”号。

运行下面两种格式可以看出它们的区别:a=[1 2 3;4 5 6;7 8 9] a=[1 2 3;4 5 6;7 8 9];a=1 2 3 不显示结果4 5 67 8 91.1.8 当输入语句过长需要换行时,应加上“…”后再回车,则可续行输入。

1.2 文件管理常用命令1.2.1 帮助(HELP)命令MATLAB有很多命令,因此很不容易记忆。

使用HELP命令可以得到有关命令的屏幕帮助信息。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验一 Matlab软件的使用1、实验目的和要求㈠掌握Matlab软件的使用方法;㈡常用离散时间信号的产生、显示和运算。

2、实验内容㈠Matlab软件常用命令和工具箱的使用,学会简单的编程;㈡编程实现常用离散时间信号;㈢离散时间信号的叠加、移位、线性卷积等基本运算。

3、实验原理利用软件生成数字信号处理系统中所涉及的信号及信号的加、移位、卷积等运算。

4、主要仪器设备计算机,Matlab软件。

实验步骤:程序1:n=-5:30;x=[zeros(1,5),1,zeros(1,30)];stem(n,x,'fill');grid on程序2:n=-5:30;x=[zeros(1,5),ones(1,8),zeros(1,23)];stem(n,x,'fill');grid on;程序3:n=-15:20;x=[zeros(1,15),1,2,3,4,5,zeros(1,16)];stem(n,x,'fill');grid on;程序4:N=10;M=10;L=N+M-1;x=[1,2,3,4,5,zeros(1,5)];h=[1,2,1,2,zeros(1,6)];y=conv(x,h);n=0:L-1;stem(n,y,'*k');grid on ;实验结果:1、2、3、4、实验二 DTFT和Z变换的Matlab实现1、实验目的和要求㈠掌握DTFT和Z变换的Matlab实现的方法;㈡进一步了解离散时域系统的时频特性。

2、实验内容㈠ DTFT和z变换的计算和基本性质;㈡离散时域系统的仿真;㈢利用函数实现离散时域系统的频率响应,分析零、极点的分布,进一步了解离散时域系统的时域特性。

3、实验原理利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统频率响应等的图像并于笔算结果进行比较,找出异同。

编译合适程序能计算取值范围不同的离散卷积。

4、主要仪器设备计算机,Matlab软件。

实验步骤:程序1:n=0:30;x=zeros(1,length(n));h=zeros(1,length(n));x([find((n>=0)&(n<=4))])=1;h([find((n>=0)&(n<=8))])=0.5; figure(1)subplot(3,1,1); stem(n,x);axis([0,30,0,2]); title('输入序列'); xlabel('n'); ylabel('x(n)'); subplot(3,1,2); stem(n,h);axis([0,30,0,2]);title('冲激响应序列'); xlabel('n');ylabel('h(n)');y=conv(x,h);subplot(3,1,3);n=0:length(y)-1; stem(n,y);title('输出响应'); xlabel('n');ylabel('y(n)');程序2:b=[0.0181,0.0543,0.0543,0.0181];a=[1.000,-1.76,1.1829,-0.2781];w=pi*freqspace(500); H=freqz(b,a,w);MH=abs(H);AH=angle(H);subplot(2,1,1);plot(w/pi,MH);grid;axis([0,1,0,1]);xlabel('w(pi)');ylabel('|H|');title('幅度、相位响应'); subplot(2,1,2);plot(w/pi,AH);grid;xlabel('w(pi)'); ylabel('angle(H)');实验结果:1、2、实验三 DFS和DFT的Matlab实现1、实验目的和要求㈠加深DFS和DFT的算法原理和基本性质的理解;㈡熟悉利用FFT对典型数字信号进行频谱分析。

2、实验内容㈠使用Matlab实现程序周期序列的DFS;㈡利用FFT和IFFT编程实现两个有限长序列的循环卷积;3、实验原理对有限长序列使用离散Fouier变换(DFT)可以很好的反映序列的频谱特性,而且易于用快速算法在计算机上实现。

4、主要仪器设备计算机,Matlab软件。

实验步骤:程序1:function shiyan22()F=50;N=64;T=0.000625;n=1:N;x=cos(2*pi*F*n*T);subplot(2,1,1);plot(n,x);hold ontitle('x(n)');xlabel('n');X=fft(x);subplot(2,1,2);plot(n,X);gridtitle('FFT|X|');xlabel('f(pi)');程序2:n=0:30;%输入x(n)和冲激响应h(n)x=zeros(1,length(n));h=zeros(1,length(n));x([find((n>=0)&(n<=4))])=1;h([find((n>=0)&(n<=8))]) =0.5;subplot(3,1,1);stem(x);title('x(n)');axis([0,30,0,2]);subplot(3,1,2);stem(h);title('h(n)');axis([0,30,0,2]);X=fft(x);H=fft(h);Y=X.*H;y=ifft(Y);subplot(3,1,3);实验结果:1、2、实验四 IIR数字滤波器的设计和实现1、实验目的和要求㈠掌握冲击响应不变法和双线形变换法设计IIR数字滤波器的原理和方法;㈡通过观察对实际信号的滤波,获得对数字滤波的感性认识。

2、实验内容㈠巴特沃思滤波器设计实现模拟低通滤波器;㈡利用冲击响应不变法和双线形变换法对模拟滤波器进行A/D转换;3、实验原理用数字滤波器的单位脉冲响应序列)(nh模仿模拟滤波器的冲激响应)(tha,让)(nh正好等于)(tha的采样值,即)()(nThnha=。

4、主要仪器设备计算机,Matlab软件。

实验步骤:设计一个巴特沃思数字低通滤波器,设计指标如下:通带内πωω2.0=<p幅度衰减不大于1dB;阻带πωω35.0=>s幅度衰减不小于15dB。

编制实现该数字滤波器程序并且实现数字滤波。

function shiyan3()fp=200*pi;fs=300*pi;wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1;Fs=1000;[N,fc]=buttord(fp,fs,Rp,As,'s'); [b,a]=butter(N,fc,'s');w=[0:1000*2*pi];[hf,w]=freqs(b,a,1000);[d,c]=impinvar(b,a,Fs);wd=[0:512]*pi/512;hw1=freqz(d,c,wd);[f,e]=bilinear(b,a,Fs);hw2=freqz(f,e,wd);OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);ep=sqrt(10^(Rp/10)-1);Ripple=sqrt(1/(1+ep.^2));Attn=1/10^(As/20);subplot(3,1,1);plot(wd/pi,abs(hw1)/abs(hw1(1)));title('幅度响应(冲击响应不变法)');xlabel('w(pi)');ylabel('H');ylabel('H');axis([0,1,0,1.1]);set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[0,Attn,Ripple,1]); grid;subplot(3,1,2);plot(wd/pi,20*log10((abs(hw2)+eps)/abs(hw2(1))));title('幅度响应(双线性变换法(dB))');xlabel('w(pi)');ylabel('H');axis([0,1,-40,5]);set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[-50,-15,-1,0]); grid;subplot(3,1,3);plot(wd/pi,-angle(hw2));title('相位响应');xlabel('w(pi)');ylabel('pi unit');%axis([0,1,0,1.1]);set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[-1,0,1]);grid;实验结果:实验五 FIR数字滤波器的设计和实现1、实验目的和要求㈠掌握用窗函数法设计FIR数字滤波器的原理和方法;㈡熟悉线性相位FIR滤波器的幅频特性和相频特性;㈢了解不同窗函数对滤波器性能的影响。

2、实验内容㈠用Hanning窗设计线性相位带通滤波器;㈡分别改用矩形窗和Blackman窗设计线性带通滤波器;3、实验原理确定数字滤波器的性能要求,根据性能要求,合理选择单位脉冲响应)(nh的奇偶对称性,从而确定理想频率响应)(ωjdeH的幅频特性和相频特性;求理想单位脉冲响应)(nhd,选择适当的窗函数)(nw,根据)()()(nwnhnhd=求所需设计的FIR滤波器单位脉冲响应;4、主要仪器设备计算机,Matlab软件。

相关文档
最新文档