MATLAB程序设计与应用 刘卫国第二章作业
第二章作业
3.写出完成下列操作的命令。
（1）将矩阵A第2~5行中第1，3，5列元素赋给矩阵B。MA TLAB程序：
（2）删除矩阵A的第七号元素。
MA TLAB程序：
(3).将矩阵A每个元素值加30
MA TLAB程序;
MA TLAB程序：
MA TLAB程序：
MA TLAB程序;
MA TLAB程序;
MA TLAB程序;
MA TLAB程序：
5.已知矩阵，完成下列操作。
（1）。取出A的前3列构成矩阵B，前两列构成矩阵C，右下角3乘2子矩阵构成矩阵D，B与C的乘积构成矩阵E
MA TLAB程序：
MA TLAB程序;
matlab课后习题解答第二章
第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
MATLAB第一章作业答案
第一章 M A T L A B 概况与基本操作 1.选择题（每题2分，共20分）： (1)最初的MATLAB 核心程序是采用D 语言编写的。 A.PASCAL B.C C.BASIC D.FORTRAN (2)即将于2011年9月发布的MATLAB 新版本的编号为C 。 A.MATLAB 2011Ra B.MATLAB 2011Rb C.MATLAB R2011b D.MATLAB R2011a (3)在默认设置中，MATLAB 中的注释语句显示的颜色是B 。 A.黑色 B.绿色 C.红色 D.蓝色 (4)如果要以科学计数法显示15位有效数字，使用的命令是A 。 A.format long e B.format long C.format long g D.format long d (5)在命令窗口新建变量a 、b ，如果只查看变量a 的详细信息，使用的命令为A 。 A.whos a B.who a C.who D.whos (6)如果要清除工作空间的所有变量，使用的命令为C 。 A.clear B.clear all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上B 。 A.冒号 B.分号 C.空格 D.逗号 (8)如果要重新执行以前输入的命令，可以使用D 键。 A.下箭头↓ B.右箭头→ C.左箭头← D.上箭头↑ (9)如果要查询函数det 的功能和用法，并显示在命令窗口，应使用命令C 。 A.doc B.lookfor C.help D.三者均可 (10)如果要启动Notebook 文档，下列D 操作是可行的。 A.在命令窗口输入notebook 命令 B.在命令窗口输入notebook filename 命令 C.在Word 中启动M-book 文档 D.三者均可 2.填空题（每空1分，共20分）： (1)MATLAB 是matrix 和laboratory 两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler 和Jack Little 。 (2)在MATLAB 的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为褐红色，注释显示的字体为绿色，错误信息显示的字体为红色。 (3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format loose ，各行之间不添加空行的命令为format compact 。 (4)在MATLAB 中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组每行各个元素，逗号的作用是分隔数组每行各个元素或函数的各个输入参数，分号的作用是作为不显示命令结果的命令行的结尾或分隔数组各列，冒号的作用是生成一维数组或表示数组全部元素，百分号的作用是引导一行注释，…的作用是连接相邻两行，感叹号的作用是调用操作系统命令。 3.程序设计题（每题10分，共40分） (1)以25m/s 的初速度向正上方投球（g=9.8m/s 2），计算到达最高点的时间tp 以及球从出发点到 最高点的距离hp 。 解：根据物理学知识，物体上抛运动的速度与经过的时间之间的关系为0p p v v gt =-，因此所需要的时间为0p p v v t g -=。而到达最高点时的速度0p v =，因此可根据此公式求出tp ： v0=25;g=9.8;vp=0; tp=(v0-vp)/g tp = 2.5510
MATLAB练习作业(已做)2012
说明:如在测试时请将程序中的中文字符改为西文(部分程序为手写,并未经matlab验证) 2012年春《MATLAB基础》第一次上机练习与作业练习１：脚本文件 建立一个文件名为helloWorld的脚本文件，运行该文件时显示出如下文字： Hello World! I am going to learn MATLAB ?提示: 用disp显示字符串，将要显示的字符号串用单引号引起来。如'This is a string' 练习２：变量 获取并保存当前的日期与时间 ●用函数clock生成一个变量，变量名为start ●用size查看start的维数，它是一个行向量还是一个列向量？ ●start包含什么内容？用help clock查看 ●用函数datestr将向量start转换成字符串，得到新的变量,名为startString ●将start与startString保存为mat文件，文件名为startTime ●在练习１建立的脚本文件helloWorld.m文件中，用load函数导入变量startTime, 并显示如下文字： I started learning MATLAB on * start date and time* 练习3：标量 你将要以指数增长的速度来学习MATLAB, 将如下内容添加到helloWorld.m文件中 ●假设你的学习时间是一个常量,为1.5 days，将此时间用秒表示，赋给变量tau ●假设课程持续时间为5 days. 将这个时间单位转换为秒，保存在变量endofClass中？ ●将学到的知识描述为t的函数，函数方程为： ●用函数datestr将向量start转换成字符串，得到新的变量startString ●在课程结束时间endofClass，你将学到多少知识？用变量knowledgeAtEnd表示(指数函数exp)
第二次作业《解释结构模型应用》
大连海事大学 实验报告 《系统工程》 2014～2015学年第一学期 实验名称：基于解释模型在大学生睡眠质量问题的研究学号姓名：马洁茹姚有琳 指导教师：贾红雨 报告时间： 2014年9月24日
《系统工程》课程上机实验要求 实验一解释结构模型在大学生睡眠质量问题中的研究 实验名称：基于MATLAB软件或C/Java/其他语言ISM算法程序设计(一) 实验目的 系统工程课程介绍了系统结构建模与分析方法——解释结构模型法（Inter pretative Structural Modeling ·ISM）是现代系统工程中广泛应用的一种分析方法，能够利用系统要素之间已知的零乱关系，用于分析复杂系统要素间关联结构，揭示出系统内部结构。ISM方法具有在矩阵的基础上再进一步运算、推导来解释系统结构的特点，对于高维多阶矩阵的运算依靠手工运算速度慢、易错，甚至几乎不可能。 本次实验的目的是应用计算机应用软件或者是基于某种语言的程序设计快速实现解释结构模型(ISM)方法的算法，使学生对系统工程解决社会经济等复杂性、系统性问题需要计算机的支持获得深刻的理解。学会运用ISM分析实际问题。 (二) 实验要求与内容： 1.问题的选择 根据对解释结构模型ISM知识的掌握，以及参考所给的教学案例论文，决定选择与我们生活有关的——大学生睡眠质量问题。 2．问题背景 睡眠与我们的生活息息相关，当每天的身体机制在不断运行的过程中身体负荷不断变大，到了夜间就需要休息。但是同一寝室的同学大多休息时段不同，有些习惯早睡，有些会由于许多原因晚睡。有些睡眠较沉不会轻易被打扰，有些睡眠较轻容易被鼾声或者其他声响惊醒。学习得知，解释系统模型是通过对表面分离、凌乱关系的研究，揭示系统内部结构的方法。因此，我想尝试通过解释模型来对该问题进行研究分析。 3.用画框图的形式画出ISM的建模步骤。
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
matlab作业第10章
10-2建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠标左键之后显示出Left Button Pressed 字样。 >> hf=figure('color','r','windowbuttondownfcn','axis off;text(0.1,0.5,''Left Button Pressed'',''fontsize'',30)') hf = 1 10-3利用图形对象绘制下列曲线，要求先利用缺省属性绘制曲线，然后通过图形句柄操作来改变曲线的颜色、线型和线宽，并利用文字对象给曲线添加文字标注。 （1）4 2 11x x y ++= x=0:0.1:10; y=((1+x.^2)/(1+x.^4)); hf=plot(x,y); >> set(hf,'color','r','linestyle',':','marker','*','linewidth',1);
（4）?????==325t y t x >> t=0:2:100; >> x=t.*t; >> y=5*t.^3; >> hf=plot(x,y); >>set(hf,'color','b','linestyle',':','marker','p','linewidth',0.3);
10-4利用图形对象绘制下列三维图形，要求与上题相同。 （1）?? ???===t z t y t x sin cos >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t; >> hf=plot3(x,y,z);
第二次数学建模作业
4. 根据表1.14 的数据，完成下列数据拟合问题： 表 1.14 美国人口统计数据(百万人) 年份1790 1800 1810 1820 1830 1840 1850 1860 人口 3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 年份1870 1880 1890 1900 1910 1920 1930 1940 人口38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 年份1950 1960 1970 1980 1990 2000 人口150.7 179.3 204.0 226.5 251.4 281.4 解答：（1）： （i）执行程序： t=1790:10:2000; x=[3.9,5.3,7.2,9.6,12.9,17.1,23.2,31.4,38.6,50.2,62.9,76.2,92.0,106.5,123.2,131.7,150.7,179.3,204 .0,226.5,251.4,281.4]; f=@(r,t)3.9.*exp(r(1).*(t-1790)); r=nlinfit(t,x,f,0.036) sse=sum((x-f(r,t)).^2) plot(t,x,'k+',1790:10:2000,f(r,1790:10:2000),'k') axis([1790,2000,0,300]),legend('测量值','理论值') xlabel('美国人口/（百万）'),ylabel('年份') title('美国人口指数增长模型图II') 运行结果： >> Untitled r = 0.0212 sse = 1.7433e+004 即，拟合效果：r =0.0212；误差平方和为：1.7433e+004. 拟合效果图（i）：
MATLAB第一章作业答案
第一章M A T L A B概况与基本操作 1.选择题： (1)最初的MATLAB核心程序是采用A语言编写的。 A.FORTRAN B.C C.BASIC D.PASCAL (2)即将于2011年9月发布的MATLAB新版本的编号为D。 A.MATLAB 2011Ra B.MATLAB 2011Rb C.MATLAB R2011a D.MATLAB R2011b (3)在默认设置中，MATLAB中的注释语句显示的颜色是D。 A.黑色 B.蓝色 C.红色 D.绿色 (4)如果要以科学计数法显示15位有效数字，使用的命令是B。 A.format long B.format long e C.format long g D.format long d (5)在命令窗口新建变量a、b，如果只查看变量a的详细信息，使用的命令为B。 A.who a B.whos a C.who D.whos (6)如果要清除工作空间的所有变量，使用的命令为 C 。 A.clear B.clear all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上D。 A.冒号 B.逗号 C.空格 D.分号 (8)如果要重新执行以前输入的命令，可以使用B键。 A.下箭头↓ B.上箭头↑ C.左箭头← D.右箭头→ (9)如果要查询函数inv的相关信息，并显示在命令窗口，应使用命令A。 A.help B.lookfor C.doc 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。 备注：本题布置给大家时有一点小错误，现在予以更正。 (4)在MATLAB中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组中每一行的各个元素，逗号的作用是分隔数组中每一行的各个元素或不同的命令，分号的作用是分隔数组中的各行或控制命令执行结果是否在命令窗口显示，冒号的作用是生成一维数组或显示全部元素，百分号的作用是注释行的开头，…的作用是把相邻两行的语句连接为一个命令，感叹号的作用是执行操作系统命令。 3.先建立自己的工作目录，再将自己的工作目录设置到MATLAB搜索路径下。请写出操作步骤或用Matlab命令实现。用help命令能查询到自己的工作目录吗？ 解：操作步骤： (1)在Windows环境中建立一个工作目录，如：c:\mywork； (2)启动MATLAB，在File菜单中选择Set Path…命令，显示出如下图的对话框：
西安交通大学数字图像处理第二次作业
数字图像处理的基本数学 工具的使用 摘要 本报告主要介绍了运用编程软件MATLAB对图像灰度级进行变换、求取图像均值与方差、采用不同的内插方法对图像进行缩放及利用仿射变换对图像进行空间变换处理的方法。同时，对最近邻内插法、双线性内插法、双三次内插法进行图像处理的效果进行了详细的对比，并对出现差异的原因做出了简要分析。 姓名： X X X 班级： 学号： 提交日期：年月日
2_1. 把lena 512*512图像灰度级逐级递减8-1显示； (1) 问题分析： 所要实现的功能是：在不改变图像大小的前提下，使得整幅图像的灰度级逐级递减并将图像显示出来；即对所有像素点的灰度依次进行除2操作； (2) 实验过程： 工具：MATLAB 软件； 利用imread()函数将图像读入MATLAB ，利用imshow()对原图像进行显示，再利用循环体对整幅图像的灰度级逐级递减并一一进行显示。 源代码附于本报告最后一部分。 (3) 实验结果： 备注：在大小为512X512的途中观察更为方便，但此处为了便于排版以及将结果进行对比对所有图像做了一定的缩小。 a b c d
e f g h 图2_1 (a)大小为512X512的256灰度级图像；(b)~(h)保持图像大小不变的同 时以灰度级128,62,32,16,8,4,2显示的图像。 (4) 结果分析： 对图2_1中的(a)~(h)图像进行对比可知，256级、128级、64级以及32级灰度的图像几乎没有太大的区别；然而在灰度级为16的图(e)中出现了较为明显的伪轮廓，这种效果是由数字图像的平滑区域中的灰度级数不足引起的。（说明：此分析为本人肉眼的观察结果，对细节的观察难免存在疏漏之处，还请批评指正。） 2_2. 计算lena 图像的均值方差； (1) 问题分析： 所要实现的功能是：计算图像‘lena.bmp ’的均值与方差； (2) 实验过程： 工具：MATLAB 软件； 利用imread()函数将图像读入MATLAB ，由于二维数字图像使用二维阵列表示的，因而可以直接利用MATLAB 中的mean2()及std2()分别求整幅图像的均值于方差； 源代码附于本报告最后一部分。 (3) 实验结果：均值 m =99.0512 方差 =52.8776。 2_3. 把lena 图像用近邻、双线性和双三次插值法zoom 到2048*2048；
MATLAB上机答案
一熟悉Matlab工作环境 1、熟悉Matlab的5个基本窗口 思考题： （1）变量如何声明，变量名须遵守什么规则、是否区分大小写。 答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则： 变量名必须以字母开头，只能由字母、数字或下划线组成。 变量名区分大小写。 变量名不能超过63个字符。 关键字不能作为变量名。 最好不要用特殊常量作为变量名。 （2）试说明分号、逗号、冒号的用法。 分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。 逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。 冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。 （3）linspace（）称为“线性等分”函数，说明它的用法。 LINSPACE Linearly spaced vector.线性等分函数 LINSPACE(X1,X2)generates a row vector of100linearly equally spaced points between X1and X2. 以X1为首元素，X2为末元素平均生成100个元素的行向量。 LINSPACE(X1,X2,N)generates N points between X1and X2. For N<2,LINSPACE returns X2. 以X1为首元素，X2为末元素平均生成n个元素的行向量。如果n＜2，返回X2。 Class support for inputs X1,X2: float:double,single 数据类型：单精度、双精度浮点型。 （4）说明函数ones（）、zeros（）、eye（）的用法。 ones（）生成全1矩阵。 zeros（）生成全0矩阵。 eye（）生成单位矩阵。 2、Matlab的数值显示格式
北科,北京科技大学,数学实验,MATLAB第二次作业
《数学实验》报告 实验名称 MATLAB绘图 学院 专业班级 姓名 学号 年月
一、 【实验目的】 了解并学习绘制MATLAB 二维曲线和三维曲线的图形。 二、 【实验任务】 1.绘制)π4,0(),3sin(3 ∈=x x e y x 的图像，要求用蓝色的星号画图，并且画出其包络线3 x e y ±=的图像，用红色的点划线画图。 3.在同一图形窗口画三个子图，要求使用指令gtext ，axis ，legend ，title ，xlabel ，ylabel ： （1）)ππ,(,cos -∈=x x x y （2）)π,4π(,sin 1 tan 3∈=x x x x y （3）]8,1[,sin 1∈=x x e y x 5.绘制圆锥螺线的图像并加各种备注，圆锥螺线的参数方程为： π]20,0[,t 2z 6π sin 6π cos { ∈===t t t y t t x 三、 【实验程序】 1. x=0:pi/50:4*pi; y1=exp(x/3).*sin(3*x); y2=exp(x/3); y3=-exp(x/3); plot(x,y1,'b*') hold on plot(x,y2,'r-.') hold on plot(x,y3,'r-.') 3. x1=-pi:pi/50:pi; y1=x1.*cos(x1); x2=pi:pi/50:4*pi; y2=x2.*tan(1./x2).*sin(x2.^3); x3=1:0.01:8; y3=exp(1./x3).*sin(x3); subplot(221),plot(x1,y1,'r-'),grid on axis tight xlabel('x 轴'),ylabel('y 轴') title('y=xcosx')
计算方法及其MATLAB实现第二章作业
作者：夏云木子 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)-…
matlab上机作业
第二次 上机作业 1、 求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹。 （1）????? ???????--=901511250324153211A （2）??????-=2149.824343.0B 1. A=[1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9] D=diag(A) C=triu(A) B=tril(A) E=rank(A) F=trace(A) a1=norm(A,1) a2=norm(A,inf) a3=norm(A,inf) c1=cond(A) c1=cond(A,1) c2=cond(A,2) c3=cond(A,inf) 2. B=[0.43,43,2;-8.9,4,21] D=diag(B) C=triu(B) B=tril(B) E=rank(B) F=trace(B) a1=norm(B,1) a2=norm(B,inf) a3=norm(B,inf) c1=cond(B) c1=cond(B,1) c2=cond(B,2) c3=cond(B,inf) 2、 求矩阵A 的特征值和相应的特征向量。
???? ??????=225.05.025.0115.011A A=[1,1,0.5;1,1,0.25;0.5,0.25,2] [V ,D]=eig(A) 3、 下面是一个线性方程组： ???? ??????=????????????????????52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解。 （2） 将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相 对变化。 （3） 计算系数矩阵A 的条件数并分析结论。 A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6] B=[0.95,0.67,0.52] X=inv(A)*b c1=cond(A,1) c2=cond(A,2) c3=cond(A,inf) 4、 利用Matlab 提供的randn 函数生成符合正态分布的10×5随机矩阵A,进行 如下操作： (1)A 各列元素的均值和标准方差 (2)A 的最大元素和最小元素 (3)求A 每行元素的和以及全部元素之和 (4)分别对A 的每列元素按升序、每行按降序排列 X=randn(10,5) M=mean(X) D=std(X) m=max(X) n=min(X) P=sum(X,2) sum(p)
matlab第2章作业
深圳大学实验报告课程名称：可视化计算语言 实验项目名称：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 《信号与系统》课程研究性学习手册 姓名吕航(第十七组) 学号09213044 同组成员马琳 刘煜晨 朱志杰 王琛珏 指导教师魏杰 时间2010/11 专题研讨二、信号与系统的频域分析 研讨题目： M3-2 题目分析： (1)根据电路知识解的方程为：y’’+2y’+6y=6x （1）用impulse函数求解单位冲激响应（2）用step函数求单位阶跃响应 仿真程序： untiled21.m ts=0;te=5;dt=0.01; sys=tf([6],[1 2 6]); t=ts:dt:te; y=impulse(sys,t); plot(t,y); xlabel('Time(sec)') ylabel('h(t)') untiled21.m ts=0;te=5;dt=0.01; sys=tf([6],[1 2 6]); t=ts:dt:te; y=step(sys,t); plot(t,y); xlabel('Time(sec)') ylabel('h(t)') 仿真结果： 结果分析： 解得的方程一阶导数前的系数很小所以这个系统的的震荡并不明显，由MA TLAB得画图就可以明显的看出来。 自主学习内容： 陈后金信号与系统P117 阅读文献： [1]杜晶晶金学波信号与系统实训指导[M]西安电子科技大学出版社2009 [2] 陈后金信号与系统[M] 高教出版社2007 发现问题(专题研讨或相关知识点学习中发现的问题)： 在画图是对时间的其实和终点的选择很重要，特别是画二阶的动态方程。要不断地调他的终点时间才可以把整个震荡图形画进去。 问题探究： 只要我的方程没列错结果就应该是正确的。 研讨题目： M3-3 题目分析： 第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 第十章 二项式期权定价（满分70分） 一、单期二项式期权定价 1.复制、定价 例10.1设某股票当前价格为100元，一年后可能上涨50%（即期末价格为150元），也可能下跌25%（即期末价格为75元），无风险债券的利率为5%（国库券年利率），债券的当前价格为1元，现在有1份执行价格为125元、有效期为1年的该股票看涨期权，试计算该看涨期权的价格。 Step1.复制，即1*Call=X*Shares+Y*Bonds ???=+=+0 05.1752505.1150Y X Y X 解出X 和Y 的值。 Step2. 定价，即期末现金流相等，则期初的价格相等: C=X*S0+Y*A0. 输入命令： A=[150,1.05;75,1.05]; B=[25;0]; X=A\B S 0=100;B 0=1; C=X'*[S 0;B 0] 2、风险中性定价 (1) 以离散计息的方式进行折现,则单期二项式期权定价公式应为： ) )(1()0,)1(max()()0,)1(max()(d u r X d S r u X u S d r C ++---+-++= 由于S 、u 、d 、r 、X 的值事先已知，故可以确定出C 的值。 (2) 以连续复利计息的方式进行折现，则单期二项式期权定价公式应为： ) ()0,)1(max()1()0,)1(max()1(d u e X d S e u X u S d e C r r r +---++-+-+= (3) 若考虑股票派发的红利。当已知股票派发红利时，如果再已知股票的波动率，可以直接用Cox-Ross-Rubinstein二项式期权定价模型来为期权进行定价，其函数命令为： [AssetPrice, OptionValue] = binprice(Price, Strike, Rate, Time, Increment,… Volatility, Flag, DividendRate, Dividend, ExDiv) 其中，Price、Strike、Rate、Time、Increment、Volatility、Flag、DividendRate、Dividend、ExDiv分别表示股票的当前价格、期权的执行价格、无风险利率、期权的有效期、时间间隔的大小、股票的波动率、期权的看涨看跌类别（看涨取值1，看跌取值0）、红利率、红利、红利支付日期（最后三个变量是可以任选的，默认值均为0）。 例10.2仍用例10.1中的数据。 输入命令： r=0.05;u=0.50;d=0.25;S=100;X=125; C1=max(S*(1+u)-X,0);C2=max(S*(1-d)-X,0); p=(r+d)/(u+d); C=(p*C1+(1-p)*C2)/(1+r) 例10.3设一看涨期权标的资产当前时刻的价格为52，期权的执行价格为50，无风险利率为10%，期权的有效期5个月，波动性为40%，在3-1/2 月时支付2.06的红利。试计算该看涨期权的价格。 输入命令： >>[Price, Option] = binprice(52, 50, 0.1, 5/12, 1/12, 0.4, 1, 0,2.06, 3.5) 则输出结果给出了标的资产的价格和二叉树的每个节点上的期权价格。 二、多期二项式期权定价 1、两期 例10.5仍用【例10.1】中的数据。假设其它条件相同，当第一期的股价上涨至150元或下跌至75元时，第二期的股价分别针对150元而上涨至225元或下跌至112.5元，而针对75元第二期股价又分别上涨至112.5元或下跌至56.25元。 第一次上机作业： 计算 clc;clear y=((2*sin(0、3*pi)*exp(3))/(1+sqrt(5)))*log10(6) 1.生成等差矩阵a=[1,3,5,7…47,49];在经度[122,128]间生成70个等差数列。 a=(1:2:50) %生成等差矩阵a=[1,3,5,7…47,49] b=linspace(122,128,70) %在经度[122,128]间生成70个等差数列生成34均匀随机矩阵b；提取b得第3行第2列元素赋值为5，求b 中最小值、平均值；提取b得第三列；b右边增添1列为 2 3 4,；删除b 得第二列。 b=rand(3,4) %生成3×4均匀随机矩阵b b(3,2)=5 %提取b得第3行第2列元素赋值为5 m=min(b(:)) %求b中得最小值 n=mean(b(:)) %求b中得平均值 a=b(:,3) %提取b得第三列 b(:,5)=[2 3 4] %在b右边增添1列为2 3 4 b(:,2)=[] %删除b得第二列 在区间x=[0,2]内做sin(x)图；按x、sin(x)、cos(x)存成3列到data、txt 文档；再读取文件data、txt，用读取得数据做cos(x)图。 x=0:0、1*pi:2*pi m=sin(x) n=cos(x) plot(x,m,'r--*') %做sin(x)图像 a=[x;sin(x);cos(x)]' save -ascii data、txt a %按x、sin(x)、cos(x)存成3列到data、txt文档 q=load('data、txt') %读取文件data、txt plot(q(:,1),q(:,3)) %用读取得data、txt做cos(x)图 第二次上机作业： 生成34，[-1,1]间得随机矩阵A；求A矩阵中得元素总个数，将A中小于0得元素等于1；再取出A中大于0、2并小于0、7得元素。 clc;clear A= -1+(1+1)*rand(3,4) %生成3×4，[-1,1]间得随机矩阵A length(A(:)) %求A矩阵中得元素总个数 A(A<0)=1 %将A中小于0得元素等于1 A(A>0、2&A<0、7) %取出A中大于0、2并小于0、7得元素 2.X在[0:0、1: 2]内，求值；并作y>0就是x-y折线图。 clc;clear x=0:0、1:2*pi %X在[0:0、1: 2π]内 y=cos(x、^2) %求y=cos(x^2 )值 a=find(y>0) %y>0 m=x(a); n=y(a); plot(m,n,'r--*') %作x-y折线图 首先由计算机随机产生一个[1,100]之间得一个整数，然后由用户猜测所产生得这个数。根据用户猜测得情况给出不同得提示，如果猜测得数大于产生得数，则显示“High”，小于则显示“Low”，等于则显示“You won！”，同时退出这个游戏，如果用户7次还没有猜出，则显示“You lost！”，并 第10章 Matlab 在最优化问题中的应用 优化理论是一门实践性很强的学科，广泛应用于生产管理、军事指挥和科学试验等各种领域，Matlab 优化工具箱提供了对各种优化问题的一个完整的解决方案。 在数学上，所谓优化问题，就是求解如下形式的最优解： Min fun (x) Sub. to [C.E.] [B.C.] 其中fun (x)称为目标函数，“Sub. to ”为“subject to ”的缩写，由其引导的部分称为约束条件。[C.E.]表示Condition Equations ，即条件方程，可为等式方程，也可为不等式方程。[B.C.]表示Boundary Conditions ，即边界条件，用来约束自变量的求解域，以lb ≤x ≤ub 的形式给出。当[C.E.]为空时，此优化问题称为自由优化或无约束优化问题；当[C.E.]不空时，称为有约束优化或强约束优化问题。 在优化问题中，根据变量、目标函数和约束函数的不同，可以将问题大致分为： ·线性优化 目标函数和约束函数均为线性函数。 ·二次优化 目标函数为二次函数，而约束条件为线性方程。线性优化和二次优化 统称为简单优化。 ·非线性优化 目标函数为非二次的非线性函数，或约束条件为非线性方程。 ·多目标优化 目标函数并非一个时，称为多目标优化问题。 本章将对以上几类优化问题在Matlab 中的实现作比较详细的讲解。另外还将介绍两个利用优化方法解非线性方程的函数。 通过本章的介绍，用户可以不必掌握艰涩的各种优化算法而轻易地解决一些常用的最优化问题了。 10.1 线性规划问题 线性规划问题即目标函数和约束条件均为线性函数的问题。 其标准形式为： min C ’x sub. To Ax = b x ≥0 其中C, b, 0∈R n ，A ∈R m ?n ，均为数值矩阵，x ∈R n 。 若目标函数为：max C ’x ，则转换成：min –C ’x 。 标准形式的线性规划问题简称为LP （Linear Programming ）问题。其它形式的线性规划问题经过适当的变换均可以化为此种标准形。线性规划问题虽然简单，但在工农业及其他生产部门中应用十分广泛。 在Matlab 中，线性规划问题由linprog 函数求解。 函数：linprog %求解如下形式的线性规划问题： x f T x min such that b x A ≤? beq x Aeq =? ub x lb ≤≤ 其中f, x, b, beq, lb, ub 为向量，A, Aeq 为矩阵。 格式：x = linprog(f,A,b) x = linprog(f,A,b,Aeq,beq) x = linprog(f,A,b,Aeq,beq,lb,ub)第二次作业
第2章--MATLAB数据及其运算-习题答案教学内容
matlab第十章
MATLAB习题
第10章__Matlab在最优化问题中的应用