matlab作业第10章
第10章 LTI系统的MATLAB辅助
(10.2)
3 .创建状态空间模型
状态空间模型是采用线性微分或差分方程来描述 系统的动态行为。 连续时间系统具有如下的一般形式
dx Ax Bu dt
y Cx Du
(10.3)
使用ss命令创建系统的状态空间模型的调用格式为
sys = ss(A, B, C, D)
例 10.1 在MATLAB中创建下面系统的状态空间模型：
–Feedin:sys1的输入向量,指定哪些sys1的输入与 反馈环相连 –Feedout:sys1的输出向量,指定sys1的哪些输出 端用于反馈
10.3 系统分析工具
控制系统工具箱为用户提供了一整套用于LTI模 型的时域和频域分析工具。 这些函数大都支持所有类
型的系统, 包括连续和离散系统、 SISO或MIMO系统甚
y=step(num,den,t)，y=step(sys, t) num和den为系统传递函数的分子和分母多项式系数，t为 仿真时间向量，一般可以由t=0:step:end等步长地产生。 该函数返回值y为系统输出。 [y,t,x]=step(num,den), [y,t,x]=step(sys) 此时时间向量t由系统模型的特性自动生成, 状态变量x 返回为空矩阵。 [y,x,t]=step(A,B,C,D,iu)： A,B,C,D为系统的状态空间描述，iu用来指明输入变量 的序号。x为系统返回的状态轨迹。
– 闭环特征根：roots(p) , p:闭环特征多项式 – [z, p, k]=ss2zp(a,b,c,d) – [z, p, k]=tf2zp(num,den)
已知某系统的模型如下所示：
1 2 x 4 7 y 2 2 1 6 7 2 5 2 1 0 3 0 x u 8 5 0 1 6 1 6 1x 7u
MATLAB基础教程(第四版) 第10章 MATLAB GUI设计
2019年9月27日
第4页
MATLAB R2018b 基础教程
GUI 简介（1/3）
GUI概述
MATLAB 中的 GUI 程序为事件驱动的程序。事件包括按下按钮，鼠标单击等。GUI 中的每个 控件与用户定义的语句相关。当在界面上执行某项操作时，则开始执行相关的语句。 MATLAB 提供了两种创建图形用户接口的方法：通过 GUI 向导创建的方法和编程创建 GUI 的方法。 MATLAB R2018b提供了另一种创建图形用户界面的方法：通过AppDesigner面向 对象的方式创建GUI。用户可以通过需要，选择适当的方法创建图形用户接口。通常可以参考 下面的建议。
2019年9月27日
第9页
MATLAB R2018b 基础教程
通过向导创建GUI界面（3/3）
创建菜单
1. 创建菜单栏 “快捷键” “在此菜单项上方放置分隔线” “在此菜单项前添加选中标记” “启用此项” “MenuSelectedFcn” “更多属性” 2. 创建右键菜单
2019年9月27日
第10页
MATLAB R2018b 基础教程
编写GUI代码（1/4）
GUI文件 通常情况下，一个GUI包含两个文件，一个FIG文件和一个
M文件。 M文件通常包含一个与文件同名的主函数，各个控件对应
的响应函数，这些响应函数为主函数的子函数。
2019年9月27日
第11页
MATLAB R2018b 基础教程
第13页
MATLAB R2018b 基础教程
编写GUI代码（4/4）
通过GUIDE创建GUI实例
1. 创建GUI界面 2. 编写响应函数代码
2019年9月27日
第14页
MATLAB7.0使用详解-第10章 SIMULINK高级仿真技术
10.3.1
使能子系统（Enable Subsystem）
• 使能子系统（Enable Subsystem）将控制信号分 为允许（enable）和禁止（disenable）两种。当 控制信号为正时，即控制信号为允许状态，系统 可以执行使能子系统中的模块；否则就禁止模块 功能。
10.3.2 触发子系统（Triggered Subsystem）
10.2 子系统的封装
• 使用子系统技术可以很好的优化系统模型得界面，是系统 模型的可读性更强。在对系统进行仿真时，首先要对系统 模块参数进行设置，对子系统也是如此，需要对子系统所 有模块进行合适的参数设置。 • 在第8章的例子中，子系统中模块的参数是逐一设置的，这 会给用户带来极大的不便。子系统的封装技术则是对子系 统的包装，使其成为一个真正意义上的SIMULINK模块，通 过封装，可以简化模型，用户使用一个动态设置的对话框 代替多个静态对话框，省去了对子系统内部结构的了解， 为用户的直接调用提供了方便。 • 一个已经创建的子系统的封装方法主要是通过使用下图所 示的封装编辑器来实现的，在封装编辑器中用户可以设置 封装参数设置、图标、初始以及文本。归纳起来，即设定 子系统初始值及特性、生成模块图标以及创建模块帮助文 档和描述信息。
10.4 SIMULINK的模型调试
• 如同诸多系统设置平台一样，SIMULINK也具有界 面友好、功能强大的调试功能。SIMULINK提供了 一个图形化的调试界面，简化了调试操作。 simulink调试器是一个定位和诊断模型错误的工 具．通过SIMULINK调试工具，用户可以采用多种 方式对模型进行调试，发现其中可能存在的问题 并加以修改，从而使得用户的模型设计、仿真、 分析更加快速、便捷。
ቤተ መጻሕፍቲ ባይዱ
(完整版)Matlab第十讲教案
西南科技大学本科生课程备课教案计算机技术在安全工程中的应用——Matlab入门及应用授课教师：徐中慧班级：专业：安全技术及工程第十章 MATLAB自定义函数课型：新授课教具：多媒体教学设备，matlab教学软件一、目标与要求✧通过解说与实例练习，掌握matlab创建函数M文件的方法✧掌握matlab中全局变量与局部变量的定义与用法✧通过解说与实例练习，掌握在matlab主函数M文件中创建子函数✧在实例练习过程中，回顾利用伪码编写简单程序的方法✧掌握通过创建matlab函数M文件解决生活中的计算问题二、教学重点与难点本堂课教学的重点在于引导学生掌握matlab中函数M文件的创建及应用。
本堂课的难点在于理解matlab中函数M文件主函数与子函数的区别及调用，局部变量与全局变量的定义与应用范围的区别。
三、教学方法本课程主要通过讲授法、演示法、练习法等相结合的方法来引导学生掌控本堂课的学习内容。
1)通过讲授法向学生讲述创建matlab函数M文件的基本方法、全局变量与局部变量的定义及用法等。
2)通过运用多媒体设备现场演示matlab创建函数M文件的应用实例。
3)在掌握创建matlab函数M文件基本方法的基础上，采用练习法引导学生创建函数M文件解决实际问题。
四、教学内容课后习题五（1）拉力测试装置在测试过程中，被测样本受均匀外力的作用产生形变。
下图中显示的是一组拉力测试数据。
根据以下公式计算应力与形变：00l l F A l σε-=和＝ 其中，σ是产生的应力，单位为lbf/in 2(psi);F 为施加的外力，单位为lbf;A 为样本的截面积，单位为in 2；ε为产生的形变，单位为in/in ；l 为样本的长度；0l 为样本的原始长度。
（a ）测试样本是直径为0.505in 的金属杆，根据直径可以计算出金属杆的截面积，进一步利用所提供的数据计算金属杆的应力和形变。
（b ）以形变为x 轴，应力为y 轴，作x-y 线图。
matlab章课后作业
matlab章课后作业MATLAB基础教程1~8章作业Matlab第⼀章1.阐述Matlab的功能Matlab作为⼀种⾼级计算软件，是进⾏算法开发、数据可视化、数据分析以及数值计算的交互式应⽤开发环境，已被⼴泛应⽤于不同领域。
Matlab的基本功能包括：数学计算功能、图形化显⽰功能、M 语⾔编程功能、编译功能、图形⽤户界⾯开发功能、Simulink建模仿真功能、⾃动代码⽣成功能。
Matlab第⼆章1．创建double的变量，并进⾏计算。
(1)a=87，b=190，计算a+b、a-b、a*b。
(2)创建uint8 类型的变量，数值与(1)中相同，进⾏相同的计算。
>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u) ans = 15.4189>> sqrt(u-3*v)/(u*v) ans =0 + 0.4410i 4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i >> sin(2-8*i) 1.3553e+003 +6.2026e+002i5.判断下⾯语句的运算结果。
数字图像处理与机器视觉-基于MATLAB实现 第10章 图像识别基础
• 以数据聚类的监督学习方法； • 以统计分类的无监督学习方法； • 通过对基本单元判断是否符合某种规则的结构模式识别方法； • 可同时用于监督或者非监督学习的神经网络分类法。 1.线性判用一条直线来划分已有的学 习集的数据，然后根据待测点在直线的那一边决定的分类。如下图可以做出一条直线来 划分两种数据的分类。但是一般情况下的特征数很多，想降低特征数维度。可以通过投 影的方式进行计算。然而使得一个多维度的特征数变换到一条直线上进行计算。可以减 少计算工作的复杂度。
10.2 模式识别方法
c.对称连接网络 对称连接网络有点像循环网络，但是单元之间的连接是对称的（它们在两个方向上 权重相同）。比起循环网络，对称连接网络更容易分析。这个网络中有更多的限制，因 为它们遵守能量函数定律。没有隐藏单元的对称连接网络被称为“Hopfield 网络”。有 隐藏单元的对称连接的网络被称为“玻尔兹曼机” 。 神经网络可以看成是从输入空间到输出空间的一个非线性映射，它通过调整权重和 阈值来“学习”或发现变量间的关系，实现对事物的分类。由于神经网络是一种对数据 分布无任何要求的非线性技术，它能有效解决非正态分布和非线性的评价问题，因而受 到广泛的应用。由于神经网络具有信息的分布存储，并行处理及自学习能力等特点，它 在泛化处理能力上显示出较高的优势。
模式识别是指对表征事物或现象的各种形式的(数值的、文字的和逻辑关系的)信息进行 处理和分析，以对事物或现象进行描述、辨认、分类和解释的过程，是信息科学和人工智 能的重要组成部分。
基于监督学习的模式识别系统由4大部分组成，即待识别对象、预处理、特征提取和分 类识别，如图10-1所示。
图10-1 模式识别流程图
matlab课后习题答案(1-9章)要点
1 数字1.5e2，1.5e3 中的哪个与1500相同吗？1.5e32 请指出如下5个变量名中，哪些是合法的？abcd-2xyz_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))^2w1 = -2.0000 + 3.4641i ;w2 = 4.0000 ;w3 =-2.0000 + 3.4641i 5 指令clear, clf, clc 各有什么用处？clear 清除工作空间中所有的变量。
clf 清除当前图形。
clc 清除命令窗口中所有显示。
第二章1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？3/7+0.1双; sym(3/7+0.1)符; sym('3/7+0.1') 符;; vpa(sym(3/7+0.1)) 符;2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') symvar(sym('sin(w*t)'),1) w a z3 （1）试写出求三阶方程05.443=-x 正实根的程序。
注意：只要正实根，不要出现其他根。
（2）试求二阶方程022=+-a ax x 在0>a 时的根。
（1）reset(symengine)syms x positive solve(x^3-44.5) ans =(2^(2/3)*89^(1/3))/2（2）求五阶方程022=+-a ax x 的实根 syms a positive %注意：关于x 的假设没有去除 solve(x^2-a*x+a^2)Warning: Explicit solution could not be found. > In solve at 83 ans =[ empty sym ]syms x clear syms a positivesolve(x^2-a*x+a^2) ans =a/2 + (3^(1/2)*a*i)/2 a/2 - (3^(1/2)*a*i)/24 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。
第10章MATLAB特征值与特征向量的计算
第10章MATLAB特征值与特征向量的计算MATLAB是一种非常强大的数值分析和科学计算工具，可以用于计算矩阵的特征值和特征向量。
特征值和特征向量是矩阵分析中非常重要的概念，具有广泛的应用。
在MATLAB中，计算矩阵的特征值和特征向量有多种方法。
下面将介绍其中两种常用的方法。
一种方法是使用eig函数。
该函数可以计算矩阵的特征值和特征向量。
用法如下：```MATLAB[A, D] = eig(B)```其中，B是要计算特征值和特征向量的矩阵，A是一个由特征向量组成的矩阵，D是一个由特征值组成的对角矩阵。
特征值和特征向量的对应关系是D的对角线上的值与A的列是一一对应的。
另一种方法是使用eigs函数。
该函数可以计算矩阵的部分特征值和特征向量。
用法如下：```MATLAB[V, E] = eigs(A, k)```其中，A是要计算特征值和特征向量的矩阵，k是要计算的特征值和特征向量的数量。
这两种方法的计算结果是一样的，只是使用的方法略有不同。
eig函数计算所有的特征值和特征向量，而eigs函数只计算部分特征值和特征向量。
计算矩阵的特征值和特征向量在很多领域中都有重要的应用。
例如，在图像处理中，特征值和特征向量可以用于图像压缩和识别；在控制系统中，特征值和特征向量可以用于分析系统的稳定性和性能等。
虽然MATLAB提供了方便的函数来计算特征值和特征向量，但在实际应用中，需要注意以下几点：1.确保输入的矩阵是方阵，否则将无法计算特征值和特征向量。
2.需要注意计算所得的特征值和特征向量的精度。
MATLAB中的计算结果可能存在误差，因此在应用中需要进行适当的误差处理。
3.特征值和特征向量的计算量较大，如果矩阵的维度较大，可能会导致计算时间过长。
对于大规模问题，可以考虑使用近似方法来计算特征值和特征向量。
总之，MATLAB提供了方便的函数来计算矩阵的特征值和特征向量，这对于很多领域的数值分析和科学计算都非常有用。
在使用这些函数时，需要注意输入矩阵的维度和计算结果的精度，并选择适当的方法来计算特征值和特征向量。
matlab10-7
电子信息工程学院
25
第4章
求解非线性方程f(x)=0
电子信息工程学院
26
第4章
求解非线性方程f(x)=0
二分法简单方便，但收敛速度慢；
迭代法虽然收敛速度稍微快点，但需要判断 能否收敛； 只要初值选取得当，切线法具有恒收敛且收 敛速度快的优点，但需要求出函数的导数；
弦截法不需要求导数，特别是我们介绍的快 速弦截法，收敛速度很快，但是需要知道两 个近似的初始根值才能作出弦，要求的初始 条件较多。
电子信息工程学院 30
第4章
求解非线性方程f(x)=0
4.4 求解非线性方程组数值解的迭代法 对于非线性方程组（以二元方程组为 例，其它可以类推）
电子信息工程学院
，
f (x) ，如此继续。 0
19
第4章
求解非线性方程f(x)=0
得到牛顿迭代公式：
x k 1 f ( xk ) xk f ( x k )
由此得出迭代序列： x1 , x2 , xk ,,
条件下收敛于方程的根x*
在一定
matlab课后答案完整版
ones表示1矩阵zeros表示0矩阵ones(4)表示4x4的1矩阵zeros(4)表示4x4的0矩阵zeros(4,5)表示4x5的矩阵eye(10,10)表示10x10的单位矩阵rand(4,5)表示4x5的伴随矩阵det(a)表示计算a的行列式inv(a)表示计算a的逆矩阵Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩[v,d]=eig(a)对角矩阵b=a’表示求a矩阵的转置矩阵sqrt表示求平方根exp表示自然指数函数log自然对数函数abs绝对值第一章一、5（1）b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12];>> d=[65 5;-23 -59;54 7];>> e=b*ce =5271 11574-11336 6641978 3112（2）a=50:1:100二、1 、x=-74;y=-27;z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z =2、a=-3.0:0.1:3.0;>> b=exp(-0.3*a).*sin(a+0.3)3、x=[2 4;-0.45 5];y=log(x+sqrt(1+x.^2))/2y = 4、a*b表示a矩阵和b矩阵相乘a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列A(m,:)表示取a矩阵第m行的全部元素A(:,n)表示取a矩阵的第n列全部元素A./B表示a矩阵除以b矩阵的对应元素，B.\A等价于A./BA.^B表示两个矩阵对应元素进行乘方运算A.^2表示a中的每个元素的平方A^2表示A*A例：x=[1,2,3];y=[4,5,6];z=x.^yz=1 32 729指数可以是标量（如y=2）.底数也可以是标量（如x=2）5、a=1+2i;>> b=3+4i;>> c=exp((pi*i)/6)c =d=c+a*b/(a+b)d =第二章二、4、（1）y=0;k=0;>> while y<3k=k+1;y=y+1/(2*k-1);end>> display([k-1,y-1/(2*k-1)])ans =第三章二1（1）x=0:pi/10:2*pi;>> y=x-x.^3/6;>> plot(x,y)（2）x=0:pi/10:2*pi; y=(exp(-x.^2/2))/2*pi;plot(x,y)（3）x=-8:0.01:8;y=sqrt((64-x.^2)/2);plot(x,y)（4）t=0:0.1:8*pi; >> x=t.*sin(t);>> y=t.*cos(t); >> plot(x,y)x=0:pi/100:2*pi;y1=exp(-0.5*x);y2=exp(-0.5*x).*sin(2*x);plot(x,y1,x,y2)>> title('x from 0 to 2{\pi} ');>> xlabel('variable x');>> ylabel('variable y');>> text(1.5,0.5,'曲线y1=e^(-0.5x)');>> text(3,0.1,'曲线y2=cos(4{\pi}x)e^{-0.5x}');>> legend('y1','y2')variable xvariabley2、(1)y1=2*x-0.5;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);>> k=find(abs(y-x)<1e-2);>> t1=t(k)t1 =>> z=sin(3.*(t1)).*cos(t1) z =>> plot(t,x,t,y,'k:',t1,z,'bp');（2）subplot(1,2,1);>> scatter(x1,y1,10)；>> title('y=2x-0.5')；>> subplot(1,2,2);>> scatter(x,y,10)3、subplot(1,2,1);x=0:0.01:pi;y=sin(1./x);plot(x,y)subplot(1,2,2);fplot('sin(1./x)',[1,100])4、t=0:pi:2*pi;y=1./(1+exp(-t));subplot(2,2,1);%图形窗口的分割bar(t,'group'); %绘制柱形图（分组）subplot(2,2,2);barh(t,'stack'); %绘制柱形图（堆积）subplot(2,2,3);loglog(t,y); %函数使用全对数坐标，x,y均采用常用对数刻度subplot(2,2,4);semilogy(t,y); %函数使用半对数坐标，y轴为常用对数刻度，x轴仍为线性刻度101010101010105、（1）theta=linspace(-pi,pi,100);ro=5.*cos(theta)+4;polar(theta,ro);（2）x=linspace(0,2*pi,100);a=1>> r=a.*(1+cos(x));polar(x,r);2706、（1）t=0:pi/10:2*pi;>> x=exp((-t)/20).*cos(t);>> y=exp((-t)/20).*sin(t);>> z=t;>> plot3(x,y,z);（2）t=0:0.01:1;x=t;>> y=t.^2;>> z=t.^3;>> plot3(x,y,z);7、x=-30:0.1:0;>> y=0:0.1:30;>> [x,y]=meshgrid(x,y);>>z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2 );>> meshc(x,y,z);绘制曲面图和等高线8、x=linspace(-3,3,100);>> y=linspace(-3,3,100);>> [x y]=meshgrid(x,y); %可以将向量转化为矩阵>> fxy=-5./(1+x.^2+y.^2);>> i=find(abs(x)<=0.8 & abs(y)<=0.5); >> fxy(i)=NaN;>>surf(x,y,fxy) %绘制三维曲面图9、u=linspace(1,10,100);v=linspace(-pi,pi,100);[u v]=meshgrid(u,v); x=3.*u.*sin(v); y=2.*u.*cos(v); z=4*u.^2; surf(x,y,z); shading interp;第五章二1、a=rand(1,30000);mean(a) %求平均数 ans = >>b=std(a) %求标准差 b =>> c=max(a) c =>> d=min(a) d =size(find(a>0.5))/size(a)ans =2、h=[466,715,950,1422,1635]; >> w=[7.04,4.28,3.40,2.52,2.13]; >> hh=[500,900,1500];>> ww=interp1(h,w,hh,'spline')ww =3、x=linspace(1,10,50); y=log(x);f=polyfit(x,y,5); %求曲线的拟合 >> yy=polyval(f,x); >> plot(x,y,'r-',x,yy,'g.')5、（1）、（2） p1=[1,2,0,7]; p2=[1,-2]; p3=[1,0,5,1]; p12=conv(p1,p2); >>p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3]; >> roots(p) ans =-3.46561.2400（3） a=[-1,4,3;2,1,5;0,5,6]; >>y1=polyval(p,a) %以矩阵a 的每一个元素为自变量 y1 =-29 291 95 19 -3 697 -13 697 1427 >>y2=polyvalm(p,a) %以矩阵a 为自变量 y2 =391 2084 3273502 2693 4207720 3775 5892 6、（1）z=fzero('3*x-sin(x)+1',0) %求x=0时附近的根z =第八章二、2t=0:pi/20:2*pi;x=sin(t);y=cos(t);x1=sin(7*t);y1=cos(7*t);h=plot(x,y,x1,y1);set(h,'marker','x','linewidth',2);set(gca,'xtick',-1:0.1:1);title('篮筐')3、x=0:pi/10:5*pi;y=exp(-0.2*x).*cos(x)+2;h=plot(x,y);set(gca,'color','red','linestyle','-','linewidth',3) ;text(5,2.4,'y=exp(-0.2*x).*cos(x)+2'); 4、t=-pi:pi/100:pi;x=cos(t);y=sin(t);z=t;h=plot(t,x,t,y,t,z);set(h,'linestyle','-','linewidth',3);字符串例ch='Welcome to Beijing';subch=ch(12:18) 选12~18个字符串（空格也算）subch =Beijing>> k=find(ch>='a'&ch<='z'); 找到所有的小写字母的位置>> ch(k)=ch(k)-('a'-'A'); 将小写字母变成大写字母>> char(ch)ans =WELCOME TO BEIJING>> length(k) 统计小写字母的个数ans =14例：已知y=1-1/2+1/3-1/4.........-1/100求y 的值y=0;>> n=100;>> for i=1:100;y=y+(-1)^(i-1)/i;end>> disp(y)绘制二维曲线图x=0:pi/100:2*pi;>> y1=0.2*exp(-0.5*x).*cos(4*pi*x);>> y2=1.5*exp(-0.5*x).*cos(pi*x);>> plotyy(x,y1,x,y2);绘制三维图像例：x=sint+tcosty=cost-tsintz=tt=0:pi/10:10*pi;x=sin(t)+t.*cos(t);y=cos(t)-t.*sin(t);z=t;plot3(x,y,z);axis([-30 30 -30 30 0 35]); 坐标轴的最大值与最小值title('line in 3-D space'); 图形的题目>> xlabel('x');ylabel('y');zlabel('z'); 标注坐标>> grid on; 加网格线xz三维例]2/,0[],,0[,cossin22ππ∈∈+=yxyxz[x,y]=meshgrid(0:pi/100:pi,0:pi/100:pi/ 2);>> z=sin(x.^2)+cos(y.^2);>> mesh(x,y,z);>> axis([0 4 0 1.8 -1.5 1.5]);t=0:pi/20:2*pi;subplot(1,2,1);[x,y,z]=cylinder(sin(t),30);surf(x,y,z); 绘制三维曲面图subplot(1,2,2);>> [x,y,z]=peaks(100);>> mesh(x,y,z); 绘制三维网格图多项式求导例：f(x)=1/x^2+5p=[1];>> q=[1,0,5];>> [p,q]=polyder(p,q)注：c=conv(a,b) 表示a多项式与b 多项式乘积[p,r]=deconv(a,b) 表示a多项式与b 多项式相除其中p为商向量r为余数向量p=polyder(p) 表示求p的导数p=poleder(p,q) 表示求p乘以q的导数[p,q]=poleder(p,q) 表示p除以q的导数多项式求根例：f(x)=2x^4-12x^3+3x^2+5p=[2,-12,3,0,5];>> x=roots(p);>> p=[2,-12,3,0,5];x=roots(p) 求方程f(x)=0的根x =5.72460.8997>> g=poly(x) 已知多项式的根求多项式g =符号求导x=a(t-tsint)y=b(1-cost) 求y对x的一阶导数syms x y a b t;>> f21=a*(t-sin(t));>> f22=b*(1-cos(t));>> diff(f22)/diff(f21) 求y对x的一阶导数ans =b*sin(t)/a/(1-cos(t))注：diff(f1,x,2) 表示f1对x的二阶导数diff（f3,x）表示z对x的偏导diff（f3,y）表示z对y的偏导求不定积分int(f) 求f的不定积分f1=int(f,a,b) 求f在a，b之间的定积分eval（f1）计算积分值符号求极限syms x h>> f=(sin(x+h)-sin(x))/h;>> limit(f,h,0) h趋向于0ans =cos(x) 例2f=sym('(1+t/x)^x');limit(f,inf) f 趋向于无穷 ans = exp(t) 例3f=sym('x*(sqrt(x^2+1)-x)');limit(f,sym('x'),inf,'left') x 趋向于正无穷ans = 1/2 大小写ch='Welcome to Beijing'; subch=ch(12:18)subch =Beijing>> k=find(ch>='A'&ch<='Z'); ch(k)=ch(k)-('A'-'a'); >> char(ch)ans =welcome to beijing>> length(k)ans =2课堂习题2一、求方程组⎪⎩⎪⎨⎧=-+=++=-+234326454510243z y x z y x z y x 的解，用三种方法。
MATLAB-第六次实验
x、y、z 都是矩阵时，plot3 命令的使用
[X,Y]=meshgrid(-pi:0.1:pi); Z=sin(X)+cos(Y); plot3(X,Y,Z)
MATLAB基础
5.4 三维图形绘制与复数的图形绘制
g.绘制三维圆柱曲面命令cylinder
cylinder命令中，柱面的轴线定义为 z 轴，只要给出母线的描述就可完成一个柱面。 调用格式为： [X,Y,Z] = cylinder(R,N)； [X,Y,Z] = cylinder(R)：缺省值 N=20； [X,Y,Z] = cylinder：缺省值 N=20，R=[1，1]。 R：是一描述柱面母线的向量； N：是旋转柱面上的分割线条数；
MATLAB基础
已知系统传递函数
n Y(S) 100 G(S) 2 2 X(S) S 10 S 100 S 2 n S n 2
2
，其中：
0.5 n 10rad / s
在输入信号x(t)=1（即单位阶跃信号），要求建立M文件，画出输出 信号y(t)（即单位阶跃响应）的图象， MATLAB程序如下： G=tf([100],[1,10,100]); %得到传递函数G(S) t=0:0.05:3; %时间t设在0—3秒间，运算间隔0.05秒 y=step(G,t); %得到单位阶跃响应y(t) plot(t,y,'-k'); %绘制出y(t) 图象 xlabel('t'); %横坐标轴标号为t ylabel('y(t)'); %纵坐标轴标号为y(t) grid %绘制出栅格以便于读取数值
MATLAB基础
计算积分
7
int(f,v,a,b): 计算定积分 int(f,a,b): 计算关于默认变量的定积分 int(f,v): 计算不定积分 int(f): 计算关于默认变量的不定积分
第10章 MATLAB图形句柄
Root 根屏幕
Figure 图形窗
Uicontrol 界面控件
Axes 轴
Uimenu 界面菜单
Line 线
Surface 面
Rectan gle 方
Patch 块
Image 图像
Text 文字
Light 光
句柄图形体系的对象树结构
2．图形对象句柄 ． MATLAB在创建每一个图形对象时，都为该 在创建每一个图形对象时， 在创建每一个图形对象时 对象分配唯一的一个值， 对象分配唯一的一个值，称其为图形对象 句柄(Handle)。句柄是图形对象的唯一标识 句柄 。 不同对象的句柄不可能重复和混淆。 符，不同对象的句柄不可能重复和混淆。 计算机屏幕作为根对象由系统自动建立， 计算机屏幕作为根对象由系统自动建立，其 句柄值为0， 句柄值为 ，而图形窗口对象的句柄值为一 正整数，并显示在该窗口的标题栏， 正整数，并显示在该窗口的标题栏，其他 图形对象的句柄为浮点数 浮点数。 图形对象的句柄为浮点数。MATLAB提供 提供 了若干个函数用于获取已有图形对象的句 柄。
get(1) Alphamap = [ (1 by 64) double array] BackingStore = on CloseRequestFcn = closereq Color = [0.8 0.8 0.8] Colormap = [ (64 by 3) double array] CurrentAxes = [] CurrentCharacter = s CurrentObject = [] CurrentPoint = [0 0] Dithermap = [ (64 by 3) double array] DithermapMode = manual DoubleBuffer = off FileName = FixedColors = [ (3 by 3) double array] IntegerHandle = on InvertHardcopy = on KeyPressFcn = MenuBar = figure MinColormap = [64] Name =
MATLAB习题答案(清华大学)
高等应用数学问题MATLAB 求解习题参考解答（薛定宇著）目录第1 章计算机数学语言概述2第2 章MATLAB 语言程序设计基础5第3 章微积分问题的计算机求解17第4 章线性代数问题的计算机求解29第5 章积分变换与复变函数问题的计算机求解43第6 章代数方程与最优化问题的计算机求解53第7 章微分方程问题的计算机求解71第8 章数据插值、函数逼近问题的计算机求解93第9 章概率论与数理统计问题的计算机求解114第10 章数学问题的非传统解法127第A章自由数学语言Scilab 简介136第1 章计算机数学语言概述1 在你的机器上安装MATLAB 语言环境，并键入demo 命令，由给出的菜单系统和对话框原型演示程序，领略MATLAB 语言在求解数学问题方面的能力与方法。
【求解】在MATLAB 提示符>> 下键入demo 命令，则将打开如图1-1 所示的窗口，窗口左侧的列表框可以选择各种不同组合的演示内容。
图1-1 MATLAB 演示程序界面例如，用户选择MATLAB ! Graphics ! Volume Vlsulization 演示，则将得出如图1-2 所示的演示说明，单击其中的Run this demo 栏目，则将得出如图1-3 所示的演示界面。
用户可以在该界面下按按钮，逐步演示相关内容，而实现这样演示的语句将在该程序界面的下部窗口中给出。
2 作者用MATLAB 语言编写了给出例子的源程序，读者可以自己用type 语句阅读一下源程序，对照数学问题初步理解语句的含义，编写的源程序说明由下表列出。
第1 章计算机数学语言概述3图1-2 MATLAB 演示程序界面举例序号文件名程序说明例1.1 c1ex1.m 利用MATLAB 的符号运算工具箱求解微分问题例1.2 c1ex2.m 分别利用MATLAB 的符号运算工具箱和数值运算功能求解多项式方程，其中用数值方法得出的结果有误差例1.3 c1ex3.m 分别利用MATLAB 的符号运算工具箱和数值运算功能计算Hilbert 矩阵的行列式，其中用数值方法得出的结果有很大误差例1.4 c1ex4.m 令x1 = y; x2 = y_，则可以将原来的二阶微分方程转换成一阶微分方程组，然后就可以求解微分方程的数值解了，原方程是非线性微分方程，故不存在解析解。
matlab习题第十章
第十章实验指导1. screen=get(0,'screensize');W=screen(3);>> H=screen(4);>> figure('position',[0.2*H,0.2*H,0.5*W,0.3*H],'Name','设计菜单','numbertitle','off','menubar','none');aa=['new item'];bb=[‘open item’];cc=[‘save item’];dd=[‘help item’];hfile=uimenu(gcf,'label','&File');)uimenu(hfile,'label','new','call',msgbox(aa,'Title'),'set(save,''enab le'',''off'')');uimenu(hfile,’label’,’open’,’call’,msgbox(bb,’Title’),set(sav e,’’enable’’,’’off’’)’);uimenu(hfile,’lable’,’save’,’call’,msgbox(cc,’Title’)); uimenu(gcf,’lable’,’&Exit’,’call’,close(gcf));hhelp=uimenu(gcf,’label’,’&Help’);uimenu(hhelp,’label’,’about,’’call’,msgbox(dd,’Title’),’set( save,’’enable’’,’’on’’))2、pdstart=uicontrol(gcf,'style','push','position',...[50,5,60,25],'string','start plot',...'callback','t=-pi:pi/20:pi;plot(t,t.*t)');ptgrid=uicontrol(gcf,'style','toggle','position',...[150,5,60,25],'string','grid','callback','grid');htxt=uicontrol(gcf,'style','text','string',...'color options','position',[200,130,150,20]);hr=uicontrol(gcf,'style','radio','string',...'red','position',[200,100,150,25],'value',1,...'callback',['set(hr,''value'',1);','set(hb,''value'',0);',...'set(hy,''value'',0);','set(plot(t,t.*t),''color'',''R'')']);hb=uicontrol(gcf,'style','radio','string',...'blue','position',[200,75,150,25],'value',1,...'callback',['set(hb,''value'',1);','set(hr,''value'',0);',...'set(hy,''value'',0);','set(plot(t,t.*t),''color'',''B'')']);hy=uicontrol(gcf,'style','radio','string',...'yellow','position',[200,50,150,25],'value',1,...'callback',['set(hy,''value'',1);','set(hr,''value'',0);',...'set(hb,''value'',0);','set(plot(t,t.*t),''color'',''Y'')']);3、function varargout = duoxiangshiqiugen(varargin)% DUOXIANGSHIQIUGEN MATLAB code for duoxiangshiqiugen.fig% DUOXIANGSHIQIUGEN, by itself, creates a new DUOXIANGSHIQIUGEN or raises the existing% singleton*.%% H = DUOXIANGSHIQIUGEN returns the handle to a new DUOXIANGSHIQIUGEN or the handle to% the existing singleton*.%% DUOXIANGSHIQIUGEN('CALLBACK',hObject,eventData,handles,...) calls the local% function named CALLBACK in DUOXIANGSHIQIUGEN.M with the given input arguments.%% DUOXIANGSHIQIUGEN('Property','Value',...) creates a new DUOXIANGSHIQIUGEN or raises the% existing singleton*. Starting from the left, property value pairs are % applied to the GUI before duoxiangshiqiugen_OpeningFcn gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to duoxiangshiqiugen_OpeningFcn via varargin. %% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one% instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help duoxiangshiqiugen% Last Modified by GUIDE v2.5 28-May-2014 09:58:02% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @duoxiangshiqiugen_OpeningFcn, ...'gui_OutputFcn', @duoxiangshiqiugen_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before duoxiangshiqiugen is made visible.function duoxiangshiqiugen_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to duoxiangshiqiugen (see VARARGIN)% Choose default command line output for duoxiangshiqiugenhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes duoxiangshiqiugen wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.function varargout = duoxiangshiqiugen_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;function edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit1 as text% str2double(get(hObject,'String')) returns contents of edit1 as a double% --- Executes during object creation, after setting all properties.function edit1_CreateFcn(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');end% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)com1=get(handles.edit1,'string');A=str2mat(com1);a=str2num(A)X=roots(a)set(handles.edit2,'string',num2str(X));function edit2_Callback(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text% str2double(get(hObject,'String')) returns contents of edit2 as a double% --- Executes during object creation, after setting all properties.function edit2_CreateFcn(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction varargout = jizuobiaoquxian(varargin)% JIZUOBIAOQUXIAN MATLAB code for jizuobiaoquxian.fig% JIZUOBIAOQUXIAN, by itself, creates a new JIZUOBIAOQUXIAN or raises the existing% singleton*.%% H = JIZUOBIAOQUXIAN returns the handle to a new JIZUOBIAOQUXIAN or the handle to% the existing singleton*.%% JIZUOBIAOQUXIAN('CALLBACK',hObject,eventData,handles,...) calls the local% function named CALLBACK in JIZUOBIAOQUXIAN.M with the given input arguments.%% JIZUOBIAOQUXIAN('Property','Value',...) creates a new JIZUOBIAOQUXIAN or raises the% existing singleton*. Starting from the left, property value pairs are % applied to the GUI before jizuobiaoquxian_OpeningFcn gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to jizuobiaoquxian_OpeningFcn via varargin. %% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one % instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help jizuobiaoquxian% Last Modified by GUIDE v2.5 28-May-2014 09:55:56% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @jizuobiaoquxian_OpeningFcn, ...'gui_OutputFcn', @jizuobiaoquxian_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before jizuobiaoquxian is made visible.function jizuobiaoquxian_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to jizuobiaoquxian (see VARARGIN)% Choose default command line output for jizuobiaoquxianhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes jizuobiaoquxian wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.function varargout = jizuobiaoquxian_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)com1=get(handles.edit1,'string');com2=get(handles.edit2,'string');com3=get(handles.edit3,'string');A=str2mat(com1);B=str2mat(com2);C=str2mat(com3);a=str2num(A);b=str2num(B);n=str2num(C);th=0:0.01:pi;polar(th,a*cos(b+th*n))axis([-5,5,-20,20]);hold onfunction edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit1 as text% str2double(get(hObject,'String')) returns contents of edit1 as a double% --- Executes during object creation, after setting all properties. function edit1_CreateFcn(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit2_Callback(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text% str2double(get(hObject,'String')) returns contents of edit2 as a double% --- Executes during object creation, after setting all properties. function edit2_CreateFcn(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit3_Callback(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit3 as text% str2double(get(hObject,'String')) returns contents of edit3 as a double% --- Executes during object creation, after setting all properties. function edit3_CreateFcn(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');end5、思考练习1.图形用户界面，就是把所需要的M文件，编程一个可以和用户交互的界面程序，这样看起来更直观的一种呈现方式。
matlab数学实验答案第十章
1.>> f=[-3,4,-2,5];A=[1 1 3 -1;2 -3 1 -2];>> b=[14;-2];>> Aeq=[4 -1 2 -1];Beq=-2;lb=zeros(3,1);>> [x,feval]=linprog(f,A,b,Aeq,Beq,lb) Optimization terminated.x =0.00008.00000.0000-6.0000feval =2.00002.>> f=[5;4;8];A=[2 -1 0;5 3 0];b=[4;15];>> Aeq=[1 2 1];Beq=6;lb=zeros(3,1);>> [x,feval]=linprog(f,A,b,Aeq,Beq,lb) Optimization terminated.x =0.00003.00000.0000feval =12.00003.M文件>> A=[1 1 0 0 0 0 0 0 0 0 0;1 -0.06 1 1 0 0 0 0 0 0 0 ;-0.15 -0.06 1 1 1 1 1 1 0 0 0 ;-0.15 -0.06 -0.15 1 -0.06 1 1 -0.06 1 1 0 ;-0.15 -0.06 -0.15 1 -0.06 1 1 -0.06 1 -0.06 1;0 0 0 1 -0.06 -0.15 1 0 0 0 0;0 0 0 1 0 0 1 0 0 0 0];>> b=[10;10;10;10;10;3;4];>> f=[0;0;0;-1.40;0;0;-1.25;0;-1.14;0;-1.06];>> lb=[0;0;0;0;0;0;0;0;0;0;0];>> ub=[10;10;10;10;10;10;10;10;10;10;10];>> [x,fval,exitflag,output,lambda]=linprog(f,A,b,[],[],lb,ub) Optimization terminated.x =6.73163.26840.00003.46457.74140.00000.00000.00008.20580.00000.0000fval =-14.2049五年末拥有的资金金额总额最大为 14.20494.>> f=[-30;-40];A=[20 100;40 80;90 60];>> b=[200;200;360];>> [x,feval]=linprog(f,A,b,[],[],[0;0])Optimization terminated.x =3.50000.7500feval =-135.00005.>> f=[-7000;-12000];A=[9 5;4 5;3 10];b=[360;200;300]; >> [x,feval]=linprog(f,A,b)Optimization terminated.x =20.000024.0000feval =-4.2800e+0056.>> f=[10,5,6,4,8,12];>> A=[-1 -1 -1 0 0 0;0 0 0 -1 -1 -1; -1 0 0 -1 0 0 ; 0 -1 0 0 -1 0 ;0 0 -1 0 0 -1]; >> b=[-60;-100;-50;-70;-40];>> Aeq=[];>> beq=[]; vlb=zeros(6,1);vub=[];>> [x,fval]=linprog(f,A,b,Aeq,beq,vlb,vub)Optimization terminated.x =0.000020.000040.000050.000050.00000.0000fval =940.00007.>> f=[8400;7500;6700];>> A=[60 65 80;-3800 -3500 -2500];>> b=[70 ;-2900];>> Aeq=[1 1 1];>> beq=1;>> lb=[0 0 0];>> ub=[inf inf inf];>> [x,fval]=linprog(f,A,b,Aeq,beq,lb,ub)Optimization terminated.x =0.00000.66670.3333fval =7.2333e+0038.>> f=[-1;-2;1;-1;-4;2];>> A=[1 1 1 1 1 1;2 1 -2 1 0 0;0 0 1 1 2 1];>> b=[6;4;4];>> lb=zeros(6,1);>> [x,fval]=linprog(f,A,b,[],[],lb)Optimization terminated.x =0.00004.00000.00000.00002.00000.0000fval =-16.000010.建立M函数fun44.mfunction f=fun44(x)f=-(sqrt(x(1))+sqrt(x(2))+sqrt(x(3))+sqrt(x(4))); 建立M函数myconl.mfunction [g,ceq]=mycon1(x)g(1)=x(1)-400;g(2)=1.1*x(1)+x(2)-440;g(3)=1.21*x(1)+1.1*x(2)+x(3)-484;g(4)=1.331*x(1)+1.21*x(2)+1.1*x(3)+x(4)-532.4;ceq=0command窗口>> x0=[1;1;1;1];>> vlb=[0;0;0;0];vub=[];>> A=[];b=[];Aeq=[];beq=[];>>[x,fval]=fmincon('fun44',x0,A,b,Aeq,beq,vlb,vub,'mycon1') 12.function f=fmb(x)f=-x(1)*x(2)*x(3);>> A=[1 -2 -2;1 2 2];b=[0;72];>> lb=[-inf;10;-inf];ub=[inf;20;inf];>> Aeq=[1 -1 0];Beq=10;>> [x,feval]=fmincon(@fmb,[20,10,0],A,b,Aeq,Beq,lb,ub,[])x =22.5850 12.5850 12.1225feval =-3.4456e+003。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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);
z
（4）3y
>> x=0:1:50;
>> y=0:1:50;
>> z=y.^3;
>> hf=plot3(x,y,z);。