前言
之前完成老师任务的时候写的,里面有一个错误,就是误差分析的时候有个公式错了,自己注意下。
代码
%作者:山科_xxin
%时间:2016.10
clear;clc;
[NUM]=xlsread('G:\xuanfuwu\nihe.xlsx');%读取文件
[NUM2]=xlsread('G:\xuanfuwu\jingdu.xlsx');
[h,l]=size(NUM);%获得行列数
[h2,l2]=size(NUM2);
num=NUM(1,2:l);
for i=1:h-1
a{i}=corr2(NUM(i+1,2:l),num)
end
[ma,I]=max(abs(cell2mat(a)))%求绝对值的最大值
B=NUM(I+1,1);
disp(['相关性最大波段;',num2str(B),'nm','相关系数;',num2str(ma)])%显示
x=NUM(I+1,2:l);
y=num;
scatter(x,y,'b');
hold on
x2=NUM2(I+1,2:l2);
y4=NUM2(1,2:l2);
scatter(x2,y4,'r');
hold on
y2=polyfit(x,y,1);%拟合曲线
yi=polyval(y2,NUM2(I+1,2:l2));%精度评价