交大MATLAB第四次上机实验

第四讲上机实验

实验内容

1．用命令或图标进入simulink 模块库，仔细观察各个模块组并分别打开，熟悉其模块的构成。回忆连续模块组中有哪些常用模块？Sum 模块在哪个模块组中？

2．用3种不同的方法，分别打开一空白模型窗口，熟悉有关的菜单及选项，并向其中拖曳一模块，如Step 模块，对其进行大小的调整、方向的旋转、标签名字的改变及隐藏、增加阴影、参数改变等操作，熟悉右键的使用。

3．在Simulink 中构建模型：对一个正弦信号进行积分处理，然后将原始正弦信号和积分后的信号送到示波器中同时显示出来。

4．一系统可由微分方程322()y y y u t ++=来描述，系统初始状态为零，求取系统响应。

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实验报告(1-4)

信号与系统MATLAB第一次实验报告一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。

举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：

MATLAB上机实验(答案)

MATLAB工具软件实验(1) （1）生成一个4×4的随机矩阵，求该矩阵的特征值和特征向量。程序： A=rand(4) [L,D]=eig(A) 结果： A = 0.9501 0.8913 0.8214 0.9218 0.2311 0.7621 0.4447 0.7382 0.6068 0.4565 0.6154 0.1763 0.4860 0.0185 0.7919 0.4057 L = -0.7412 -0.2729 - 0.1338i -0.2729 + 0.1338i -0.5413 -0.3955 -0.2609 - 0.4421i -0.2609 + 0.4421i 0.5416 -0.4062 -0.0833 + 0.4672i -0.0833 - 0.4672i 0.4276 -0.3595 0.6472 0.6472 -0.4804 D = 2.3230 0 0 0 0 0.0914 + 0.4586i 0 0 0 0 0.0914 - 0.4586i 0 0 0 0 0.2275 （2）给出一系列的a值，采用函数 22 22 1 25 x y a a += - 画一组椭圆。程序： a=0.5:0.5:4.5; % a的绝对值不能大于5 t=[0:pi/50:2*pi]'; % 用参数t表示椭圆方程 X=cos(t)*a; Y=sin(t)*sqrt(25-a.^2); plot(X,Y) 结果：（3）X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2], （a）写出计算其负元素个数的程序。程序： X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2]; L=X<0; A=sum(L) 结果： A =

Matlab上机实验

Matlab 上机实验一、实验目的 1、掌握绘制MATLAB 二维、三维和特殊图形的常用函数； 2、熟悉并掌握图像输入、输出及其常用处理的函数。二、实验内容 1 绘制函数的网格图和等高线图。42 2cos cos y x ye x z +-= 其中x 的21个值均匀分布在[-5,5]范围，y 的31个值均匀分布在[0,10]，要求将产生的网格图和等高线图画在同一个图形窗口上。 2 绘制三维曲面图，使用纯铜色调色图阵进行着色，并进行插值着色处理。?????===s z t s y t s x sin sin cos cos cos 230,20ππ≤≤≤≤t s 3 已知 ???????>++≤+=0),1ln(210,22x x x x e x y π 在-5<=x<=5区间绘制函数曲线。 4 已知y1=x2，y2=cos(2x)，y3＝y1*y2，其中x 为取值-2π~2π的等差数列(每次增加0.02π)，完成下列操作： a) 在同一坐标系下用不同的颜色和线型绘制三条曲线，给三条曲线添加图例； b) 以子图形式，分别用条形图、阶梯图、杆图绘制三条曲线，并分别给三个图形添加标题“y1=x^2”，“y2=cos(2x)”和“y3＝

y1*y2”。 5 在xy 平面内选择区域[][],,-?-8888 ，绘制函数 z =的三种三维曲面图。 6 在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注 “y=sin(x)”, “y=cos(x)” ,x 轴,y 轴,标题为“正弦余弦函数图象”. 7 分别用线框图和曲面图表现函数z=cos(x)sin(y)/y ，其中x 的取值为 [-1.5pi,1.5pi]，y=x ，要求：要有标题、坐标轴标签 8 有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线，并加入标题和图列框（用代码形式生成） 9 2 2y x xe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图 10 x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。 11 用sphere 函数产生球表面坐标，绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。 12 以自己的个人画像或照片(JPG)为对象，读入该图像并了解图像的信息，同时利用所学函数对其进行灰度、二值、旋转及缩放等处理，并以PNG 形式输出。

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 （3）矩阵的转置及共轭转置

matlab上机实验指导书

MATLAB应用基础实验指导书

