matlab第6章MATLAB图形绘制
MATLAB基础教程(第四版) 第6章 MATLAB 编程

命令
说明
disp(A)
显p(‘text’) Format
显示单引号内部的文本串 控制屏幕输出的显示格式
Fprintf x = input('text')
执行格式化的写入到屏幕或者写入到一个文件 显示单引号中的文本,等待用户的键盘输入,并且将输入值存储在x中
x = input('text','s') k=menu('title','option1','option2',...
有变量,就需要考虑使用函数文件。这将避免用变量名“弄乱”工作空间,并将 同时减少内存需求。
(6) 在不使用文本编辑器打开M-文件的情况下,用户可以使用type命令来查看文件的 内容。
2019年9月27日
第7页
MATLAB R2018b 基础教程
脚本文件和编辑器(4/11)
有效使用命令窗口和编辑器/调试器窗口
脚本文件和编辑器(8/11)
使用脚本文件存储数据 可能会有一些应用程序要求用户频繁地访问同一组数据。
如果是这样的话,用户可以将数组中的数据存储在一个脚 本文件之中。
2019年9月27日
第12页
MATLAB R2018b 基础教程
脚本文件和编辑器(9/11)
控制输入和输出
MATLAB为从用户获取输入以及为格式化输出(执行MATLAB命令而获得的结果)提供了许多有用的命令。
2019年9月27日
第5页
MATLAB R2018b 基础教程
脚本文件和编辑器(2/11)
创建和使用脚本文件
2019年9月27日
第6页
MATLAB R2018b 基础教程
数字图像处理及应用(MATLAB)第6章习题答案

第六章习题与思考题参考答案1. 简述飞机红外图像识别中用到的五个红外特征量各自的作用?1)长宽比:反应了目标的几何形状;2)复杂度:反应了红外目标轮廓的情况;3)紧凑度:反应了红外目标在其所在最小外接矩形中所占比重;4)均值对比度:反映了目标的物理特性与背景的物理特性之间的关系;5)部分最亮像素点数与目标总像素数的比值:反应了目标的明暗变化情况;2. 可视密码共享中,如果实现(4 , 5)门限的可视密码分享,程序将做如何编写?clear allclose allM=imread('0.jpg');ss=rgb2gray(M);figureimshow(ss);[m n]=size(ss);for i=1:m*nif ss(i)>250ss(i)=250;endends=double(ss)+1;x=[1 2 3 4 5];g1=zeros(m,n);g2=zeros(m,n);g3=zeros(m,n);g4=zeros(m,n);g5=zeros(m,n);yy1=zeros(m,n);yy2=zeros(m,n);yy3=zeros(m,n);yy4=zeros(m,n);yy5=zeros(m,n);y1=zeros(m,n);y2=zeros(m,n);y3=zeros(m,n);y4=zeros(m,n);y5=zeros(m,n);for j=1:m*na1=mod(2*j,251);a2=mod(3*j,251);a3=mod(5*j,251);f=[a1 a2 a3 s(j)];g1(j)=polyval(f,x(1));yy1(j)=mod(g1(j),251);g2(j)=polyval(f,x(2));yy2(j)=mod(g2(j),251);g3(j)=polyval(f,x(3));yy3(j)=mod(g3(j),251);g4(j)=polyval(f,x(4));yy4(j)=mod(g4(j),251);g5(j)=polyval(f,x(5));yy5(j)=mod(g5(j),251);endy1=uint8(yy1-1)y2=uint8(yy2-1);y3=uint8(yy3-1);y4=uint8(yy3-1);y5=uint8(yy3-1);figure,imshow(y1);figure,imshow(y2)figure,imshow(y3);figure,imshow(y4);figure,imshow(y5);l1=(x(2)*x(3)*x(4)*x(5))/[(x(1)-x(2))*(x(1)-x(3))*(x(1)-x(4))*(x(1)-x(5))];l2=(x(1)*x(3)*x(4)*x(5))/[(x(2)-x(1))*(x(2)-x(3))*(x(2)-x(4))*(x(2)-x(5))];l3=(x(1)*x(2)*x(4)*x(5))/[(x(3)-x(1))*(x(3)-x(2))*(x(3)-x(4))*(x(3)-x(5))];l4=(x(1)*x(2)*x(3)*x(5))/[(x(4)-x(1))*(x(4)-x(2))*(x(4)-x(3))*(x(4)-x(5))];l5=(x(1)*x(2)*x(3)*x(4))/[(x(5)-x(1))*(x(5)-x(2))*(x(5)-x(3))*(x(5)-x(4))];rr1=zeros(m,n);r=zeros(m,n);for j=1:m*nrr1(j)=mod(yy1(j)*l1+yy2(j)*l2+yy3(j)*l3+yy4(j)*l4+yy5(j)*l5,251);endr=uint8(rr1-1);figure,imshow(r);3. 已知图像⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=805020016821024015010017018013070901262068M 。
第06章_MATLAB数值计算_参考解答

