MATLAB作业1参考答案_2016_
《MATLAB仿真技术》实验指导书2016附答案分析
实验项目及学时安排实验一 MATLAB环境的熟悉与基本运算 2学时实验二 MATLAB数值计算实验 2学时实验三 MATLAB数组应用实验 2学时实验四 MATLAB符号计算实验 2学时实验五 MATLAB的图形绘制实验 2学时实验六 MATLAB的程序设计实验 2学时实验七 MATLAB工具箱Simulink的应用实验 2学时实验八 MATLAB图形用户接口GUI的应用实验 2学时实验一 MATLAB环境的熟悉与基本运算一、实验目的1.熟悉MATLAB开发环境2.掌握矩阵、变量、表达式的各种基本运算二、实验基本知识1.熟悉MATLAB环境:MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下:(1)变量名可以由英语字母、数字和下划线组成(2)变量名应以英文字母开头(3)长度不大于31个(4)区分大小写MATLAB中设置了一些特殊的变量与常量,列于下表。
MATLAB运算符,通过下面几个表来说明MATLAB的各种常用运算符4.MATLAB的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式5.MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、学习使用help命令,例如在命令窗口输入help eye,然后根据帮助说明,学习使用指令eye(其它不会用的指令,依照此方法类推)2、学习使用clc、clear,观察command window、command history和workspace等窗口的变化结果。
3、初步程序的编写练习,新建M-file,保存(自己设定文件名,例如exerc1、exerc2、exerc3……),学习使用MA TLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
MATLAB课后习题集附标准答案
第2章MATLAB概论1、与其他计算机语言相比较,MA TLAB语言突出的特点是什么?答:起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富,可扩展性强.2、MA TLAB系统由那些部分组成?答:开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装?答:在安装MATLAB时,安装内容由选择组件窗口中各复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上?答:在MATLAB操作桌面上有五个窗口,在每个窗口的右下角有两个小按钮,一个是关闭窗口的Close 按钮,一个是可以使窗口称为独立的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面称为独立窗口,在独立窗口的view菜单中选择Dock,菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器?答:在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗?如何操作?答:存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外,还有什么用途?答:命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别?答:当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径?答:(1)帮助浏览器:选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.(2)help命令:在命令窗口键入“help”命令可以列出帮助主题,键入“help 函数名”可以得到指定函数的在线帮助信息.(3)lookfor命令:在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.(4)模糊查询:输入命令的前几个字母,然后按Tab键,就可以列出所有以这几个字母开始的命令和函数.注意:lookfor和模糊查询查到的不是详细信息,通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375,并将其赋予变量a ? 答:在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法?各有什么优点? 答:(1)直接输入法,如a=[1 2 3 4],优点是输入方法方便简捷;(2)通过M 文件建立矩阵,该方法适用于建立尺寸较大的矩阵,并且易于修改; (3)由函数建立,如y=sin(x),可以由MATLAB 的内部函数建立一些特殊矩阵; (4)通过数据文件建立,该方法可以调用由其他软件产生数据.3.3在进行算术运算时,数组运算和矩阵运算各有什么要求?答:进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则,如矩阵a 与b 相乘(a*b )时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。
MATLAB作业1 参考答案-推荐下载
3、试用符号元素工具箱支持的方式表达多项式 f (x) x5 3x4 4x3 2x2 3x 6 ,并 Nhomakorabea令x
s 1
,将f(x)
s 1
替换成s
的函数。
【求解】可以先定义出f 函数,则由subs() 函数将x 替换成s 的函数
>> syms s x;f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;F=subs(f,x,(s-1)/(s+1))
end 其中,前者语句结构简单,但适用范围更广,允许使用矩阵型x,后者只能使用向量型的 x,但不能处理矩阵问题。
63
7、用数值方法可以求出 S 2i 1 2 4 8 262 263 ,试不采用循环的形式求 i0
出 和式的数值解。由于数值方法采用double 形式进行计算的,难以保证有效位数字,所以结 果不一定精确。试采用符号运算的方法求该和式的精确值。 【求解】用符号运算的方式可以采用下面语句 >> sum(sym(2).^[1:63]) ans = 18446744073709551614 由于结果有19 位数值,所以用double 型不能精确表示结果,该数据类型最多表示16 位有 效数字。其实用符号运算方式可以任意保留有效数字,例如可以求200 项的和或1000 项的 和可以由下面语句立即得出。 >> sum(sym(2).^[1:200]) ans = 3213876088517980551083924184682325205044405987565585670602750 >> sum(sym(2).^[1:1000]) 第一种方法 i=0; s=0; for i=0:63
MATLAB上机习题一-2016-解答
4. 在命令窗口输入 ispc,按键盘的 F1 键会弹出帮助,阅读并写出该内置变量的 含义。 (也可通过主菜单中[Help]->[Product Help]菜单打开帮助,并搜索 ispc 来完 成) 解答:
表示 MATLAB 的版本是否是 Windows 系统平台下的版本
5. 观察 MATLAB 命令窗口中关键字、字符串、注释、错误提示的显示颜色。提 示:输入相关表达式来观察,例如字符串可以这样输: ss=’abc’ 解答: 关键字:蓝色 字符串:紫色 注释:绿色需要 MATLAB7.0 以上版本,请确认系统装安装有。第 1、2 机房需要在开机时 选择【APPLY】系统。 截图方法:按键盘上右上角的[Print Screen]键,在作业文档中按 Ctrl+V 键粘贴即可。
1. 打开 MATLAB,观察界面中包含哪些子窗口,写出窗口名字 解答: (不需截图) 命令窗口( Command Window ) 、工作区窗口( Workspace ) 、历史命令窗口 (Command History) 、当前文件夹窗口(Current Folder) 、...
u
(e 2)
+ v)
2
v2 − u
3) 解答:
u = 2; v = 3;
u − 3v uv
%1) 4*u*v/log(v) %2) (exp(u)+v)^2/(v^2-u) %3) sqrt(u-3*v)/(u*v) 运行屏幕截图:
9. 计算如下表达式: 1) ( 3 − 5i )( 4 + 2i ) 2) sin ( 2 − 8i ) 解答:
2. 使用 path 命令将计算机的 C 盘添加到系统路径,使用 cd 命令将 C 盘设置为 当前工作路径 解答:(不需截图,仅需要写出相应的指令表达式即可) path('C:\', path) cd('C:\') %或者 path(path, 'C:\')
(完整版)MATLAB)课后实验答案[1]
1 + e2 (2) z = 1 ln( x + 1 + x 2 ) ,其中 x = ⎡⎢ 2⎣-0.45 ⎦2 2 ⎪t 2 - 2t + 1 2 ≤ t <3 ⎨实验一MATLAB 运算基础1. 先求下列表达式的值,然后显示 MATLAB 工作空间的使用情况并保存全部变量。
(1) z = 2sin 8501221 + 2i ⎤5 ⎥(3) z = e 0.3a - e -0.3asin(a + 0.3) + ln 0.3 + a ,a = -3.0, - 2.9, L , 2.9, 3.03⎧t 2 0 ≤ t < 1 (4) z = ⎪t 2 - 11 ≤ t <2 ,其中 t=0:0.5:2.5 4⎩解:M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(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=⎣O2⨯3⎥,其中E、R、O、S分别为单位矩阵、随机矩阵、零矩S⎦阵和对角阵,试通过数值计算验证A=⎢⎥。
matlab课后习题答案第一章
第 1 章基础准备及入门习题1及解答1 数字1.5e2,1.5e3 中的哪个与1500相同吗?〖解答〗1.5e32 请指出如下5个变量名中,哪些是合法的?abcd-2 xyz_3 3chan a变量 ABCDefgh〖解答〗2、5是合法的。
3 在MATLAB环境中,比1大的最小数是多少?〖解答〗1+eps4 设 a = -8 , 运行以下三条指令,问运行结果相同吗?为什么?w1=a^(2/3)w2=(a^2)^(1/3)w3=(a^(1/3))^2〖解答〗(1)不同。
具体如下w1=a^(2/3) %仅求出主根w2=(a^2)^(1/3) %求出(-8)^2的主根w3=(a^(1/3))^2 %求出(-8)主根后再平方w1 =-2.0000 + 3.4641iw2 =4.0000w3 =-2.0000 + 3.4641i(2)复数的多方根的,下面是求取全部方根的两种方法:(A)根据复数方根定义a=-8;n=2;m=3;ma=abs(a);aa=angle(a);for k=1:m %m决定循环次数sa(k)=(aa+2*pi*(k-1))*n/m; %计算各根的相角endresult=(ma^(2/3)).*exp(j*sa) %计算各根result =-2.0000 + 3.4641i 4.0000 - 0.0000i -2.0000 - 3.4641i(B )利用多项式023=-a r 求根p=[1,0,0,-a^2];r=roots(p)r =-2.0000 + 3.4641i-2.0000 - 3.4641i4.00005 指令clear, clf, clc 各有什么用处?〖解答〗clear清除工作空间中所有的变量。
clf清除当前图形。
clc清除命令窗口中所有显示。
6 以下两种说法对吗?(1)“MATLAB 进行数值的表达精度与其指令窗中的数据显示精度相同。
” (2)MATLAB 指令窗中显示的数值有效位数不超过7位。
matlab2016年课后习题及答案详解【精选】
2012年12月7日星期五第2章 MATLAB 矩阵运算基础2.1 在MATLAB 中如何建立矩阵,并将其赋予变量a ?⎥⎦⎤⎢⎣⎡194375>> a=[5 7 3;4 9 1]2.2 有几种建立矩阵的方法?各有什么优点?可以用四种方法建立矩阵:①直接输入法,如a=[2 5 7 3],优点是输入方法方便简捷;②通过M 文件建立矩阵,该方法适用于建立尺寸较大的矩阵,并且易于修改;③由函数建立,如y=sin(x),可以由MATLAB 的内部函数建立一些特殊矩阵;④通过数据文件建立,该方法可以调用由其他软件产生数据。
2.3 在进行算术运算时,数组运算和矩阵运算各有什么要求?进行数组运算的两个数组必须有相同的尺寸。
进行矩阵运算的两个矩阵必须满足矩阵运算规则,如矩阵a 与b 相乘(a*b )时必须满足a 的列数等于b 的行数。
2.4 数组运算和矩阵运算的运算符有什么区别?在加、减运算时数组运算与矩阵运算的运算符相同,乘、除和乘方运算时,在矩阵运算的运算符前加一个点即为数组运算,如a*b 为矩阵乘,a.*b 为数组乘。
2.5 计算矩阵与之和。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6];>> a+bans =777914131512142.6 求的共轭转置。
⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x >> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i 2.7 计算与的数组乘积。
MATLAB第一章作业答案
M A T L A B第一章作业答案(共5页)--本页仅作为文档封面,使用时请直接删除即可----内页可以根据需求调整合适字体及大小--第一章M A T L A B概况与基本操作1.选择题:(1)最初的MATLAB核心程序是采用A语言编写的。
(2)即将于2011年9月发布的MATLAB新版本的编号为D。
2011Ra 2011Rb R2011a R2011b(3)在默认设置中,MATLAB中的注释语句显示的颜色是D。
A.黑色B.蓝色C.红色D.绿色(4)如果要以科学计数法显示15位有效数字,使用的命令是B。
long long e long g long d(5)在命令窗口新建变量a、b,如果只查看变量a的详细信息,使用的命令为B。
a a(6)如果要清除工作空间的所有变量,使用的命令为 C 。
all C.两者都可 D.两者都不可(7)在创建变量时,如果不想立即在命令窗口中输出结果,可以在命令后加上D。
A.冒号B.逗号C.空格D.分号(8)如果要重新执行以前输入的命令,可以使用B键。
A.下箭头↓B.上箭头↑C.左箭头←D.右箭头→(9)如果要查询函数inv的相关信息,并显示在命令窗口,应使用命令A。
D.三者均可(10)如果要启动Notebook文档,下列D操作是可行的。
A.在命令窗口输入notebook命令B.在命令窗口输入notebook filename命令C.在Word中启动M-book文档D.三者均可2.填空题:(1)MATLAB是MATrix和LABoratory两个单词前三个字母的组合,意为“矩阵实验室”,它的创始人是Cleve Moler和Jack Little。
(2)在MATLAB的默认设置中,关键字显示的字体为蓝色,命令、表达式、计算结果显示的字体为黑色,字符串显示的字体为紫色,注释显示的字体为绿色,错误信息显示的字体为红色。
(3)在命令窗口中,输出结果显示为各行之间添加空行的命令为format compact,各行之间不添加空行的命令为format compact。
matlab实验:运算基础并且附有答案 (1)
实验二、MATLAB运算基础一、实验目的掌握MATLAB各种表达式的书写规则及常用函数的使用。
掌握MATLAB中字符串、元胞数组和结构的常用函数的使用。
二、实验内容及步骤1、设有矩阵A和B,A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 1718 19 20;21 22 23 24 25],B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 1311]1)求它们的乘积C >>C=A*B2)将矩阵C的右下角3x2子矩阵赋给D >>I=[3 4 5];J=[2 3];D=C(I,J)也可以用>>D=C([3 4 5],[2 3])D =520 397705 557890 7172、完成下列操作1)求[100,999]之间能被61整除的数及其个数(提示:先利用冒号表达式,再利用find和length函数。
)>> a=100:999;find(rem(a,61)==0)ans =23 84 145 206 267 328 389 450 511572 633 694 755 816 877>> b=a(ans)b =122 183 244 305 366 427 488 549 610 671 732 793 854 915 976>> length(b)ans =152)建立一个字符串向量,删除其中的大写字母(提示:利用find函数和空矩阵。
)a=’I am maying’;a( find(a>’A’&a<’Z’))=[]3、已知A=[23 10 -78 0;41 -45 65 5;32 5 0 32;6 -54 92 14],取出其前3行构成矩阵B,其前两列构成矩阵C,其左下角3x2子矩阵构成矩阵D,B与C的乘积构成矩阵E,分别求E<D、E&D、E|D、~E|~D。
matlab练习一、二答案
练习一答案1、一个三位整数各位数字的立方和等于该数本身,则称该数为水仙花数。
试编写程序找到,并输出全部水仙花数。
for n=100:999a=floor(n/100);b=floor(rem(n,100)/10);c=rem(n,10);if n==a^3+b^3+c^3s=nendend2、若一个数等于它的各个真因子之和,则称该数为完数,如6=1+2+3,所以6 是完数。
求[1 500]之间的全部完数。
for i=2:500s=1;for j=2:(i-1)if (rem(i,j)==0)s=s+j;endendif s==isendend3、写一个函数pifun 计算下列级数:F(n)=4*(1-1/3+1/5-1/7+...)并使用tic 和toc 指令来测量pifun(1000000)的计算时间。
n=input('请输入n值:');tic;f=0;for i=1:nf=f+4*((-1)^(i+1))/(2*i-1) ;endftoc4、输入一个百分制成绩,要求输出成绩等级为A~E,即90~100 为A,80~89 为B,70~79为C,60~69 为D,60 以下为E。
m=input('请输入一个百分制成绩:');m=floor(m/10);switch mcase {10,9}disp('A');case 8disp('B');case 7disp('C');case 6disp('D');disp('E');otherwisedisp('E');end5、输入x,y 的值,并将它们的值互换后输出。
x = input('x=');y = input('y=');disp ([x,y]);t=x;x=y;y=t;disp ([x,y]);6、某商场对顾客所购买的商品实行打折销售,标准如下(商品价格用price 来表示):price<200 没有折扣;200≤price<500 3%折扣;500≤price<1000 5%折扣;1000≤price<2500 8%折扣;2500≤price<5000 10%折扣;5000≤price 14%折扣输入所售商品的价格,求其实际销售价格。
matlab课后习题答案1到6章
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”?答:因为向量可以看成是仅有一行或一列的矩阵,单个数据(标量)可以看成是仅含一个元素的矩阵,故向量和单个数据都可以作为矩阵的特例来处理。
因此,矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵,问:(1)A*B和A.*B的值是否相等?答:不相等。
(2) A./B和B.\A的值是否相等?答:相等。
(3)A/B和B\A的值是否相等?答:不相等。
(4)A/B和B\A所代表的数学含义是什么?答:A/B等效于B的逆右乘A矩阵,即A*inv(B),而B\A等效于B矩阵的逆左乘A矩阵,即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答:B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答:A(7)=[](3)将矩阵A的每个元素值加30。
答:A=A+30;(4)求矩阵A的大小和维数。
答:size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答:t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34⨯矩阵。
答:reshape(x,3,4);(7)求一个字符串的ASCII码。
答:abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答:char(49);4.下列命令执行后,L1、L2、L3、L4的值分别是多少?A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答:L1的值为[0, 0, 0, 0, 1, 0,0, 0, 0]L2的值为[1, 1, 1, 1,1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1,1, 1, 0, 0, 0]L4的值为[4, 5, 6]5.已知23100.7780414565532503269.5454 3.14A-⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦完成下列操作:(1)取出A的前3行构成矩阵B,前两列构成矩阵C,右下角32⨯子矩阵构成矩阵D,B与C的乘积构成矩阵E。
MatLab练习册1(答案)(1)
MatLab 练习题(1) 参考答案编题: 周德强 (长江大学数学学院)1, 基本运算符输入练习>> (5*2+1.3-0.8)*10^2/25 答案:422,计算y =23log )10sin(4)43.0(2+⋅+-e π 答案:2.2618>>sin(10*pi)*exp(-0.3+4^2)+log(23)/log(4) %将23log 4化成自然对数的形式 3,换行输入练习: >>z = 10*sin(pi/3)* ...cos(pi/3) 答案:4.33014,函数求值运算 >>x= sin(223/3); >>y = x^2; >>z = y*10;>>x+2y-5z 答案:-37.60255,ie z i z i z ⋅==+=π23215,2,72,计算3221z z z z z +=答案:-2.1379 + 1.6552i>>z1=2+7*i >>z2=2*i>>z3=5*exp(2*pi*i) >>z=(z1*z2)/(z2+z3)6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x 答案:(填写程序语句) >>x=(3:5.5:44) 或 >>x=[3:5.5:44] 或 >>x=3:5.5:44 7, 建立等差一维数组x :首项为0,末项为π,项数为15 答案:(填写程序语句) >>x=linspace(0,pi,15)8, 计算行列式的值187624323=A 答案:-8>>A = [3,2,3; 4,2,6; 7,8,1] %注意n 行n 列的矩阵输入方式>>det(A)9, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=187624323A ,矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积.答案:>>A = [3,2,3; 4,2,6; 7,8,1] >>B = [1,1,1; 2,2,2; 3,3,3] >> A*B16 16 16 26 26 2626 26 26>> A.*B %体会点运算 3 2 3 8 4 12 21 24 310,写出计算 Sin(30o )的程序语句.答案:(填写程序语句) >> sin(pi*30/180) 或 sin(pi/6)11,,计算f(2), f(1)f(2)+f 2(3) 答案:>> x=2>> f2=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14) f2 =-4.4792 >> x=1>> f1=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14) >> x=3>> f3=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14) >> f1*f2+f3^2-12.602312,,计算f 在x=[2 3 4 6]各处的函数值; 答案:-4.47920.9484 2.2741 4.3335方法1: >> x=[2 3 4 6];>> f=(x.^3 - 2*x.^2 + x - 6.3)./(x.^2 + 0.05*x - 3.14) %注意用点运算方法2: >> syms x %定义x为符号变量>> f=(x^3 - 2*x^2 + x - 6.3)/(x^2 + 0.05*x - 3.14);%写出f的符号表达式>> f1=subs(f,x,[2 3 4 6]) %计算f在x=[2 3 4 6]各处的函数值MATLAB提供基本的算术运算有:加 (+)、减 (-)、乘 (*)、除 (/)、幂次方 (^),范例为:5+3, 5-3, 5*3, 5/3, 5^3MATLAB常用数学函数●三角函数和双曲函数●指数函数●复数函数●圆整函数和求余函数● 矩阵变换函数● 其他函数例:>>y = sin(10)*exp(-0.3*4^2) y = -0.0045【例】复数ie z i z i z 63212,21,43π=+=+=表达,及计算321z z z z =。
(完整版)MATLAB)课后实验答案[1]
实验一 MATLAB 运算基础1. 先求下列表达式的值,然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =,其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0:0.5:2.5 解:4. 完成下列操作:(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量,删除其中的大写字母。
解:(1) 结果:(2). 建立一个字符串向量例如:ch='ABC123d4e56Fg9';则要求结果是:实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦,其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵,试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下;5. 下面是一个线性方程组:1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解,并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解: M 文件如下: 123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现,分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
Matlab程序设计(2016大作业)
Matlab程序设计课程大作业题目名称:_________________________________ 班级:_________________________________ 姓名:_________________________________ 学号:_________________________________ 课程教师:温海骏学期:2015-2016学年第2学期完成时间:MATLAB 优化应用§1 线性规划模型一、线性规划问题:问题1:生产计划问题假设某厂计划生产甲、乙两种产品,现库存主要材料有A 类3600公斤,B 类2000公斤,C 类3000公斤。
每件甲产品需用材料A 类9公斤,B 类4公斤,C 类3公斤。
每件乙产品,需用材料A 类4公斤,B 类5公斤,C 类10公斤。
甲单位产品的利润70元,乙单位产品的利润120元。
问如何安排生产,才能使该厂所获的利润最大。
问题2:投资问题某公司有一批资金用于4个工程项目的投资,其投资各项目时所得的净收益(投入资金百分比)如下表:工程项目收益表由于某种原因,决定用于项目A 的投资不大于其他各项投资之和而用于项目B 和C 的投资要大于项目D 的投资。
试确定该公司收益最大的投资分配方案。
问题3:运输问题有A 、B 、C 三个食品加工厂,负责供给甲、乙、丙、丁四个市场。
三个厂每天生产食品箱数上限如下表:四个市场每天的需求量如下表:从各厂运到各市场的运输费(元/每箱)由下表给出:求在基本满足供需平衡的约束条件下使总运输费用最小。
§2 多目标规划模型多目标规划定义为在一组约束下,多个不同的目标函数进行优化设计。
数学模型:12min ()()().()0,1,2,,m j f x f x f x st g x j k⎡⎤⎣⎦≤=其中x=(x 1 ,x 2 , … ,x n )为一个n 维向量;f i (x)为目标函数,i=1,2, … ,m; g j (x)为系统约束, j=1, 2, … ,k 。
matlab课后练习习题及答案详解
matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较,MATLAB语言突出的特色是什么?MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。
MATLAB系统由那些部分构成?MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。
安装MATLAB时,在选择组件窗口中哪些部分一定勾选,没有勾选的部分此后怎样补安装?在安装MATLAB时,安装内容由选择组件窗口中个复选框能否被勾选来决定,能够依据自己的需要选择安装内容,但基本平台(即MATLAB选项)一定安装。
第一次安装没有选择的内容在补安装时只要依据安装的过程进行,不过在选择组件时只勾选要补装的组件或工具箱即可。
MATLAB操作桌面有几个窗口?怎样使某个窗口离开桌面成为独立窗口?又怎样将离开出去的窗口从头搁置到桌面上?在MATLAB操作桌面上有五个窗口,在每个窗口的右上角有两个小按钮,一个是封闭窗口的Close按钮,一个是能够使窗口成为独立窗口的Undock 按钮,点击Undock按钮就能够使该窗口离开桌面成为独立窗口,在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。
怎样启动M文件编写/调试器?在操作桌面上选择“成立新文件”或“翻开文件”操作时,M文件编写/调试器将被启动。
在命令窗口中键入edit命令时也能够启动M文件编写/调试器。
储存在工作空间中的数组能编写吗?怎样操作?储存在工作空间的数组能够经过数组编写器进行编写:在工作空间阅读器中双击要编写的数组名翻开数组编写器,再选中要改正的数据单元,输入改正内容即可。
命令历史窗口除了能够察看前方键入的命令外,还有什么用途?命令历史窗口除了用于查问从前键入的命令外,还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
怎样设置目前目录和搜寻路径,在目前目录上的文件和在搜寻路径上的文件有什么差别?目前目录能够在目前目录阅读器窗口左上方的输入栏中设置,搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MATLAB 作业1参考答案
1 、试用符号元素工具箱支持的方式表达多项式5432()34236f x x x x x x =+++++,并令11
s x s -=+,将f(x) 替换成s 的函数。
【求解】可以先定义出f 函数,则由subs() 函数将x 替换成s 的函数
>> syms s x;f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;F=subs(f,x,(s-1)/(s+1)) F =
(s-1)^5/(s+1)^5+3*(s-1)^4/(s+1)^4+4*(s-1)^3/(s+1)^3+
2*(s-1)^2/(s+1)^2+3*(s-1)/(s+1)+6
2、用数值方法可以求出63626302
124822i i S ===++++++∑,试不采用循环的形式求出
和式的数值解。
由于数值方法采用double 形式进行计算的,难以保证有效位数字,所以结果不一定精确。
试采用符号运算的方法求该和式的精确值。
【求解】用符号运算的方式可以采用下面语句
>> sum(2.^[1:63])
ans =
1.844674407370955e+019
>> sum(sym(2).^[1:63])
ans =
18446744073709551614
由于结果有19 位数值,所以用double 型不能精确表示结果,该数据类型最多表示16 位有效数字。
其实用符号运算方式可以任意保留有效数字,例如可以求200 项的和或1000 项的和可以由下面语句立即得出。
>> sum(sym(2).^[1:200])
ans = 3213876088517980551083924184682325205044405987565585670602750
>> sum(sym(2).^[1:1000])
ans = 2143017214372534641896850098120003621122809623411067214887500776740 70 2102249872244986396757631391716255189345835106293650374290571384628 08 7196915514939714960786913554964846197084214921012474228375590836430 60 9294996716388253479753511833108789215412582914239295537308433532085 96
63305248773674411336138750
3、编写一个矩阵相加函数mat_add() ,使其具体的调用格式为A=mat_add(A1 ,A2 ,A3 ,…) ,要求该函数能接受任意多个矩阵进行加法运算。
(注:varargin 变量的应用)
【求解】可以编写下面的函数,用varargin 变量来表示可变输入变量
function A=mat_add(varargin)
A=0;
for i=1:length(varargin), A=A+varargin{i}; end
如果想得到合适的错误显示,则可以试用try, catch 结构。
function A=mat_add(varargin)
try
A=0;
for i=1:length(varargin), A=A+varargin{i}; end
catch, error(lasterr); end
4、下面给出了一个迭代模型
211
1 1.40.3k k k k k x y x y x ++⎧=+-⎨=⎩ 写出求解该模型的M-函数(M-脚本文件),如果取迭代初值为000,0x y == ,那么请进行30000 次迭代求出一组x 和y 向量,然后在所有的k x 和k y 坐标处点亮一个点( 注意不要连线) ,最后绘制出所需的图形。
(提示这样绘制出的图形又称为Henon 引力线图,它将迭代出来的随机点吸引到一起,最后得出貌似连贯的引力线图。
)
【求解】用循环形式解决此问题,可以得出所示的Henon 引力线图。
>> x=0; y=0;
for i=1:29999
x(i+1)=1+y(i)-1.4*x(i)^2;
y(i+1)=0.3*x(i);
end
plot(x,y,'.')
上述的算法由于动态定义x 和y ,所以每循环一步需要重新定维,这样做是很消耗时间的,所以为加快速度,可以考虑预先定义这两个变量,如给出x=zeros(1,30000)。
5、对合适的θ范围选取分别绘制出下列极坐标图形(注:要求把图形窗口分为4块,每块绘一个图)
①21.0013ρθ=,②cos(7/2)ρθ=,③sin()/ρθθ=,④31cos (7)ρθ=-
【求解】绘制极坐标曲线的方法很简单,用polar( ) 即可以绘制出极坐标图。
注意绘制图形时的点运算:
>> t=0:0.01:2*pi; subplot(221), polar(t,1.0013*t.^2),% (a)。