代码改变世界

MATLAB绘制函数图

2013-09-04 10:25  chloe_zhou  阅读(10133)  评论(0编辑  收藏  举报

 序言

   Matlab可以根据用户给出的数据绘制相应的函数图.对于单个2D函数图,需要给出一个行向量x作为函数图上离散点集的横坐标,以及一个与x列数一样的横坐标y作为函数图上点集的纵坐标. 向量xy的取值既可以用离散数值直接列出,也可以用取值范围(冒号)和点运算(在幂运算、除法等算术运算符之前加个点号)得出相应的向量.

 正文

1. 定义函数, 绘制初步图形

 绘制函数y=f(x)的图形分为三步:定义x的取值范围和取值间隔/取值个数;通过定义yx的关系f(x)来给出y向量;用plot(x,y)命令绘制函数图形.

  (1) 定义x的取值范围和取值间隔用x=(a:k:b) 的形式, 表示x属于范围[a,b], 且间隔为k,即x的第一个值为a, 第二个值为a+k, 第三个值为a+2k, 以此类推, 最后一个值为b. 例如:

  

  其中k值越小则同一定义域范围上的点数越多,绘出的函数曲线越光滑.

      另外还可以用linspace函数生成x的值, x=linspace(a,b,n)是为x生成ab之间均匀的n个值, 结果是一个1*n的行向量(这n个数之间的差是相等的). 如果不指定个数n, 仅仅用linspace(a,b)则默认生成ab之间的100个值.

  (2) 定义yx之间的关系直接利用Matlab的算术运算符和已经定义的函数来给出表达式,例如y=sin(x)、y=5*x-1等等. 注意,如果涉及xn次幂, 除法或指数函数,则不能直接用x^n, x/n, a^x而要用x.^nx./n或者a.^x来表示,指对x的分量进行运算得到y的对应的分量.

  (3) 完成前面两步后,最后直接用plot(x,y)就能绘制出函数图形,其中x作为横轴,y作为纵轴. 例如函数在[0, 1]之间的图形绘制代码和结果如下:

s=linspace(0,1,50);
y=-(s-2+sqrt(5*s.^2-8*s+4))/2;
plot(s,y);

  

  如果我们需要在一个图形浮动窗口中同时绘出多个相关的函数以进行对比,则给出的y要一个矩阵,它的列数与x的列数相等,行数等于函数个数,其第i行表示第i个函数上离散点的纵坐标取值.

  例如,若我们想绘出三个函数 y1=sin(x)/√(x+1), y2=sin(x/2)/√(x+1), y3=sin(x/3)/√(x+1)的对比图,我们可以用如下代码:

x=[0:0.2:20];
y(1,:)=sin(x)./sqrt(x+1);
y(2,:)= sin(x/2)./sqrt(x+1);
y(3,:)= sin(x/3)./sqrt(x+1);
plot(x,y);

   运行结果如下:

  

 2. 代码设置绘图属性

 在绘图时,可以用代码直接设置曲线颜色,图例,默认坐标轴的名称等等,可参考如下代码

h1=plot(x,y(1,:),':b');%无记忆行走,蓝色虚线条
hold on;%保证不覆盖前面已绘的图形
h2=plot(x,y(2,:),'--g');%一步记忆
hold on;
h3=plot(x,y(3,:),'-r');%两步记忆
hold on;
legend([h1,h2,h3],'QW0M','QW1M','QW2M');%为各个曲线分别添加图例
hold on;
xlabel ('Position');ylabel('Probability');%设置横纵坐标的标签名称

 3. 单独绘制坐标轴

 单独绘制轴线的思路是在轴线位置画直线,然后在两直线的端点处画三角形,然后填充三角形作为箭头。

%设置函数定义域和值域的绘图范围
rangeXmin=-3.5;
rangeXmax=3.5;
rangeYmin=-1.5;
rangeYmax=1.5;
%axis off;关闭默认的轴显示
hold on;%将后续绘图保持在本图形窗口

x=linspace(rangeXmin,rangeXmax,50);
y1=sin(x);y2=x;%设置两个函数
plot(x,y1,'r',x,y2,'g');%y1函数曲线为红色r,y2函数曲线为绿色g
axis([rangeXmin*1.1 rangeXmax*1.1 rangeYmin*1.1 rangeYmax*1.1]);%设置坐标范围
line([rangeXmin,rangeXmax],[0,0],'color','black');%绘制x轴直线,黑色
line([0,0],[rangeYmin,rangeYmax],'color','black');%绘制y轴直线
%绘制x轴箭头
arrowXx=[rangeXmax 0.97*rangeXmax 0.97*rangeXmax rangeXmax];
arrowXy=[0 0.03*rangeYmax -0.03*rangeYmax 0];
plot(arrowXx,arrowXy,'black');
fill(arrowXx,arrowXy,'black');%填充三角形
%绘制y轴箭头
arrowYx=[0 0.03*rangeXmax -0.03*rangeXmax 0];
arrowYy=[rangeYmax 0.97*rangeYmax 0.97*rangeYmax rangeYmax];
plot(arrowYx,arrowYy,'black');
fill(arrowYx,arrowYy,'black');
%添加坐标轴名称
text(0.95*rangeXmax,0.1*rangeYmax,'x');
text(0.05*rangeXmax,rangeYmax,'y');
%添加重点坐标以及刻度
text(0.05*rangeXmax,-0.05*rangeYmax,num2str(0),'Color','black');%原点坐标
text(pi/2,-0.1*rangeYmax,'\pi/2');%默认文本颜色为黑色,π字符用\pi输出
text(pi/2,0,'|');%刻度
text(-pi/2,-0.1*rangeYmax,'-\pi/2');text(-pi/2,0,'|');
text(pi,-0.1*rangeYmax,'\pi');text(pi,0,'|');
text(-pi,-0.1*rangeYmax,'-\pi');text(-pi,0,'|');
text(0.05*rangeXmax,1,'1');text(0,1,'-');
text(0.05*rangeXmax,-1,'-1');text(0,-1,'-');

 绘图结果如下

 

注意,其中line(X,Y,'propertyname','propertyvalue')函数的参数X为行向量,给出直线起点和终点的横坐标,Y也为行向量,给出直线起点和终点的纵坐标;text(x,y,'textvalue')是在坐标(x,y)处添加文本'textvalue'。

 4. 手动设置函数图属性

 在绘制出函数曲线之后, 我们可以通过图形窗口上的菜单和按钮来对函数图进行编辑, 如添加图例说明、设置坐标标签、更改曲线颜色、显示某些点的坐标值等等。

  (1) 添加图例. 利用图例(Legend),我们可以明确地为每一条函数曲线取一个名字. 方法是点击Figure窗口上的Insert Legend按钮, 此时图形的右上角会出现图例框, 其默认的曲线名称为Data1, Data2, Data3. 双击曲线名称文本, 就可以重新命名曲线了. 右击图例框, 利用弹出的快捷菜单还可以设置图例文本字体、文本颜色、图例外框颜色等等外观属性。

  

  (2) 设置坐标标签. 对于我们前面的函数, 横坐标是x, 纵坐标是y1,y2,y3, 我们可以通过Figure窗口的菜单Insert-X Label, Insert-Y Label来添加坐标标签. 如果要删除标签, 右击标签选择Delete即可.

  (3) 更改曲线颜色粗细形状等属性. 点击Figure窗口上的Edit Plot按钮, 然后双击某一条曲线, 在弹出的属性设置窗口中设置参数即可.

  (4) 显示曲线上某些点的坐标值. 如果曲线上的某个点比较重要, 需要额外标记出其坐标值, 则单击Figure窗口上的Data Cursor按钮, 然后单击曲线上需要标记的点即可. 如果想额外再多标记几个点, 则先右击Figure区域, 选择Create New Datatip, 然后再单击待标记的点. 如果想删除某个标记点, 右击那个点, 选择Delete Current Datatip.