当前位置：文档库 › 2016_SHI_06_数学实验-MatLab图形界面编程_V2.0

2016_SHI_06_数学实验-MatLab图形界面编程_V2.0

石振锋

2016 年春季学期

哈工大

计算数学研究所

MatLab 之图形用户界面设计

1 图形用户界面编程3

§Matlab UI 编程◆图形用户界面编程4

§Matlab UI 编程

计算机屏幕

root

图形窗口对象

figure

右键菜单对象uicontextmenu 菜单对象uimenu

控件对象uicontrol 坐标轴对象axes

线对象line 文字对象text 图像对象image 光源对象light

兄

5§Matlab UI 编程

findobj

获取具有指定的属性值的对象的句柄

1. 图形对象句柄

§Matlab UI 编程

【例1】绘制正弦曲线并查看有关对象的句柄

§Matlab UI 编程

2. 图形对象属性及其设定8

§Matlab UI 编程

（2）图形对象的公共属性

§Matlab UI 编程

（3）属性的设置

当创建一个对象时，必须给对象的各种属性赋予必要的属将【例1】窗口位置设为向量p_vect 所指定的值

§Matlab UI 编程

get 函数的调用格式为：V=get(句柄,属性名)

figure('position',[1,1,0.5*V(3),0.5*v(4)])

§Matlab UI 编程

【例2】在同一坐标下绘制红、绿两根不同曲线，希望获得绿

色曲线的句柄，并对其线型重新设置

§Matlab UI 编程

3. 图形对象的创建

§Matlab UI 编程

MATLAB 为每个图形窗口提供了很多属性15

§Matlab UI 编程

6）Position 属性——设定对象的位置和大小，为前两个值为窗口左下角横纵坐标值，后两个元素值为窗口的宽为(1.0,1.0)。该属性将影响一切定义的属性项。

§Matlab UI 编程

8）KeyPressFcn 属性、WindowButtonDownFcn 意义意义

(17

§Matlab UI 编程

【例3】建立一个图形窗口。该图形窗口没有菜单条，标题名称为18

§Matlab UI 编程

（2）坐标轴对象

axes(H) %这里引用的H 必须存在

§Matlab UI 编程

MATLAB 为每个坐标轴对象提供了很多属性外，常见坐标轴对象属性。

实现图形窗口的任意分割。SubWindowByAxis.m

§Matlab UI 编程

（3）曲线对象

建立曲线对象使用line 函数，调用它的命令形式为：DrawCurveWithCurveObj.m

21§Matlab UI 编程

（4）文字对象

使用text 函数可以根据指定位置和属性值添加文字说明，并保存【例7】利用曲线对象绘制曲线并利用文字对象完成标注。

请参考：DrawCurveWithCurveObjAndTextLabeling.m

22§Matlab UI 编程

（5）曲面对象

建立曲面对象使用surface 函数，调用它的命令形式为：【例8】利用曲面对象绘制三维曲面z=sin(y)cos(x)。

请参考：Draw3DSurfWithCurveObj.m

§Matlab UI 编程

（6）区域块对象

区域块对象是指一个任意多边形区域。在DrawPatch.m

24§Matlab UI 编程

（7）方框对象

在MATLAB 中，矩形、圆角矩形、椭圆、圆及对应的区域椭圆和圆各一个，并使用了不同的线型和线宽。请参考：

DrawMultiShape.m

§Matlab UI 编程

1.2 菜单设计

灰色

§Matlab UI 编程

3）Callback 属性——MATLAB 回调字符串，选择菜单项时，回调串传给函数eval 作为对菜单的响应；初始值为空矩阵定义热键，它由Alt-字符激活

§Matlab UI 编程

8）Position 属性——菜单对象的相对位置。顶层菜单从左到右编号，子菜单从上至下编号

yes ：回调串可中断

§Matlab UI 编程

14）Parent 属性——父对象的句柄；如果菜单对象是顶层】建立如下图所示的菜单。

§Matlab UI 编程

3. 快捷菜单

快捷菜单是用鼠标右键单击某对象时在屏幕上弹出的菜函数将该快捷菜单和某图形对象联系起来。以控制曲线的颜色。请参考：

Demo_MenuForDraw.m

§Matlab UI 编程

1.3 界面设计

§Matlab UI 编程

1）BackgroundColor ——控件背景色。时，可编辑文本框是多行文本。

2. 控件对象的属性

§Matlab UI 编程

（6）Position ——位置向量[left bottom width height]height]表示相对于图形对象左下角的控件的左下角位置。（10）FontName 、FontSize 、FontAngle 等用来控制字体的属性

