码迷,mamicode.com
首页 > 其他好文 > 详细

Matlab之数据处理

时间:2015-08-16 10:40:32      阅读:245      评论:0      收藏:0      [点我收藏+]

标签:

写在前面的,软件不太强大,每次保存都需要生成rec和dark的文件,在处理是只需要一个就行了,所有网上查看了下运用批处理的命令去掉多余的文件:

解决办法:windows命令模式下CMD进入文件的目录,输入命令:del /q *rec*.csv(*rec*就是需要删除的包含rec的文件)

昨天下午搞了一个下午的数据提取,想过了用多种方法包括java,delphi等,最后还是想到用Matlab比较方便,选对工具真的很重要啊

我的数据格式:有多个csv文件,对我有用的只有B15:B270的数据,如何提取,看下面代码:

wave = csvread(.\dark\dark.csv,14,0,[14,0,269,0]);%看参数,具体参数的功能参考相应的Matlab书籍
dark =csvread(.\dark\dark.csv,14,1,[14,1,269,1]);%获取背景参数,当前目录的下级目录
ref = csvread(.\ref\ref.csv,14,1,[14,1,269,1]);%获取参比
B = ref-dark;
files = dir(*.csv);%获取当前目录的文件数,注意读取的是csv文件
l = length(files);
C = [];
for i = 1:l
    filename = files(i).name;%获得文件名
    A=csvread(filename,14,1,[14,1,269,1]);%读入指定列的数据,详细参数功能参考对应的代码
    D = B./(A-dark);%矩阵对应的数相除
    C = [C D];%矩阵列合并
end
E = [wave log10(C)];
%dlmwrite(result.txt,C,delimiter,\t);%以文本的形式输出,注意指定精确度precision,%f32
xlswrite(result.xlsx,E,Sheet1);%xls的形式输出
msgbox([complete!  num2str(l)  xls files.]);

后来跑数据不好,分析原因是没有去三次的平均值,我是每一个位置光谱扫了三次,但是对应的化学只有一个值(一开始我的办法是三次都是取同一个化学值,但是这样的做法是不对的),所有在原有的基础上,每隔3列取一次平均值,算法如下(代码比较简单):

A = xlsread(result.xlsx,Sheet1,B3:PE258);
[m,n]=size(A);
B=[];
for i=3:3:n%第二个参数是步长
    temp = A(:,i-2:i);
    C = mean(temp,2);
    B=[B C];
end
xlswrite(mean_3.xlsx,B,Sheet1);%xls的形式输出
msgbox([complete! ]);

以上就提取完成了,但是化学值是必须自己手动写的,那个真叫一个辛苦,420个数据还好吗,真的还好吗?

方法是通用的,以后读取xls的文件可以对应的用xlsread和xlswrite方法来写,类似能读的文本还有txt、dat......

关于常见的光谱数据的处理,对于同一个点的三次的取平均值的matlab的算法:

B=[];%产生一个空矩阵,
for i=0:14
data=zeros(396,1);%生成一个行*列为396*1的零矩阵     
for j=1:3
    s=strcat(num2str(i),-,num2str(j),.txt);%-是取得i和j之间的链接符。如果是1.2那就是.。
    k=importdata(s);
    t=k(:,2);%如果原数据的第一列是光的波数,第二列才是吸光度,那么这里写2
    data=data+t;%注意矩阵的相加
end
    data=data/3;%举证的除法运算,其中每个点都除以3
    B=[B,data];
end

 关于近红外光谱的吸光度的算法,偶然在一个java程序中看见的计算方法:

result[i] = Math.abs(Math.round(((resultQ[i] - resultDark[i]) / (double)(resultLight[i] - resultDark[i])) * 10000) / 10000.0);//此处得到的是透射率T%
result[i] = -Math.log10(result[i]);//得到吸光度abs

 Matlab实现曲线拟合

function string = curveFitting(a,x,y)
%spfit  拟合两组数据的线性关系
%a  拟合的次数
%x  x轴数据
%y  y轴数据
fx=polyfit(x,y,a)%结果返回a次拟合多项式系数,从高次到低次依次存放到变量fx中; 
yvalue=polyval(fx,x);%求得多项式在x0处的值y0;
ymean=mean(y);
Rsquare=(norm(yvalue-ymean)/norm(y-ymean))^2 %求得相关系数R的平方进行输出
SC=[];
Temp=[];
for i =1:a
    if fx(i) > 0
        if i >1
            string = [+,num2str(fx(i))];
        else
            string = num2str(fx(i));
        end
    else
        string = num2str(fx(i));
    end
    Temp = [string,[x^,num2str(a+1-i)]];
    SC = [SC,Temp];
end
if fx(a+1) > 0
    lastString = [+,num2str(fx(a+1))];
    SC =[SC,lastString];
else
    SC = [SC,num2str(fx(a+1))];
end
%画图程序(x,y)对应的点,x1y1对应的是带入拟合后的点
x1=x(1):0.2:x(length(x));%0.2是步长,不是特别要求
y1=polyval(fx,x1);
plot(x,y,b*,x1,y1,g);
xlabel(strcat(拟合方程:y=,SC));

技术分享

Matlab批量到数据(二维荧光数据处理)

B=[];%产生一个空矩阵
for i=0:14
data=zeros(396,1);%生成一个行*列为396*1的零矩阵     
for j=1:3
    s=strcat(num2str(i),-,num2str(j),.txt);%-是取得i和j之间的链接符。如果是1.2那就是.。
    k=importdata(s);
    t=k(:,2);%如果原数据的第一列是光的波数,第二列才是吸光度,那么这里写2
    data=data+t;
end
    data=data/3;
    B=[B,data];
end

 

Matlab之数据处理

标签:

原文地址:http://www.cnblogs.com/sowhat4999/p/4439862.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!