如何用Matlab绘制曲线图.docx

各位同学：
在写论文和报告时，为了很好地表达你研究和开发的结果，不仅要用文字详细地描述你方法、步骤和结果，还必须配以各种图来说明问题。下面是我们实验室张媛媛老师申请博士学位论文中的部分曲线图、硬件框图、软件流程图和实验装置原理框图。她将在部分曲线图下面给出绘制图形的Matlab程序和相关步骤，供大家学习和参考。
例一：

图2-3-6动态线性环节的输入输出信号图2-3-7模型输出和消噪后实验时数据比较
1,输入信号u(k)；2,输出信号y(k)1,实验数据；2,模型输出
绘图程序如下：
figure(1)
plot(t,y,'k',t,x,'k','LineWidth',1.4)
xlabel('Time(s)','fontname','宋体','Fontsize',9);％绘制横坐标
ylabel('Voltage(v)','fontname','宋体','Fontsize',9);％绘制纵坐标
%xlabel('时间(s)','fontname','宋体','Fontsize',9);
%ylabel('电压(v)','fontname','宋体','Fontsize',9);
％设置合适的图框大小.可将下面四句变为子程序，以便调用。
set(gcf,'color',[1,1,1]);
set(gca,'xcolor',[0,0,0],'ycolor',[0,0,0]);
set(gcf,'units','centimeters','position',[5,10,6.8,5.2]);
set(gca,'box','on','fontname','宋体','Fontsize',9);
％设置指向线的位置
annotation1=annotation(figure(1),'line',[0.55850.6038],[0.72250.6459]);
annotation1=annotation(figure(1),'line',[0.47550.4453],[0.71290.6651]);
％标注数字“1”“2”
annotation1=annotation(...
figure(1),'textbox',...
'Position',[0.38490.54860.33960.1404],...
'LineStyle','none',...
'FontSize',8,...
'String',{'1'},...
'FitHeightToText','on');
annotation1=annotation(...
figure(1),'textbox',...
'Position',[0.59740.53820.33960.1404],...
'LineStyle','none',...
'FontSize',8,...
'String',{'2'},...
'FitHeightToText','on');
例二：

图2-4-3的静态特性
绘图程序如下：
figure(1)
plot(t,y,'k',t,x,'k','LineWidth',1.4)
xlabel('时间(s)','fontname','宋体','Fontsize',9);
ylabel('质量流量(g/s)','fontname','宋体','Fontsize',9);
legend('1,实验数据','2,拟合曲线')
※当根据程序画出的线段位置不理想时，可用下面的步骤来做。
一：选中“箭头”，使图形处于可编辑状态，如图A所示。将线段拖到所需的位置。
二：右键打开下菜单，双击“showM-code”，如右图所示。
三：则打开此时线段位置所对应的m文件，例如下面。
annotation1=annotation(figure1,'line',[0.59690.6265],[0.77270.7259]);
四：将程序拷贝到自己的原程序中，注意修改figure1为figure(i)，对应该图。

图A图B
例三：
例3.1

(e)77.30g/s
图2-5-4小波消噪处理后的实验数据的动态性能指标

(a)线性补偿后(b)非线性校正后
图3-2-10阶跃输入17.35g/s时的动态非线性校正结果
此为利用程序stepspecs.m在图上标注调节时间，两条横虚线为误差带。
figure(1)
plot(x,y,’k’)
[os,ts,tr]=stepspecs(t,y,5)％输出超调量，调节时间，上升时间，并绘制上面图形。
标准图框用第一页程序。
例3.2

输入幅值112.151g/s
图3-3-4Wiener模型校正系统的在线非线性校正结果
4条虚线是给出4组数据点。画出虚线。

1,传感器模型输出；2,等效系统输出；3,动态线性补偿后
图3-2-2动态线性补