§Matlab UI 编程

（11）Style ——定义控件对象的类型，取值可以为：

pushbutton ：（默认）命令按钮（13）Value ——控件的当前值。文本对象和按钮不设置该属性。

§Matlab UI 编程

（14）ButtonDownFcn ——当控件被选择时，传给函数eval 。初始值为空矩阵

§Matlab UI 编程

1.4 标准对话框函数及其调用

图标’)

Custom ，缺省（默认）36

§Matlab UI 编程

2. 询问对话框

§Matlab UI 编程

3. 文件选择对话框

例如：[f,p]=uigetfile('*.m;*.txt','请选择一个文件38

§Matlab UI 编程

FontAngle: 'italic'

4. 字体设置对话框

39§Matlab UI 编程

5. 颜色设置对话框

，则在图中指向初始颜色向量所定义的颜色位置，且按下“取消”按钮时返回初始颜色向量的值。')

§Matlab UI 编程

6. 警告、错误与帮助信息对话框

h=helpdlg('请选择目标,再进行相关操作！','帮助')

§Matlab UI 编程

2 图形用户界面设计工具GUIDE

为了更方面简洁地进行界面设计，MATLAB 提供了一个可视化图形界面开发环境GUIDE （Graphical User Interface 即可得到下图：

§Matlab UI 编程

单击OK 后，便可以进入GUIDE 设计主窗口。

MATLAB 为GUI 设计一共准备了4种模板，分别是GUI with Uicontrols (带控件对象的GUI with Axes and Menu (带坐标轴与菜单的GUI

43§Matlab UI 编程

Guide 工具栏

界面编辑客户区

1.窗口组成

Edit 、View 、Layout 、Tools 和Help 6对象浏览

运行

客户区大小调整

§Matlab UI 编程

面板控件

§Matlab UI 编程

种函数和代码）。

2.基本操作

§Matlab UI 编程

2.2 GUIDE 常用工具：对选中的控件进行剪切操作。：对选中的控件进行复制并粘贴。

：对选中的控件打开属性编辑窗口对象浏览窗口。

：按下鼠标时控件回调的函数。 DeleteFcn ：定义在对象的删除阶段执行的回调例程。

§Matlab UI 编程

2. 属性编辑器(Property Inspector)48

§Matlab UI 编程

（2）对象的常规信息

§Matlab UI 编程

（3）控件回调函数的执行

BusyAction ：处理回调函数的中断。有两种选项：即 Value ：控件的当前值。

§Matlab UI 编程

3. 对象浏览器(Object Browser)

51§Matlab UI 编程

4. Tab 顺序编辑器(Tab Order Editor)

§Matlab UI 编程

5. 对齐工具(Alignment Tool)

进行调整。

垂直

底—顶

关闭

间距水平

底—顶

右—左

左—左

53§Matlab UI 编程

6. 菜单编辑器(Menu Editor)

弹出式菜单（Context Menu ）

移动菜单项

删除菜单项

弹出式菜单

§Matlab UI 编程

【例16】编写一个3D 立体图的演示程序。它有一个下拉菜单“3D 演示”包含子菜单“三维曲线“三维曲面”，用于显示相应图形。

function mnPlot3_Callback(hObject, eventdata, handles)

[x,y,z]=peaks;plot3(x,y,z);

function mnMesh_Callback(hobject,eventdata,handels)mesh(peaks);

function mnSurf_Callback(hobject,eventdata,handels)surf(peaks);

编写菜单

55§Matlab UI 编程

7. 网格标尺设置编辑器(Grid and Rulers)

界面中添加网格与标尺。

§Matlab UI 编程

8. GUI 应用属性设置编辑器(GUI Options)

固定、比例、自

命令行适应性57

§Matlab UI 编程

2.3 GUIDE 应用示例

使用Guide 进行GUI 的界面设计，就像在纸上绘制界面例如：handles.figure1

§Matlab UI 编程

如果需要在函数间传递数据，可先给建立一个字段，然后调用guidata 综合示例——温度转换

根据公式：

，其中将输入的华氏温度转换成摄氏温度。

§Matlab UI 编程

本程序相应函数中插入的代码：

end

§Matlab UI 编程

【例17】综合示例——初等函数绘制

请参考：myfirst.m 和myfirst.fig

Myfirst.m

&myfirst.fig

§Matlab UI编程

MATLAB数学实验练习题

————————————————————————————————作者：————————————————————————————————日期：

“MATLAB”练习题要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) > ezplot('exp(x)-3*x^2') >> grid on ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2）1sin 02 x x - =至少三个根

>> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2sin cos 0x x x -= 所有根

>> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3,x,0) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x)

