matlab第六章

编号：3109002964

《信号与系统仿真实验》

实验报告

实验时间：2011年5月6日

实验班级：09通信3班

实验报告总份（片）数： 1 份（片）实验教师：梁秀玲老师

信息工程学院工学一号馆404实验室

广东工业大学

广东工业大学实验报告

信息工程学院通信工程专业 3 班成绩评定_______

学号3109002964 姓名叶连响(合作者____号____) 教师签名_______

实验六题目信号的抽样和恢复第_10周星期_五_第_3-4_节

一、实验目的

（1）验证抽样定理；

（2）熟悉信号的抽样与恢复过程；

（3）通过实验，观察欠采样时信号频谱的混跌现象；

（4）掌握采样前后的信号频谱的变化，加深对采样定理的理解；

（5）掌握采样频率的确定方法。

二、程序设计实验

设计一模拟信号：x(t)=3sin(2pi*f*t),采样频率fs为5120Hz,取信号频率f=150Hz (正常取样)，和f=3000Hz（欠采样）两种情况进行采样分析。

1、t=0:0.006:1

f=150;

xa=3*sin(2*pi*f*t);

subplot(2,1,1);

plot(t,xa);grid

xlabel('时间,mseec');ylabel('幅值'); title('连续时间信号x_{a}(t)'); axis([0 1 -3.5 3.5]) subplot(2,1,2); Ts=1/5120; n=0:Ts:0.01; xs=3*sin(2*pi*f*n); k=0:length(n)-1; stem(k,xs);grid;

xlabel('时间,mseec');ylabel('幅值'); title('离散时间信号x[n]'); axis([0 length(n)-1 -3.5 3.5])

0.1

0.2

0.3

0.40.50.60.7

0.8

0.9

1

时间,mseec 幅值

连续时间信号x a (t)

5

10

15

20

253035

40

45

50

时间,mseec

幅值

离散时间信号x[n]

2、t=0:0.00001:1;

f=3000;

fs=5120;

x=3*sin(2*pi*f*t);

subplot(2,1,1)

plot(t,x);

xlabel('时间,mseec');ylabel('幅值'); title('连续时间信号');

axis([0 0.001 -3.2 3.2])

subplot(2,1,2);

T=1/fs;

n=0:T:1;

xs=3*sin(2*pi*f*n);

k=0:49;

stem(k,xs(1:50));

xlabel('时间,mseec');ylabel('幅值'); title('输出序列');

axis([0 50 -3.2 3.2])

0.1

0.2

0.3

0.40.50.60.7

0.8

0.9

1

x 10

-3

-20

2

时间,mseec

幅值

连续时间信号

5

10

15

20253035

40

45

50

时间,mseec

幅值

输出序列

三、实验总结

通过了学习本次实验更进一步了解matlab 的应用，利用matlab 来显示采样前后的信号频谱的变化，让人们通过图像对采样定理有更深的理解和认识。 四、实验思考问题

在域抽样定理中,为什么要求被抽样信号必须是带限信号？如果是频带无限的应如何处理？

带限信号即是带宽有限的信号，采样定理要求采样频率必须是信号最高频率的2倍以上，否则会出现频率混叠。如果是非带限（无限带宽），最高采样频率理论上要求无穷大，这是做不到的。典型的无限带宽信号有白噪声。一般需要在采样前面加抗混叠滤波器。加低通滤波器，调节滤波器的截止频率等于你希望保留的的信号的最高频率即可。

信号处理引论第五章及matlab作业答案

a=input('type in the first sequence ='); b=input('type in the second sequence ='); c=conv(a,b); M=length(c)-1; n=0:1:M; disp('output sequence =');disp(c) stem(n,c) xlabel('Time index n');ylabel('Amplitude'); type in the first sequence =[2 4 6 4 2 0 0 0] type in the second sequence =[3 -1 2 1] output sequence = Columns 1 through 9 6 10 18 16 18 12 8 2 0 Columns 10 through 11 0 0 ??? Undefined function or variable 'ylable'. Error in ==> E:\Matlab6p5FULL\bin\win32\Untitled.m On line 8 ==> xlabel('Time index n');ylabel('Amplitude'); type in the first sequence =[2 4 6 4 2 0 0 0] type in the second sequence =[3 -1 2 1] output sequence = Columns 1 through 9 6 10 18 16 18 12 8 2 0 Columns 10 through 11

第一章 基本MATLAB基本命令

