MATLAB从文件得出数据并计算吸收光谱
这一系列就是科研用的真实程序了，也是对自己的一个备忘录
真的收购每次都重写了
但真的文件太多了找不到啊！！！
好吧是我废物
废话不多说，这就开始
基础的清理：
clear
clc
close all读取文件中的数据（来自示波器直接保存）
data=csvread('1\DC\Tek000_ALL.csv',13,0);%csvread只能读取纯数据，跳过前面的字符
%此时data中应该储存了对应的五列数据，我提取了两个通道，中间有一个空列
%读取另一个文件
data2=csvread('1\TDC\Tek000_ALL.csv',13,0);
%将data分成两部分
ch4 = data(:,2);%4通道结果
MATH1 = data(:,5);
T_MATH1 = data2(:,5);
%分别得到两列数据
详细的注释写在每一行了
csvread函数读取文件中的数值，但是只能读取数字，需要汉字会读取失败
这里从第十三行开始读取，结果是一个矩阵
首先绘制原始数据的图像：
%绘制原始图像：
figure(1)
x = 1:length(ch4);%两信号的x轴理论上是相同的
[AX,H1,H2] = plotyy(x,ch4,x,MATH1,'plot');
set(get(AX(1),'Ylabel'),'String','ch4');%设置左边的y轴标签
set(get(AX(2),'Ylabel'),'String','MATH1');%设置右边的y轴标签
xlabel('Time');
set(H1,'Linestyle','--');%设置曲线1的线型
set(H2,'Linestyle',':');%设置曲线2的线型这里用到了在同一张图中绘制两个y轴的数据
在我的MATLAB基础操作专栏中有详细解释：
https://blog.csdn.net/fendoudecheng/article/details/130218116
实验数据，当然是不会放效果图的，自己想象~
使用移动平均函数平滑曲线
主要是测量的时候噪音较大：
%可以看到还是有很多高频的震荡，这里再次取平均
%直接使用移动均值函数
%M = movmean(A,k) 返回由局部 k 个数据点的均值组成的数组
ave_ch4 = movmean(ch4,16);%4通道结果
ave_MATH1 = movmean(MATH1,16);%1通道相邻16个点取平均后的结果
ave_T_MATH1 = movmean(T_MATH1,16);%1通道相邻16个点取平均后的结果
这里分别取十六个采样点进行平均
移动平均就是一个数组相邻的十六个取平均
原数组长度：200000，改变后的长度：200000
注意这可不是变成十六分之一哦
绘制平滑后的曲线：
% 另一种方式绘制：
figure(2)
yyaxis left
plot(x,ave_ch4)
yyaxis right
plot(x,ave_MATH1)
set(H1,'Linestyle','--');%设置曲线1的线型
set(H2,'Linestyle',':');%设置曲线2的线型同样是一张图两个y轴，注意这里用了一种新的方法
按照MATLAB官方的说法，这个更提倡
画图：
%比较直接透射和经过铷泡
figure(3)
plot(x,10*ave_MATH1,x,5*ave_T_MATH1);
hold on
plot(x,5*ave_T_MATH1-10*ave_MATH1);我这里直接的投射是单调递减的，用来标定一个周期：
%得到真正的谱
[m1,index1] = max(ave_T_MATH1);%取数组的最大值和对应的序号
[m2,index2] = min(ave_T_MATH1);%取数组的最小值和对应的序号
x = index1:index2;
figure(4)
plot(x,5*ave_T_MATH1(index1:index2)-10*ave_MATH1(index1:index2));最后是归一化：
%归一化
% T = 5*ave_T_MATH1(index1:index2)-10*ave_MATH1(index1:index2);
% T = T/max(T);
ave_MATH1 = ave_MATH1(index1:index2);
ave_MATH1 = max(ave_T_MATH1)-ave_MATH1;
ave_MATH1 = ave_MATH1/max(ave_MATH1);
figure(5);
plot(x,ave_MATH1);注意这里，求归一化的过程，也就是归一化怎么进行：
1、用整张图的最大值减去整个图，得到吸收谱，归一化
2、用透射谱减，得到吸收谱，归一化
这里不太精确，用第一种就行
最后是保存在文件内：
saveas(gca,'my2.png')结束
放个图做封面吧：
（图片来自于网络）
上图搜索关键词：做实验哪有不疯的