陈后金信号与系统第二版 第二章matlab

52页例2-13

>> %program2-8

>> t=-3:0.001:3;

>> ft=tripuls(2*t,4,0.5);

>> ft1=tripuls(2*t,4,0.5);

>> subplot(2,1,1)

>> plot(t,ft1)

>> title('f(2t)')

>> ft2=tripuls((2-2*t),4,0.5);

>> subplot(2,1,2)

>> plot(t,ft2)

>> title('f(2-2*t)')

53页例2-14

>> %program2_9 the energy of exponential sequence >> k=0:10;

>> A=1;a=-0.6;

>> fk=A*a.^k;

>> W=sum(abs(fk).^2)

W =

1.5625

54页例2-15

function yt = f2_2(t)

yt=tripuls(t,4,0.5);

%program2_10 differentiation

h=0.001;t=-3:h:3;

y1=diff(f2_2(t))*1/h;

plot(t(1:length(t)-1),y1)

title('df(t)/dt')

-3-2-10123

-1.2-1

-0.8

-0.6

-0.4

-0.2

0.2

0.4

df(t)/dt

%program2_11 integration

t=-3:0.1:3;

for x =1:length(t)

y2(x)=quad('f2_2',-3,t(x)); end

plot(t,y2)

title('integral of f(t)')

48页图2-2

>> %program2_1 decaying exponential signal

>> A=1;a=-0.4;

>> t=0:0.01:10;

>> ft=A*exp(a*t);

>> plot(t,ft)

01234567891000.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

48页图2-3

>> %program2_2 sinusoidal signal

>> A=1;

>> w0=2*pi;

>> phi=pi/6;

>> t=0:0.001:8;

>> ft=A*sin(w0*t+phi);

>> plot(t,ft)

>>

012345678-1-0.8

-0.6

-0.4

-0.2

0.2

0.4

0.6

0.8

1

>> %program2_3sample function

>> t=-3*pi:pi/100:3*pi;

>> ft=sinc(t/pi);

>> plot(t,ft)

-10-8-6-4-20246810-0.4-0.2

0.2

0.4

0.6

0.8

1

Matlab上机实验答案

Matlab上机实验答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/180)/(1+exp(2)) z1 = >> x=[2 1+2i; 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = - + + -

>> a=::; >> z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) (>> z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)可以验证z3==z33，是否都为1） z3 = Columns 1 through 5 + + + + + Columns 6 through 10 + + + + + Columns 11 through 15 + + + + + Columns 16 through 20 + + + + +

Columns 21 through 25 + + + + + Columns 26 through 30 + + + + + Columns 31 through 35 + + + + + Columns 36 through 40 + + + + + Columns 41 through 45 + + + + + Columns 46 through 50

+ + + + + Columns 51 through 55 + + + + + Columns 56 through 60 + + + + + Column 61 + (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:: >> t=0::; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^ 2-2.*t+1) z4 =

MATLAB全部实验及答案

MATLAB全部实验及答案 实验一、MATLAB基本操作 实验内容及步骤 4、有关向量、矩阵或数组的一些运算 （1）设A=15；B=20；求C=A+B与c=a+b？ （2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与 A.*B？ A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？ （4）设a=[1 -2 3;4 5 -4;5 -6 7] 请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全 下标的形式），并将其单下标转换成全下标。 clear,clc a=[1 -2 3;4 5 -4;5 -6 7]; [x,y]=find(a<0); c=[]; for i=1:length(x) c(i,1)=a(x(i),y(i)); c(i,2)=x(i); c(i,3)=y(i); c(i,4)=(y(i)-1)*size(a,2)+x(i); end c

（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？前面那 个是虚数矩阵，后面那个出错 （6）请写出完成下列计算的指令： a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？ a^2= 22 16 16 25 26 23 26 24 28 a.^2= 1 4 9 9 16 4 25 4 9 （7）有一段指令如下，请思考并说明运行结果及其原因 clear X=[1 2;8 9;3 6]; X( : ) 转化为列向量 （8）使用三元组方法，创建下列稀疏矩阵 2 0 8 0 0 0 0 1 0 4 0 0 6 0 0 0 方法一： clear,clc