MATLAB学习指南 第1章．基本MATLAB指令 1.1.基本数量运算 首先，我们来谈一谈怎样向程序中添加注释（例如此行）。注释就是我们向程序中添加的文字，用来解释我们在做什么。所以，如果我们或其他人以后读到此代码，就能很容易地理解代码在执行什么指令。在一个MATLAB文件中，如果百分号%出现在一行文字中，百分号后面的所有文字都是注释，而不是MATLAB想要进行解释的指令。首先，我们向屏幕输入信息，告诉计算机我们开始运行1.1节。 指令disp('字符串')在屏幕上展示了文本字符串。 接下来，我们令一个变量等于1。 如果x没有被声明，这一指令先为变量x在内存中分配一个空间，然后又将x的值1存储在所分配的空间中。同时，这一指令会将"x = 1"显示在屏幕上。通常，我们不希望像这样的输出结果把屏幕弄得杂乱，所以我们在指令的最后加上一个分号就可以使指令变得“不可见”。例如，我们使用下列指令把x的值“不可见地”变为2，然后将结果显示在屏幕上。x=2；x的值发生了改变，但是却没有显示在屏幕上 disp('我们已经改变了x的值')； 然后，通过输入“x”，不带分号，我们显示x的值。 现在，我们来看如何声明其它变量。 y=2*x; 对y的值进行初始化，使其等于x的值的2倍。 x=x+1;使x的值增加1。 z=2*x; 声明了另一变量z。 因为在声明变量z时x的值已经发生了变化，所以z不等于y。 接下来，我们想看存储在内存中的变量的列表。我们使用指令“who”来实现。 通过使用“whos”我们能得到更多的信息 我们也可使用这些指令来得到有关仅有的某些变量的信息。

我们想要去掉变量“差”。 使用指令“clear”来实现。 接下来，我们想要去掉变量x和y。 我们再次使用指令“clear”。 一般来说，好的程序设计模式都要求每行只编写一个指令；但是，MATLAB却允许将多个指令放在一行。 更普遍的情况是，由于语句结构的长度使得我们希望一个指令能够分成多行。这可通过使用3个点来实现。 最后，当使用指令“clear all”时，我们可一次去掉所有变量。 1.2.基本向量运算 声明一个变量最简单的，但不值得推荐的方法就是逐一输入其分量。 x显示x的值 一般来说，最好一下子声明一个变量，因为这样的话MATLAB从一开始就知道它需要分配多少内存。对于大型向量，这会更加有效。 y=[1 4 6] 与上述代码作用相同。 注意到这次声明了一个行向量。为得到一个列向量，我们或者使用转置（复合x的伴随矩阵）算符xT=x’;把实行向量x进行转置，或从一开始就将其变成列向量。 为了看到行向量和列向量在大小上的不同，使用指令“size”便可返回向量或矩阵的大小。

MATLAB语言及控制系统仿真_参考答案解析_第6章

6.4 控制系统频域分析MATLAB 仿真实训 6.4.1实训目的 1. 学会利用MATLAB 绘制开环系统的伯德图； 2. 学会利用MATLAB 绘制开环系统的极坐标图； 3. 掌握通过编程或相关命令求取系统稳定裕度的方法； 4. 通过仿真进一步理解掌握系统频域分析的有关知识。 6.4.2实训内容 1.已知单位负反馈系统的开环传递函数为 ()()() 10.5s 10.2s 1s 1.0k s G +++= ）（ 要求编程绘制50=k 时的极坐标图，确定曲线与负实轴的交点坐标及频率值； n=50; d=conv([0.1,1],conv([0.2,1],[0.5,1])); sys=tf(n,d); nyquist(sys) 曲线与负实轴的交点坐标为-3.76； 曲线与负实轴的交点频率值9.2； 2.绘制下列系统的伯德图，并要求在图上显示出幅值裕度、相角裕度等信息。 （1）()()() 18s 12s 2.6 s G ++= >> n=2.6; >> d=conv([2,1],[8,1]); >> sys=tf(n,d) Transfer function: 2.6

---------------------- 16 s^2 + 10 s + 1 >> margin(sys) 从图上信息可知，幅值裕度为无穷，相角裕度为88.2度。 （2）()()() 110s 1s 10 s G ++= s >> n=10; >> d=conv([1,0],conv([1,1],[10,1])); >> sys=tf(n,d) Transfer function: 10 ------------------------ 10 s^3 + 11 s^2 + s >> margin(sys)

