欧美极品高清xxxxhd,国产日产欧美最新,无码AV国产东京热AV无码,国产精品人与动性XXX,国产传媒亚洲综合一区二区,四库影院永久国产精品,毛片免费免费高清视频,福利所导航夜趣136

標(biāo)題: Matlab串口波形刷新和數(shù)據(jù)分析m文件 [打印本頁(yè)]
作者: CHENYIXIAO    時(shí)間: 2021-9-25 11:58
標(biāo)題: Matlab串口波形刷新和數(shù)據(jù)分析m文件
V7-203_Matlab串口波形刷新和數(shù)據(jù)分析m文件
  1. % ********************************************************************************************************************
  2. % ********************************************************************************************************************
  3. close all
  4. clear all

  6. %刪除所有已經(jīng)打開(kāi)的串口，這條很重要，防止之前運(yùn)行沒(méi)有關(guān)閉串口
  7. delete(instrfindall);  

  9. %打開(kāi)串口COM1，波特率115200，8位數(shù)據(jù)位，1位停止位，無(wú)奇偶校驗(yàn)，無(wú)流控制
  10. s = serial('COM1', 'BaudRate', 115200, 'DataBits', 8, 'StopBits', 1, 'Parity', 'none', 'FlowControl', 'none');
  11. s.ReadAsyncMode = 'continuous';
  12. fopen(s);

  14. fig = figure(1);

  16. % ********************************************************************************************************************
  17. AxisMax =  65536;    %坐標(biāo)軸最大值
  18. AxisMin = -65536;    %坐標(biāo)軸最小值
  19. window_width = 800;  %窗口寬度

  21. g_Count =0;          %接收到的數(shù)據(jù)計(jì)數(shù)
  22. SOF = 0;             %同步幀標(biāo)志
  23. AxisValue = 1;       %坐標(biāo)值
  24. RecDataDisp = zeros(1,100000); %開(kāi)辟100000個(gè)數(shù)據(jù)單元，用于存儲(chǔ)接收到的數(shù)據(jù)。
  25. RecData = zeros(1,100);        %開(kāi)辟100個(gè)數(shù)據(jù)單元，用于數(shù)據(jù)處理。
  26. Axis = zeros(1,100000);        %開(kāi)辟100000個(gè)數(shù)據(jù)單元，用于X軸。

  28. window = window_width * (-0.9); %窗口X軸起始坐標(biāo)
  29. axis([window, window + window_width, AxisMin, AxisMax]); %設(shè)置窗口坐標(biāo)范圍

  31. %子圖1顯示串口上傳的數(shù)據(jù)
  32. subplot(2,1,1);
  33. grid on;
  34. title('串口數(shù)據(jù)接收');
  35. xlabel('時(shí)間');
  36. ylabel('數(shù)據(jù)');

  38. %子圖2顯示波形的幅頻響應(yīng)
  39. subplot(2,1,2);
  40. grid on;
  41. title( 'FFT');
  42. xlabel( '頻率');
  43. ylabel( '幅度');

  45. Fs = 100;        % 采樣率
  46. N = 50;         % 采樣點(diǎn)數(shù)
  47. n = 0:N-1;      % 采樣序列
  48. f = n * Fs / N; %真實(shí)的頻率

  50. % ********************************************************************************************************************

  52. while ishandle(fig)
  53.    
  54.     %設(shè)置同步信號(hào)標(biāo)志， = 1表示接收到下位機(jī)發(fā)送的同步幀
  55.     SOF = 0;  
  56.    
  57.     %發(fā)送同步幀
  58.     fwrite(s, 13);
  59.    
  60.     %獲取是否有數(shù)據(jù)
  61.     bytes = get(s, 'BytesAvailable');
  62.     if bytes == 0
  63.          bytes = 1;
  64.     end
  65.    
  66.     %讀取下位機(jī)返回的所有數(shù)據(jù)
  67.     RecData = fread(s, bytes, 'uint8');
  68.    
  69.     %檢索下位機(jī)返回的數(shù)據(jù)中是否有字符$
  70.     StartData = find(RecData == 13);
  71.    
  72.     %如果檢索到$，讀取10個(gè)字節(jié)的數(shù)據(jù)，也就是5個(gè)uint16的數(shù)據(jù)
  73.     if(StartData >= 1)
  74.         RecData = fread(s, 5, 'uint16');
  75.         SOF =1;
  76.         StartData = 0;
  77.     end
  78.    
  79. %更新接收到的數(shù)據(jù)波形
  80. if(SOF == 1)
  81.         %更新數(shù)據(jù)
  82.         RecDataDisp(AxisValue) =  RecData(1);
  83.         RecDataDisp(AxisValue + 1) =  RecData(2);
  84.         RecDataDisp(AxisValue + 2) =  RecData(3);
  85.         RecDataDisp(AxisValue + 3) =  RecData(4);
  86.         RecDataDisp(AxisValue + 4) =  RecData(5);
  87.        
  88.         %更新X軸
  89.         Axis(AxisValue) = AxisValue;
  90.         Axis(AxisValue + 1) = AxisValue + 1;
  91.         Axis(AxisValue + 2) = AxisValue + 2;
  92.         Axis(AxisValue + 3) = AxisValue + 3;
  93.         Axis(AxisValue + 4) = AxisValue + 4;

  95.         %更新變量
  96.         AxisValue = AxisValue + 5;
  97.         g_Count = g_Count + 5;
  98.        
  99.         %繪制波形
  100.         subplot(2,1,1);
  101.         plot(Axis(1:AxisValue-1),  RecDataDisp(1:AxisValue-1), 'r');
  102.         window = window + 5;
  103.         axis([window, window + window_width, AxisMin, AxisMax]);
  104.         grid on;
  105.         title('串口數(shù)據(jù)接收');
  106.         xlabel('時(shí)間');
  107.         ylabel('數(shù)據(jù)');
  108.         drawnow
  109. end
  110.    
  111. if(g_Count== 50)
  112.    subplot(2,1,2);
  113.    %對(duì)原始信號(hào)做 FFT 變換
  114.    y = fft(RecDataDisp(AxisValue-50:AxisValue-1), 50);
  115.    
  116.    %求 FFT 轉(zhuǎn)換結(jié)果的模值
  117.    Mag = abs(y)*2/N;  
  118.    
  119.    %繪制幅頻相應(yīng)曲線(xiàn)
  120.    plot(f, Mag, 'r');
  121.    grid on;
  122.    title( 'FFT');
  123.    xlabel( '頻率');
  124.    ylabel( '幅度');
  125.    g_Count = 0;
  126.    drawnow
  127. end

  129. end

  131. fclose(s);
  132. delete(s);

  134. % ********************************************************************************************************************
復(fù)制代碼






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