%可参见教材第 157 页例 6.19
运行结果: dx =
000
%当 x=pi/2 时单独计算 x=pi/2; f=inline('sin(x).^2+cos(x).^2'); dx=diff(f([x, pi]))/(pi/2)
dx = 0
(2) 程序设计: clear all; close all; clc; x=1:3; f=inline('sqrt(x.^2+1)'); dx=diff(f([x, 4]))
运行结果:
U= 1.0e-004 * -0.0675
0.1715
fmin = 1.9920e-010
(2) 程序设计: clear all; close all; clc; f=inline('-sin(x)-cos(x.^2)'); fmax=fminbnd(f, 0, pi)
%用内联函数,求负的最小值 %注意函数名 f 不加单引号'
高教社刘卫国《MATLAB 程序设计与应用》(第二版)习题参考解答
第 6 章:MATLAB 数值计算
教材 P189 习题六
第 6 章 MATLAB 数值计算
1. 利用 MATLAB 提供的 randn 函数生成符合正态分布的 10×5 随机矩阵 A,进行如下操 作:
(1) A 各列元素的均值和标准方差。 (2) A 的最大元素和最小元素。 (3) 求 A 每行元素的和以及全部元素之和。 (4) 分别对 A 的每列元素按升序、每行元素按降序排列。 答:
运行结果: P=
0
15.0000
7.0000 -3.5000
0.5000 -2.0000 -2.0000
自动控制原理习题全解及MATLAB实验 第6章习题解答

系统开环传递函数为 G0 s
s0.1s
K
10.2s
1
,要求:
(1)系统响应斜坡信号 r(t)=t 时,稳态误差 ess 0.01 ;
(2) 系统相位裕量 ' 40 。
试用分析法设计一个串联滞后-超前校正装置。
解:(1)系统为Ⅰ型系统,在单位斜坡信号下
分稳态误差为
essr
1 k
令 essr
稳态性能与动态性能? 答:PID 兼有 PI、PD 控制的特点,它相当于提供了一个积分环节与两个一阶微分环节。
积分环节改善稳态性能,两个一阶微分环节改善动态性能。 试分别叙述利用比例负反馈和微分负反馈包围振荡环节所起到的作用。
答:二阶振荡环节的频率特性为
1
T 2S 2 2 S 1
用比例负反馈 H(s)=h
0.2s 1 0.0143s 1
(5) Gc (s)
s 1 14s 1
0.2s 1 0.0143s
(6)
G
k
(s)
s(14s
100(s 1) 1)(0.1s 1)(0.0143s
1)
' 180 [90 arctan 7 arctan(14 7) arctan(0.1 7) arctan(0.0143 7)] 41.9 40
10lg( 12)
6dB
,
最后得出 c' m 4.47rad/s>4.4rad/s
(4) 确定校正装置的转折频率
1 m
2.2rad/s ,2 m
8.8rad/s ,T 1 0.45s , 1
G(s)=
s
2.2 s
1 1
0.45s 0.11s
1 1
第六章MATLAB数据分析与功能函数