数学软件MATLAB实验作业

数学软件与数学实验作业一.《数学软件》练习题（任选12题，其中19－24题至少选2题）： 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0)

matlab数学实验复习题(有标准答案)

复习题 1、写出3 2、i ｎv(A)表示Ａ的逆矩阵； 3、在命令窗口健入 clc,４、在命令窗口健入clea 5、在命令窗口健入６、x=-1:０.2:17、det(A)表示计算Ａ的行列式的值；8、三种插值方法：拉格朗日多项式插值,分段线性插值，三次样条插值。 9、若A=123456789?? ????????,则fliplr （Ａ）＝321654987?????????? Ａ-３=210123456--??????????A ．＾２=149162536496481?????????? tril(A)=100450789?????????? tri ｕ（Ａ，－1)=123456089??????????ｄiag(A ）＝100050009?????????? Ａ(:,2),＝2 58A(3,:)=369 10、nor ｍcd ｆ（1，1，2)=0.5%正态分布mu=1,s ｉgm ａ＝２,x ＝1处的概率ｅ45(@ｆ,［ａ,b ］,x0),中参数的涵义是@fun 是求解方程的函数M 文件,[a,b ］是输入向量即自变量的范围a 为初值,ｘ0为函数的初值，t 为输出指定的[ａ,b],x 为函数值１5、写出下列命令的功能：te ｘt （1,2,‘ｙ=s ｉｎ(x)’

hold on 16fun ｃｔion 开头; 17 ,4) ３，４） 21、设x 是一向量，则）的功能是作出将Ｘ十等分的直方图２2、interp １(［1,2,3]，[3,４,5]，2.5) Ans=4.５ 23、建立一阶微分方程组? ??+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来）二、写出运行结果: １、>＞ｅy ｅ（3，4)=1000 01000010 2、>>s ｉze([1,2,3])=1;3 3、设b=ro ｕnd （unifrnd(-５，5,１，4)),则=3 ５ 2 -5 >>[x,m]=ｍin(b)；x ＝－5；m=4 ，［x,n ］=sort(b ）－５ 2 3 ５ 4 3 １ 2 ｍea ｎ(ｂ）=１.2５，m ｅdian(b)＝2.5,range(ｂ）=10 ４、向量ｂ如上题,则＞>ａn ｙ(b),all(ｂ<2),ａｌl(b<6） Ans ＝１ 0 １５、＞>[5 6;7 8]>[7 8；5 6]=00 11 6、若1234B ??=???? ,则 7、>>ｄｉag(d ｉag （B ）)=10 04 8、＞>[4:-2:１].＊[-1,６]=-4 12 9、>>acos(0.5),a ｔan(1） ans ＝ 1.６5９8 ans=

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“”练习题要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;

>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =

0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1

6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5

8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =

matlab数学实验练习题

Matlab 数学实验实验一插值与拟合实验内容：预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。（1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s 10 ≤≤-=x x y （4）、22),exp(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(t e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面一组t ，V 数据确定0V 和τ。实验二常微分方程数值解试验实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。实验内容：

实验三地图问题 1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量：以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土的近似面积，并与它的精确值41288km2比较。

数学实验练习题(MATLAB)

注意:在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上）. 第一次练习题 1.求解下列各题： 1)30sin lim x mx mx x ->- 2)(4)cos ,1000.0=x mx y e y 求 3）21/2 0mx e dx ?(求近似值，可以先用inline 定义被积函数，然后用quad 命令) 4）4 224x dx m x +? 5 0x =展开(最高次幂为8). 2.对矩阵21102041A m -?? ?= ? ?-?? ，分别求逆矩阵，特征值，特征向量，行列式，并求矩阵,P D （D 是对角矩阵），使得1A PDP -=。 3. 已知2 1(),()2f x e x μσ=--分别在下列条件下画出)(x f 的图形： (1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). 4.画（1）sin 020cos 02100x u t t y u t u t z m ??=≤≤?=?≤≤??=?

(2) sin()03,03z mxy x y =≤≤≤≤ （3）sin()(/100cos )02cos()(/100cos )02sin x t m u t y t m u u z u π π=+?≤≤?=+?≤≤?=? 的图（第4题只要写出程序）. 5．对于方程50.10200 m x x --=，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会. 第二次练习题判断迭代收敛速度的程序 x0=1;stopc=1;eps=10^(-8);a=1;c=1;b=2*c;d=a;k=0; f=inline('(a*x+b)/(c*x+d)'); kmax=100; while stopc>eps&k

MATLAB数学实验100例题解

一元函数微分学实验1 一元函数的图形（基础实验）实验目的通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧. 初等函数的图形 2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势. 解：程序代码： >> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps); plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象：程序代码： >> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps); plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象： cot(x)

4在区间]1,1[-画出函数x y 1 sin =的图形. 解：程序代码： >> x=linspace(-1,1,10000); y=sin(1./x); plot(x,y); axis([-1,1,-2,2]) 图象：二维参数方程作图 6画出参数方程? ??==t t t y t t t x 3cos sin )(5cos cos )(的图形：解：程序代码： >> t=linspace(0,2*pi,100); plot(cos(t).*cos(5*t),sin(t).*cos(3*t)); 图象：

极坐标方程作图 8 作出极坐标方程为10/t e r =的对数螺线的图形. 解：程序代码： >> t=0:0.01:2*pi; r=exp(t/10); polar(log(t+eps),log(r+eps)); 图象： 90270 分段函数作图 10 作出符号函数x y sgn =的图形. 解：程序代码： >> x=linspace(-100,100,10000); y=sign(x); plot(x,y); axis([-100 100 -2 2]);

浅析Matlab数学实验报告

数学实验报告姓名: 班级: 学号: 第一次实验任务过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)

cos()sin(2=-=++y x y x y x ,其中、计算的图形。下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);

matlab与数学实验大作业

《数学实验与MATLAB》 ——综合实验报告实验名称：不同温度下PDLC薄膜的通透性与驱动电压的具体关系式的研究学院：计算机与通信工程学院专业班级：姓名：学号：同组同学： 2014年 6月10日

一、问题引入聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。二、问题分析想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。三、实验数据选择10、20、30摄氏度三个不同温度，其他条件一致。

（1）、10摄氏度实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：

南邮MATLAB数学实验答案(全)

第一次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案

“MATLAB”练习题要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)