MATLAB)课后实验答案[1]

实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221ln(1)2 z x x =++，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 2242011 122123t t z t t t t t ?≤=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)

4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果： m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 (2). 建立一个字符串向量例如： ch='ABC123d4e56Fg9';则要求结果是： ch='ABC123d4e56Fg9'; k=find(ch>='A'&ch<='Z'); ch(k)=[] ch =

matlab所有实验及答案

实验二 习题 1、 矩阵Y= ???? ? ???? ???3472123100451150425 ，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。取出子矩阵?? ? ? ??21301，并求该矩阵的维数。 解：命令为： Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10) sub2ind([4 3],2,3) [i,j]=ind2sub([4 3],10) find(Y==100) sub2ind([4 3],3,2) B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y) 2、 建立一个数值范围为0—100内4*5的整数随机矩阵。 求出大于50的元素的位置。 解：命令为： G=int8(100*rand(4,5)) find(G>50) 3、 已知矩阵A=[1 0 -1 ；2 4 1； -2 0 5]，B=[0 -1 0；2 1 3；1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ，A/B 、A\B 解：命令为： A=[1 0 -1 ；2 4 1； -2 0 5] B=[0 -1 0；2 1 3；1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B 4、 利用函数产生3*4阶单位矩阵和全部元素都为8的 4*4阶矩阵，并计算两者的乘积。 解：命令为： A=eye(3,4) B=8*ones(4) C=A*B 5、 创建矩阵a=????? ???????------7023021.5003.120498601 ,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b 构成矩阵d ，计算a*b 、c

matlab实验四及其答案

实验四GUI 设计专业 学号姓名成绩电气工程及其自 动化201409140305 杨诚1.创建GUI 绘制方程c bx ax y ++=2图形，需要显示绘图结果的坐标系窗口，还能够输入 a ， b ， c 的值和x 取值范围（最大值和最小值）。 回调函数的编写： a=str2num(get(handles.edit1,'String'));b=str2num(get(handles.edit2,'String'));c=str2num(get(handles.edit3,'String'));xmin=str2num(get(handles.edit4,'String'));xmax=str2num(get(handles.edit5,'String'));x=xmin:0.1:xmax;y=a*x.^2+b*x+c;plot(x,y); 设计的运行界面截图：

2.设计一个GUI，完成画出y=sin(x)、y=cos(x)和y=x 3.的波形图。回调函数的编写： x=-5:0.1:5 plot(x,sin(x)) x=-5:0.1:5 plot(x,cos(x))

x=-5:0.1:5 y=x.^3 plot(x,y)

设计的运行界面截图： 3.创建一个GUI，含有下拉菜单，下拉菜单中有背景颜色选择。回调函数的编写： yanse=get(handles.popupmenu1,'value'); switch yanse case1 set(gcf,'color','r'); case2 set(gcf,'color','y'); case3 set(gcf,'color','g'); end 设计的运行界面截图：

MATLAB实验题目及答案

实验二一维二维数组的创建和寻访 一、实验目的 1、掌握一维数组、二维数组创建和寻访的几种方法。 2、区别数组运算和矩阵运算的差别。 3、熟悉执行数组运算的常用数组操作函数。 4、掌握数组运算中的关系和逻辑操作及常用的关系、逻辑函数。 5、掌握“非数”、“空”数组在MA TLAB中的应用。 二、实验主要仪器与设备 装配有MA TLAB7.6软件的计算机 三、预习要求 做实验前必须认真复习第三章MATLAB的数值数组及向量化运算功能。 四、实验内容及实验步骤 1、一维数组的创建方法有哪几种？举例说明。 答：一维数组的创建方法有： ①递增/递减型一维数组的创建：冒号生成法：x=a：inc：b 线性（或对数）定点法：x=linspace(a,b,n)，x=logspace(a,b,n) ②逐个元素输入法：如x=[0.1，sin(pi/5)，-exp(-3)，-2*pi] ③运用MA TLAB函数生成法：例ones，rand等。 2、输入以下指令，并写出运行结果。本例演示：数组元素及子数组的各种标识和寻访格式；冒号的使用；end的作用。 A=zeros(2,6) %创建(2×6)的全零数组 A(:)=1:12 %赋值号左边：单下标寻访(2×6) 数组A的全部12个元素 %赋值号右边：拥有12个元素的一维数组 A(2,4) %双下标：A数组的第2行第4列元素 A(8) %单下标：数组A的第8个元素 A(: , [1,3]) %双下标：显示A的“第1列和第3列上全部行的元素” A([1, 2, 5, 6]') %单下标：把A数组第1,2,5,6个元素排成列向量 A(: , 4:end) %双下标：显示A的“从第4起到最后一列上全部行的元素” %在此end用于“列标识”，它表示“最后一列” A(2,1:2:5)=[-1, -3, -5] %把右边的3个数分别赋向A数组第2行的第1,3,5个元素位置 B=A([1, 2, 2, 2], [1, 3, 5]) %取A数组的1,3,5列的第1行元素作为B的第1行 %取A数组的1,3,5列的第2行分别作为B的第2,3,4行 L=A<3 %产生与A维数相同的“0，1”逻辑数组 A(L)=NaN %把逻辑1标识的位置上的元素赋为“非数” 运行结果： A = 0 0 0 0 0 0 0 0 0 0 0 0

Matlab实验指导书(含答案)汇总

实验一：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/2 3．两个函数之间的操作 求和 (1) sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5 乘积 (1) exp(-x)*sin(x)

matlab实验指导答案详解(非常详细正确)

实验一 MATLAB 工作环境熟悉及简单命令的执行 一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。 二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单 （Menu ）、工具条 （Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作 空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、启动MATLAB ，熟悉MATLAB 的桌面。 2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。 （1）（365-52?2-70）÷3 >>(365-52*2-70)/3 ans = 63.6667 （2）>>area=pi*2.5^2 area = 19.6350 （3）已知x=3，y=4，在MATLAB 中求z ： () 2 3 2y x y x z -= >>x=3 >>y=4 >>z = x ^2 * y ^3 / (x - y) ^2 z = 576 （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。 m1=???? ? ???? ???11514412679810115133216 执行以下命令 >>m1 =[16 2 3 13 ; 5 11 10 8 ; 9 7 6 12 ; 4 14 15 1 ] >>m1( 2 , 3 ) ans = 10 >>m1( 11 ) ans = 6 >>m1( : , 3 ) ans =3 10 6 15 >>m1( 2 : 3 , 1 : 3 ) ans =5 11 10 9 7 6 >>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) ans = 34 （5）执行命令>>help abs

Matlab数学实验一2015(标准答案版)

Matｌab数学实验一——matlab初体验 一、实验目的及意义 ［1] 熟悉MAＴＬＡB软件的用户环境； ［2] 了解MAＴLAB软件的一般目的命令； ［３] 掌握MATLＡB数组操作与运算函数； 通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。 二、实验内容 １.认识ｍatｌab的界面和基本操作 2.了解mａtlab的数据输出方式（ｆormat) 3. MATLAＢ软件的数组(矩阵）操作及运算练习; 三、实验任务 根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式和要求填写实验报告 1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。 iｊeps inf nａn pi realｍaｘrealｍin 2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。 3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。 4. 写出在命令窗口中的计算步骤和运行结果。 （1)计算 1.22 10 (ln log) 81 e ππ +- ； >＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns = ０．０3４8 (2) >> x=2;y＝４； >> z=x^２+exｐ(x+y)-y*lｏg(x)-3 z ＝ 40１.65６2 （3)输入变量 13 5.3, 25 a b ?? ==?? ?? ，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0 １.ｍat”文件。测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3 ａ＝

