matlab 一般的实验数据展示

experData=[0.00614	164.40122
0.00706	166.68591];


%% 原函数 
d=0.001;
e=[0:d:0.2];
Q=366.70;
A=163.25;
e0=0.03425;
p=1.03529;
sigma=Q-(Q-A)./(1+power((e-0.006)./e0,p));

e=e(:);
sigma=abs(sigma(:));
% sigma=abs(sigma(:));  % 带绝对值

% close  all
figure
plot(e,sigma,'r','linewidth',1)
axis([0 0.2 100 400])
grid on
hold on
plot(experData(:,1),experData(:,2),'.b','markersize',7);
legend('Theoretical Data','Experimental Data')

%% 求导
mean_e=e(1:end-1)+d/2;
mean_e=mean_e(:);
d_sigma=diff(sigma)/d;
d_sigma=d_sigma(:);
mean_exper=mean( [experData(1:end-1,1),experData(2:end,1)],2);
d_exper1=diff(experData(:,1));
d_exper2=diff(experData(:,2));
diff_exper=d_exper2./d_exper1;

figure
plot(mean_e,d_sigma,'r','linewidth',1)
grid on
hold on
plot(mean_exper,diff_exper,'.-b','markersize',7);
legend('Derivative of Theoretical Data','Difference of Experimental Data')