画图

plot(1,1,'r.','markersize',50)

21.81.61.41.210.8

例二：

0.60.40.2000.20.40.60.811.21.41.61.82 y=sin( x )曲线图10.8x=0:0.05:4*pi; y=sin(x); plot(x,y) grid on, title(' y=sin( x )曲线图' ) xlabel(' x = 0 : 0.05 : 4Pi ')

0.60.40.20-0.2-0.4-0.6-0.8-102468 x = 0 : 0.05 : 4Pi 101214 title xlabel ylabel text grid hold 图形标题 x坐标轴标注 y坐标轴标注 标注数据点 给图形加上网格 保持图形窗口的图形 表4.1.1.1 Matlab图形命令 §4.1.2多重线 例三： x=0:0.1:9; y1=2*x; y2=x+1; y3=x.^2; plot(x,y1,x,y2,x,y3) 908070605040302010001234567 例四：

x=0:0.1:9; y1=2*x; y2=x+1; y3=x.^2; plot(x,y1); hold on plot(x,y2); plot(x,y3);

线型 - : -. -- 实线 虚线 点划线 间断线 . o x + * s d ^ v > < p h 908070605040302010001234567点标记 点 小圆圈 叉子符 加号 星号 方格 菱形 朝上三角 朝下三角 朝右三角 朝左三角 五角星 六角星 y m c r g b w k 颜色 黄 棕色 青色 红色 绿色 蓝色 白色 黑色 表4.1.3.1线型和颜色控制符 例一1例4：fplot('sin',[0 4*pi]) 0.80.60.40.20-0.2-0.4-0.6-0.8-1024681012

例5：fplot('sin(1 ./ x)', [0.01 0.1])

10.80.60.40.20-0.2-0.4-0.6-0.8-10.010.020.030.040.050.060.070.080.090.1109例6：fplot('abs(exp(-j*x*(0:9))*ones(10,1))',[0 2*pi],'-o')

87632100123456

例7：fplot('[sin(x), cos(x) , tan(x)]',[-2*pi 2*pi -2*pi 2*pi])

20

例5：半对数坐标绘图 t=0.001:0.002:20; y=5 + log(t) + t; semilogx(t,y, 'b') hold on

semilogx(t,t+5, 'r')

例6：极坐标绘图 t=0:0.01:2*pi; polar(t,sin(6*t))

例7：a=normrnd(10,2,10000,1);

histfit(a)

例8：比较正态分布与平均分布的分布图： yn=randn(30000,1); %%正态分布 x=min(yn) : 0.2 : max(yn); subplot(121) hist(yn, x)

yu=rand(30000,1); %%平均分布 subplot(122) hist(yu, 25)

例10：用函数patch绘制填充图 fplot('humps',[0,2],'b') hold on

patch([0.5 0.5:0.02:1 1],[0 humps(0.5:0.02:1) 0],'r'); hold off

title('A region under an interesting function.') grid （humps（一个Matlab演示函数））

-2-4-6-6-4-20246302520151050-510-310-210-110010110290 112060 0.8 0.615030 0.4 0.2180021033024030027035030025020015010050002468101214161825001400120020001000150080010006004005002000-505000.51A region under an interesting function.100806040200-2000.20.40.60.811.21.41.61.82 例9：绘制子图 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(2*x),'-.*'); title('sin(2x)'); subplot(2,2,4);

plot(x,cos(3*x),':d') title('cos(3x)')

例11：用函数fill绘制填充图 x=0:pi/60:2*pi; y=sin(x); x1=0:pi/60:1; y1=sin(x1); plot(x,y,'r'); hold on

fill([x1 1],[y1 0],'g')

sin(x)10.50-0.5-1024sin(2x)10.50-0.5-10246810.50-0.5-1026810.50-0.5-102cos(x)4cos(3x)6846810.80.60.40.20-0.2-0.4-0.6-0.8-101234567例11：绘制sin(r)/r函数的图形．建立图形用以下方法： x=-8:.5:8; y=x';

x=ones(size(y))*x; y=y*ones(size(y))'; R=sqrt(x.^2+y.^2)+eps; z=sin(R)./R; mesh(x,y,z)

(各语句的意义是：首先建立行向量x，列向量y；然后按向量的长度建立1-矩阵；用向量乘以产生的1-矩阵，生成网格矩阵，它们的值对应于x-y坐标平面；接下来计算各网格点的半径；最后计算函数值矩阵Z．用mesh函数即可以得到图形．上述命令系列中的前4行可用以下一条命令替代：[x, y]=meshgrid(-8:0.5:8))

