最近在做数据采集工作,需要将labview采集回来的lvm格式文件用matlab进行处理,因此做了如下的处理:
1.读取采集回来的数据
2.对数据做离散傅里叶变换
于是做了如下程序,需要示例数据请留言,可提供此格式的样例数据
本实验涉及几组对照实验,用的时候只需将注释取消即可
clc;
clear all;
N=256;
Fs=25600;
row_num_first=100;
% [s_0]=textread('0.lvm','%s');
% figure(1)
% s_0=str2double(s_0(row_num_first:24000,1));
% plot(s_0)
% xlabel('simple');ylabel('amplitude');title('----')
% [s_10,s_20]=textread('10_20.lvm','%s%s');
% row_num_first=100;
% s_10=str2double(s_10(row_num_first:24000,1));
% s_20=str2double(s_20(row_num_first:24000,1));
% figure(2)
% plot(s_10)
% xlabel('simple');ylabel('amplitude');title('-----')
% figure(3)
% plot(s_20)
% xlabel('simple');ylabel('amplitude');title('-------')
% Y=fft(s_0,N);
% Y1=fftshift(Y);
% f=(-N/2+1:N/2)*Fs/N;
% Mag=abs(Y1);
% figure(4);
% plot(f,10*log10(Mag));
%
% Y_10=fft(s_10,N);
% Y1_10=fftshift(Y_10);
% f=(-N/2+1:N/2)*Fs/N;
% Mag1=abs(Y1_10);
% figure(5);
% plot(f,10*log10(Mag1));
%
% Y_20=fft(s_20,N);
% Y1_20=fftshift(Y_20);
% f=(-N/2+1:N/2)*Fs/N;
% Mag2=abs(Y1_20);
% figure(6);
% plot(f,10*log10(Mag2));
%% %
% [s4_0,s4_1,s4_2,s4_3]=textread('4-1-0.lvm','%s%s%s%s');
% figure(1)
% s4_0=str2double(s4_0(row_num_first:24000,1));
% plot(s4_0)
% xlabel('simple');ylabel('amplitude').
% [s4_10,s4_20,s4_30]=textread('4-1-10-20.lvm','%s%s%s');
% s4_10=str2double(s4_10(row_num_first:24000,1));
% s4_20=str2double(s4_20(row_num_first:24000,1));
% figure(2)
% plot(s4_10)
% xlabel('simple');ylabel('amplitude');
% figure(3)
% plot(s4_20)
% xlabel('simple');ylabel('amplitude');
%
% Y=fft(s4_0,N);
% Y1=fftshift(Y);
% f=(-N/2+1:N/2)*Fs/N;
% Mag=abs(Y1);
% figure(4);
% plot(f,10*log10(Mag));
%
% Y_10=fft(s4_10,N);
% Y1_10=fftshift(Y_10);
% f=(-N/2+1:N/2)*Fs/N;
% Mag1=abs(Y1_10);
% figure(5);
% plot(f,10*log10(Mag1));
%
% Y_20=fft(s4_20,N);
% Y1_20=fftshift(Y_20);
% f=(-N/2+1:N/2)*Fs/N;
% Mag2=abs(Y1_20);
% figure(6);
% plot(f,10*log10(Mag2));title
%%
% [s4_0,s4_1,s4_2,s4_3]=textread('w-0.lvm','%s%s%s%s');
% row_num_first=100;
% figure(1)
% s4_0=str2double(s4_0(row_num_first:24000,1));
% plot(s4_0)
% xlabel('simple');ylabel('amplitude');
% [s4_10,s4_20,s4_30]=textread('w-10-20.lvm','%s%s%s');
% row_num_first=100;
% s4_10=str2double(s4_10(row_num_first:24000,1));
% s4_20=str2double(s4_20(row_num_first:24000,1));
% figure(2)
% plot(s4_10)
% xlabel('simple');ylabel('amplitude');
% figure(3)
% plot(s4_20)
% xlabel('simple');ylabel('amplitude');
%
% Y=fft(s4_0,N);
% Y1=fftshift(Y);
% f=(-N/2+1:N/2)*Fs/N;
% Mag=abs(Y1);
% figure(4);
% plot(f,10*log10(Mag));
%
% Y_10=fft(s4_10,N);
% Y1_10=fftshift(Y_10);
% f=(-N/2+1:N/2)*Fs/N;
% Mag1=abs(Y1_10);
% figure(5);
% plot(f,10*log10(Mag1));
%
% Y_20=fft(s4_20,N);
% Y1_20=fftshift(Y_20);
% f=(-N/2+1:N/2)*Fs/N;
% Mag2=abs(Y1_20);
% figure(6);
% plot(f,10*log10(Mag2));
%% test
% [s4_0,s4_1,s4_2,s4_3]=textread('10.lvm','%s%s%s%s');
% row_num_first=31;
% s4_0=str2double(s4_0(row_num_first:end,1));
% s4_1=str2double(s4_1(row_num_first:end,1));
% s4_2=str2double(s4_2(row_num_first:end,1));
% s4_3=str2double(s4_3(row_num_first:end,1));
% [s4_0,s4_1,s4_2]=textread('10.lvm','%s%s%s');
% row_num_first=31;
% s4_0=str2double(s4_0(row_num_first:end,1));
% s4_1=str2double(s4_1(row_num_first:end,1));
% s4_2=str2double(s4_2(row_num_first:end,1));
%%
%
% [s_0,s_1]=textread('10_20.lvm','%s%s');
% row_num_first=33;
% s_0=str2double(s_0(row_num_first:24000,1));
% s_1=str2double(s_1(row_num_first:24000,1));
%
% N=256;
% Fs=25600;
% Y=fft(s_0,N);
% Y1=fftshift(Y);
% f=(-N/2+1:N/2)*Fs/N;
% Mag=abs(Y1);
% figure(1);
% plot(f,10*log10(Mag));
% Y_2=fft(s_1,N);
% Y1_2=fftshift(Y_2);
% f=(-N/2+1:N/2)*Fs/N;
% Mag2=abs(Y1_2);
% figure(2);
% plot(f,10*log10(Mag2));