中南大学材料学院MATLAB题库答案3.0版

合集下载

MATLAB课后习题集附标准答案

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操作题集答案

中南大学材料院matlab操作题集答案

的491操作题集操作题1 《Matlab与科学计算》P28 例2.1要求计算水在温度为0℃,20℃,40℃,60℃,80℃时的黏度,已知水的黏度随/(1+at+bt2)温度的变化公式为μ=μ为0℃水的黏度,值为1.758*10-3。

其中μ要求分别以长短格式显示计算结果。

操作题2 《Matlab与科学计算》P31 例2.2数值数组【1:5】和字符串转换示例。

操作题3《Matlab与科学计算》P44 例2.9比较用左除和右除法分别求解恰定方程的解。

操作题4 《Matlab与科学计算》 P48 例2.14计算矩阵magic(3)的指数,并比较不同函数的结果。

指数运算expm(a)操作题5 《Matlab与科学计算》 P50 例2.18计算下面矩阵[-149 -50 -154;537 180 546;-27 -9 -25]的特征值条件数。

特征值函数condeig(A)操作题6《Matlab与科学计算》P62 例2.29对矩阵a = 1 1 1 11 2 3 41 3 6 101 4 10 20 进行三角抽取使之分别变成0 0 0 0 1 1 1 01 0 0 0 及 123 41 3 0 0 1 3 6 101 4 10 0 1 4 10 20操作题7《Matlab与科学计算》P62 例2.30用3种方法建立向量v 为[1 2 6 20]多项式的伴随矩阵。

操作题8《Matlab与科学计算》P64 例2.31已知a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5];求该矩阵的3次方及数组的3次方。

操作题9《Matlab与科学计算》P66 例2.32 数组逻辑运算演示。

已知a=[1:3;4:6;7:9],b=[0 1 0;1 0 1;0 0 1], x=5, y=ones (3)*5,求x<=a,ab=a&b,~b操作题10《Matlab与科学计算》P71 例2.41对多项式p =[2 -5 6 -1 9]求根,计算其在x=3的值,进行微分计算并写出其表达式。

MATLAB习题及参考答案经典.doc

MATLAB习题及参考答案经典.doc

习题:1, 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。

2, 对于B AX =,如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ,求解X 。

3, 已知:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ,分别计算a 的数组平方和矩阵平方,并观察其结果。

4, 角度[]604530=x ,求x 的正弦、余弦、正切和余切。

(应用sin,cos,tan.cot)5, 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵: (1)组合成一个4⨯3的矩阵,第一列为按列顺序排列的a 矩阵元素,第二列为按列顺序排列的b 矩阵元素,第三列为按列顺序排列的c 矩阵元素,即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574(2)按照a 、b 、c 的列顺序组合成一个行矢量,即 []2965318772546, 将(x -6)(x -3)(x -8)展开为系数多项式的形式。

(应用poly,polyvalm)7, 求解多项式x 3-7x 2+2x +40的根。

(应用roots)8, 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。

(应用poly,polyvalm)9, 计算多项式9514124234++--x x x x 的微分和积分。

(应用polyder,polyint ,poly2sym)10, 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。

(应用x=a\b)11, 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。

(应用pinv) 12, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ,计算a 的行列式和逆矩阵。

(应用det,inv)13, y =sin(x ),x 从0到2π,∆x =0.02π,求y 的最大值、最小值、均值和标准差。

matlab考试题目及答案

matlab考试题目及答案

matlab考试题目及答案1. 题目:编写一个MATLAB函数,实现计算并返回一个向量中所有元素的平方和。

答案:以下是一个简单的MATLAB函数,用于计算并返回一个向量中所有元素的平方和。

```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目:给定一个3x3的矩阵A,使用MATLAB编写代码,求出矩阵A 的转置。

答案:可以通过简单的转置操作来求得矩阵A的转置。

```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目:编写一个MATLAB脚本,实现对一个二维数组进行排序,并返回排序后的数组。

答案:以下是一个MATLAB脚本,用于对一个二维数组进行排序,并返回排序后的数组。

```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目:给定一个向量x,使用MATLAB编写代码,计算并返回向量x的元素个数。

答案:可以通过内置函数`numel`来计算向量x的元素个数。

```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目:编写一个MATLAB函数,实现计算并返回两个向量元素的点积。

答案:以下是一个简单的MATLAB函数,用于计算两个向量的点积。