第一章 MATLAB及其工作环境介绍 (1) 1.1 MATLAB简介 (1) 1.2 MATLAB的工作环境介绍 (1) 1.3 MATLAB的基本管理命令 (4) 第二章 MATLAB的数值计算功能 (5) 2.1 变量与赋值语句 (5) 2.2 MATLAB矩阵 (5) 2.3 MATLAB表达式 (10) 2.4 MATLAB常用数学函数 (11) 2.5 矩阵的基本运算 (12) 2.6 数组运算 (16) 2.7 多项式及其运算 (17) 第三章 MATLAB程序设计入门 (19) 3.1 M文件 (19) 3.2 数据的输入输出 (21) 3.3 全局变量和局部变量 (23) 3.4 程序流程控制 (23) 第四章 MATLAB的符号运算功能 (28) 4.1 建立符号对象 (28) 4.2 符号算术运算 (29) 4.3 符号微积分运算 (32) 4.4 符号函数的可视化 (34) 第五章 MATLAB的可视化功能 (37) 5.1 二维图形 (37) 5.2绘制三维图形 (42) 5.3 特殊坐标图形 (44) 5.4 图形句柄 (45)

第一章 MATLAB及其工作环境介绍 1.1 MATLAB简介 MATLAB是matrix和laboratory前三个字母的缩写，意思是实验室矩阵。MATLAB 语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，自1984年由美国MathWorks公司推向市场以来，经过十多年的发展与完善，MATLAB已发展成为由MATLAB语言、MATLAB工作环境、MATLAB图象处理系统、MATLAB数学函数库和MATLAB 应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的体系。MATLAB由“主包”和三十多个扩展功能和应用学科性的工具箱组成。 MATLAB具有以下基本功能： ●数值计算功能 ●符号计算功能 ●图形处理及可视化功能 ●可视化建模及动态仿真功能 MATLAB语言是以矩阵计算为基础的程序设计语言，语法规则简单易学。其指令格式与数学表达式非常相近，用MATLAB编写程序犹如在便笺上列写公式和求解，因而被称为“便笺式”的编程语言。另外，MATLAB还具有功能丰富和完备的数学函数库及工具箱，大量繁杂的数学运算和分析可通过调用MATLAB函数直接求解，大大提高效率，其程序编译和执行速度远远超过了传统的C和FORTRAN语言，因而用MATLAB 编写程序，往往可以达到事半功倍的效果。在图形处理方面，MATLAB可以给数据以二维、三维乃至四维的直观表现，并在图形色彩、视角、品性等方面具有较强的渲染和控制能力，使技术人员对大量原始数据的分析变得轻松和得心应手。 MATLAB的上述特点，使它深受工程技术人员及科技专家的欢迎，并成为应用学科计算机辅助分析、设计、仿真、教学等领域不可缺少的基础软件。目前MATLAB已成为国际上公认的最优秀的科技应用软件。 1.2 MATLAB的工作环境介绍一、MATLAB的工作环境

实验二 MATLAB程序设计含实验报告

实验二 MATLAB 程序设计一、实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。二、实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写：启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正三、实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因 c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。（提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。要求：（1）用switch 语句实现。（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。（提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 的值，调用该函数后，

河南城建学院MATLAB上机实验答案

一熟悉Matlab工作环境 1、熟悉Matlab的5个基本窗口思考题：（1）变量如何声明，变量名须遵守什么规则、是否区分大小写。答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则：?变量名必须以字母开头，只能由字母、数字或下划线组成。 ?变量名区分大小写。 ?变量名不能超过63个字符。 ?关键字不能作为变量名。 ?最好不要用特殊常量作为变量名。（2）试说明分号、逗号、冒号的用法。分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。（3）linspace（）称为“线性等分”函数，说明它的用法。 LINSPACE Linearly spaced vector. 线性等分函数 LINSPACE(X1, X2) generates a row vector of 100 linearly equally spaced points between X1 and X2. 以X1为首元素，X2为末元素平均生成100个元素的行向量。 LINSPACE(X1, X2, N) generates N points between X1 and 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上机指导书 ************************ 昆明理工大学机电学院彭用新 2015年3月