Matlab实验第一次实验答案

实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找 sqrt（开方）函数的使用方法； 解：sqrt Square root Syntax B = sqrt(X) Description B = 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. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i

0 + 1.0000i 1.0000 1.4142 2、矩阵运算 （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

Matlab实验四答案

实验四采用SIMULINK的系统仿真 一、实验目的： 1、熟悉SIMULINK 工作环境及特点 2、掌握线性系统仿真常用基本模块的用法 3、掌握SIMULINK 的建模与仿真方法 4、子系统的创建和封装设计 二、实验基本知识： 1．了解SIMULINK模块库中各子模块基本功能

2. SIMULINK 的建模与仿真方法 （1）打开模块库，找出相应的模块。鼠标左键点击相应模块，拖拽到模型窗口中即可。（2）创建子系统：当模型大而复杂时，可创建子系统。 （3）设置仿真控制参数。 三、实验容

练习A 1．SIMULINK 仿真实际应用 （1） 双环调速的电流环系统的方框图模型为： 图中参数设为 Ks=44；Ts=0.00167；Ta=0.017；R=1；Tm=0.075；Ce=0.1925；Kt=0.01178；T1=0.049；T2=0.088， (1)在Simulink 集成环境下建立模型，在给定信号作用点处输入单位给定阶跃响应信号，0.3秒后在扰动信号点输入单位阶跃响应信号。并绘制相应的响应曲线 (2)计算仿真结果的超调量、上升时间、调节时间、稳态误差。 (3)设计PID 调节器替代图中的比例积分调节器，调节P K ， i T ，d T 用使系统满足超调量15%，上升时间0.3s ，调节时间0.4s 的要求。 (4)要求对加入的PID 控制器封装成一个模块使用 Gain3