matlab课后答案完整版

ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值 第一章 一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z = 2、a=::; >> b=exp*a).*sin(a+ 3、x=[2 4; 5]; y=log(x+sqrt(1+x.^2))/2 y =4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素 A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = + d=c+a*b/(a+b) d = + 第二章 二、4、（1） y=0;k=0; >> while y<3 k=k+1; y=y+1/(2*k-1); end >> display([k-1,y-1/(2*k-1)]) ans = 第三章 二1（1）

第五章MATLAB的绘画与图形处理

第五章 MATLAB 的绘画与图形处理 MATLAB 具有非常强大的二维和三维绘图功能，尤其擅长于各种科学运算结果的优秀可视化。 5.1二维曲线的绘制 5.1.1基本绘图命令plot plot 命令是MATLAB 中最简单而且使用最广泛的一个绘图命令，用来绘制二维曲线。 1. plot(x,y) 语法： plot(x,y) %绘制以x 为横坐标y 为纵坐标的二维曲线 【例5.1】绘制正弦曲线y=sin(x)和方波曲线，如图5.1所示。 >> x1=0:0.1:2*pi; >> y1=sin(x1); %y1为x1的正弦函数 >> plot(x1,y1); >> x2=[0 1 1 2 2 3 ]; >> y2=[1 1 0 0 1 1 ]; >> plot(x2,y2); >> axis([0 4,0 2]) %将坐标轴范围设定为0-4和0-2 2. plot(x1,y1,x2,y2,…)绘制多条曲线 plot 命令还可以同时绘制多条曲线，用多个矩阵对为参数，MATLAB 自动以不同的颜色绘制不同曲线。 【例5.2】绘制三条曲线，如图5.2所示。 图5.1 (a) 正弦曲线 (b) 方波曲线

>> x=0:0.1:2*pi; >> plot(x,sin(x),x,cos(x),x,sin(3*x)) %画三条曲线 图5.2 三条曲线 5.1.2绘制曲线的一般步骤 表5.1为绘制二维、三维图形一般步骤的归纳。 表5.1绘制二维、三维图形的一般步骤 内容 曲线数据准备： 对于二维曲线，横坐标和纵坐标数据变量； 说明： ?步骤1和3是最基本的绘图步骤，如果利用MA TLAB的默认设置通常只需要这两个基本步骤就可以基本绘制出图形，而其他步骤并不完全必需。 ?步骤2一般在图形较多的情况下，需要指定图形窗口、子图时使用。

matlab答案第五章

第五章 1.绘制下列曲线。 （1）2 1100 x y += （2）22 21x e y -=π （3）12 2 =+y x （4）? ??==3 2 5t y t x 答: (1) x=-10::10; y=100./(1+x.^2); plot(x,y) (2) x=-10::10; y=1/(2*pi)*exp(-x.^2/2); plot(x,y)

(3) ezplot('x^2+y^2=1')

(4) t=-10::10; x=t.^2; y=5*t.^3; plot(x,y) 2. 绘制下列三维图形。 （1）?? ? ??===t z t y t x sin cos （2）?? ? ??=+=+=u z v u y v u x sin sin )cos 1(cos )cos 1( （3）5=z （4）半径为10的球面（提示：用sphere 函数 答: (1) t=0:pi/100:2*pi; x=cos(t); y=sin(t);

z=t; plot3(x,y,z) (2) u=0:pi/100:2*pi; v=0:pi/100:2*pi; x=(1+cos(u)).*cos(v); y=(1+cos(u)).*sin(v); z=sin(u); plot3(x,y,z)

(3) [x,y]=meshgrid(-100:100); z=5*ones(size(x)); mesh(x,y,z) (4) [x,y,z] = sphere; %等价于sphere(20) 20是画出来的球面的经纬分面数...20的话就是20个经度, 20个纬度,或者说“就是划出20╳20面的单位球体，把球的表面分成400个小网格” surf(x,y,z); %绘制单位球面 x = 10*x; y = 10*y; z = 10*z; surf(x,y,z); %绘制半径为10的球面 axis equal

最新MATLAB应用基础第一章matlab基础

M A T L A B应用基础第一章m a t l a b基础

MATLAB应用 赵国瑞 天津大学电子信息工程学院计算机基础教学部 2000.3 制作

概述 MATLAB是世界流行的优秀科技应用软件之一。具有功能强大(数值计算、符号计算、图形生成、文本处理及多种专业工具箱)、界面友好，可二次开发等特点。 自1984年由美国MathWorks公司推向市场以来，先后发布了多个版本，1993年发布4.0版，1996年发布5.0版，1999年发布5.3版。目前发布的为6.5版。 MATLAB有专业和学生版之分。二者功能相同，但计算规模和计算难度有差别。 在国内外，已有许多高等院校把MATLAB列为本科生、研究生必须掌握的基本技能。我校自1999年列为研究生选修课程。而且有很多教师、研究生把它作为进行科研的重要工具。 国内关于MATLAB的书籍很多，如： 《精通MATLAB 5.3》张志涌等编著北京航空航天大学出版社，2000.8

《科学计算语言MATLAB简明教程》杜藏等编著南开大学出版社，1999.6 《精通MATLAB 5》张宜华编写清华大学出版社，1999.6 《精通MATLAB--综合辅导与指南》 Duane Hanselman、Bruce Littlefield编著李人厚等译较西安交通大学出版社，1998.1 等等 本课程主要介绍MATLAB 5.3的基本功能和基础知识。至于其包含的多种工具箱，如仿真工具箱、解非线性方程（组）工具箱、优化工具箱等，应通过本学习后，结合各专业自己进一步学习和使用。 第1章MATLAB基础 1.1 源文件（M－文件） 分为两类：函数文件和非函数文件。 都用扩展名.M 1.1.1函数文件

控制系统仿真matlab第六章习题答案

控制系统仿真第六章课后题 作业6.1 在图6.1中，已知单位负反馈系统被控对象的传递函数为) 1001.0)(11.0()(++=s s s K s G 试编写matlab 程序，设计系统的超前矫正器Gc(s),要求： 1）在斜坡信号r （t)=2t 作用下，系统的稳态误差ess<=0.002; 2)校正后系统的相位裕度Pm 范围为：45~55； 3）绘制系统校正后的bode 图和阶跃响应曲线。 程序： >> s=tf('s'); >> G=1000/(s*(0.1*s+1)*(0.001*s+1)); >> margin(G) % 绘制校正前的bode 图 >> figure(2) >> sys=feedback(G,1); >> step(sys) %绘制校正前的单位阶跃响应曲线

