MATLAB绘图常用示例语句

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’：最近邻点插值，直接完成计算；