To Work s pace Warning: Unable to reduce the step size without violating minimum step size of 1.015e-015 at time 0.0001. Continuing simulation with the step size restricted to 1.015e-015 and using an effective relative error tolerance of 0., which is greater than the specified relative error tolerance of 0.001. >> wentai=y(end) chaotiao=100*(max(y)-wentai)/wentai INDEXtm=find(max(y)==y); Tm=t(INDEXtm) ; TT=t(find((abs(y-wentai)/wentai)>0.05)); Ts=max(TT) m=length(y);

MATLAB实验题答案

1、求以下变量的值,并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ;a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算就是否合法,为什么？如合法, 结果就是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree、 >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a、*b result6 = 2 8 -3 4 1 5 30 >> result7=a、/b result7 = 0、5000 0、5000 -3、0000 4、0000 1、6667 1、2000 >> result8=a、c Attempt to reference field of non-structure array、 >> result9=a、\b result9 = 2、0000 2、0000 -0、3333 0、2500 0、6000 0、8333 >> result10=a、^2 result10 = 1 4 9 16 25 36 >> result11=2、^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。

MATLAB实验题答案

1、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A

Matlab数学实验一2015答案版

Matlab数学实验一——matlab初体验 一、实验目的及意义 [1] 熟悉MATLAB软件的用户环境; [2] 了解MATLAB软件的一般目的命令; [3] 掌握MATLAB数组操作与运算函数; 通过该实验的学习,使学生能熟悉matlab的基础应用,初步应用MATLAB软件解决一些简单问题。 二、实验内容 1.认识matlab的界面与基本操作 2.了解matlab的数据输出方式(format) 3、 MATLAB软件的数组(矩阵)操作及运算练习; 三、实验任务 根据实验内容与步骤,完成以下具体实验,要求写出实验报告(实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式与要求填写实验报告 1.在command window中分别输入如下值,瞧它们的值等于多少,并用matlab的help中查询这些缺省预定义变量的含义,用中文写出它们的意义。 i j eps inf nan pi realmax realmin 2.分别输入一个分数、整数、小数等,(如:a=1/9),观察显示结果,并使用format 函数控制数据的显示格式,如:分别输入format short、format long、format short e、format long g、format bank、format hex等,然后再在命令窗口中输入a,显示a的值的不同形式,并理解这些格式的含义。 3.测试函数clear、clc的含义及所带参数的含义(利用matlab的help功能)。 4、写出在命令窗口中的计算步骤与运行结果。 (1)计算 1.22 10 (ln log) 81 e ππ +- ; >> (log(pi)+log(pi)/log(10)-exp(1、2))^2/81 >>ans = 0、0348 (2) >> x=2;y=4; >> z=x^2+exp(x+y)-y*log(x)-3 z = 401、6562 (3)输入变量 13 5.3, 25 a b ?? ==?? ?? ,在工作空间中使用who,whos,并用save命令将变量存入”D:\exe01、 mat”文件。测试clear命令,然后用load命令将保存的”D:\exe01、mat”文件载入>> a=5、3 a =

陈后金信号与系统第二版 第二章matlab

52页例2-13 >> %program2-8 >> t=-3:0.001:3; >> ft=tripuls(2*t,4,0.5); >> ft1=tripuls(2*t,4,0.5); >> subplot(2,1,1) >> plot(t,ft1) >> title('f(2t)') >> ft2=tripuls((2-2*t),4,0.5); >> subplot(2,1,2) >> plot(t,ft2) >> title('f(2-2*t)') 53页例2-14 >> %program2_9 the energy of exponential sequence >> k=0:10; >> A=1;a=-0.6; >> fk=A*a.^k; >> W=sum(abs(fk).^2) W = 1.5625

54页例2-15 function yt = f2_2(t) yt=tripuls(t,4,0.5); %program2_10 differentiation h=0.001;t=-3:h:3; y1=diff(f2_2(t))*1/h; plot(t(1:length(t)-1),y1) title('df(t)/dt') -3-2-10123 -1.2-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 df(t)/dt

%program2_11 integration t=-3:0.1:3; for x =1:length(t) y2(x)=quad('f2_2',-3,t(x)); end plot(t,y2) title('integral of f(t)')

2013MATLAB数学实验答案(全)

第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

MATLAB上机实验1答案

实验1 Matlab 初步 一、 问题 已知矩阵A 、B 、b 如下： ??????????????????-------------=031948 1187638126542 86174116470561091143A ???????? ??????????------=503642237253619129113281510 551201187851697236421B []1187531=b 应用Matlab 软件进行矩阵输入及各种基本运算。 二、 实验目的 学会使用Matlab 软件构作已知矩阵对应的行（列）向量组、子矩阵及扩展矩阵，实施矩阵的初等变换及线性无关向量组的正交规范化，确定线性相关相关向量组的一个极大线性无关向量组，且将其余向量用极大线性无关向量组线性表示，并能编辑M 文件来完成所有的实验目的。 三、 预备知识 1、 线性代数中的矩阵及其初等变换、向量组的线性相关性等知识。 2、 Matlab 软件的相关命令提示如下； （1） 选择A 的第i 行做一个行向量：ai=A(i,:)； （2） 选择A 的第j 行做一个列向量：ai=A(j,:)； （3） 选择A 的某几行、某几列上的交叉元素做A 的子矩阵：A([行号],[列号])； （4） n 阶单位阵：eye(n)；n 阶零矩阵：zeros(n)； （5） 做一个n 维以0或1为元素的索引向量L ，然后取A(:,L)，L 中值为1的对 应的列将被取到。 （6） 将非奇异矩阵A 正交规范化，orth(A) ；验证矩阵A 是否为正交阵，只需做 A*A'看是否得到单位阵E 。 （7） 两个行向量a1和a2的内积：a1*a2'。 （8） 让A 的第i 行与第j 列互换可用赋值语句：A([i,j],:)=A([j,i],:)；

MATLAB全部实验及答案

实验项目序号1 MATLAB基本操作一、实验目的 熟悉MATLAB软件环境，掌握命令窗口的使用。 二、实验基本知识 1、熟悉MATLAB的桌面工具和环境 2、Start菜单 三、主要仪器及耗材

计算机 四、实验内容和步骤 1、MATLAB的启动及操作界面认识 说明有哪几种启动的方法；操作界面有哪些窗口；对常用的窗口作些简单的介绍。 ?使用Windows的“开始”菜单。 ?运行MATLAB系统的启动程序matlab.exe。 ?利用桌面上的快捷方式。 ?MATLAB主窗口 ?命令窗口(Command Window) ?当前目录窗口(Current Directory) ?工作空间窗口(Workspace) ?命令历史窗口(Command History) 2、下列符号中哪些是MATLAB的合法变量名？用给变量赋值的方法在机器上验证你的 答案（3vars、global、help、My_exp、sin、X+Y、_input、E-4、AbCd、AB_C_。）help，My_exp，sin，AbCd，AB_C_； 3、命令窗口的简单使用 1）简单矩阵的输入 1 2 3 A = 4 5 6 7 8 9 2）求[12+2×(7-4)]÷32的算术运算结果 4、有关向量、矩阵或数组的一些运算 1）设A=15；B=20；求C=A+B与c=a+b？ 2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与A.*B？ A*B遵循矩阵运算，A.*B遵循每个对应元素之间进行相乘

3）设A=10，B=20；求C=A/B与D=A\B？

Matlab实验题目与答案

（温馨提示：实验课结束后，请将所有作业（题目、代码、结果）利用word 整 理成一个完整的实验报告，加上封面，打印，纸质档于18周周一交） 第一次上机作业 目的： 1. 掌握MATLAB 各种表达式的书写规则 2. 运行课堂上讲过的例子，熟悉矩阵、表达式的基本操作和运算。 作业： 1. 熟悉matlab 集成环境界面。回答以下问题，并操作相关的指令： （1） 分别写出清除命令窗口和清除变量的指令。 答： clc 和clear （2）在命令行输入命令后，matlab 的搜索过程是怎样的？ 答： （1）检查该命令对象是不是一个变量。 （2）检查该命令对象是不是一个内部函数。 （3）检查该命令对象是否为当前目录下的程序文件。 （4）检查该命令对象是否为MATLAB 搜索路径中其他目录下的M 文件。 （3）什么是matlab 的当前工作目录？写出两种设置当前工作目录的方法？ 答： 就是matlab 当前文件读取和存储的默认路径 （1）在当前目录窗口中更改 （2）在MATLAB 桌面工具栏中更改 （3）使用cd 命令:cd c:\mydir---将c ：\mydir 设置为当前目录 （4）什么是matlab 的搜索路径？写出两种设置搜索路径的方法？ 答： 指Matlab 运行文件时进行搜索的目录。 （1）用path 命令设置： （2）用Set Path 对话框设置 （5）help 命令和doc 命令有什么作用，它们有什么区别？ 答： help 命令：最基本的帮助命令，查询信息直接显示在命令窗口。 doc 命令：在帮助窗口中显示HTML 帮助文档，显示函数的详细用法及 例子，比help 命令更详细。 2. 在matlab 中输入下列表达式，并求各表达式的值，显示MATLAB 工作空间的使用情况 并用两种方式保存全部变量，变量保存的文件名必须包含自己的学号后四位数： （1）)1034245.01(26-?+?=w w=sqrt(2)*(1+0.34245*10^-6) w = 1.4142 （2）,)tan(22a c b e abc c b a x ++-+++ =ππ 其中a=3.5,b=5,c=9.8。 a=3.5;b=5;c=9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) x =

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

