標題: STM32F103ZET6進行AD采集并通過串口將數(shù)據(jù)傳至MATLAB動態(tài)接收并進行處理 [打印本頁]
作者: HLHLHLHLHLHL    時間: 2020-2-20 11:07
標題: STM32F103ZET6進行AD采集并通過串口將數(shù)據(jù)傳至MATLAB動態(tài)接收并進行處理
串口接收程序：function Serial()   %創(chuàng)建函數(shù)
delete(instrfindall);   %先關閉串口，否則可能導致出錯

global x            %全局變量，供串口中斷函數(shù)使用

global t;                        %全局變量，這里根據(jù)需要繪制圖形的個人需要而設
global m;
global i;

t = [0];      %時間軸
m = [0];      %數(shù)據(jù)軸
i = 0;        %用于計數(shù)

p = plot(t,m);
axis([i-50 i+50 -1 10]);
grid on;

%%----------------------------配置部分----------------------------

%創(chuàng)建串口對象，直接在命令框輸入serial('com1')可查看串口參數(shù)
x = serial('com4');
%設置或修改串口通信參數(shù)
set(x,'BaudRate',115200);

%設置中斷觸發(fā)方式
set(x,'BytesAvailableFcnMode','Terminator') %ASCII觸發(fā)，字符觸發(fā)
set(x,'Terminator','CR/LF')         %接收到\r\n后觸發(fā)中斷

% set(x,'BytesAvailableFcnMode','byte')     %數(shù)據(jù)觸發(fā)，根據(jù)發(fā)送來的數(shù)據(jù)長度來觸發(fā)中斷
% set(x,'BytesAvailableFcnCount',1)


x.BytesAvailableFcn = @Callback  %定義中斷響應函數(shù)對象，類似于中斷函數(shù)名

fopen(x);                %打開串口，類似于開關

fwrite(x,255);           %發(fā)送握手信號0xff，stm32接收到后開始回傳數(shù)據(jù)
pause                 %需要暫停，來接收回傳數(shù)據(jù)

fclose(x);

end


串口中斷服務子程序：
function Callback(obj,event)   %創(chuàng)建中斷服務函數(shù)，繪制圖像

global t; %時間
global m; %縱坐標
global i; %時間變化值
%
if i<3000
    out = fscanf(obj);
    data = str2num(out)        %將接收到的字符轉(zhuǎn)換為數(shù)值
    save('a.mat','data','-ascii','-append');

%     x=load('a.mat','-ascii');

%       out=fread(obj,5,'uint32');
%       data = str2num(out)
     %%----------------------以下根據(jù)需要自行編寫-------------------------------
    t = [t i];                                                       
    m = [m data];
    plot(t,m)
    xlabel('t');
    ylabel('data');
    axis([i-50 i+50 0 3.3]);
    grid on;

    i=i+1;
end

end


對數(shù)據(jù)進行FFT處理：
% clear all
% close all
% clc
%% Parameter Interface
SampleFre = 1000; %采樣頻率
SignalData=load('a.mat','-ascii');
% SignalLen = 20000;  %數(shù)據(jù)長度
SignalLen=length(SignalData);
%% Main
%--------------------------------------------------------------------------------------
t = (0:SignalLen-1)/SampleFre; %時間
% SignalData = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); %采樣頻率>2倍信號最大頻率
NoiseData = 2*randn(size(t));
SignalAddNoiseData = SignalData + NoiseData;
plot(1000*t(1:50),SignalAddNoiseData(1:50))
title('信號添加隨機噪聲')
xlabel('t(毫秒)')
%--------------------------------------------------------------------------------------
FFTData = fft(SignalAddNoiseData);
FFTAmplitude0 = abs(FFTData)/SignalLen;  %頻譜關于SignalLen/2對稱
FFTAmplitude1 = FFTAmplitude0(1:SignalLen/2);
FFTAmplitude0(2:end) = 2*FFTAmplitude0(2:end); %復數(shù)域的幅值需要進行變換才能與時域的幅值對應
FFTAmplitude1(2:end) = 2*FFTAmplitude1(2:end);
Frequence0 = SampleFre*(0:(SignalLen-1))/SignalLen;
Frequence1 = SampleFre*((0:(SignalLen/2)-1))/SignalLen;
subplot(2,1,1)
plot(Frequence0,FFTAmplitude0)
title('雙邊頻譜')
xlabel('f(Hz)')
ylabel('FFTAmplitude')
subplot(2,1,2)
plot(Frequence1,FFTAmplitude1)
title('單邊頻譜')
xlabel('f(Hz)')
ylabel('FFTAmplitude')


mat.zip

11.55 KB, 下載次數(shù): 46, 下載積分: 黑幣 -5

內(nèi)含三個matlab子程序



歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1