第六章MATLAB数据分析与功能函数在MATLAB中,有很多强大的数据分析和功能函数,可以帮助我们对数据进行处理、分析和可视化。
本章将介绍一些常用的数据分析和功能函数,并使用实例来介绍它们的用法。
1.统计函数MATLAB提供了丰富的统计函数,可以进行各种统计计算,例如平均值、标准差、中位数等。
- mean函数用于计算数组的平均值,例如:mean([1, 2, 3])的结果是2- std函数用于计算数组的标准差,例如:std([1, 2, 3])的结果是0.8165- median函数用于计算数组的中位数,例如:median([1, 2, 3])的结果是22.数据拟合函数MATLAB提供了拟合函数,可以用来对数据进行曲线拟合,从而得到数据的数学模型。
- polyfit函数用于进行多项式拟合,例如:x = [1, 2, 3, 4, 5]; y = [1, 4, 9, 16, 25]; p = polyfit(x, y, 2);表示对数据点进行二次多项式拟合,并返回拟合的系数。
3.数据聚类函数MATLAB提供了强大的数据聚类函数,可以将数据集分成多个不同的类别。
- kmeans函数用于进行k-means聚类,例如:x = [1, 1.5, 3, 3.5, 4, 6]; idx = kmeans(x, 2);表示对数据进行2类聚类,并返回每个数据点所属的类别。
- gmdistribution函数用于进行高斯混合模型聚类,例如:x = [1, 1.5, 3, 3.5, 4, 6]; gm = gmdistribution.fit(x', 2);表示对数据进行高斯混合模型聚类,并返回聚类结果。
4.数据可视化函数MATLAB提供了各种数据可视化函数,可以将数据以图表的形式展示出来,便于数据分析和理解。
- plot函数用于绘制折线图,例如:x = [1, 2, 3, 4, 5]; y = [1, 4, 9, 16, 25]; plot(x, y);表示将x和y的数据点用折线连接起来。
第六章 (2)MATLAB高级图形设计

。此时,Value数据就已经存储在句柄值为 ui_handle的对象内。
② 在执行的过程中若要取回该变量,可以通过以下 方式在任一Callback中获取该数据值:
Value=get(‘ui_handle’,’UserData’);
用的有hObject、 eventdata和 handles。
hObject:当前回调函数的图形对象句柄;
eventdata:预留的输入参数;
handles:存放图形窗口中所有图形对象句柄的结 构体,存储了所有在图形界面中的控件、菜单、 坐标轴对象的句柄,可以用于在function之间传 递数据。例如:handles.pushbutton1就是按钮 pushbuton1的句柄,handles.axes1就是axes1 的句柄。
global A;
25
多个回调函数之间变量的传递
方法二:直接通过对象的UserData属性来进行各个
C注al意lb:ack此之方间法的在数操据作存上取简操单作且。方便,但每个对
①象次于首设为仅变复先该V能量杂a必对存时的lu须象取,操e,一先作先的则个前会将句输变的造数柄入量 变 成据值以值 量 很存为下, 值 多储ui程因 就 不_到h序此 会 便一an即被当 。个d可覆同l特e盖一:,定掉对需的,象要对这存存象样储储中对两的,值假
➢string属性是按钮上显示的文字,默认的string属 性与Tag相同。
3.滚动条(Slider) (1)常用属性 滚动条的常用回调函数是callback,当每次单击滚动条
改变当前值时调用该函数。 (2)回调函数 常用属性有value、max、min和SliderStep。
第6章 MATLAB高级图形设计

%h_obj是图形对象的句柄; %funname是函数名,每个创建的图形对象函数名与对象
名相同; %'PropertyName'是属性名,属性名是字符串; %PropertyValue是属性值。
例: >> h_f=figure('Position',[200 300 500 400]) >> h_a1=axes('position',[0.1,0.05,.85,.85]) %创建条形图的坐标轴 练习
gcf %查找符合指定属性值的对象句柄
例
第6章 MATLAB高级图形设计
3. 句柄图形对象属性的获取和设置 a=set(h_obj,'PropertyName',PropertyValue,...)
%设置图形对象的属性值
a=get(h_obj,'PropertyName')
Uicontrol 控件对象
Uicontextmenu 上下文菜单
Uimenu 菜单
Uipanel 面板
Uitoolbar 工具栏
第6章 MATLAB高级图形设计
6.1.2 句柄图形对象的操作
1. 创建句柄图形对象 h_obj=
funname('PropertyName',PropertyValue,……)
%获取图形对象的属性值PropertyValue
第6章 MATLAB高级图形设计
例6-2 使用句柄图形对象绘制正弦曲线
>> x=0:0.1:10; >> y=sin(x); %创建无标题窗口 >> h_f=figure('Position',[200 300 300 300],'menubar','none'); >> h_a1=axes('position',[0.1,0.1,.8,.8]); >> h_t=title(h_a1,'正弦曲线'); %创建标题 >> h_l=line(x,y); %设置坐标轴刻度及刻度标注 >> set(gca,'xtick',[0 pi/2 pi 3*pi/2 2*pi 5*pi/2 3*pi]) >> set(gca,'xticklabel',{'0','pi/2','pi','3*pi/2','2pi','5*pi/2','3pi'}) >> set(gca,'xgrid','on','ygrid','on'); %设置坐标轴属性 >> set(h_l,'linewidth',2) %设置线属性 >> set(get(h_t,'parent'),'color','y') %设置标题的父对象属性 %创建矩形框 >> h_ann0=annotation(gcf,'rectangle',[0.1 0.5 .8 0.4],... 'FaceAlpha',.7,'FaceColor','red');
MATLAB图形绘制

18
【例6-4】 在图形中添加文本字符串。 x=0:0.1:10; y=sin(x); plot(x,y) xlabel('x') ylabel('y=sinx') text(0,sin(0),‘\leftarrowsin(x)=0’) % 在指定位置添
加左箭头及字符串
text(3*pi/4,sin(3*pi/4),'\rightarrowsin(x)=0.707') text(7*pi/4,sin(7*pi/4),'\leftarrowsin(x)=-0.707')
x=x+(x==0)*eps; % 用一个“机器0”小数代替0
y=sin(x)./x; % 用可逻辑运算的sin(esp)/esp近似代替sin(0)/0的极限
plot(x,y)
xlabel('x')
% 在x轴上标注x
ylabel('y=sinx/x')
% 在y轴上标注y=sinx/x
title('门函数的频谱') % 在图形上方添加标题
9
④plotyy指令的常用调用格式 plotyy(x1,y1,x2,y2) plotyy(x1,y1,x2,y2,'f') plotyy(x1,y1,x2,y2,'f1','f2')
指令中出现的参数f、f1、f2等代表绘制数据的 方式,可选择plot、semilogx、semilogy、 loglog等不同的形式。
28
⑤二维区域图的绘制 区域图的绘制使用area指令,该指令用于
在图形窗口中显示一段曲线,该曲线可由一 个矢量生成,也可由矩阵中的列生成(其实 在MATLAB中,矢量是矩阵的一种特殊形式 ,即列数为1的矩阵就是矢量)。如果矩阵的 列数大于1,则area指令将矩阵中每一列的 值都绘制为独立的曲线,并且对曲线之间和 曲线与x轴之间的区域进行填充。这种图形 在MATLAB中就称为区域图。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
【例6-9】 垂直条形图的绘制。 x=[1 2 3]; %定义条形的位置 y=[3 5 2; 4 6 8; 7 5 3]; %定义条形的高度 bar(x,y) 【例6-10】 绘制一个二维水平且堆叠的条形图。 x=[1 2 3]; y=[3 5 2; 4 6 8; 7 5 3]; barh(x,y)
3
6.1 MATLAB图形绘制基础
1. MATLAB图形绘制的基本步骤
在MATLAB中,一般按照下述的几个步骤绘制 图形。 ①准备需绘制的数据或函数,常用典型指令如下。
x=0:0.1:10; y1=bessel(1,x); y2=bessel(2,x); y3=bessel(3,x);
4
②选择图形输出的窗口及位置,常用典型指令如下。 figure(1) subplot(m,n,k) ③调用基本的绘图函数,常用典型指令如下。 plot(x,y1,x,y2,x,y3) plot3(x,y,z,'r :') ④设置坐标轴的范围、标记号和网格线,常用典型指 令如下。 axis([0,10,-3,3]) axis([x1,x2,y1,y2,z1,z2]) grid on
13
MATLAB虽然会自动为每条曲线赋予不同的颜 色以示区别,但有时却很难判断曲线和函数的对 应关系,可以通过两种方法来解决这个问题。第 一种方法,把这些曲线在同一个绘图窗口的不同 区域分别显示,把例6-2程序修改如下。
x=0:0.01:10; y1=sin(x); y2=x.*sin(x); y3=exp(2*cos(x)); subplot(2,2,1),plot(x,y1) % 在第1个子图中显示y1 subplot(2,2,2),plot(x,y2) % 在第2个子图中显示y2 subplot(2,2,3),plot(x,y3) % 在第3个子图中显示y3
17
2. 二维图形的标注 表6-3对常用的一些标注指令进行了归纳,我 们通过一些例子来说明表6-3中各指令的功能。 【例6-3】 坐标轴及标题的标注。 x=-10:0.1:10; % 定义图形的横轴坐标范围及采样步长 x=x+(x==0)*eps; % 用一个“机器0”小数代替0 y=sin(x)./x; % 用可逻辑运算的sin(esp)/esp近似代替sin(0)/0的极限 plot(x,y) xlabel('x') % 在x轴上标注x ylabel('y=sinx/x') % 在y轴上标注y=sinx/x title('门函数的频谱') % 在图形上方添加标题
28
⑤二维区域图的绘制 区域图的绘制使用area指令,该指令用于 在图形窗口中显示一段曲线,该曲线可由一 个矢量生成,也可由矩阵中的列生成(其实 在MATLAB中,矢量是矩阵的一种特殊形式 ,即列数为1的矩阵就是矢量)。如果矩阵的 列数大于1,则area指令将矩阵中每一列的 值都绘制为独立的曲线,并且对曲线之间和 曲线与x轴之间的区域进行填充。这种图形 在MATLAB中就称为区域图。
29
【例6-11】 根据矩阵数据来绘制区域图。 A=[1 2 3 4 2468 3573 7532 6 3 2 1]; area(A) %绘制区域 set(gca,'xtick',1:5) %设定x轴的标示 grid on %显示网格 set(gca,'layer','top') %将网格显示在图形之上
14
运 行 后 结 果 如 图 6.6 所 示 。 程 序 中 subplot(2, 2, 3),plot(x, y3)的含义是把绘图窗 口划分成2行2列共4个区域(可同时显示4个子 图),把y3显示在第2行第1列,即第3个子图的 位置。此时,可以方便地区分y1、y2、y3并观 察它们的形状。
15
3. 典型二维图形的绘制
①对数、半对数坐标轴图形的绘制 在有些时侯,需要的函数可能在两个坐标轴或 某个坐标轴上有较大的取值范围,这时可以通过 loglog、semilogx、semilogy等指令在x轴和(或)y轴 按对数比例绘制二维图形。
22
【例6-6】 对数、半对数坐标轴图形的绘制。 x=0:0.1:10; y=exp(x); subplot(1,3,1) % 显示在第1个子图上 plot(x,y) subplot(1,3,2) loglog(x,y) % 在x轴和y轴都按对数比例绘制图形 subplot(1,3,3) semilogy(x,y) % 在x轴按线性比例、y轴按对数比例绘
7
2. MATLAB基本绘图命令
MATLAB提供了大量的指令用于将矢量数据
以曲线图形的方式进行显示以及这些曲线图形的
注释和打印。 详见表6-1。
①plot指令的常用调用格式
plot(y,'s') plot(x,y,'s') plot(x1,y1,'s1',x2,y2,'s2') h=plot(…)
11
有时为了便于观察,可以在图形上加上网格, 此时只需在上例程序后加上grid on即可。 x=0:0.01:10; y=sin(x); plot(x,y) grid on 运行后结果如图6.4所示。
12
【例6-2】在一个图形窗口中绘制多条函数曲线。 x=0:0.01:10; y1=sin(x); y2=x.*sin(x); %y2=xsinx 2cosx y3=exp(2*cos(x)); % y3 = e plot(x,y1,x,y2,x,y3) 运行后结果如图6.5所示。
25
③极坐标图的绘制 极坐标也是一种常用的坐标形式,在有些场合 使用起来非常方便。极坐标图的绘制使用的指令是 polar,其调用格式为polar(theta,rho,linespec), 即用极角theta和极径rho画出极坐标图形,参量 linespec则可以指定极坐标图中线条的线型、标记 符号和颜色等。 【例6-8】 极坐标图的绘制。 x=0:0.01:2*pi; polar(x,sin(2*x).*cos(2*x),'r:')
制二维图形
23
②双y轴图形的绘制 利用MATLAB的plotyy指令可以同时绘制两条 函数曲线,这两条曲线共用一个x轴,而y轴则为 两个,分别位于图形的左边和右边。这时,可以 将具有不同取值范围的两条函数曲线放到一个图 形中,以便进行分析和比较。
24
【例6-7】 双y轴图形的绘制。 x=0:1000; a=1000;b=0.01;c=0.01; y1=a*exp(-b*x); y2=cos(c*x); plotyy(x,y1,x,y2,'semilogy','plot')
title('八瓣玫瑰图')
26
④二维条形图的绘制
在MATLAB中,用指令bar和barh来绘制二维条形 图,其中指令bar用来绘制垂直条形图,barh用来绘 制水平条形图。指令的调用格式为 bar(x, y, width, ‘style’, linespec) 或barh(x, y, width, ‘style’, linespec), 其中的参数width代表条形的宽度,默认值为0.8, 当width的值大于1时,条形将会出现交叠;参数style 用来定义条形的类型,可选值为group或stack,其默 认值为group,如选stack,则对mn矩阵只绘制n组条 形,每组一个条形,且条形的高度为这一列中所有元 素的和;参数linespec用来定义条形的颜色。
18
【例6-4】 在图形中添加文本字符串。 x=0:0.1:10; y=sin(x); plot(x,y) xlabel('x') ylabel('y=sinx') text(0,sin(0),‘\leftarrowsin(x)=0’) % 在指定位置添
加左箭头及字符串
text(3*pi/4,sin(3*pi/4),'\rightarrowsin(x)=0.707') text(7*pi/4,sin(7*pi/4),'\leftarrowsin(x)=-0.707')
另外一种方法,可以通过自定义曲线的颜色、 线型等来区别不同的曲线。对例6-2程序的最后 一句修改如下 plot(x,y1,'r:',x,y2,'g--',x,y3,'b-.') 运行后结果如图6.7所示 在图6.7中,用红色的虚线(在程序中用r:表示) 表示函数y1,用绿色的双画线(在程序中用g--表 示)表示函数y2,用蓝色的点画线(在程序中用b-. 表示)表示y3。这样就能方便区分同一窗口中不 同的曲线。
19
如果在不要求精确定位的情况下对图形进行 标注,还可以使用gtext指令实现以交互的方式将 标注字符串放置在图形中。例如,在图6.8中的 正弦曲线上执行下面的指令 gtext('第一个零点') gtext('第二个零点') gtext('第三个零点') 按回车键后打开图形窗口,当光标进入图形 窗口时,会变成一个大十字,表明系统正在等待 用户的动作。单击想要加入标注的地方即可。
16
此外还可以在不同函数曲线上标注不同的数据 点型以观察数据点。比如对例6-2程序的第一句 及最后一句修改如下 x=0:0.2:10; plot(x,y1,'r:+',x,y2,'g--d',x,y3,'b-.o') 修改第一句的目的是增加数据取值步长,以便 于观察数据点。 运行后结果如图6.8所示 。 在图6.8中,y1上的数据点用加号表示,y2上 数据点用菱形表示,y3上的数据点用圆形表示。
20
【例6-5】 图形的图例标注。 x=linspace(0,2*pi,50); y1=cos(x); y2=sin(x); plot(x,y1,x,y2) xlabel('x的取值范围') ylabel('y1和y2的值') legend(' y1=cos(x)', ' y2=sin(x)')