《MATLAB 语言及其用》 实验指导书 目录

实验一Matlab 使用方法和程序设计........................ 实验二控制系统的模型及其转换............................. 实验三控制系统的时域、频域和根轨迹分析........... 实验四动态仿真集成环境-Simulink......................... 实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找 sqrt（开方）函数的使用方法； 在 CommandWindowL里输入help，接在在search里输入sqr即可。 sqrt Square root Syntax B = sqrt(X) Description B = 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.

Tips See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i 1.0000 1.4142 See Also nthroot | realsqrt | sqrtm 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B A=[1 2;3 4];B=[5 5;7 8]; C=A^2*B >> format compact C = 105 115 229 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/B A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; C=A\B,D=A/B C = 1.0e+016 * 0.3152 -1.2609 0.9457

Matlab实验五答案

实验五 控制系统的频域与时域分析 一、实验目的： 1、掌握控制系统数学模型的基本描述方法和相互转化 2、了解控制系统的稳定性分析方法 3、掌握控制系统频域与时域分析基本方法 二、实验基本知识： 1. 系统数学模型的几种表示方法 方案一 函数转换表 [z ,p ,k]=tf2zp(num , den); [num , den ]=zp2tf(z ,p ,k) [num , den]=ss2tf(A ,B ,C ,D , iu ) %iu —表示输入的序号（对多输入系统） [z, p ,k]=ss2zp(A ,B ,C ,D , iu ) %[z ,p ,k]— 表示对第iu 个输入信号的传递函数的零极点。 [A , B ,C ,D]= tf2ss (num , den) [A , B ,C ,D]=zp2ss(z ,p ,k) 2、在MATLAB 工具箱中，提供了子系统的连接处理函数： 1） series()函数：系统串联实现 格式：sys=series(1sys ,2sys ) [A,B,C,D]=series(A 1,B 1,C 1,D 1, A 2,B 2,C 2,D 2) 2） parallel()函数：系统并联实现。 格式：sys=parallel(1sys ,2sys ) 3） feedback()函数：系统反馈连接

