MATLAB绘图常用示例语句
clear;%初始化数据
clc;%清屏
close all;%关闭所有绘图
x=[];
y=[];
plot(x,y,'k.');%绘制曲线，单引号内为线型与颜色
plot(x1,y1,'k+-',x2,y2,'k--',x3,y3,'k-.','linewidth',2.5);%在一幅图上同时绘制多条曲线，及线宽设置
%字母颜色标点线型
% y 黄色? 点线
% m 粉红○圈线
% c 亮蓝××线
% r 大红＋＋字线
% g 绿色－实线
% b 蓝色* 星形线
% w 白色：虚线
% k 黑色(--) 点划线
xlabel('频率(Hz)');%横坐标说明
ylabel('流量(L/min)');%纵坐标说明
title('');%图形标题
legend('','','');%线型说明
axis([xmin xmax ymin ymax])%设置x轴，y轴范围
grid on;%网格
hold on;%在原图上继续绘图
figure;%重新生成新的绘图
fid1=fopen('F:\实验数据\D6_100V_1Hz_2.txt','r'); %读取.txt 数据
title1=fgetl(fid1);
title2=fgetl(fid1);
title3=fgetl(fid1);
title4=fgetl(fid1);%去除txt文件中的前四行
fz=50000;%设置采样频率
dianshu=200000;%所取点数
[A,N]=fscanf(fid1,'%f',[2,dianshu+100]);%从文件中读出数据[数据数组，总数据个数]（数据存放变量，'数据类型'，[列数，行数]）
x=A(1,1:dianshu);%读取数据第1列
x1=x';%行列向量转置
y=A(2,1:dianshu);%读取数据第二列
%以下为滤波程序（主要调节参数c的大小）
c=100;
wn=0.1;
fs=50000; %采样频率；
b=fir1(c,wn/(fs/2),hamming(c+1));
y1=filtfilt(b,1,y);%对y滤波
pj=mean(y1);%对y1求均值
plot(x,y1);
cz=interp1(x,y1,a);%求图形中横坐标为a，对应的纵坐标的值figure;
[ax,h1,h2]=plotyy(t1,y1,t1,y2); %双坐标函数，各参数意义，参考帮助文件
set(ax(1),'XLim',[0,10000],'XTick',[0:1000:10000]);%设置第一x 坐标范围及间距
set(ax(1),'YLim',[0,3000],'YTick',[0:300:3000],'ycolor','k'); %设置第一y坐标范围及间距，及颜色
set(ax(2),'XLim',[0,10000],'XTick',[0:1000:10000]);%应同ax(1),’XLim’, 否则会错位
set(ax(2),'YLim',[0,50],'YTick',[0:5:50],'ycolor','k');%两纵坐标的间距数应一致，否则也会错位
set(h1,'linestyle','-','color','k');%设置h1曲线线型、颜色
set(h2,'linestyle',':','color','k'); %设置h2曲线线型、颜色xlabel('频率（Hz）');%横坐标标识
set(get(ax(1),'Ylabel'),'String','位移（um）'); %第1纵坐标标识set(get(ax(2),'Ylabel'),'String','流量/（L/min)'); %第2纵坐标标识
legend([h1 h2],'阀芯位移','出口流量'); %添加曲线说明
set(h1,'linewidth',2);%设线型宽度
set(h2,'linewidth',2);
grid on;%显示网格
hold on;%保持窗口不变，以便添加曲线
%直线拟合(a,b）
a=[0 0.5 2 3.5 5];
b=[0.00411 0.02030 0.06768 0.115 0.16343];
c=1;
p=polyfit(a,b,c);%参数c为拟合曲线次数，p为降幂系数
y=polyval(p,a);
plot(a,y,a,b);
%曲线求导,求速度与加速度
v=(diff(s)./diff(t));
plot(t1,v)%时间点数比位移曲线对应时间点数少一点
F=A(:,1);%读取excel文件数据，Inport Date时应把变量名改成A,F为第一列所有点；A(m,n)表示第m行第n列的点
max(A);%A阵列中的最大值
[y,i]=max(A);%y为最大值，i为最大值所在位置
syms a b;
ezplot(a^2+b^2-1);%作方程'a^2+b^2=1'的图形;
fc=var(A);%求向量A的方差
bzc=std(A);%求向量A的标准差
%求积分
syms x;
int(x^2+x,x,1,2);%函数x^2+x对x在区间[1,2]上求积分
%解方程1
syms I;
f=-pi*((0.25e-4)/S+(8e-7)*I*S*((S^2+1.6e-5)/(S^2-1.6e-5)))-B; B=0.1;
S=0.2;
solve(f,'I')%直接求出数值解
%解方程2
B=0.1;
S=0.2;
X=solve('B=-pi*((0.25e-4)/S+(8e-7)*I*S*((S^2+1.6e-5)/(S^2-1.6e -5)))','I')
subs(X)%将B,S值代入到X中
%求插值
yi=interp1(X,Y,xi,'method');%xi为所求插值的横坐标
%‘nearest’：最近邻点插值，直接完成计算；
%‘spline’：三次样条函数插值；
%‘linear’：线性插值（缺省方式），直接完成计算；%‘cubic’：三次函数插值；
s=xlsread('A.xlsx',1,'A1:B2');%读取EXCL文件数据。