10.50-0.51051005-50-5-10-10

例13：（附带等高线的图） [x,y]=meshgrid([-4:.5:4]); z=sqrt(x.^2+y.^2); meshc(x,y,z)

63210420-2-4-4-2204

例14：增加了边界面屏蔽的图 [x,y]=meshgrid([-4:.5:4]); z=sqrt(x.^2+y.^2); meshz(z)

632102015105005151020例15：球体图形 sphere(30); axis square;

（若只输入sphere画图，则是默认了n=20的情况．）

10.50-0.5-110.50-0.5-1-1-0.50.501:例16：绘制地球表面的气温分布示意图(着色的三维表面图)． [a,b,c]=sphere(40); t=abs(c); %求绝对值 surf(a,b,c,t); axis equal

colormap('hot')

@@@坐标轴的控制函数axis，调用格式如下： axis([xmin,xmax,ymin,ymax,zmin,zmax]) 用此命令可以控制坐标轴的范围． 与axis相关的几条常用命令还有：

axis auto 自动模式，使得图形的坐标范围满足图中一切图元素 axis equal 严格控制各坐标的分度使其相等 axis square 使绘图区为正方形

axis on 恢复对坐标轴的一切设置 axis off 取消对坐标轴的一切设置

axis manual 以当前的坐标图形的绘制 grid on 在图形中绘制坐标网格． grid off 取消坐标网格．

10.50-0.5-110.50-0.5-110.50-0.5-1

例17：

[x, y]=meshgrid(-4:.2:4); R=sqrt(x.^2+y.^2); z=-cos(R); mesh(x,y,z)

xlabel('x\\in[-4,4]','fontweight','bold'); ylabel('y\\in[-4,4]','fontweight','bold');

zlabel('z=-cos(sqrt(x^2+y^2))','fontweight','bold');

title('旋转曲面','fontsize',15,'fontweight','bold','fontname','隶书');

122z=-cos(sqrt(x+y))旋转曲面0.50-0.5-14242-2-4-40-2x[-4,4]0y[-4,4]例18：已知如下点列，求其回归直线，并计算最小误差平方和． X 0.1 Y 42 0.11 43.5 .12 45 .13 45.5 .14 45 .15 47.5 .16 49 .17 53 .18 50 .2 55 .21 55 .23 60 参考的程序如下：

x=[0.1 0.11 .12 .13 .14 .15 .16 .17 .18 .2 .21 .23]; y=[42 43.5 45 45.5 45 47.5 49 53 50 55 55 60]; 60n=length(x); 5856xb=mean(x);

yb=mean(y);

52x2b=sum(x.^2)/n;

50xyb=x*y'/n;

48b=(xb*yb-xyb)/(xb^2-x2b);

46a=yb-b*xb;

44y1=a+b.*x;

42plot(x,y,'*',x,y1);

400.1serror=sum((y-y1).^2)

例19：多项式的曲线拟合 x=[1 2 3 4 5];

y=[5.6 40 150 250 498.9];

p=polyfit(x,y,n) %数据的n次多项式拟合 x2=1:0.1:5; %%n取1时，即为最小二乘法 y2=polyval(p,x2); %计算多项式的值 %（polyvalm计算矩阵多项式） plot(x,y,'*',x2,y2);grid on

(poly：矩阵的特征多项式、根集对应的多项式) [图为n=3时的拟合曲线]

linewidth markeredgecolor makerfacecolor markersize 0.120.140.160.180.20.220.2450045040035030025020015010050011.522.533.4.55指定线条的粗细 指定标记的边缘色 指定标记表面的颜色 指定标记的大小

例20：对图线添加标注

x = 0:.2:12;

plot(x,bessel(1,x),x,bessel(2,x),x,bessel(3,x));

legend('First','Second','Third','Location','NorthEastOutside')

0.60.50.40.30.20.10-0.1-0.2-0.3-0.4024681012FirstSecondThird

