% Ts=1/30e6;
% fid = fopen('C:\Users\Administrator\Desktop\AD_test\2017-5-27-1M-2Vpp.txt');
% test = fscanf(fid, '%s\n');
% tt = fi(0,1,16,15);
% x1 = zeros(1,length(test)/4);
% for i = 1:length(test)/4
% tt.hex = test((4*i-3):4*i);
% x1(1,i)=tt.data;
% end
%
% fs=30e6;%设定采样频率
% N=1024;
% n=0:N-1;
% x=x1(1:n); %n为付立叶变换的点数,是2的幂,一般取512或1024
%
% %***************1.正弦波****************%
%
% t=n/fs;
% figure(1);
% subplot(211);
% plot(1:length(x),x,'b');%作输入信号的时域波形
% %stem(x);
% hold on;
% xlabel('t');
% ylabel('x');
% title('输入信号x时域波形');
% grid on;
%
%
% %进行FFT变换并做频谱图
% y=fft(x,N);%进行fft变换
% mag=abs(y)*2/N;%求幅值
% f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
% figure(1);
% subplot(212);
% plot(f,mag,'b');%做频谱图
% hold on;
% xlabel('频率(Hz)');
% ylabel('幅值');
% title('输入信号x幅频谱图N=1024');
% grid;
a = load('C:\Users\Administrator\Desktop\AD_test\2017-5-27-1M-2Vpp-1.txt');%读数据
% t=a(:,1)';%时间
% x=a(:,2)';%数据
% %dt=t(2)-t(1);%用diff(t)看了下,你的采样时间间隔一样,就只算了一个差值
% fs=32e6;%设定采样频率%傅里叶变换后的坐标轴的最大值就是采样频率fs,与dt是倒数关系
% ax=abs(fft(x));%fft是matlab中傅里叶变换的命令,abs是求模,我画的是幅值谱
% f=(0:length(t)-1)*fs/length(t);%频谱图横轴,总共49个点,最大值是fs,每个点就是fs/49
% figure();%画在一起了
% subplot(211);plot(t,x);%上:时域图
% subplot(212);plot(f(1:25),ax(1:25));%下