1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408

3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans =

数学实验(MATLAB)课后习题答案

数学实验练习2.1 画出下列常见曲线的图形。（其中a=1,b=2,c=3）1、立方抛物线3x y= 解：x=-5:0.1:0;y=(-x).^(1/3); y=-y; x=0:0.1:5; y=[y,x.^(1/3)]; x=[-5:0.1:0,0:0.1:5]; plot(x,y) 2、高斯曲线2x e = y- 解：fplot('exp(-x.^2)',[-5,5])

3、笛卡儿曲线)3(13,13332 2 2 axy y x t at y t at x =++=+= 解：ezplot('x.^3+y.^3-3*x*y',[-5,5]) x y x.3+y.3-3 x y = 0 或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)

4、蔓叶线)(1,13 2 23 2 2x a x y t at y t at x -=+=+= 解：ezplot('y.^2-x.^3/(1-x)',[-5,5]) x y y.2-x.3/(1-x) = 0 或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)

5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi; x=t-sin(t); y=2*(1-cos(t)); plot(x,y) 6、星形线)(sin ,cos 3 23 23 233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3;

(完整版)MATLAB数学实验第二版答案(胡良剑)

数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点 >> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500

数学实验与数学软件(Mathmaticandmatlab)

数学软件与数学实验2013-2014学年度秋季学期期末试卷专业：统计学班级：11级2班学号：20110723 姓名：晏静

一、按要求计算出下列表达式的值 (1)318, 3 162 53 ?? + ? ?? , 21 eπ+, 2.5 tg, 2 log15； (2)给出π的9位和e的10位近似值； (3)求658和4102的最大公约数及35和25的最小公倍数； (4)产生10个0与10之间随机数的一个表; (5)求虚数1453 i i i i +- -的实部,虚部,模,共轭,辐角。 (6)自己运用Table建立两个表，并进行表运算，如连接、并集、交、排序等操作。

二、因式分解 22212321332112322 1 22(1)()()()4;(2)21;x x x x x x x x x x x x x x x +++++---- 解：三、解方程(组) 1234234124234-2+344-+-3(1)+31-73+3 x x x x x x x x x x x x x -=??=? ? +=??+=-? 65432(2)5232002000.x x x x x x -+--++= 四、求极限 () 20 (1)1sin ;(2);(3)56! ctg x n x n n n Lim x Lim n n →→∞ →∞++

(1) (2) (3) 五、求导数 32 22(1)()=ln(x+1+);(2)()=cos 2,; (3)=log (),Z . x f x x f f x e y x y Z xy x y y ???求的导数已知求求关于的二阶导 (1) (2) (3) 六、求下列定积分与不定积分： ()()()12201+sin ln 1+(1);(2);(3)sin (1+cos ) +1(1+)(2+-) x x dx dx x x x x x x ? ? ?2 2-(4)=0,=1,==.y D D x y y x I x e d σ??设是由直线围成的区域，计算的值 (1) (2)

matlab实验报告

Matlab实验报告 ——定积分的近似计算学生姓名：学号：专业：数学与应用数学专业

数学实验报告实验序号:1001114030 日期：2012年10月20日班级应一姓名陈璐学号1001114030 实验名称：定积分的近似运算问题背景描述：利用牛顿—莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适合于被积分函数的原函数能用初等函数表达出来的情形。如果这点办不到或不容易办到，这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只能应用近似方法去计算相应的定积分。实验目的：本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线发。对于定积分的近似数值计算，Matlab有专门函数可用。实验原理与数学模型： 1.sum（a）：求数组a的和。 2.format long：长格式，即屏幕显示15位有效数字。 3.double（）：若输入的是字符则转化为相应的ASCII码；若输入的是整型数之则转化为相应的实型数值。 4.quad（）：抛物线法求数值积分。格式：quad（fun，a，b）。此处的fun是函数，并且

为数值形式，所以使用*、/、^等运算时要在其前加上小数点。 5.trapz():梯形法求数值积分。格式：trapz（x，y）。其中x为带有步长的积分区间;y为数值形式的运算。 6.fprintf（文件地址，格式，写入的变量）：把数据写入指定文件。 7.syms 变量1变量2……：定义变量为符号。 8.sym（'表达式'）：将表达式定义为符号。 9.int（f，v，a，b）：求f关于v积分，积分区间由a到b。 10.subs（f，'x'，a）：将a的值赋给符号表达式f中的x，并计算出值。若简单地使用subs （f），则将f的所有符号变量用可能的数值代入，并计算出值。实验所用软件及版本：Matlab 7.0.1

Matlab数学实验一2015(标准答案版)

Matｌab数学实验一——matlab初体验一、实验目的及意义［1] 熟悉MAＴＬＡB软件的用户环境；［2] 了解MAＴLAB软件的一般目的命令；［３] 掌握MATLＡB数组操作与运算函数；通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。二、实验内容１.认识ｍatｌab的界面和基本操作 2.了解mａtlab的数据输出方式（ｆormat) 3. MATLAＢ软件的数组(矩阵）操作及运算练习; 三、实验任务根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式和要求填写实验报告 1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。 iｊeps inf nａn pi realｍaｘrealｍin 2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。 3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。 4. 写出在命令窗口中的计算步骤和运行结果。（1)计算 1.22 10 (ln log) 81 e ππ +- ； >＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns = ０．０3４8 (2) >> x=2;y＝４； >> z=x^２+exｐ(x+y)-y*lｏg(x)-3 z ＝ 40１.65６2 （3)输入变量 13 5.3, 25 a b ?? ==?? ?? ，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0 １.ｍat”文件。测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3 ａ＝

