MATLAB作业（八章）
8.1某随动系统的微分方程为
0.07523dt d ?+0.7522dt
d ?+dt d ?
+K ?=K ψ
已知：?(0)=??(0)=?
??(0)=0，ψ(t)=1(t)。试分别就K=2.5、5.0、12.5三种情况对系统进行仿真，考察?(t)的动态性能。
8.2已知系统模型
?????+--=+=?
?
pu
x x x pu
x x 212212, x 1(0)=x 2(0)=0 当p=1、2、10时，对以下几种情况进行仿真，并比较不同输入幅值下的系统输出响应。 (1)
u （t ）=1； （2）u （t ）=t ； （3）u （t ）=sin （t ）； （4）
u （t ）=1+sin （t ）； （5）u(t)=1+t+sin(t).
(1)
(2)
(3)
(4)
(5)
8.3 对如图所示的系统进行仿真。 （1）输入信号为r （t ）=1（t ）； （2）输入信号为如下分段函数 s
5.0t t 1s
5.0t t 2{
t r ><==）（）（
(1)
(2).
f=-1的阶跃扰动，试就系统在该扰动作用下IV环节的动态过程。
8.8非线性控制系统如图所示，求系统在有饱和非线性和无饱和非线性两种情况下的动态响应，分析饱和非线性对系统动态性能的影响。
(1).
(2)
第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1)
ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t
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）
第8章M文件函数 使用MATLAB函数时，例如inv, abs, angle和sqrt，MATLAB获取传递给它的变量，利用所给的输入，计算所要求的结果。然后，把这些结果返回。由函数执行的命令，以及由这些命令所创建的中间变量，都是隐含的。所有可见的东西是输入和输出，也就是说函数是一个黑箱。 这些属性使得函数成为强有力的工具，用以计算命令。这些命令包括在求解一些大的问题时，经常出现的有用的数学函数或命令序列。由于这个强大的功能，MATLAB提供了一个创建用户函数的结构，并以M文件的文本形式存储在计算机上。MATLAB函数fliplr 是一个M文件函数良好的例子。 function y = fliplr(x) % FLIPLR Flip matrix in the left/right direction. % FLIPLR(X) returns X with row preserved and columns flipped % in the left/right direction. % % X = 1 2 3 becomes 3 2 1 % 4 5 6 6 5 4 % % See also FLIPUD, ROT90. % Copyright (c) 1984-94 by The MathWorks, Inc. [m, n] = size(x); y = x(: , n : -1 : 1); 一个函数M文件与脚本文件类似之处在于它们都是一个有.m扩展名的文本文件。如同脚本M文件一样，函数M文件不进入命令窗口，而是由文本编辑器所创建的外部文本文件。一个函数的M文件与脚本文件在通信方面是不同的。函数与MATLAB工作空间之间的通信，只通过传递给它的变量和通过它所创建的输出变量。在函数内中间变量不出现在MATLAB工作空间，或与MATLAB工作空间不交互。正如上面的例子所看到的，一个函数的M文件的第一行把M文件定义为一个函数，并指定它的名字。它与文件名相同，但没有.m扩展名。它也定义了它的输入和输出变量。接下来的注释行是所展示的文本，它与帮助命令：? help fliplr相对应。第一行帮助行称为H1 行，是由lookfor命令所搜索的行。最后，M文件的其余部分包含了MATLAB创建输出变量的命令。 8.1 规则和属性
《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书 班级：T1243-7 姓名：柏元强 学号：20120430724 总评成绩： 汽车工程学院 电测与汽车数字应用中心
目录 实验04051001 MATLAB语言基础 (1) 实验04051002 MATLAB科学计算及绘图 (18) 实验04051003 MATLAB综合实例编程 (31)
实验04051001 MATLAB语言基础 1实验目的 1)熟悉MATLAB的运行环境 2)掌握MATLAB的矩阵和数组的运算 3)掌握MATLAB符号表达式的创建 4)熟悉符号方程的求解 2实验内容 第二章 1.创建double的变量，并进行计算。 (1)a=87，b=190，计算 a+b、a-b、a*b。 clear,clc a=double(87); b=double(190); a+b,a-b,a*b (2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。 clear,clc a=uint8(87); b=uint8(190); a+b,a-b,a*b 2．计算：
(1) () sin 60 (2) e3 (3) 3cos 4??π ??? clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4) 3．设2u =，3v =，计算： (1) 4 log uv v (2) () 2 2 e u v v u +- (3) clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v) 4．计算如下表达式： (1) ()() 3542i i -+ (2) () sin 28i - clear,clc (3-5*i)*(4+2*i) sin(2-8*i)
第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？
第一大题： （1） a = 7/3 b = sym(7/3) c = sym(7/3,'d') d = sym('7/3') v1=vpa(abs(a-d)) v2=vpa(abs(b-d)) v3=vpa(abs(c-d)) a = 2.3333 b = 7/3 c = 2.3333333333333334813630699500209 d = 7/3 v1 = 0.0 v2 = 0.0 v3 = 0.00000000000000014802973661668756666666667788716（2） a = pi/3 b = sym(pi/3) c = sym(pi/3,'d') d = sym('pi/3') v1=vpa(abs(a-d)) v2=vpa(abs(b-d)) v3=vpa(abs(c-d)) a = 1.0472 b = pi/3 c = 1.047197551196597631317786181171 d = pi/3 v1 = 0.0 v2 = 0.0 v3 = 0.00000000000000011483642827992216762806615818554 (3) a = pi*3^(1/3) b = sym(pi*3^(1/3)) c = sym(pi*3^(1/3),'d') d = sym('pi*3^(1/3)') v1=vpa(abs(a-d)) v2=vpa(abs(b-d)) v3=vpa(abs(c-d))
a = 4.5310 b = 1275352044764433/281474976710656 c = 4.5309606547207899041040946030989 d = pi*3^(1/3) v1 = 0.00000000000000026601114166290944374842393221638 v2 = 0.00000000000000026601114166290944374842393221638 v3 = 0.0000000000000002660111416629094726767991785515 第二大题： （1） c1=3/7+0.1 c1 = 0.5286 双精度 （2） c2=sym(3/7+0.1) c2 = 37/70 符号 （3） c3=vpa(sym(3/7+0.1)) c3 = 0.52857142857142857142857142857143 完整显示精度 第三大题： （1） findsym(sym('sin(w*t)'),1) ans = w (2) findsym(sym('a*exp(-X)' ) ,1) ans = a (3) findsym(sym('z*exp(j*theta)'),1) ans = z
第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分 diff或数值梯度gradient指令计算)(t y'曲线 y'，然后把)(t y和)(t 绘制在同一图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1);
yc=gradient(y)/dt; plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =
美赛数学建模专用- 第八章MATLAB的GUI 程序设计 Chapter 8: Design of MATLAB of GUI program GUI(Graphical User Interfaces):由各种图形对象组成的用户界面，在这种用户界面下，用户的命令和对程序的控制是通过“选择”各种图形对象来实现的。 基本图形对象分为控件对象和用户界面菜单对象，简称控件和菜单。 一.控件对象及属性(Object and its attributes of controller)) 1. GUI控件对象类型(The mode of controller object) 控件对象是事件响应的图形界面对象。当某一事件发生时，应用程序会做出响应并执行某些预定的功能子程序（Callback）. 控件对象及其功能：（表7—1） 2.控件对象的描述(Description of controller object) MATLAB中的控件大致可分为两种，一种为动作控件，鼠标点击这些控件时会产生相应的响应。一种为静态控件，是一种不产生响应的控件，如文本框等。 每种控件都有一些可以设置的参数，用于表现控件的外形、功能及效果，既属性。属性由两部分组成：属性名和属性值，它们必须是成对出现的。 （1）按钮(Push Buttons)：执行某种预定的功能或操作； （2）开关按钮(Toggle Button)：产生一个动作并指示一个二进制状态（开或关），当鼠点击它时按钮将下陷，并执行 callback（回调函数）中指定的内容，再次点击，按钮复 原，并再次执行callback 中的内容； （3）单选框(Radio Button)：单个的单选框用来在两种状态之间切换，多个单选框组成一个单选框组时，用户只能在 一组状态中选择单一的状态，或称为单选项； （4）复选框(Check Boxes)：单个的复选框用来在两种状态之
第4章 MATLAB程序流程控制 习题4 一、选择题 1．下列关于脚本文件和函数文件的描述中不正确的是（）。 A A．函数文件可以在命令行窗口直接运行 B．去掉函数文件第一行的定义行可转变成脚本文件 C．脚本文件可以调用函数文件 D．函数文件中的第一行必须以function开始 2．下列程序的输出结果是（）。D y=10; if y==10 y=20; elseif y>0 y=30 end disp(y) A．1 B．30 C．10 D．20 3．有以下语句： a=eye(5); for n=a(2:end,:) for循环的循环次数是（）。C A．3 B．4 C．5 D．10 4．设有程序段 k=10; while k k=k-1 end 则下面描述中正确的是（）。A A．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段： x=reshape(1:12,3,4); m=0;
n=0; for k=1:4 if x(:,k)<=6 m=m+1; else n=n+1; end end 则m和n的值分别是（）。C A．6 6 B．2 1 C．2 2 D．1 2 6．调用函数时，如果函数文件名与函数名不一致，则使用（）。A A．函数文件名 B．函数名 C．函数文件名或函数名均可 D．@函数名 7．如果有函数声明行为“function [x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是（）。B A．x=f1(a,b,c) B．[x,y,z,w]=f1(a,b,c) C．[x,b,z]=f1(a,y,c) D．[a,b]=f1(x,y,z) 8．执行语句“fn=@(x) 10*x;”，则 fn是（）。A A．匿名函数 B．函数句柄 C．字符串 D．普通函数 9．执行下列语句后，变量A的值是（）。D >> f=@(x,y) log(exp(x+y)); >> A=f(22,3); A．22,3 B．22 C．3 D．25 10．程序调试时用于设置断点的函数是（）。A A．dbstop B．dbclear C．dbcont D．dbstack 二、填空题 1．将有关MATLAB命令编成程序存储在一个扩展名为.m的文件中，该文件称为。M文件 2．有语句“for k=[12;34]”引导的循环结构，其循环体执行的次数为。1 3．MATLAB中用于控制不确定重复次数的循环语句为，若在循环执行过程中需要终止该循环时采用的语句为。while…end，break 4．函数文件由语句引导。在函数定义时，函数的输入输出参数称为参数，简称。在调用函数时，输入输出参数称为参数，简称。 function，形式，形参，实际，实参 5．在MATLAB中，函数文件中的变量是变量。定义变量是函数间传递信息的一种手段，可以用命令定义。局部，全局，global 6．应用程序的错误有两类，一类是错误，另一类是运行时的错误，即错
作者：夏云木子 1、 >> syms re(x) re(y) re(z) >> input('计算相对误差：'),re(x)=10/1991,re(y)=0.0001/1.991,re(y)=0.0000001/0.0001991 所以可知re(y)最小，即y精度最高 2、 >> format short,A=sqrt(2) >> format short e,B=sqrt(2) >> format short g,C=sqrt(2)
>> format long,D=sqrt(2) >> format long e,E=sqrt(2) >> format long g,F=sqrt(2) >> format bank,H=sqrt(2) >> format hex,I=sqrt(2) >> format +,J=sqrt(2) >> format,K=sqrt(2)
3、 >> syms A >> A=[sqrt(3) exp(7);sin(5) log(4)];vpa(pi*A,6) 4、1/6251-1/6252=1/6251*6252 5、（1）1/（1+3x）-(1-x)/(1+x)=x*(3*x-1)/[(1+3*x)*(1+x)] (2) sqrt(x+1/x)-sqrt(x-1/x)=2/x/[sqrt(x-1/x)+sqrt(x+1/x)] (3) log10(x1)-log(x2)=log10(x1/x2) (4) [1-cos(2*x)]/x =x^2/factorial(2)-x^4/factorial(4)+x^6/factorial(6)-…
第八章 实验指导 1、 >> syms x y; >> A=x^4-y^4; >> factor(A) ans = (x - y)*(x + y)*(x^2 + y^2) >> a=5135; >> factor(a) ans = 5 13 79 2、(1) >> syms x >> f=(x-2)/(x^2-4); >> limit(f,x,2)
ans = 1/4 （2） >> syms x f=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1); limit(f,x,-1,'right') ans = -Inf 3、 （1） >> syms x y >> f=sin(1/x); >> diff(f) ans = -cos(1/x)/x^2 >> diff(f,x,2)
ans = (2*cos(1/x))/x^3 - sin(1/x)/x^4 >> f=(1-cos(2x))/x; (2) >> f=(1-cos(2*x))/x f = -(cos(2*x) - 1)/x >> diff(f) ans = (2*sin(2*x))/x + (cos(2*x) - 1)/x^2 >> diff(f,x,2) ans =
(4*cos(2*x))/x - (4*sin(2*x))/x^2 - (2*(cos(2*x) - 1))/x^3 4、 (1)、 >> x=sym('x'); >> f=sqrt(exp(x)+1); >> int(f) ans = atan((exp(x) + 1)^(1/2)*i)*2*i + 2*(exp(x) + 1)^(1/2) (2) >> x=sym('x'); y=sym('y'); f=x/(x+y); >> int(f,y) ans = x*log(x + y) (3) >> syms x;
深圳大学实验报告课程名称：可视化计算语言 实验项目名称：Matlab矩阵及其运算 学院：医学院 专业：生物医学工程 指导教师： 报告人：学号：班级： 实验时间：2012-9-26 实验报告提交时间： 教务部制
实验目的与要求： 1．掌握MATLAB 数据对象的特点以及数据的运算规则。 2．掌握MATLAB 中建立矩阵的方法以及矩阵处理的方法。 3．掌握MATLAB 分析的方法。 实验过程及内容： 1．求下列表达式的值。 （1）w=2×(1+0.34245×10-6) （2）a c b abc c b a x ++-+++ =)tan(e 22ππ，其中a=3.5，b=5，c=-9.8。 （3），])4 8333.0()41[(22απβπαy ---=π其中α=3.32，β=-7.9。 （4）)1ln(2122t t e z t ++=，其中t =??????--65.05i 312。 2．已知 A = ??????????--7613870451， B =???? ??????--023352138 求下列表达式的值： （1）A+6B 和A 2-B+I （其中I 为单位矩阵）。 （2）A*B 、A.*B 和B*A 。 （3）A/B 及B\A 。 （4）[A,B]和[A([1,3],:);B^2]。 3．已知 ????????????---=14.35454 .9632053256545410778.01023Α 完成下列操作： （1）输出A 在[10，25]范围内的全部元素。 （2）取出A 前3行构成矩阵B ，前两列构成矩阵C ，右下角3×2子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。 （3）分别求表达式E matlab 第四章作业 4.2编写函数文件求半径为r 的圆的面积和周长 步骤：1.首先建立函数文件，命名为Abb.m 程序如下： function [s,c]=Abb(r) s=pi*r*r;%求圆的面积 c=2*pi*r;%求圆的周长 end 在工作空间中调用Abb 求圆的面积和周长（将半径设为3） 程序结果如下 [s c]=Abb(3) s = 28.2743 c = 18.8496 4.3已知y=2 2221312111n ++++ 。当n=100时，求y 的值 步骤：建立函数，命名为Cbb.m 求平方倒数的和。其程序如下 function F=Cbb(n) F=0;k=1; while k<=n F=F+(1/k)/k; k=k+1; end end 在工作空间中求y 值 y= Cbb(100) y= 1.6350 4.5编制程序，计算1+2+ +n<2000时的最大n 值 步骤：建立函数文件Dbb.m,求累加和小于n 的最大数，使n=2000程序如下 function k=Dbb(n) k=1; f=0; while f k=k-2; end 将2000作为输入值，求n.在工作空间中输入如下： Dbb(2000) ans = 62 4.6分别用for 和while 循环结构编写程序，求出K=∑=6302 i i 分别用for 和while 循环结构编写命令文件Ebb,Fbb.并在工作空间中求K Ebb 程序如下 f=1; J(1)=1; for i=1:63 J(i+1)=J(i)*2; f=f+J(i+1); End 求K 程序： Ebb.m f f = 1.8447e+019 Fbb 程序如下 I=1; H=1; i=1; while i<=63 H=H*2; I=I+H; i=i+1; end 在工作空间中求K Fbb.m I I = 1.8447e+019 第2章--M A T L A B数据及其运算-习题答案 第2章 MATLAB数据及其运算 习题2 一、选择题 1．下列可作为MATLAB合法变量名的是（）。D A．合计 B．123 C．@h D．xyz_2a 2．下列数值数据表示中错误的是（）。C A．+10 B．1.2e-5 C．2e D．2i 3．使用语句t=0:7生成的是（）个元素的向量。A A．8 B．7 C．6 D．5 4．执行语句A=[1,2,3;4,5,6]后，A(3)的值是（）。B A．1 B．2 C．3 D．4 5．已知a为3×3矩阵，则a(:,end)是指（）。D A．所有元素 B．第一行元素 C．第三行元素 D．第三列元素 6．已知a为3×3矩阵，则运行a (1)=[]后（）。A A．a变成行向量 B．a变为2行2列 C．a变为3行2列 D．a变为2行3列 7．在命令行窗口输入下列命令后，x的值是（）。B >> clear >> x=i*j A．不确定 B．-1 C．1 D．i*j 8．fix(354/100)+mod(354,10)*10的值是（）。D A．34 B．354 C．453 D．43 9．下列语句中错误的是（）。B A．x==y==3 B．x=y=3 C．x=y==3 D．y=3,x=y 10．find(1:2:20>15)的结果是（）。C A．19 20 B．17 19 C．9 10 D．8 9 11．输入字符串时，要用（）将字符括起来。C A．[ ] B．{ } C．' ' D．" " 12．已知s='显示"hello"'，则s的元素个数是（）。A A．9 B．11 C．7 D．18 第八章 层次分析法 层次分析法（Analytic Hierarchy Process ，简称AHP ）是对一些较为复杂、较为模糊的问题作出决策的简易方法，它特别适用于那些难于完全定量分析的问题。它是美国运筹学家T. L. Saaty 教授于70年代初期提出的一种简便、灵活而又实用的多准则决策方法。 §1 层次分析法的基本原理与步骤 人们在进行社会的、经济的以及科学管理领域问题的系统分析中，面临的常常是一个由相互关联、相互制约的众多因素构成的复杂而往往缺少定量数据的系统。层次分析法为这类问题的决策和排序提供了一种新的、简洁而实用的建模方法。 运用层次分析法建模，大体上可按下面四个步骤进行： （i ）建立递阶层次结构模型； （ii ）构造出各层次中的所有判断矩阵； （iii ）层次单排序及一致性检验； （iv ）层次总排序及一致性检验。 下面分别说明这四个步骤的实现过程。 1.1 递阶层次结构的建立与特点 应用AHP 分析决策问题时，首先要把问题条理化、层次化，构造出一个有层次的结构模型。在这个模型下，复杂问题被分解为元素的组成部分。这些元素又按其属性及关系形成若干层次。上一层次的元素作为准则对下一层次有关元素起支配作用。这些层次可以分为三类： （i ）最高层：这一层次中只有一个元素，一般它是分析问题的预定目标或理想结果，因此也称为目标层。 （ii ）中间层：这一层次中包含了为实现目标所涉及的中间环节，它可以由若干个层次组成，包括所需考虑的准则、子准则，因此也称为准则层。 （iii ）最底层：这一层次包括了为实现目标可供选择的各种措施、决策方案等，因此也称为措施层或方案层。 递阶层次结构中的层次数与问题的复杂程度及需要分析的详尽程度有关，一般地层次数不受限制。每一层次中各元素所支配的元素一般不要超过9个。这是因为支配的元素过多会给两两比较判断带来困难。 下面结合一个实例来说明递阶层次结构的建立。 例1 假期旅游有1P 、2P 、3P 3个旅游胜地供你选择，试确定一个最佳地点。 在此问题中，你会根据诸如景色、费用、居住、饮食和旅途条件等一些准则去反复比较3个侯选地点。可以建立如下的层次结构模型。 目标层O 选择旅游地 准则层C 景色 费用 居住 饮食 旅途 措施层P 1P 2P 3P 1.2 构造判断矩阵 第八章 exercize 2 x=[480 500 380 1100 1100 490 230 250 300 510 1300]; >> y=[180 150 170 350 460 240 60 90 110 250 200]; >> plot(x,y,'*') >> c=corrcoef(x,y)； exercize 7(1) a=1/sqrt(2*pi)% 被积函数的最大值 a = 0.3989 n=10000;x=rand(n,1);y=unifrnd(0,0.4,n,1); c1=y>0;c2=y 第九章 exercize 3(1) >>x=linspace(0,2*pi,10); >> y=sin(x); >> xi=linspace(0,2*pi,50);yi=sin(xi); >> yi1=interp1(x,y,xi); >> yi2=interp1(x,y,xi,'spline'); >> plot(xi,yi1,'r',xi,yi2,'b',x,y,'*') e1=norm(yi1-yi),e2=norm(yi2-yi) %均方误差 e1 = 0.3097 e2 = 0.0158 exercize 3(4) x=linspace(-2,2,10);y=exp(-x.^2); >> xi=linspace(-2,2,80); >> yi1=interp1(x,y,xi);yi2=interp1(x,y,xi,'spline'); >> yi=exp(-xi.^2); >> plot(xi,yi1,'r',xi,yi2,'g',x,y,'*') >> e1=norm(yi1-yi),e2=norm(yi2-yi) e1 = 0.1518 e2 = 0.0069 exercize 12 y=[100 75 80 70 50 65 90 100 110 60]'; >> x1=[1000 600 1200 500 300 400 1300 1100 1300 300]'; >> x2=[5 7 6 6 8 7 5 4 3 9]'; >> plot(x1,y,'*') >> plot(x2,y,'o') >> [b,bint,r,rint,stats]=regress(y,[ones(10,1),x1,x2]); 4. (1) syms x positive assumeAlso(x,'integer') f1=x^4-5.1*x^3+58.04*x^2-264.384*x+321.408; r1=solve(f1) (2) syms x clear syms a positive solve(x^2-a*x+a^2) 5. A=sym('[a11,a12,a13;a21,a22,a23;a31,a32,a33]') 或 A=sym('a%d%d',[3,3]) Dt=det(A) Iv=inv(A) [Ivs,d]=subexpr(Iv,'d') 6. syms x k f=x^k; s=symsum(f,k,0,inf) subs(s,x,{sym('-1/3'),sym('1/pi'),sym('3')}) 7. syms k syms x positive f=2/(2*k+1)*((x-1)/(x+1))^(2*k+1); s=simple(symsum(f,k,0,inf)) 8. syms t y=abs(sin(t)) d=diff(y) d0=limit(d,t,0,'left') d1=limit(d,t,pi/2) (or subs(d,t,sym('pi/2'))) 9. syms x clear syms x f=exp(-abs(x))*abs(sin(x)) vpa(int(f,x,-5*pi,1.7*pi),64) 10. syms x y int(int(f,y,1,x^2),x,1,2) 11. syms t x f=sin(t)/t y=int(f,t,0,x) y5=subs(y,x,sym('4.5')) ezplot(y,[0,2*pi]) 12. syms x syms n positive f=sin(x)^n y=int(f,x,0,pi/2) ys=vpa(subs(y,n,sym('1/3'))) yd=vpa(subs(y,n,1/3)) 13. syms a b k n x=a^k; h=b^k; f=subs(h,k,n)*subs(x,k,k-n) y=symsum(f,n,0,k) 15. syms A t w syms a positive ft=A*exp(-a*abs(t)); Fw=fourier(ft,t,w) 16. syms t A w syms tao positive ft=A*((1+t/tao)*(heaviside(t+tao)-heaviside(t))+(1-t/tao)*(heaviside(t)-heaviside(t-tao))); Fw=simple(fourier(ft,t,w)) Fw1=simple(subs(Fw,[A,tao],[2,2])) ezplot(Fw1) 17. syms s t Fs=(s+3)/(s^3+3*s^2+6*s+4) ft=ilaplace(Fs, s,t) 18. syms t s y=sym('f(t)') 习题三 1.选择题 (1)已知a=4，b='4'，下面说法错误的为C 。 A. 变量a比b占用的存储空间大 B. 变量a和b可以进行加、减、乘、除运算 C. 变量a和b的数据类型相同 D. 变量b可以用eval命令执行 (2)已知s='显示"hello"'，则s的元素个数是B 。 A. 12 B. 9 C. 7 D. 18 (3)运行命令>>a=sym('pi','d')则对于变量a的描述D是正确的。 A. a是符号变量 B. a显示为10位的数值 C. a显示为32位的数值 D. a不存在 (4)对于符号表达式g=sym('sin(a*z)+cos(w*v)')，如果运行命令diff(g)，则求导自变量是D。 A. a B. z C. w D. v (5)运行命令>>a=double(sym('sin(pi/2)'))，则变量a是 C。 A. 符号变量 B. 字符串'1' C. double型的数值1 D. 出错 (6)运行命令>>y=dsolve('x*D2y-3*Dy=x^2','t')求解微分方程，则D。 A. Dy是指dy/dx B. 得出的y是通解有一个常数符C1 C. D2y是指d2y/dx D. 得出的y是通解有两个常数符C1和C2 (7)有一个2行2列的元胞数组c，则c(2)是指C 。 A. 第1行第2列的元素内容 B. 第2行第1列的元素内容 C. 第1行第2列的元素 D. 第2行第1列的元素 (8)对于一个2行2列的元胞数组c，如果要删除第二列的元素，应该使用的命令是B 。 A.c{:,2}=[] B.c(:,2)=[] C.A、B两种方式均可 D.A、B两种方式均不可 (9)有一个3行10列的构架数组student，每个构架有name（姓名）、scores两个字段，其中scores是以1×5数组表示的五门课的成绩，那么要查看第4个学生的第2门课成绩，采用的命令中达不到要求的是。 A.student(1,2).scores(2) B.student(4).scores(2) C.两种方式均可 D.两种方式均不可 (10)构架数组结构同上，如果增添一个字段ave用来存放每个同学的平均成绩，则要查看所有同学的平均成绩，应该采用的命令是D。 A.[student.ave] B.{student.ave} C.两种方式均可 D.两种方式均不可 (11)运行以下命令： >>x=[1 2 3;4 5 6]; >>y=x+x*i >>plot(y) 则在图形窗口中绘制A条曲线。 A. 3 B. 2 C. 6 D. 4 (12)运行以下命令： >>x=[1 2 3;4 5 6]; >>plot(x,x,x,2*x) 则在图形窗口中绘制B条曲线。 A. 4 B. 6 C. 3 D. 5 (13)subplot(2,1,1)是指A的子图。 A. 两行一列的上图 B. 两行一列的下图 C. 两列一行的左图 D. 两列一行的右图 (14)运行命令>>figure(3)，则执行D 。 A. 打开三个图形窗口 B. 打开一个图形窗口 C. 打开图形文件名为“3.fig” D. 打开图形文件名为“figure 3.fig” (15)运行以下命令：matlab第四章作业
第2章--MATLAB数据及其运算-习题答案教学内容
Matlab建模教程-第八章-层次分析法
matlab数学实验第八章~第十章部分习题答案
MATLAB2014b第二章部分答案
MATLAB第三章和第四章作业
相关主题
文本预览