```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目:给定一个矩阵B,使用MATLAB编写代码,求出矩阵B的行列式。

中南大学matlab课后习题(10)

中南大学matlab课后习题(10)

中南大学matlab课后习题(10)Unit 1实验内容1.答:用help命令可以查询到自己的工作目录。

输入help命令:help <函数名>2.答:MATLAB的主要优点:通过例1-1至例1-4的验证,MATLAB的优点是MATLAB以矩阵作为数据操作的基本单位,使得矩阵运算变得非常简捷,方便,高效。

还提供了丰富的数值计算函数。

MATLAB绘图十分方便,只需输入绘图命令,MATLAB便可自动绘出图形。

3.答:INV(X) is the inverse of the square matrix X。

A warning message is printed if X is badly scaled or nearly singular. PLOT(X,Y) plots vector Y versus vector X. If X or Y is a matrix, then the vector is plotted versus the rows or columns of the matrix, whichever line up. If X is a scalar and Y is a vector, length(Y) disconnected points are plotted. PLOT(Y) plots the columns of Y versus their index. If Y is complex, PLOT(Y) is equivalent to PLOT(real(Y),imag(Y)).In all other uses of PLOT, the imaginary part is ignored. For vectors, MAX(X) is the largest element in X. For matrices,MAX(X) is a row vector containing the maximum element from each column. For N-D arrays, MAX(X) operates along the first non-singleton dimension. [Y,I] = MAX(X) returns the indices of the maximum values in vector I. If the values along the first non-singleton dimension contain more than one maximal element, the index of the first one is returned. ROUND(X) rounds the elements of X to the nearest integers. MAX(X,Y) returns an array the same size as X and Y with the largest elements taken from X or Y. Either one can be a scalar。

Matlab编程与应用习题和一些参考答案

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二1.安装Matlab 软件。

2.验证所学内容和教材上的例子。

3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x>> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59944.设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ,求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。

>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'C1 =19 -82 3012 27 3-38 54 29>> C2=A'*BC2 =-15 16 -24 3663 -17 93 -10522 6 117 -6019 46 84 -10>> C3=A.*BC3 =5 16 24 -26-18 -12 -15 72-2 -21 108 -56>> inv(C1)ans =0.0062 0.0400 -0.0106-0.0046 0.0169 0.00300.0168 0.0209 0.0150>> inv(C2)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.e-019.ans =1.0e+015 *-0.9553 -0.2391 -0.1997 0.27000.9667 0.2420 0.2021 -0.2732-0.4473 -0.1120 -0.0935 0.1264-1.1259 -0.2818 -0.2353 0.3182>> inv(C3)??? Error using ==> invMatrix must be square.5.设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ,把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线。

Matlab考试题库+答案

Matlab考试题库+答案

367
25.已知 A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写______; A ./ B 的运行结果_______。 ans=[0 1 1;1 1 0] ans= [0 2.25 2;1 0.6 NaN] 26. 产生 4 阶全 0 方阵的命令为_______; 产生 3 阶全 1 方阵的命令为。 zeros(4) ones(3) 27. A=rand(2,5);b=size(A);c=length(A);则 b 和 c 的值分别 为_______。 2 5和 5 28. 用 if 判断语句判断 80≥ x >60,在 MATLAB 中 if 语言后的 判断应写为_______。 if x<=80 & x>60 29.为了使两个 pl ot 的图形在同一个坐标显示,可以使用 _______命令进行图形保持;可以使用_______命令为图形添加网 格。 hold on grid on 30. MATLAB 的工作空间中有三个变量 v1, v2, v3,写出把它们 保存到文件 my_data.mat 中的指令_______;
容是不同的。 15. 产生四维单位矩阵的语句为(). A.ones(4) B.eye(4) C.zeros(4) D.rand(4) >> eye(4) ans = 1000 0100 0010 0001 >> ones(4) ans = 1111 1111 1111 1111 >> zeros(4) ans = 0000 0000
输入 12.已知 A=[1 2 3;4 5 0;7 8 9];B=[1 0 3;1 5 0;0 1 2]; 写出下列各指令运行的结果。 A+B; A.*B; A==B ; ans= [2,2,6;5,10,0;7,9,11] ans= [1,0,9;4,25,0;0,8,18] ans= [1,0,1;0,1,1;0,0,0] 13.已知 A 是矩阵,求 A 的对角矩阵函数是_______, 求 A 的下三角矩阵函数是_______。 diag tril 14.MATLAB 的程序文件的扩展名分别是_______。 .m 15.MATLAB 最基本的绘图函数为_______。 plot() 16. A=[1,2,3;4,5,6]; A(:,[1,3])=[];A=__________________ [2;5] 17. fix(-1.5)=________, round(-1.5)____________. -1 -2 18.已知 A=[1 3;4 6];C=[A,eye(size(A));A ,ones(size(A))]

Matlab 题与答案

Matlab 题与答案

MATLAB作业11 、在你的机器上安装MATLAB 语言环境,并键入demo 命令,由给出的菜单系统和对话框原型演示程序,领略MATLAB 语言在求解数学问题方面的能力与方法。

>> demo>>左击MATLAB前‘+’,Programming前‘+’,左击Manipulating Multidimensional ArraysMATLAB具备高性能的数值计算和可视化的科学工程计算功能,支持解释性语言输入,编程实现简单,具有丰富的数学函数功能支持。

MATLAB允许与C、Fortran语言接口,其部件Simulink可以采用图形输入的方式来搭构所研究的系统。

包含丰富的工具包,在系统仿真、数字信号处理、图形图像分析、数理统计、通信及自动控制领域得到广泛的应用。

总的来说,该软件有三大特点。

一是功能强大。

具有数值计算和符号计算、计算结果和编程可视化、数学和文字统一处理、离线和在线计算等功能;二是界面友善、语言自然。

MATLAB以复数处理作为计算单元,指令表达与标准教科书的数学表达式相近;三是开放性强。

该公司本身就推出了30多个应用工具箱,而世界上超过200家公司开发出与MATLAB兼容的第三方产品,这些产品向用户提供更多的工具箱、模块集、与其他商业产品的接口等。

2、启动MATLAB 环境,并给出语句tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),toc,试运行该语句,观察得出的结果,并利用help 命令对你不熟悉的语句进行帮助信息查询,逐条给出上述程序段与结果的解释。

>> tic, A=rand(500); B=inv(A); norm(A*B-eye(500)),tocans =1.2333e-012Elapsed time is 0. seconds.tic, A=rand(500);B=inv(A);norm(A*B-eye(500)),toc3、试用符号元素工具箱支持的方式表达多项式5432()34236f x x x x x x =+++++,并令11s x s -=+,将f(x) 替换成s 的函数。

Matlab试题和答案

Matlab试题和答案

M a t l a b试题和答案(总7页)--本页仅作为文档封面,使用时请直接删除即可----内页可以根据需求调整合适字体及大小--Matlab 上机考试试题考试要求:1、从10道题目中随机抽取3道独立完成,时间1小时。

(输入randperm (10),取前三个数)2、每个题目兴建一个.m的文件,命名方式ks+N.m(N为题号),然后将所选三个题目放入一个文件夹,文件名为学号+姓名。

考试完成后将文件夹通过FTP 提交。

3、考试完成后要写一份报告,内容包括以下:(建一个.Doc的文档,文件名为学号+姓名)(1)题号,题目;(2)运行结果及其分析;(3)图也要粘贴在文档中。

4、查阅资料写一篇2000字左右的关于matlab在电子信息中的应用的小论文或综述,也可以具体的写matlab在电子信息中某一个方面或某一个点的应用。

(打印或手写都可,打印版要交电子文档)5、所有要交的东西在1月3号之前必须交齐。

(由学习委员统一收齐交给我,电子文档也拷到学习委员处,统一拷给我)。

所交项目包括:考试报告打印版,小论文打印版(两个装订在一起,考试报告在上,小论文在下,最好做一个统一的封皮),考试报告doc文档,小论文doc文档。

Matlab 上机考试试题1.求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3X +8z- 5w=9-6x+5y-2z+10w=-8(1)求系数矩阵的秩; (2)求出方程组的解。

2.在[-10,10;-10,10]范围内画出函数2222 sinyx yx z++=的三维图形。

3.试画出系统321()221H s s s s =+++的零极点分布图,判断系统是否稳定,同时求其单位冲激响应和频率响应(幅频特性和相频特性)。

4. 将一个屏幕分4幅,选择合适的步长在右上幅与左下幅绘制出下列函数的图形。

(1)]22[)cos(ππ,,-∈x x (曲线图); (2)4)y 2,-4x (-242),(2222≤≤≤≤+=;y x y x f (曲面图)。

中南大学材料学院MATLAB题库答案3.0版

中南大学材料学院MATLAB题库答案3.0版

P.28 ex 2.1 数值变量的运算>> format short e>> clear>> muw0=1.785e-3;>> a=0.03368;>> b=0.000221;>> t=0:20:80;>> muw=muw0./(1+a*t+b*t.^2)muw =1.7850e-003 1.0131e-003 6.6092e-004 4.6772e-004 3.4940e-004p.31 ex 2.2 数值驻足和字符串转换>> a=[1:5];>> b=num2str(a);>> aa =1 2 3 4 5>> bb =1 2 3 4 5>> b*2ans =98 64 64 100 64 64 102 64 64 104 64 64 106P.44 ex 2.9 比较用左除和右除分别求解恰定方程(线性方程组如果方程数等于未知数个数,叫做恰定方程组,如果方程多于未知数,叫做超定方程组,反之称为欠定。

换个角度说,系数矩阵如果是方阵,就是恰定方程组)的解见课本P.48 ex 2.14 计算矩阵的指数并比较不同函数的结果>> b=magic(3);>> expm(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897 >> expmdemo2(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897>> expm1(b)ans =1.0e+003 *2.9800 0.0017 0.40240.0191 0.1474 1.09560.0536 8.1021 0.0064>> expmdemo3(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897P50 ex 2.18 计算矩阵的特征值条件数>> a=rand(3)a =0.9649 0.9572 0.14190.1576 0.4854 0.42180.9706 0.8003 0.9157>> [V,D,s]=condeig(a)V =0.4913 0.6696 0.66960.3158 -0.4476 + 0.2831i -0.4476 - 0.2831i0.8117 -0.2332 - 0.4655i -0.2332 + 0.4655iD =1.8146 0 00 0.2757 + 0.3061i 00 0 0.2757 - 0.3061is =1.17921.23661.2366P62 ex 2.29矩阵的抽取、三角抽取>> a=pascal(4)a =1 1 1 11 2 3 41 3 6 101 4 10 20>> diag(a,-2)ans =14>> v=diag(diag(a))v =1 0 0 00 2 0 00 0 6 00 0 0 20%diag简单来说就是抽取矩阵各对角线上的元素,如上是抽取主对角线以下第二条对角线之元素,其另一功能是建立对角矩阵>> tril(a)ans =1 0 0 01 2 0 01 3 6 01 4 10 20>> triu(a)ans =1 1 1 10 2 3 40 0 6 100 0 0 20%triu&tril用法与diag非常类似,用途是提取下、上三角矩阵P62 ex2.30 建立多项式之伴随矩阵这道题有点凌乱了…….求解释P64 ex 2.31 数组的幂运算>> a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5];>> a^3ans =32 -28 -101 3499 -12 -151 239-1 49 93 851 -17 -98 139>> a.^3ans =8 1 -27 -127 1 0 343-1 8 64 -81 0 -1 125P66 ex 2.32 数组之逻辑运算>> a=[1:3;4:6;7:9];>> b=[0 1 0;1 0 1;0 0 1];>> x=5;y=ones(3)*5;>> x<=aans =0 0 00 1 11 1 1%此处小于等于,对a中元素和x注意比较,大者为假为0,不大者为1>> ab=a&bab =0 1 01 0 10 0 1%此处为与运算,就是同真才为真(同为非零数)>> ~bans =1 0 10 1 01 1 0%逻辑非运算,即全都非,真变假假变真;还有逻辑或运算,看下面即懂:>> a|bans =1 1 11 1 11 1 1%总结多项式运算的函数:poly:Polynomial with specified roots特征多项式的生成p=poly(a)a 为n阶特征矩阵,所得一般为n阶特征多项式;poly2sym数值2符号;polyval 求多项式的值;roots 求多项式的根;conv 多项式的乘法(向量之卷积)conv(p,d);polyder 多项式微分;polyfit 多项式拟合。

MatLab考试题题库(必做题)(带答案)

MatLab考试题题库(必做题)(带答案)

MatLab考试题题库(必做题)(带答案)一,1.请登陆美国MathWork公司的网站,查看看现在大概有多少本MATLAB-baedbook(以MATLAB为基本软件,来说明各个专业领域的教科书或工具书)。

哪一个领域的MATLAB-baedbook最多中文书共有几本答:1612本,数学方面的最多,中文书共有37本。

2.请在MATLAB中直接输入下列常数,看它们的值是多少:a.ib.jc.epd.infe.nanf.pig.realma某h.realmin依次解为:an=0+an=0+an=an=Infan=NaNan=an=+308an=3.试写一函数regPolygon(n),其功能为画出一个圆心在(0,0)、半径为1的圆,并在圆内画出一个内接正n边形,其中一顶点位于(0,1)。

例如regPolygon(8)可以画出如下之正八边型:解:新建文件如下:functiony=regPolyfon(n)n=8;%要画的n边形R=1;%圆的半径t=0::2某pi;某=R某co(t);y=R某in(t);m=linpace(pi/2,5/2某pi,n+1);某z=R某co(m);yz=R某in(m);holdonplot(某,y,某z,yz);a某i'equal';4.一条参数式的曲线可由下列方程式表示:某=in(t),y=1-co(t)+t/10当t由0变化到4某pi时,请写一个MATLAB的脚本,画出此曲线在某Y平面的轨迹。

解:新建:t=linpace(0,4某pi);某=in(t);y=1-co(t)+t/10;plot(某,y,'-o');5.当一个小圆轮沿着一条曲线行进时,轮缘任一点的轨迹就会产生变化丰富的摆线。

假设小圆轮的半径r=2。

a.当小圆轮绕着一个大圆(半径R=5)的外部滚动时,请画此「圆轮摆线」或「外花瓣线」。

b.重复上小题,但改成在大圆的内部滚动,请画出此「内花瓣线」。

MATLAB习题及参考答案

MATLAB习题及参考答案

2 9 0 1310,解方程组3 4 11 x 6。

(应用x=a\b)2 2 6 611,求欠定方程组294 73 54x68的最小范数解。

(应用pinv)5习题:6 9 3 2 4 1与b 的数组乘积。

2 7 5 4 6 81,计算a2, 对于AX B,如果A 3726,求解X。

283,已知:a 1 2 34 5 6,分别计算a的数组平方和矩阵平方,并观察其结果。

7 8 94,角度x 30 45 60,求x的正弦、余弦、正切和余切。

(应用sin,cos,4 2 7 15 95,将矩阵a 、b 和c 组合成两个新矩阵:5 7 8 36 2(1)组合成一个4 3的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的元素,第三列为按列顺序排列的c矩阵元素,即b矩阵4 7 55 8 62 1 97 3 2(2)按照a、b、c的列顺序组合成一个行矢量,即4527781356926,将(x-6)(x-3)(x-8)展开为系数多项式的形式。

(应用poly,polyvalm)7,求解多项式X3-7X2+2X+40的根。

(应用roots)8,求解在x=8 时多项式(x-1)( x-2) ( x-3)( x-4)的值。

(应用poly,polyvalm)9, 计算多项式4x412x314x25x 9的微分和积分。

(应用polyder,polyint ,poly2sym)2 2计算表达式z 10 x 3 y 5 e x y 的梯度并绘图。

(应用meshgrid, gradient, con tour, holdon, quiver)15,用符号函数法求解方程a t 2+b*t +c=0。

(应用solve )16,用符号计算验证三角等式:(应用syms,simple )用 syms,ezplot)用 plot,title,text,legend)24, x= [66 49 71 56 38] ,绘制饼图,并将第五个切块分离出来。

Matlab考试题库及答案(教师出卷参考专用)

Matlab考试题库及答案(教师出卷参考专用)

Matlab考试题库及答案(教师出卷参考专用)一、选择题1.以下哪个函数用于在Matlab中创建一个图形窗口?A. figureB. plotC. graphD. window答案:A2.在Matlab中,以下哪个选项可以用来定义一个矩阵?A. A = [1 2 3; 4 5 6]B. A = (1, 2, 3, 4, 5, 6)C. A = {1, 2, 3, 4, 5, 6}D. A = 1 2 3; 4 5 6答案:A3.以下哪个函数用于求解线性方程组Ax=b?A. solveB. linsolveC. solve(A, b)D. linsolve(A, b)答案:D4.在Matlab中,如何计算矩阵A和矩阵B的乘积?A. A BB. A \ BC. A . BD. A .\ B答案:A5.以下哪个函数用于在Matlab中绘制三维散点图?A. scatterB. scatter3C. plot3D. bar3答案:B二、填空题1.在Matlab中,要创建一个名为"myfig"的图形窗口,可以使用______函数。

答案:figure('Name', 'myfig')2.在Matlab中,要计算矩阵A的行列式,可以使用______函数。

答案:det(A)3.在Matlab中,若要计算变量x的平方,可以使用______运算符。

答案:.^24.在Matlab中,若要计算矩阵A的逆矩阵,可以使用______函数。

答案:inv(A)5.在Matlab中,要绘制一个正弦波形,可以使用______函数。

答案:plot(sin(x))三、判断题1.在Matlab中,矩阵的索引从1开始计数。

()答案:正确2.在Matlab中,可以使用逻辑运算符"&&"和"||"。

()答案:错误3.在Matlab中,矩阵乘法满足交换律。

中南大学matlab课后习题(1)

中南大学matlab课后习题(1)

第二章1·求下列表达式的值。

(1)w=sqrt(2)*(1+0.34245*10^-6)w=1.4142(2)a=3.5;b=5;c=-9.8;x=(2*pi*a+(c+b)/(pi+a*b*c)-exp(2))/(tan(b+c)+a);xx =0.9829(3)a=3.32;b=-7.9;y=2*pi*a^2*[(1-pi/4)*b-(0.8333-pi/4)*a];yy=-128.4271(4)t=[2,1-3i;5,-0.65];z=1/2*exp(2*t)*log(t+sqrt(1+t^2));zz=1.0e+004*0.0057 - 0.0007i 0.0049 - 0.0027i1.9884 -0.3696i 1.7706 - 1.0539i2,已知a,b,求下列表达式的值。

a=[-1,5,-4;0,7,8;3,61,7];b=[8,3,-1;2,5,3;-3,2,0]; (1)a+6*bans=4723-1012 3726-15 73 7a^2-b+eye(3)ans=-18 -217 1722 533 10921867 526(2)a*bans=14 14 16-10 51 21125 328 180a.*bans =-8 15 40 35 24-9 122 0b*aans =-110-157 228533 -1 28(3)a/bans =1.2234 -0.9255 2.9787-0.9468 2.3511 -0.95744.6170 3.8723 13.8936b\aans =-0.5106 -8.6170 -1.12770.7340 17.5745 1.8085-0.8830 -21.2128 0.4043(4)[a,b]ans=-1 5 -4 8 3 -10 7 8 2 5 3361 7 -3 2 0[a([1,3],:);b^2]ans =-15 -43 61 773 37 117 37 13-20 1 93.已知a,完成下列操作。

matlab课后习题参考答案

matlab课后习题参考答案

Matlab课后习题部分参考答案习题二2.student=struct('name',{'zhang','wang','li'},'age',{18,21,[] },'email',{['zh ang@','zhang@'],'',''})5、a=[1 0 0;1 1 0;0 0 1];b=[2 3 4;5 6 7;8 9 10];c=[a b]d=[a;b]c =1 0 023 41 1 0 5 6 70 0 1 8 9 10d =1 0 01 1 00 0 12 3 45 6 78 9 1010、reshape(c,2,9)1 0 1 0 1 5 3 9 71 0 0 02 8 6 4 10 reshape(d,2,9)ans =1 0 5 0 0 6 0 1 71 2 8 1 3 9 0 4 10 (注意:重新排列矩阵,是将原来的矩阵按列排序)11、a.’12、a+b a.*ba\b (即inv(a)*b,考试时,要求能自己手算出a的逆)ans =2 3 43 3 38 9 1013、c=4*ones(3);a>=cb>=c14.strcat('The picture is ','very good')a=['The picture is ','very good']18.double('very good')1.norm(A),det(A),rank(A)2.b=[1 1 1 1 1];c=b';a=[17,24,1,8,50;23,5,7,14,49;4,6,13,20,43;10,12,19,21,62;11,18,25,2,56]; x=a\c习题四1、a=14*rand(1,10)-5;for i=1:9max=i;for j=i+1:10if a(j)>a(max)max=j;endendtemp=a(i);a(i)=a(max);a(max)=temp;enda将上述脚本保存为script41.m2、function y=function42(n)3*randn(1,n)+3;for i=1:n-1max=i;for j=i+1:nif a(j)>a(max)max=j;endendtemp=a(i);a(i)=a(max);a(max)=temp;endy=a;将上述函数保存为function42.m 3、r=input(‘please input 1/2: ’);if r= =1script41;elsen=input(‘please input the length: ’);y=funtion42(n);end4.function f=function1(x,y)if y==1f=sin(x);elseif y==2f=cos(x);elsef=sin(x).*cos(x);end在命令窗口调用函数f=function1(1,3) 5.function result=function3(x,n)result=0;for i=1:nresult=result+sin(x.*i)+(-1).^i.*cos(i.*x); endfunction result=function4(n,t)result=0;A=[1 2 3;0 1 2;0 0 1];for i=1:nresult=result+A*i.*exp(A*i.*t);end习题五3.x1=-2:0.01:2;x2=-2:0.01:2;y1=x1.*sin(x2);y2=x2.*cos(x1);plot3(x1,x2,y 1,'-d',x1,x2,y2,'-d')14、x=8*rand(1,100)-2;y=reshape(x,10,10);save mydata.mat yclearclcload mydatak=y>ones(10);totel=0;for i=1:100if k(i)==1totel=totel+1;endend15、x=rand(5); %生成5*5均匀分布的随机矩阵(该处只作举例,未按题目要求)fid=fopen(‘text.txt’,’w’);count=fwrite(fid,x,’int32’); closestatus=fclose(fid);清除内存,关闭所有窗口fid=fopen(‘text.txt’,’r’);x=fread(fid,[5,5],’int32’); closestatus=fclose(fid); inv(x)16.随便生成一个矩阵,x1=-1:0.2:0.8;x2=-1:0.2:0.8;y1=x1.*sin(x2);y2=x2.*cos(x1);a=reshape([y1;y2],10,10);fid=fopen('table.txt','w'); fprintf (fid,‘%f’,a); fclose(fid) ;清除内存,关闭所有窗口a=fscanf('table.txt','%f') ; exp(a) ;fclose(fid) ;习题六1.c=ploy2str(A);B=[2,0,0,1,3,5];2.x=1:10;B=ployval(A,x);。

matlab课后练习习题及答案详解

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菜单项来达成。

Matlab考试题库答案(教师出卷参考专用)

Matlab考试题库答案(教师出卷参考专用)
ans=[0 1 1;1 1 0] ans= [0 2.25 2;1 0.6 NaN] 26. 产生 4 阶全 0 方阵的命令为_______; 产生 3 阶全 1 方阵的命令为。 zeros(4) ones(3) 27. A=rand(2,5);b=size(A);c=length(A);则 b 和 c 的值分别为_______ 。 2 5和 5 28. 用 if 判断语句判断 80≥ x >60,在 MATLAB 中 if 语言后的判断应写为 _______。 if x<=80 & x>60 29.为了使两个 plot 的图形在同一个坐标显示,可以使用_______命令进行图形 保持;可以使用_______命令为图形添加网格。
1111 1111 1111 1111
>> rand(4) ans =
0.8147 0.9058 0.1270 0.9134
0.6324 0.0975 0.2785 0.5469
0.9575 0.9649 0.1576 0.9706
0.9572 0.4854 0.8003 0.1419
CCaDb
16.求解线性方程组 Ax=b,当 det(A)≠0 时,方程的解是().
a=input('a=?');
b=input('b=?');
c=input('c=?');
d=b*b-4*a*c;
x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];
[1 2 3;1 4 9;0 7 27] 24. A=[4,15,-45,10,6;56,0,17,-45,0]; find(A>=10 & A<=20)=___ __________.

matlab基础考试题及答案

matlab基础考试题及答案

matlab基础考试题及答案1. MATLAB中,如何创建一个名为"myMatrix"的3x3矩阵,其元素为1到9的连续整数?答案:可以使用以下命令创建矩阵:```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用`magic`函数:```matlabmyMatrix = magic(3);```2. 在MATLAB中,如何计算一个向量的元素平方?答案:可以使用`.^`运算符来计算向量元素的平方:```matlabv = [1 2 3 4];v_squared = v.^2;```3. MATLAB中,如何将一个向量逆序?答案:可以使用`flip`函数或者`end`索引来逆序一个向量:```matlabv = [1 2 3 4];v_reversed = flip(v);```或者```matlabv_reversed = v(end:-1:1);```4. 在MATLAB中,如何计算矩阵的行列式?答案:可以使用`det`函数来计算矩阵的行列式:```matlabA = [1 2; 3 4];detA = det(A);```5. MATLAB中,如何使用循环结构来计算1到100的自然数之和?答案:可以使用`for`循环来计算总和:```matlabsum = 0;for i = 1:100sum = sum + i;end```或者使用`sum`函数直接计算:```matlabsum = sum(1:100);```6. 在MATLAB中,如何创建一个复数?答案:可以通过指定实部和虚部来创建一个复数:```matlabcomplexNumber = 3 + 4i;```7. MATLAB中,如何绘制一个正弦波?答案:可以使用`plot`函数和`sin`函数来绘制正弦波:```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);```8. 如何在MATLAB中计算两个矩阵的点积?答案:可以使用`dot`函数或`*`运算符来计算两个矩阵的点积:```matlabA = [1 2; 3 4];B = [5 6; 7 8];dotProduct = dot(A, B);```或者```matlabdotProduct = A * B;```9. MATLAB中,如何使用条件语句来检查一个数是否为偶数?答案:可以使用`if`语句和模运算符来检查一个数是否为偶数:```matlabnum = 10;if mod(num, 2) == 0disp('The number is even.');elsedisp('The number is odd.');end```10. 在MATLAB中,如何将一个字符串转换为小写?答案:可以使用`lower`函数来将字符串转换为小写:```matlabstr = 'Hello World';lowerStr = lower(str);```结束语:以上是MATLAB基础考试题及答案,希望这些题目能帮助你更好地理解和掌握MATLAB的基础知识。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