matlab数学实验复习题(有答案)

复习题 1、写出3个常用的绘图函数命令 2、inv （A ）表示A 的逆矩阵； 3、在命令窗口健入clc 4、在命令窗口健入 clear 5、在命令窗口健入6、x=-1：0.2：17、det （A ）表示计算A 的行列式的值；8、三种插值方法：拉格朗日多项式插值，分段线性插值，三次样条插值。 9、若A=123456789?? ???????? ，则fliplr （A ）= 321654987?? ???????? A-3=210123456--??????????A ．^2=149162536496481?????? ???? tril （A ）=100450789?? ???????? triu （A ，-1）=123456089??????????diag （A ）=100050009?? ???? ???? A(:,2),=258A(3,:)=369 10、normcdf （1，1，2）=0.5%正态分布mu=1，sigma=2，x=1处的概率 [t,x]=ode45(f,[a,b],x0),中参数的涵义是fun 是求解方程的函数M 文件，[a,b]是输入向量即自变量的围a 为初值，x0为函数的初值，t 为输出指定function 开头；17、二种数值积分的库函数名为：quad;quadl

4 3，4 21、设x ）的功能是作出将X 十等分的直方图 22、interp1([1,2,3],[3,4,5],2.5) Ans=4.5 23、建立一阶微分方程组???+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来) 二、写出运行结果： 1、>>eye(3，4)=1000 01000010 2、>>size([1,2,3])=1；3 3、设b=round （unifrnd （-5，5，1，4）），则=3 5 2 -5 >>[x,m]=min(b)；x=-5；m=4 ,[x,n]=sort(b) -5 2 3 5 4 3 1 2 mean(b)=1.25，median （b ）=2.5，range （b ）=10 4、向量b 如上题，则 >>any(b),all(b<2),all(b<6) Ans=1 0 1 5、>>[5 6;7 8]>[7 8;5 6]=00 11 6、若1234B ?? =?? ??，则 7、>>diag(diag(B))= 10 04 8、>>[4:-2:1].*[-1,6]=-4 12 9、>>acos(0.5),atan(1) ans= 1.6598 ans= 0.7448 10、>>norm([1,2,3]) Ans=3.3941 11、>>length （[1，3，-1]）=3