实验三符号计算一、操作部分：在命令窗口执行命令完成以下运算，记录运算结果。 1.findsym：帮助我们获取系统定义的自变量 f= sym('sin(a*x+b*y)'); findsym(f) 2.numden（获取分子分母）, sym2poly,（获取多项式时系数）poly2sym（根据多项式系数获得符号表达式） [n,d]=numden(sym('x*x+y')+sym('y^2')) p=sym('2*x^3+3*x^2+4'); sym2poly(p) x=[2,3,0,4]; poly2sym(x) 3. collect ：合并同类项；expand：展开多项式；horner: 分解成嵌套形式；factor：因式分解；simplify: 对表达式化简 syms x y; collect(x^2*y+y*x-x^2-2*x) collect((x+y)*(x^2+y^2+1), y) syms x y; expand((x-2)*(x-4)) syms x;horner(x^3-6*x^2+11*x-6) syms x;factor(x^3-6*x^2+11*x-6) syms x;simplify((x^2+5*x+6)/(x+2)) 4. finverse ：求得符号函数的反函数。 syms x y; finverse(1/tan(x)) f= x^2+y; finverse(f,y) finverse(f) https://www.360docs.net/doc/733614021.html,pose 求符号函数的复合函数 syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g) 6. subs :表达式替换。 syms a b;subs(a+b,a,4)

MATLAB第二次上机实验报告

电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析学生姓名：学号：指导教师：

一、实验名称实验二线性方程组求解和函数的数值逼近二、实验目的通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。实验涉及的核心知识点：病态方程求解、矩阵分解和方程组求解、Lagrange 插值。实验重点与难点：算法设计和MATLAB 编程三、实验内容 1. 对高阶多项式 ()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 2. 对2 20n =，生成对应的Hilbert 矩阵，计算矩阵的条件数；通过先确定解获得常向量 b 的方法，确定方程组 ()n H x b = 最后，用矩阵分解方法求解方程组，并分析计算结果。 3. 对函数 ()2 1 125f x x = + []1,1x ∈- 的Chebyshev 点 ()()21cos 21k k x n π ?? -= ? ?+? ? ，1,2,,1k n =+ 编程进行Lagrange 插值，并分析插值结果。四、实验数据及结果分析 1. 对高阶多项式

()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 p=[1,-1]; for i=2:20 n=[1,-i]; p=conv(p,n); % 求多项式乘积 end m=zeros(1,21); % m 的最高次幂为20，有21项 hold on x=1:20; d=[-1,0,0.1,0.5,1]; for i=1:5 delt=d(i); m(2)=delt; y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); end title('方程p(x)=0的解与扰动量delt 的关系') legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1') 2468101214161820 010 20 30 40 50 60 方程p(x)=0的解与扰动量delt 的关系 delt=-1delt=0delt=0.1delt=0.5delt=1

MATLAB上机实验练习题答案

数学建模 MATLAB上机实验练习题 1、给出一个系数矩阵A[234;541;132]，U=[123]，求出线性方程组的一个精确解。 2、给出两组数据x=[00.30.81.11.62.3]’y=[0.820.720.630.600.550.50]’，我们可以简单的认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。 3、解线性方程 4、通过测量得到一组数据： 5、已知一组测量值 6、从某一个过程中通过测量得到：分别采用多项式和指数函数进行曲线拟合。 7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近插值法，双线性插值法，二重三次方插值法） 8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。分别使用不同的颜色，线形和标识符。 9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额

占总销售额的比例，并分离第三种产品的切片。 X=19.322.151.6 34.270.382.4 61.482.990.8 50.554.959.1 29.436.347.0 10、对应时间矢量t，测得一组矢量y t00.30.8 1.1 1.6 2.3 y0.50.82 1.14 1.25 1.35 1.40 采用一个带有线性参数的指数函数进行拟合，y=a0+a1e-t+a2te-t，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。 11、请创建如图所示的结构数组（9分）姓名编号指标江明顺071023身高：176，体重：82 于越忠060134身高：168，体重：74 邓拓050839身高：182，体重：77 12、创建如图所示的元胞数组。（9分） 13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。年份19901992199419961998200020022004200620082010 产量（万吨）75.99591.972105.711123.203131.669150.697179.323203.212226.505249.633256.344 14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。 time=[2.55.07.510.013.017.020.030.040.050.060.070.0] res=[0.290.560.771.051.361.522.002.272.813.053.253.56]

Matlab上机实验答案 (1)

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

>> a=-3.0:0.1:3.0; >> z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) (>> z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)可以验证z3==z33,是否都为1） z3 = Columns 1 through 5 0.7388 + 3.2020i 0.7696 + 3.2020i 0.7871 + 3.2020i 0.7920 + 3.2020i 0.7822 + 3.2020i Columns 6 through 10 0.7602 + 3.2020i 0.7254 + 3.2020i 0.6784 + 3.2020i 0.6206 + 3.2020i 0.5496 + 3.2020i Columns 11 through 20 0.4688 + 3.2020i 0.3780 + 3.2020i 0.2775 + 3.2020i 0.2080 + 3.2020i 0.0497 + 3.2020i

安徽大学matlab第三次试验