P.28 ex 2.1 数值变量的运算>> format short e>> clear>> muw0=1.785e-3;>> a=0.03368;>> b=0.000221;>> t=0:20:80;>> muw=muw0./(1+a*t+b*t.^2)muw =1.7850e-003 1.0131e-003 6.6092e-004 4.6772e-004 3.4940e-004p.31 ex 2.2 数值驻足和字符串转换>> a=[1:5];>> b=num2str(a);>> aa =1 2 3 4 5>> bb =1 2 3 4 5>> b*2ans =98 64 64 100 64 64 102 64 64 104 64 64 106P.44 ex 2.9 比较用左除和右除分别求解恰定方程(线性方程组如果方程数等于未知数个数,叫做恰定方程组,如果方程多于未知数,叫做超定方程组,反之称为欠定。

换个角度说,系数矩阵如果是方阵,就是恰定方程组)的解见课本P.48 ex 2.14 计算矩阵的指数并比较不同函数的结果>> b=magic(3);>> expm(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897>> expmdemo2(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897>> expm1(b)ans =1.0e+003 *2.9800 0.0017 0.40240.0191 0.1474 1.09560.0536 8.1021 0.0064>> expmdemo3(b)ans =1.0e+006 *1.0898 1.0896 1.08971.0896 1.0897 1.08971.0896 1.0897 1.0897P50 ex 2.18 计算矩阵的特征值条件数>> a=rand(3)a =0.9649 0.9572 0.14190.1576 0.4854 0.42180.9706 0.8003 0.9157>> [V,D,s]=condeig(a)V =0.4913 0.6696 0.66960.3158 -0.4476 + 0.2831i -0.4476 - 0.2831i0.8117 -0.2332 - 0.4655i -0.2332 + 0.4655iD =1.8146 0 00 0.2757 + 0.3061i 00 0 0.2757 - 0.3061i s =1.17921.23661.2366P62 ex 2.29矩阵的抽取、三角抽取>> a=pascal(4)a =1 1 1 11 2 3 41 3 6 101 4 10 20>> diag(a,-2)ans =14>> v=diag(diag(a))v =1 0 0 00 2 0 00 0 6 00 0 0 20%diag简单来说就是抽取矩阵各对角线上的元素,如上是抽取主对角线以下第二条对角线之元素,其另一功能是建立对角矩阵>> tril(a)ans =1 0 0 01 2 0 01 3 6 01 4 10 20>> triu(a)ans =1 1 1 10 2 3 40 0 6 100 0 0 20%triu&tril用法与diag非常类似,用途是提取下、上三角矩阵P62 ex2.30 建立多项式之伴随矩阵这道题有点凌乱了…….求解释P64 ex 2.31 数组的幂运算>> a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5];>> a^3ans =32 -28 -101 3499 -12 -151 239-1 49 93 851 -17 -98 139>> a.^3ans =8 1 -27 -127 1 0 343-1 8 64 -81 0 -1 125P66 ex 2.32 数组之逻辑运算>> a=[1:3;4:6;7:9];>> b=[0 1 0;1 0 1;0 0 1];>> x=5;y=ones(3)*5;>> x<=aans =0 0 00 1 11 1 1%此处小于等于,对a中元素和x注意比较,大者为假为0,不大者为1>> ab=a&bab =0 1 01 0 10 0 1%此处为与运算,就是同真才为真(同为非零数)>> ~bans =1 0 10 1 01 1 0%逻辑非运算,即全都非,真变假假变真;还有逻辑或运算,看下面即懂:>> a|bans =1 1 11 1 11 1 1%总结多项式运算的函数:poly:Polynomial with specified roots特征多项式的生成p=poly(a)a 为n阶特征矩阵,所得一般为n阶特征多项式;poly2sym数值2符号;polyval 求多项式的值;roots 求多项式的根;conv 多项式的乘法(向量之卷积)conv(p,d);polyder 多项式微分;polyfit 多项式拟合。

P71 ex 2.41 多项式拟合,用5阶多项式对正弦函数进行最小二乘拟合>> x=0:pi/20:pi/2;>> y=sin(x);>> a=polyfit(x,y,5);>> x1=0:pi/30:pi*2;>> y1=sin(x1);>> y2=a(1)*x1.^5+a(2)*x1.^4+a(3)*x1.^2+a(4)*x1.^2+a(5)*x1+a(6);>> plot(x,y,'b-',x1,y2,'r*')>> legend('原曲线','拟合曲线')>> axis([0,7,-1.2,4])>> axis([0,7,-1.2,1.5]) %调整坐标轴显示范围2.符号运算命名和基本运算法则P 79 符号矩阵的运算>> A=sym('[a,b;c,d]')A =[ a, b][ c, d]>> syms a b c d>> B=[a,b;c,d]B =[ a, b][ c, d]>> B-Aans =[ 0, 0][ 0, 0]>> A/Bans =[ 1, 0][ 0, 1]>> A^2ans =[ a^2 + b*c, a*b + b*d][ a*c + c*d, d^2 + b*c]>> A.^2ans =[ a^2, b^2][ c^2, d^2]>> det(A)ans =a*d - b*c%行列式>> inv(A)ans =[ d/(a*d - b*c), -b/(a*d - b*c)] [ -c/(a*d - b*c), a/(a*d - b*c)] %逆>> rank(A)ans =2%秩%SUMMARY 矩阵的分解:奇异值分解[U,S,V]=svd(A);特征值分解[V,D]=eig(A);正交分解[Q,R]=qr(A);三角分解[L,U]=lu(A);P88 ex 3.7 利用函数gradient绘制一个矢量图>> [x,y]=meshgrid(-2:.2:2,-2:.2:2);>> z=x.*exp(-x.^2-y.^2);>> [px,py]=gradient(z,.2,.2);>> contour(z),%等高线绘制函数>> hold on>> quiver(px,py)%矢量图绘制函数%绘图SUMMARY 1.二维plot 不解释;fplot:绘制y=f(x)图形fplot(fname,lims,’s’);ezplot:绘制隐函数图形help吧;极坐标polar,对数坐标semilogx,semilogy,loglog;bar 条形图,pie 饼状图;hist直方图有些疑惑!;grid on/off:控制是否画网格线。

box on/off:控制是否加边框线。

hold on/off控制是否刷新当前轴及图形%2.三维:plot3;meshgrid函数:产生平面区域内的网格坐标矩阵;mesh 画格子;surf 面;figure(n)开窗户;subplot:割图。

3.二维图形函数运用P98 基本绘图命令>> y=rand(100,1);>> plot(y)>> x=rand(100,1);>> z=x+y.*i;>> plot(z)P101 ex4.1 绘制带有显示属性的二维属性>> x=1:0.1*pi:2*pi;>> y=sin(x);>> z=cos(x);>> plot(x,y,'--k',x,z,'-.rd')P104 ex4.5 条状图和矢量图>> x=1:10;>> y=rand(10,1);>> bar(x,y)>> x=0:0.1*pi:2*pi;>> x=0:0.1*pi:2*pi;>> y=x.*sin(x);>> feather(x,y)P104 ex4.6 函数图形绘制>> lim=[0,2*pi,-1,1];>> fplot('[sin(x),cos(x)]',lim)P105 ex 4.7 绘制饼状图>> x=[2,4,6,8];>> pie(x,{'math','english','chinese','music'})4.三维图形函数应用P107 ex4.9 绘制三维螺旋线>> x=0:pi/50:10*pi;>> y=sin(x);z=cos(x);>> plot3(x,y,z)P107 ex 4.10 绘制参数为矩阵的三维图>> [x,y]=meshgrid(-2:0.1:2,-2:0.1:2);>> z=x.*exp(-x.^2-y.^2);>> plot3(x,y,z)P109 ex 4.11 使用mesh函数绘制三维面图>> x=-8:0.5:8;y=x';>> a=ones(size(y))*x;>> b=y*ones(size(x));>> c=sqrt(a.^2+b.^2)+eps;>> z=sin(c)./c;>> mesh(z)P110 ex4.13 meshc函数绘制的三维面图>> [X,Y]=meshgrid([-4:0.5:4]);>> Z=sqrt(X.^2+Y.^2);>> meshc(Z)P111 ex 4.16 绘制三维饼状图>> clear>> x=[2,4,6,8];>> pie(x)>> pie(x,[0,0,1,0])>> pie3(x,[0,0,1,0])>> pie3(x,[0,0,1,1])P113 ex 4.19 绘制如图柱面图>> x=0:pi/20:pi*3;>> r=5+cos(x);>> [a,b,c]=cylinder(r,30);>> mesh(a,b,c)P113 ex 4.20 绘制地球表面的气温分布示意图>> [a,b,c]=sphere(40);>> t=abs(c);>> surf(a,b,c,t);>> axis('equal')>> axis('square')>> colormap('hot')5.图形控制命令P118 ex 4.24 坐标标注函数应用>> x=1:0.1*pi:2*pi;>> y=sin(x);>>plot(x,y)>> xlabel('x(0-2\pi)')>> ylabel('y=sin(x)')>> title('正弦函数','FontSize')Error using title (line 29)Incorrect number of input argumentsError in title (line 23)h = title(gca,varargin{:});>> title('正弦函数','FontSize',12)>> title('正弦函数','FontSize',12,'FontWeight',bold) %课本上不知由于版本问题还是什么,与2011b不同,最后尝试用Bold成功解决ADD 已解决,用单引号引起bold即可Undefined function or variable 'bold'.>> title('正弦函数','FontSize',12,'FontWeight',Bold)Undefined function or variable 'Bold'.>> title('正弦函数','FontSize',12,'FontWeight','Bold')>> title('正弦函数','FontSize',12,'FontWeight','Bold','FontName','隶书')P123 ex 4.30 在同一张途中绘制几个三角函数图>> x=0:0.1*pi:2*pi;>> y=sin(x);>> z=cos(x);>> plot(x,y,'-*')>> hold on>> plot (x,z,'-o')>> plot(x,y+z,'-+')>> legend('sin(x)','cos(x)','sin(x)+cos(x)',0)P124 ex 4.31 在4个子图中绘制不同的三角函数>> x=0:0.1*pi:2*pi;>> subplot(2,2,1);>> plot(x,sin(x),'-*')>> title('sin(x)')>> subplot(2,2,2)>> plot(x,cos(x),'-o')>> title('cos(x)')>> subplot(2,2,3)>> plot(x,sin(x).*cos(x),'-x')>> title('sin(x)*cos(x)')>> subplot(2,2,4)>> plot(x,sin(x)+cos(x),'-h')>> title('sin(x)+cos()')>> title('sin(x)+cos(x)')% SUMMARY interp1:1-D data interpolation (table lookup)%yi = interp1(x,Y,xi,method,'extrap') method:nearest/linear/echip(hermite)/spline,etcP222 ex 7.3 正弦曲线的插值实例>> x=0:0.05:10;>> y=sin(x);>> xi=0:.125:10;>> yi=interp1(x,y,xi);%一维插值>> plot(x,y,'*',xi,yi)P227 例7.7二次拟合>> x=[0.5:0.5:3.0];>> y=[1.75 2.45 3.81 4.80 8.00 8.60];>> a=polyfit(x,y,2)%用最小二乘法拟合a =0.4900 1.2501 0.8560>> x1=0.5:0.05:3.0;>> y1=a(1)*x1.^2+a(2)*x1+a(3);%拟合出的函数>> plot(x1,y1,'-or',x,y,'-+b')P228 例7.8拟合,用求解矩阵的方法解,求解超定方程>> xi=19:6:44xi =19 25 31 37 43>> yi=[19.0 32.3 49.0 73.3 98.8];>> format short>> a=polyfit(xi,yi,2)a =0.0635 -0.5932 7.2811>> x1=19:.1:44;>> y1=a(1)*x1.^2+a(2)*x1+a(2);%此处发生了一点错误,将a(3)误输为a(2),后面会有改正>> plot(x1,y1,'-x')>> plot(x1,y1,'-')>> x2=xi.^2;%采用求解矩阵的方法来求解此拟合问题>> x2=[ones(5,1),x2']x2 =1 3611 6251 9611 13691 1849>> ab=x2\yi'ab =-1.35220.0540>> y3=ab(1)+ab(2)*x1.^2;>> hold on;plot(x1,y3,'--g')>> y1=a(1)*x1.^2+a(2)*x1+a(3);>> plot(x1,y1,'-')%从图中可以看出,两种拟合方法较为吻合P233 例7.10 分别用矩形和梯形公式求积分>> y=inline('exp(-0.5*t).*sin(t+pi/6)');%内联函数>> d=pi/1000;>> t=0:d:3*pi;>> nt=length(t)nt =3001>> y1=y(t);>> sc=cumsum(y1)*d;>> scf=sc(nt)scf =0.9016>> format long>> scfscf =0.901618619310013%矩形方法>> z=trapz(y1)*dz =0.900840276606885%梯形方法P237 ex7.11 采用自适应Simpson公式求积分>> f=inline('x./(x.^2+4)');>> quad(f,0,1)ans =0.111571765994935P237 ex 7.12 用quadl求积分>> f=inline('exp(-x/2)');>> f1=quadl(f,1,3)f1 =0.7668>> f2=quad(f,1,3,1e-10)f2 =0.7668>> format long e>> f1f1 =7.668009991284686e-001>> f2f2 =7.668009991284349e-001P237 7.12 quadl gauss-labatto 求积分>> quadl('exp(-x/2)',1,3)ans =。

相关文档
最新文档