数学实验matlab练习题

2015-2016数学实验练习题一、选择题 1.清除Matlab工作空间（wordspace）变量的命令是（B ） A. clc B. clear C. clf D.delete 2. 清除当前屏幕上显示的所有内容，但不清除工作空间中的数据的命令是（ A ） A. clc B. clear C. clf D.delete 3. 用来清除图形的命令（ C ） A. clc B. clear C. clf D.delete 4. 在MATLAB程序中，使命令行不显示运算结果的符号是（ A ） A. ; B. % C. # D. & 5. 在MATLAB程序中，可以将某行表示为注释行的符号是（ B ） A. ; B. % C. # D. & 6.在循环结构中跳出循环，执行循环后面代码的命令为 ( B ) A. return B. break C. continue D. Keyboard 7.在循环结构中跳出循环，但继续下次循环的命令为（ C ） A. return B. break C. continue D. Keyboard 8. MATLAB中用于声明全局变量的关键字是( C ) A. inf B. syms C. global D. function 9. 用户可以通过下面哪项获得指令的使用说明（ A ） A. help B. load C. demo D. lookfor 10．在MATLAB命令窗口中键入命令S=zoros(3)；可生成一个三行三列的零矩阵，如果省略了变量名S，MATLAB表现计算结果将用下面的哪一变量名做缺省变量名（ A ） A. ans; B. pi; C. NaN; D. Eps. 11. 9/0的结果是（ B ） A. NAN； B. Inf； C. eps； D. 0 12．在MATLAB中程序或语句的执行结果都可以用不同格式显示，将数据结果显示为分数形式，用下面哪一条命令语句（ D ） A. format long； B. format long e； C. format bank； D. fromat rat 13. 下列MATLAB命令中是构造1行3列的(-1,1)均匀分布随机矩阵的命令的是（D）

MATLAB与数学实验练习