'North' 图例标识放在图顶端 'South' 图例标识放在图底端 'East' 图例标识放在图右方 'West' 图例标识放在图左方 'NorthEast' 图例标识放在图右上方（默认） 'NorthWest 图例标识放在图左上方 'SouthEast' 图例标识放在图右下角 'SouthWest' 图例标识放在图左下角 (以上几个都是将图例标识放在框图内） 'NorthOutside' 图例标识放在图框外侧上方 'SouthOutside' 图例标识放在图框外侧下方 'EastOutside' 图例标识放在图框外侧右方 'WestOutside' 图例标识放在图框外侧左方 'NorthEastOutside' 图例标识放在图框外侧右上方 'NorthWestOutside' 图例标识放在图框外侧左上方 'SouthEastOutside' 图例标识放在图框外侧右下方 'SouthWestOutside' 图例标识放在图框外侧左下方 'Best' 图标标识放在图框内不与图冲突的最佳 'BestOutside' 图标标识放在图框外使用最小空间的最佳位置

例20:对 y=cosx的数据进行一维插值，比较各种插值方法 x=[-2*pi:0.5*pi:2*pi];y=cos(x);xi=-2*pi:0.3*pi:2*pi; y_nearest=interp1(x,y, xi,'nearest'); y_linear= interp1(x,y,xi);

y_spline= interp1(x,y,xi, 'spline' ); y_cubic= interp1(x,y,xi, 'cubic' );

plot(x,y,'y',xi,y_nearest,'-m',xi,y_linear, 'c ', xi,y_spline,'r',xi,y_cubic,'g'); legend ('original data','nearest','linear','spline','cubic')

YI = interp1(X,Y,XI,'method') 'method'是插值所用的方法. 'nearest' - nearest neighbor interpolation 'linear' - linear interpolation 'spline' - piecewise cubic spline interpolation (SPLINE) 'pchip' - piecewise cubic Hermite interpolation (PCHIP) 'cubic' - same as 'pchip'

10.80.60.40.20-0.2-0.4-0.6-0.8-1-8-6-4-202468original datanearestlinearsplinecubic例21: 下表给出了美国从1950-1990年工作年限10,20,30年的工资情况，使用插值计算1975年工作15年的工资 Y S 10年 1950 1960 1970 1980 1990 150.697 179.323 203.212 226.505 249.633 20年 199.592 195.072 179.092 153.706 120.281 30年 187.625 250.287 322.767 426.730 598.243 years = 1950:10:1990; service = 10:10:30;

wage = [150.697,199.592,187.625;

179.323,195.072,250.287; 203.212,179.092,322.767; 226.505,153.706,426.730; 249.633,120.281,598.243];

w = interp2(service,years,wage,15,1975,’linear’) 结果为：w = 190.6288

例22. 在某处测得海洋不同深处水温 如下： 深度(m) 水温(℃) 446 7.04 714 4.28 950 3.40 1422 2. 8原始结果spline 插值结果71634 2.13

解：h=[446,714,950,1422,1634]; w=[7.04,4.28,3.40,2.,2.13]; hi=[500,1000,1500]; p=spline(h,w); ci=ppval(p,hi)

plot(h,w, 'or') ,hold on,plot(hi,ci,'*') legend ('原始结果','插值结果') ZI = interp2(X,Y,Z, XI ,YI ,method) method有： ‘linear’：双线性插值算法（缺省算法）. ‘nearest’：最近邻插值. ‘spline’：三次样条插值. ‘cubic’：双三次插值.

632400600800100012001400160018009例23：用matlab绘制直方图

X=[84 82 60 75 66 92 92 80 86 85 100 84 91 95 78 85 81 85 91 83

82 87 72 91 81 63 88 88 88 92 93 98 84 84]; hist(X，n)（n为直方图的等分数，图中n为10）

例24：用matlab画柱形图

y=[300 311;390 425; 312 321;

250 185; 550 535; 420 432; 410 520;]; b=bar(y);

例25：用matlab画饼图 A=[0.5,1,2,2.3,10,8.8,4] f1=A(find(A<2&A>=0)) f2=A(find(A<6&A>=2)) f3=A(find(A>=6))

ratio=[length(f1) length(f2) length(f3)]/length(A) format rat ratio

pie(ratio)

legend('0-2','2-6','>6')

例26：用matlab画折线图 years=2001:2010;

h=[165 168 172 176 178 179 180 181 182 183]; plot(years,h,'-*')

例27：用matlab画折线图

8763210626670747882869094986005004003002001000123456729%29%0-22-6>3%1841821801781761741721701681661200120022003200420052006200720082009201010y=x， 当x>0 y=sin(x), 当x <=0

86x=-10:0.01:10;

y=x .*(x>0) + sin(x) .*(x <=0); plot(x,y)

420-2-10-8-6-4-20246810