格式：sys=feedback(1sys ,2sys ) sys=feedback(1sys ,2sys ,sign) sign —定义反馈形式：正反馈，sign=+1；负反馈，sign=-1 3.控制系统根轨迹绘制 1 rlocfind():计算给定根的根轨迹增益 格式：①[k,poles]=rlocfind(sys) ②[k,poles]= rlocfind(sys,p) ① 在LTI 对象的根轨迹图中显示出十字光标，当用户选择其中一点时，其相应的增益 由k 记录，与增益相关的所有极点记录poles 中。若要使用该函数，必须首先在当前窗口上绘制系统的根轨迹。 ② 定义要得到增益的根矢量P ，即事先给出极点。除了显示出该根对应的增益以外， 还显示出该增益对应的其它根。 2 rlocus() 函数：功能为求系统根轨迹。 格式：rlocus(sys)--- 计算SISO 开环LTI 对象的根轨迹，增益自动选取。 rlocus(sys,k)---显式设置增益。 [r,k]= rlocus(sys)--- 返回系统的增益k （向量）和闭环极点r （向量），即对应于增益的闭环极点为r(i)。 4.控制系统频域分析基本方法 （1）Bode 图：bode()函数 （2）Nyquist 图：nyquist()函数 （3）稳定裕度计算：margin( )函数 综合应用： A=[0 1 –1；-6 –11 6；-6 –11 5]； B=[0 0 1]； C=[1 0 0]； D=0； ω=logspace(-1,1)； [m,p]=bode(A,B,C,D,1,ω) subplot(211),semilogx(ω,20*log10(m))； subplot(212),semilogx(ω,p)； [X,Y]=nyquist(A,B,C,D,1,ω); plot(X,Y) [m G ，m P ，cg ω，CP ω]=margin(A,B,C,D) 5.线性系统时间响应分析 step( )函数---求系统阶跃响应 impulse( )函数：求取系统的脉冲响应 lsim( )函数：求系统的任意输入下的仿真 三、实验内容 练习A ： 1、表示下列传递函数模型,并转化成其他的数学模型