>> [Gm,Pm]=margin(G); %该句值计算bode图的增益裕量Gm和相位裕量Pm >> [mag,phase,w]=bode(G); %该句只计算bode图上多个频率点w对应的幅值和相位>> QWPm=50; %取矫正后的相位为50 >> FIm=QWPm-Pm+5; >> FIm=FIm*pi/180; >> alfa=(1-sin(FIm))/(1+sin(FIm)); >> adb=20*log10(mag); >> am=10*log10(alfa); >> wc=spline(adb,w,am); >> T=1/(wc*sqrt(alfa)); >> alfat=alfa*T; >> Gc=tf([T 1],[alfat 1]) %校正器的传递函数 Transfer function: 0.01794 s + 1 ------------- 0.00179 s + 1 >> figure(3) >> margin(Gc*G) % 系统矫正后的bode图

Matlab第六章作业

Y X Z+ = 1 XY Z= 2 Matlab作业第六章 姓名：袁伟薯学号：2014301510030 习题1： 使用Combinatoinal Logic模块完成对以下函数功能的建模和仿真： AC BC AB X+ + =) )( )( (A C C B B A Y+ + + = 仿真模型实现步骤如下： 1.模型分析 根据逻辑表达式写出真值表 2.模型搭建 首先建立一个名为combinational_logic_hm1的simulink仿真文件，搭建模型。需要以下模块：3个Pulse Generator模块，2个Logical Operator模块，2个combinatorial logic模块，scope 模块，2个Mux模块以及3个Data Type Conversion模块 注：如果在没有Data Type Conversion的情况下，仿真不能运行会出现如下错误

查阅PPT和搜阅资料得知combinatorial logic模块接收boolean类型输出，输出也是boolean 类型，可以有两种方法解决：（1）将仿真参数中Optimization中的Implement logic signals as boolean data(V.S. double)去掉，避免数据类型的不匹配。（2）把数据转换模块Data Type Conversion加在conbinatorial Logic模块输出后。得到搭建的模型如下 3.各模块参数的设定如下： Pulse Generator模块参数设置 4.分析总结： （1）输出波形分析：仿真模型搭建完毕后，在Simulation/Simulation parameters内调

MATLAB第六章simulink仿真答案

实验四 SIMULINK 仿真实验 一、 实验目的 1. 学习SIMULINK 的实验环境使用。 2. 掌握SIMULINK 进行结构图仿真的方法。 二、 实验内容 1．控制系统结构图仿真 给定被控对象) 1(10 )(+= s s s G ,控制器111.0145.0)(++=s s s D ,按以下两种情 况设计SIMULINK 仿真结构图（给定信号是单位阶跃信号）。 （1） 无控制器时被控对象单位负反馈。 （2） 控制器与被控对象串连接成单位负反馈。 给定的仿真参数： （1） 信号源参数设置: 阶跃信号（Step ）的Step time 设为0秒。 （2）仿真参数设置: 仿真时间 0～10秒，求解器选定步长(Fixed-step)的ode5,仿真步长（Fixed step size ）设为0.02秒。 实验要求： (1) 在SIMULINK 中对设计的结构图进行仿真，观察输入信号，输出信号和控制信号。 (1) 记录保存两种情况下的响应波形（适当调整时间轴和纵轴坐标，使图形显示适中，同时在图中求出系统的超调和调节时间（按2％的误差带）。

123 45678910 00.20.40.60.8 1 1.2 1.4 时间 (seconds) d a t a 时序图: 0123 45678910 0.20.40.60.81 1.21.41.6 1.8时间 (seconds) d a t a 时序图:

2．动态系统微分方程仿真 在SIMULINK 中求解下列二阶微分方程代表的动态系统在阶跃信号作用下的状态响应。 给定的仿真参数： （1） 信号源参数设置：阶跃信号（Step ）的Step time 设为0秒。 （2） 仿真参数设置：仿真时间 0～8秒，求解器选变步长 (Variable-step)的ode45,最大仿真步长（Max step size ）设为0.01秒。 实验要求： （1） 据微分方程构造结构图。 （2） 结构图仿真。 a) 零状态仿真：x 1=0,x 2=0， b) 非零状态仿真：x 1=1，x 2=-1， c) 记录保存两种情况下的响应波形（适当调整时间轴和纵轴坐 标，使图形显示适中）。 1 22122110) (1,||210x y t u u x x x x x x ==++--==

自动控制原理习题全解及MATLAB实验 第6章习题解答

第6章控制系统的校正 本章主要讨论利用频率法对单输入-单输出的线性定常系统的综合和设计。在介绍控制系统校正的基本概念、控制系统的基本控制规律的基础上，介绍了各种串联校正装置（超前校正装置、滞后校正装置、滞后-超前校正装置）的特性及按分析进行相应设计的基本步骤和方法；还介绍了期望设计法的基本概念、常见的期望特性和设计步骤；另外还介绍了根轨迹法的串联校正和反馈校正的基本概念和方法；最后介绍了利用MATLAB进行控制系统校正。 教材习题同步解析 试分别说明系统的固有频率特性与系统期望频率特性的概念。 答：系统本身固有元件所具有的频率特性称为固有频率特性。设计者希望系统所能达到的频率特性称为系统期望频率特性。 试比较串联校正和反馈校正的优缺点。 答：a、校正装置和未校正系统的前向通道环节相串联，这种叫串联校正，串联校正是最常用的设计方法，设计与实现比较简单，通常将串联装置安置在前向通道的前端。 b、并联校正也叫反馈校正，它是和前向通道的部分环节按反馈方式连接构成局部反馈回路，设计相对较为复杂。并联校正一般不需要加放大器，它可以抑制系统的参数波动及非线性因素对系统性能的影。 PD控制为什么又称为超前校正？串联PD控制器进行校正为什么能提高系统的快速性和稳定性？ 答：加入PD控制相当于在系统中加入一个相位超前的串联校正装置，使之在穿越频率处有较大的相位超前角。因此，PD控制称为超前控制。PD控制的传递函数为G s Kp sτ =+，由比例控制和微分控制组合而成。增大比例系数Kp，可以展宽系统的()(1) 通频带，提高系统的快速性。微分控制反映信号的变化率的预报作用，在偏差信号变化前给出校正信号，防止系统过大地偏离期望值和出现剧烈振荡倾向，有效地增强系统的相对稳定性。 PI控制为什么又称为滞后校正？串联PI控制器进行校正为什么能提高系统的稳态性能？如何减小它对系统稳定性的影响？ 答：PI控制在低频段产生较大的相位滞后，所以滞后校正。PI控制的比例部分可以提高系统的无差度，改善系统的稳态性能。在串入系统后应使其转折频率在系统幅值穿越频率

第5章--MATLAB绘图-习题答案

第5章 MATLAB绘图 习题5 一、选择题 1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。D A．12 B．7 C．4 D．3 2．下列程序的运行结果是（）。A x=0:pi/100:2*pi; for n=1:2:10 plot(n*sin(x),n*cos(x)) hold on end axis square A．5个同心圆 B．5根平行线 C．一根正弦曲线和一根余弦曲线 D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。C A．{\alpha}+{\beta} B．{\α}+{\β} C．α+β D．\α+\β4．subplot(2,2,3)是指（）的子图。A A．两行两列的左下图 B．两行两列的右下图 C．两行两列的左上图 D．两行两列的右上图 5．要使函数y=2e x的曲线绘制成直线，应采用的绘图函数是（）。C A．polar B．semilogx C．semilogy D．loglog 6．下列程序的运行结果是（）。B [x,y]=meshgrid(1:5); surf(x,y,5*ones(size(x))); A．z=x+y平面 B．与xy平面平行的平面 C．与xy平面垂直的平面 D．z=5x平面 7．下列函数中不能用于隐函数绘图的是（）。D A．ezmesh B．ezsurf C．ezplot D．plot3 8．下列程序运行后，看到的图形（）。C t=0:pi/20:2*pi; [x,y]=meshgrid(-8:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);

第6章matlab数据分析与多项式计算_习题答案

第6章 MATLAB数据分析与多项式计算 习题6 一、选择题 1．设A=[1,2,3,4,5;3,4,5,6,7]，则min(max(A))的值是（）。B A．1 B．3 C．5 D．7 2．已知a为3×3矩阵，则运行mean(a)命令是（）。B A．计算a每行的平均值 B．计算a每列的平均值 C．a增加一行平均值 D．a增加一列平均值 3．在MATLAB命令行窗口输入下列命令： >> x=[1,2,3,4]; >> y=polyval(x,1); 则y的值为（）。 D A．5 B．8 C．24 D．10 4．设P是多项式系数向量，A为方阵，则函数polyval(P,A)与函数polyvalm(P,A)的值（）。D A．一个是标量，一个是方阵 B．都是标量 C．值相等 D．值不相等 5．在MATLAB命令行窗口输入下列命令： >> A=[1,0,-2]; >> x=roots(A); 则x(1)的值为（）。 C A．1 B．-2 C． D． 6．关于数据插值与曲线拟合，下列说法不正确的是（）。A A．3次样条方法的插值结果肯定比线性插值方法精度高。 B．插值函数是必须满足原始数据点坐标，而拟合函数则是整体最接近原始数据点，而不一定要必须经过原始数据点。 C．曲线拟合常常采用最小二乘原理，即要求拟合函数与原始数据的均方误差达到极小。 D．插值和拟合都是通过已知数据集来求取未知点的函数值。 二、填空题 1．设A=[1,2,3;10 20 30;4 5 6]，则sum(A)= ，median(A)= 。 [15 27 39]，[4 5 6[ 2．向量[2,0,-1]所代表的多项式是。2x2-1 3．为了求ax2+bx+c=0的根，相应的命令是（假定a、b、c已经赋值）。为了

MATLAB语言与控制系统仿真-参考答案-第5章

5.6 控制系统的时域响应MATLAB 仿真实训 5.6.1实训目的 1. 学会利用MATLAB 绘制系统的单位阶跃响应曲线，掌握读取系统动态性能指标的方法； 2. 学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法； 3. 掌握利用MATLAB 绘制系统的零输入响应曲线的方法； 4. 掌握利用MATLAB 绘制系统的一般输入响应曲线的方法； 5. 学会通过仿真曲线读取相关信息，并依据有关信息进行系统的时域分析。 5.6.2实训内容 1.编写程序求取下列各系统的单位阶跃响应，完成表5-5并记录相关曲线。 162.316)(21++= s s s G 16 4.216 )(22 ++=s s s G 166.116)(2 3++=s s s G 1616 )(24++=s s s G 解：>> n1=16; >> d1=[1,3.2,16]; >> sys1=tf(n1,d1); >> step(sys1) >> n2=16; >> d2=[1,2.4,16]; >> sys2=tf(n2,d2); >> step(sys2)

>> n3=16; >> d3=[1,1.6,16]; >> sys3=tf(n3,d3); >> step(sys3) >> n4=16; >> d4=[1,1,16]; >> sys4=tf(n4,d4); >> step(sys4)

w=4; cmax1=1+exp(-z1*pi/sqrt(1-z1^2)); tp1=pi/(w*sqrt(1-z1^2)); ts1=3.5/(z1*w); [cmax1,tp1,ts1] ans = 1.2538 0.8569 2.1875 >> z2=0.3; w=4; cmax2=1+exp(-z2*pi/sqrt(1-z2^2)); tp2=pi/(w*sqrt(1-z2^2)); ts2=3.5/(z2*w); [cmax2,tp2,ts2]

matlab第五章课后答案

第一题：某河床的横截面积，为了计算最大的排洪量，需要计算他的横断面积，用梯形法计算其断面积。 >> x=[0 4 10 12 15 22 28 34 40]; y=[0 1 3 6 8 9 5 3 0]; trapz(x,y) ans = 178.5000 答：用梯形积分法求的其面积为178.5平方米。 第五题：求下列积分的数值解。 （2）. e 2x 2π 0cos ?3x dx >>fun=inline('exp(2*x).*cos(x).^3'); >>quadl(fun,0,2*pi) ans = 9.7054e+004 答：所求积分的数值解为：9.7054*10^4。 （4） sinx x 10dx .>>fun=@(x)sin(x)./x; quadl(fun,1e-10,1) ans = 0.9461 答：所求积分的数值解为：0.9461. （6） dθ 2sinθ1 02π0

.>>fun=inline('sqrt(1+r.^2.*sin(th))','r','th'); dblquad(fun,0,1,0,2*pi) ans = 6.1879 答：所求积分的数值解为：6.1879. （7）1+x+y dydx,D为x2+y2≤2x .function s=dblquad2(fun,a,b,clo,dhi,n) ifnargin<6,n=100;end x=linspace(a,b,n+1); s=0; for i=1:n s=s+dblquad(fun,x(i),x(i+1),feval(clo,(x(i)+x(i+1))/2),feval(dhi,(x(i)+x(i+1))/ 2)); end >>clear; >>fun=inline('1+x+y.^2','x','y'); >>clo=inline('-sqrt(2*x-x.^2)'); >>dhi=inline('sqrt(2*x-x.^2)'); >>dblquad2(fun,0,2,clo,dhi,100) ans = 7.0705

matlab第六章部分习题解答

MATLAB实验报告 第6章习题 实验日期：2014.12.3 1.>> clear >> syms k x y a b; >> r1=symsum(x^(2*k)/sym('k!'),k,0,inf) 结果： r1 = exp(x^2) 2.>> clear >> syms x a; >> r11=a*x/sin(x); >> r12=[2*sin(x)/x,x/tan(x),x,x^2]; >> r13=(2+2/x)^x; >> r14=a*(1+1/x)^x; >> r21=limit(r11,x,0) >> r22=limit(r12,x,0) >> r23=limit(r13,x,inf)

>> r24=limit(r14,x,inf) 结果： r21 = a r22 = [ 2, 1, 0, 0] r23 = Inf r24 = exp(1)*a 3.（1）>> clear >> syms x; >> y1=taylor(cos(x),pi/2,6) >> y11=subs(y1,x,pi/2) 结果： y1 = -x+1/2*pi+1/6*(x-1/2*pi)^3-1/120*(x-1/2*pi)^5

（2）>> z1=taylor(sin(x),pi/2,4) >> z11=subs(z1,x,pi/2) 结果： z1 = 1-1/2*(x-1/2*pi)^2 z11 = 1 （3）>> p1=taylor(2/(x-2),1.5,4) >> p11=subs(p1,x,1.5) 结果： p1 = 8-8*x-16*(x-3/2)^2-32*(x-3/2)^3 p11 = -4 （4）>> q1=taylor(x/sin(x),pi/2,4) >> q11=subs(q1,x,pi/2) 结果：

matlab习题第六章

第六章 实验指导 1、(1) >> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.52]'; x=A\p x = 1.2000 0.6000 0.6000 (2)、 >> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.52]'; >> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]'; >> x=A\p x = 3.0000 -6.6000 6.6000 （3）、 >> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; p=[0.95,0.67,0.53]'; x=A\p cond(A) x = 3.0000 -6.6000 6.6000 ans =

1.3533e+003 2、（1） M文件： function fx=funx(x) fx=x^41+x^3+1; 运行; >> z=fzero(@funx,-1) z = -0.9525 (2) M文件： function fx=funx2(x) fx=x-sin(x)/x; 运行： >> x2=fzero(@funx2,0.5) x2 = 0.8767 （3） M文件： function q=fun3(p) x=p(1); y=p(2); z=p(3); q(1)=sin(x)+y^2+log(z)-7; q(2)=3*x+2^y-z^3+1; q(3)=x+y+z-5; 运行： >> options=optimset('Display','off'); x=fsolve(@fun3,[1,1,1]',options) q=fun3(x) x = 0.5991 2.3959 2.0050

第5章MATLAB绘图_习题答案

第5章MATLAB绘图 习题5 一、选择题 1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。D A．12B．7C．4D．3 2．下列程序的运行结果是（）。A x=0:pi/100:2*pi; forn=1:2:10 plot(n*sin(x),n*cos(x)) holdon end axissquare A．5个同心圆B．5根平行线 C．一根正弦曲线和一根余弦曲线D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。C A．{\alpha}+{\beta}B．αβ}C．α+βD．αβ 4．subplot(2,2,3)是指（）的子图。A A．两行两列的左下图B．两行两列的右下图 C．两行两列的左上图D．两行两列的右上图 x的曲线绘制成直线，应采用的绘图函数是（）。C 5．要使函数y=2e A．polarB．semilogxC．semilogyD．loglog 6．下列程序的运行结果是（）。B [x,y]=meshgrid(1:5); surf(x,y,5*ones(size(x))); A．z=x+y平面B．与xy平面平行的平面 C．与xy平面垂直的平面D．z=5x平面 7．下列函数中不能用于隐函数绘图的是（）。D A．ezmeshB．ezsurfC．ezplotD．plot3 8．下列程序运行后，看到的图形（）。C t=0:pi/20:2*pi; [x,y]=meshgrid(-8:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);

第一章MATLAB基础知识

第一章MATLAB基础知识 1.1 MATLAB开发环境 1.MATLAB操作桌面 MATLAB4.0以上版本都是在Windows以上环境支持下操作与运行的，因此，这里必须假定读者对Windows有一定的了解和掌握。本书以Windows98环境，MATLAB6.5版本为例介绍。书中绝大多数例子在MA TLAB5.3版本中亦能使用. 要想进入MATLAB系统，方法有二： （1）在Windows98的桌面上双击MATLAB快捷图标。 （2）单击Windows98的“开始”按钮，再依次单击：程序/MA TLAB/MA TLAB6.5； 图1-1 MATLAB6.5除保留了传统的命令窗口外，还增加了启动平台窗口、工作空间窗口、命令历史窗口、当前路径窗口等窗口，与新的M文件编辑器和新的在线帮助浏览器等共同构成了MA TLAB6.5的开发环境。 MATLAB的开发环境是MATLAB语言的基础和核心部分，MATLAB语言的全部功能都是在MA TLAB的开发环境中实现的，MA TLAB的仿真工具Simulink、MATLAB

的工具箱等其他附加功能的实现也必须使用MATLAB开发环境，因此，掌握MATLAB 的开发环境是掌握MATLAB语言的关键。 启动MA TLAB后，将显示包括命令窗口、启动平台窗口、工作空间窗口、命令历史窗口和当前路径窗口等5个窗口和主菜单组成的操作桌面（图1-1）。本节对操作桌面的各个窗口作简要介绍，部分窗口的功能和使用将在以后的章节中详细介绍。 操作桌面缺省状态下显示3个窗口，当前路径和工作空间窗口在同一位置显示，可以通过该窗口下方的箭头或窗口标签来切换，命令历史窗口和当前目录窗口在同一位置显示，可以通过该窗口下方窗口标签来切换。每个窗可以脱离操作桌面独立出来。也可以通过菜单View来选择显示哪些窗口。 MATLAB还设定了几种特定的窗口布局方式，在View菜单的Desktop Layout选项中，给定了6种布局方式： Default（缺省方式） Command Window Only（只显示命命令窗口方式） Simple（简单方式，只有命令历史窗口和命令窗口） Short History（低命令历史窗口方式） Tall History（高命令历史窗口方式） Five Panel（5个窗口平铺方式） 2.MATLAB的通用参数设置 MATLAB的通用参数和各功能窗口的参数可以通过主菜单中的file/Preferences项设置，这里先介绍通用参数的设置。 在主菜单中选择Preferences项，打开Preferences窗口（图1-2），缺省状态为通用参数设置，其选项包括： Display（显示选择）Show tooltips（显示相关信息），当单选框选择后，鼠标放在工具栏的快捷按钮上时显示相关信息。 Toolbox caching（工具箱缓冲区）对于远程使用MATLAB的用户，应选择单选框Enable Toolbox caching（打开工具箱缓冲区），建立一个高速缓冲区，以提高使用速度，对单机用户该选项作用不大。 Figure window printing（图形窗口打印）有3个选项，分别为Use printer default（按打印机缺省设置输出）、Always send as black and white（按黑白图形输出）和Always send as color（按彩色图形输出）。

第5章--MATLAB绘图-习题答案上课讲义

第5章--M A T L A B绘图-习题答案

第5章 MATLAB绘图 习题5 一、选择题 1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。D A．12 B．7 C．4 D．3 2．下列程序的运行结果是（）。A x=0:pi/100:2*pi; for n=1:2:10 plot(n*sin(x),n*cos(x)) hold on end axis square A．5个同心圆 B．5根平行线 C．一根正弦曲线和一根余弦曲线 D．5根正弦曲线和5根余弦曲线 3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。C A．{\alpha}+{\beta} B．{\α}+{\β} C．α+β D．\α+\β 4．subplot(2,2,3)是指（）的子图。A A．两行两列的左下图 B．两行两列的右下图 C．两行两列的左上图 D．两行两列的右上图 5．要使函数y=2e x的曲线绘制成直线，应采用的绘图函数是（）。C A．polar B．semilogx C．semilogy D．loglog 6．下列程序的运行结果是（）。B [x,y]=meshgrid(1:5); surf(x,y,5*ones(size(x))); A．z=x+y平面 B．与xy平面平行的平面 C．与xy平面垂直的平面 D．z=5x平面 7．下列函数中不能用于隐函数绘图的是（）。D A．ezmesh B．ezsurf C．ezplot D．plot3 8．下列程序运行后，看到的图形（）。C t=0:pi/20:2*pi; [x,y]=meshgrid(-8:0.5:8);

最新第6章 MATLAB数据分析与多项式计算_习题答案

精品好文档，推荐学习交流 第6章 MATLAB数据分析与多项式计算 习题6 一、选择题 1．设A=[1,2,3,4,5;3,4,5,6,7]，则min(max(A))的值是（）。B A．1 B．3 C．5 D．7 2．已知a为3×3矩阵，则运行mean(a)命令是（）。B A．计算a每行的平均值B．计算a每列的平均值 C．a增加一行平均值D．a增加一列平均值 3．在MA TLAB命令行窗口输入下列命令： >> x=[1,2,3,4]; >> y=polyval(x,1); 则y的值为（）。D A．5 B．8 C．24 D．10 4．设P是多项式系数向量，A为方阵，则函数polyval(P,A)与函数polyvalm(P,A)的值（）。D A．一个是标量，一个是方阵B．都是标量 C．值相等D．值不相等 5．在MA TLAB命令行窗口输入下列命令： >> A=[1,0,-2]; >> x=roots(A); 则x(1)的值为（）。C A．1 B．-2 C．1.4142 D．-1.4142 6．关于数据插值与曲线拟合，下列说法不正确的是（）。A A．3次样条方法的插值结果肯定比线性插值方法精度高。 B．插值函数是必须满足原始数据点坐标，而拟合函数则是整体最接近原始数据点，而不一定要必须经过原始数据点。 C．曲线拟合常常采用最小二乘原理，即要求拟合函数与原始数据的均方误差达到极小。 D．插值和拟合都是通过已知数据集来求取未知点的函数值。 二、填空题 1．设A=[1,2,3;10 20 30;4 5 6]，则sum(A)= ，median(A)= 。 [15 27 39]，[4 5 6[ 2．向量[2,0,-1]所代表的多项式是。2x2-1 仅供学习与交流，如有侵权请联系网站删除谢谢1