Matlab 第三次实验作业 -----------2015年11月17日星期二班级：软件工程学号： E21314003 姓名：李世 1. 设计一个程序对“北京、上海、广州、深圳、南京、武汉、天津、重庆、成都、沈阳、杭州、合肥”12个城市进行随机排序。代码：>>clear >> a={'南京','武汉','天津','重庆','成都','沈阳','北京','上海','广州','深圳', '杭州','合肥'}; a(randperm(numel(a))) ans = Columns 1 through 11 '合肥' '广州' '天津' '深圳' '重庆' '北京' '武汉' '杭州' '上海' '沈阳' '南京' Column 12 '成都 ' 2. 根据下面函数表达式编写M 函数，分别计算x 在1,0,3 1 ,1--处的函数值，并画出函数曲线图

?????????-<<≤-+ ≥=-210 0212 102 1)(2x x x x e x f x >>clear >> x=input('请输入变量值：'); if x>=0 f=0.5*exp((-0.5).*x) else if -1/2<=x&&x<0 f=x+0.5 else f=0 end end 请输入变量值：-1 f = >> 0 ans = >> 1 ans = 1 >> -1/3 ans = -0.3333

a=-3:0.01:-1/2; b=0; c=-1/2:0.01:0; d=c+1/2; e=0:0.01:3; f=1/2*exp(-e/2); plot(a,b,c,d,e,f)

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

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w 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\b 4．设???? ??????------=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' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ????? ?++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。（mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g) 或者 u=reshape(a,1,48); p1=mean(u)

第二次数学实验报告Matlab 二维曲线绘图

《数学实验》报告实验名称 Matlab 二维曲线绘图 2011年 5月

一、【实验目的】学习Matlab 绘图的运用，学会制作二维曲线，三维图形的绘画。二、【实验任务】 P79 第3，5，9题。 1，在同一图形窗口画三个子图…… 2，绘制圆锥螺线的图像并加各种标注…… 3，画三维曲面z=5-x^2-y^2与平面z=3的交线。三、【实验程序】 1. >> clear >> x=-pi:pi/50:4*pi; y1=x.*cos(x); y2=x.*tan(1./x).*sin(x.^3); y3=exp(1./x).*sin(x); subplot(3,1,1) plot(x,y1,'r*'),grid on title('y1=xcosx') xlabel('x轴'),ylabel('y轴') axis([-pi pi -pi pi]) gtext('y1=xcosx'),legend('y1=xcosx') subplot(3,1,2),plot(x,y2,'b'),grid on title('y=xtan(1/x)sin(x^3)') gtext('y=xtan(1/x)sin(x^3)') legend('y=xtan(1/x)sin(x^3)') axis([pi 4*pi -2 2]) subplot(3,1,3),plot(x,y3,'y'),grid on title('y=exp(1/x)sinx') xlabel('x轴'),ylabel('y轴') gtext('y=exp(1/x)sinx') legend('y=exp(1/x)sinx') axis([1 8 -3 3]) 2. >> clear >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2.*t; plot3(x,y,z) title('圆锥螺线') xlabel('x轴'),ylabel('y轴'),zlabel('z轴') >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t);

matlab实验内容答案

