a1=load(‘4mm06mm6x6x3e60.dat’); %打开参考数据文件
b1=load(‘4mm06mm6x6x3e60.rcs’); %%打开计算数据文件
c=max(abs(a1(:,2)));%将参考数据绝对值最大的放在c中
n=numel(a1(:,2));%所有数据的个数
d=a1(:,2)-b1(:,2);%两份数据对应项差值
d1=abs(d);%差值的绝对值
e=d1/c;%差值的绝对值除以参考数据绝对值最大值
% e1=e.*e;
e2=dot(e,e);%求每项值的平方
% e3=sum(e1);
e4=e2/n;
rms1=sqrt(e4);