一、填空 1、命令clear、clf、clc、who、whos的含义分别是clear用于清除内存中的所有变量与函数；clf用于清除图形窗口； clc用于清除命令窗口中的所有显示内容；who将内存中的当前变量以简单的形式列出；whos列出当前内存变量的名称、大小和类型等信息。 2、若矩阵 2310 , 5112 A B ???? == ? ? ???? ,则A.*B= （对应元素相乘）（2 0； 5 2） ,A*B=（矩阵的乘法）（5 6； 6 2）。 3、生成4行3列的元素全为1的矩阵的命令是 ones(4,3) ；生成4阶全零矩阵的命令是 zeros(4) 。生成对角线元素为[2 4 6]矩阵的命令是v=[2 4 6] A=diag(v,0) 。 4、假定A是一个6阶方阵，选取矩阵A第3行的指令是A(3,:) ，选取第4行第2列元素的指令是A(4,2) 。 5、生成一个从1到50，步长为10的等差数列构成的数组，可以使用的命令是 (from:step:to) 1:10:50 。 6、求x的平方根使用的命令是sqrt(x)，round(pi)的结果是（取整）3 7、命令subplot(2,2,2)的功能是把图形窗口分为2*2=4 个子图，并把第 2 个子图作为当前图形窗口。 8、matlab中，求解线性方程组时，矩阵行变换化简的命令是

reff 。 9、在用命令 p=polyfit （x ，y ，n ）对数据进行多项式拟合时，参数n 的含义为 n 次多项式。 10、求函数4y x cos7x =的5阶导数的命令是 diff （y ，5）。 11、用符号法求解微分方程2xy -3y =x ,y(1)=0,y(5)=0'''的解的指令是 dsolve(‘x*D2y-3*Dy=x^2’,’y(1)=0,y(5)=0’,’x ’) 12、Matlab 可以输入字母、汉字，但是M 文件中标点符号必须（英文）状态下输入。 14、求x e 的命令是（exp(x)），求x 的自然对数lnx 的命令（log(x )）。 15、画图时，在x 轴旁边加注文字说明的命令是 (xlabel(‘string ’) )，图名的标注命令是（title(‘string ’)），图例标注命令是（legend(‘string ’,’string ’,……)），在鼠标指定位置上标注的命令是（gtext(‘string ’)），将一个图形窗口分成多个子图的命令是（subplot(m,n,i)），画空间曲线的命令是（plot3(x,y,z)）。 16、绘制三维空间曲面图形的命令是（mesh ）和（surf ），生成格点矩阵的命令是( meshgrid ) 17、用矩形法、复合梯形公式、复合辛普生公式求12 4dx 1x +?的定积分。(详解见P150 7-28) h=0.01;x=0:h:1; y=4./(1+x.^2); format long

Matlab数学实验知识点与函数集

1.1 数学实验教学内容 1.1.1知识点（初稿）课程考核涉及函数主要为下列知识点对应的Matlab函数。知识点 Matlab函数1入门基础 1.1创建向量、矩阵（如rand，eye） 1.2常数，全局变量 1.3算术运算符 1.4关系运算符 1.5逻辑运算符 1.6数据输入、输出，输出格式 1.7绘图函数 1.7.1绘制曲线 1.7.2绘制曲面 1.7.3极坐标、参数方程 1.7.4绘图导出 1.7.5其他函数 1.8常用函数 1.9数学函数 1.10字符串操作函数 1.11文件操作函数 2控制语句 2.1分支语句 2.2循环语句 2.3其他语句、函数 3函数 3.1inline 3.2主函数 3.3子函数 4线性代数实验:,[ ], linspace, zeros, rand, randn, eye, ones, vander ans, pi, realmax, realmin, eps, inf, NaN, global +, -, *, /, .*, ./, ^, .^ <, <=, >, >=, ~= &, |, ~ load, save, format, vpa plot, plot3, ezplot, ezplot3, fplot, figure meshgrid, mesh, surf, contour polar bar, hold on, hold off, size, find, length, whos, sum, diag, class, min, max, sort, abs, input, pause, disp, cputime exp, sqrt, log, sin, cos, tan, cot, asin, acos, atan, acot, conj, real, imag, fix, floor, ceil, round, pow2, power, rem, mod, rat strcat, strvcat, str2num, num2str, sprintf fopen, fclose, fgetl, fprintf if, elseif, else, end, switch, otherwise for, while continue, break, error, warning inline function, nargin, nargout