实验报告说明： matlab 课程实验需撰写8个实验报告，每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。第一次实验内容：实验一 MATLAB 运算基础一、实验目的 1.熟悉启动和退出MA TLAB 的方法。 2.熟悉MA TLAB 命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握MA TLAB 各种表达式的书写规则以及常用函数的使用。二、实验内容 1.先求下列表达式的值，然后显示MA TLAB 工作空间的使用情况并保存全部变量。（1）2 2sin 8511z e ?= + （2 ）12ln(2 z x =+ ，其中2120.45 5i +? ? =? ?-?? （3）0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02 a a e e z a a --= +=--- 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。（4）2 2 2 01 41 1221 23 t t z t t t t t ?≤

12344347873657A -????=??????，131203327B -???? =????-?? 求下列表达式的值：（1）A+6=B 和A-B+I(其中I 为单位矩阵)。（2）A*B 和A.*B 。（3）A^3和A^.3 。（4）A/B 和B\A 。（5）[A ，B]和[A([1,3],;);B^2] 。 3.设有矩阵A 和B 12345678910111213141516171819202122232425A ????????=????????， 30 161769 23497041311B ?? ?? -?? ??=-?? ????? ? （1）求它们的乘积C 。（2）将矩阵C 的右下角3×2子矩阵赋给D （3）查看MA TLAB 工作空间使用情况。 4.完成下列操作：（1）求[100，999]之间能被21整除的数的个数。提示：先利用冒号表达式，再利用find 和length 函数。（2）建立一个字符串向量，删除其中的大写字母。提示：利用find 函数和空矩阵。第二次实验内容：实验三选择结构程序设计一、实验目的 1. 掌握建立和执行M 文件的方法。 2. 掌握利用if 语句实现选择结构的方法。 3. 掌握利用switch 语句实现多分支选择结构的方法。 4. 掌握try 语句的使用。二、实验内容

MATLAB上机习题一

MATLAB上机习题一请按以下步骤完成上机实验： 1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内； 2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面； 3）如果习题是问答题，请将答案写在题目后； 4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后； 5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”； 6）上传该文件到FTP的相关目录。 1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。 2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。 3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。 4. 通过使用帮助确定内置变量ispc 的含义。 6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。 5. 创建double类型的变量，并进行计算 1）a=87，b=190，计算a+b、a-b、a*b；

2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致。 6. 计算如下表达式： 1）()sin 60 2）3e 3）3cos 4π?? ??? 4）2 7562323336 +?-?? 7. 设u=2，v=3，计算： 1）4log uv v 2）() 2 2u e v v u +- 3 8. 计算如下表达式： 1）()()3542i i -+ 2）()sin 28i - 9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：

1）420< 2）420≤ 3）420== 4）4~20= 5）''''b B < 10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b > 2）a b < 3）&&a b b c >> 4）a d == 5）|a b c > 6）~~d 11. 判断下列变量名哪些是合法的MATLAB 变量： 1）fred 7）fred! 2）book_1 8）book-1 3）2ndplace 9）Second_Place 4）#1 10）No_1 5）vel_5 11）vel.5 6）tan 12）while 提示，可以使用isvarname 对上述变量名进行检验。 12. 编写脚本程序，命名为barbell.m ，完成如下计算：如图1所示的杠铃，

Matlab上机实验二

1. 对Runge 函数2()1/(125)R x x =+在区间[-1,1]作下列插值逼近，并和()R x 的图像进行比较，并对结果进行分析。（1）用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的20次Newton 插值多项式的图像。（2）用节点21cos()(0,1,2...,20)42 i i x i π+==，绘出它的20次Lagrange 插值多项式的图像。

（3）用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的分段线性插值函数的图像。（4）用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的三次自然样条插值函数的图像。

结果分析：等距节点的Newton插值在边界附近呈现不收敛的状态，也就是所谓的Runge 现象。这是因为被插函数2 ()1/(125) R x x =+及其任意阶导数并不是一致有界的。由Lagrange插值的余项 (n1) 1 () ()()()() (1)! n n n f R x f x L x x n ξ ω + + =-= + ，由于高阶导数值非常大，在等距节点的高次插值中，余项的值很大，亦即误差特别大，表现为Runge现象。在实际中，除非被插函数的光滑性十分好，否则等距节点的高次插值一般不使用。切比雪夫多项式零点插值没有出现Runge现象，这是因为对于Lagrange插值的余项 (n1) 1 () ()()()() (1)! n n n f R x f x L x x n ξ ω + + =-= + ，切比雪夫多项式零点插值可以将 1() n x ω +项控制在 1 2n 以下（ 1 [1,1] max() n x x ω + ∈- = 1 2n ），由此可得切比雪夫多项式零点插值的误差范围 (n1) [1,1] () max() 2(1)! n n x f R x n ξ + ∞ ∈- ≤ + ，比等距节点的高次插值小很多，因此对于被插函数的光滑性要求要比等距节点的插值也低很多。分段低次线性插值是常用的一种插值方式。线性插值的优点在于计算量小，但是从图中可以看出，其光滑性不够。三次样条插值要求其一阶导数是连续的，因此其准确性好，光滑性也较优。

matlab实验报告

MATLAB 数学实验报告指导老师：班级：小组成员：时间：201_/_/_

Matlab第二次实验报告小组成员： 1题目：实验四，MATLAB选择结构与应用实验目的：掌握if选择结构与程序流程控制，重点掌握break，return，pause语句的应用。问题：问题1：验证“哥德巴赫猜想”，即：任何一个正偶数（n>=6）均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。问题分析：由用户输入一个大于6的偶数，由input语句实现。由if判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数，再向屏幕输出两个质数即可。编程：function [z1,z2]=gede(n); n=input('please input n')

if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n-i) break end

end end 结果分析：如上图，用户输入了大于6的偶数返回两个质数5和31，通过不断试验，即可验证哥德巴赫猜想。纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰，更快的解决问题。 2题目：实验四，MATLAB选择结构与应用实验目的：用matlab联系生活实际，解决一些生活中常见的实际问题。