色哟哟视频在线观看-色哟哟视频在线-色哟哟欧美15最新在线-色哟哟免费在线观看-国产l精品国产亚洲区在线观看-国产l精品国产亚洲区久久

0
  • 聊天消息
  • 系統消息
  • 評論與回復
登錄后你可以
  • 下載海量資料
  • 學習在線課程
  • 觀看技術視頻
  • 寫文章/發帖/加入社區
會員中心
創作中心

完善資料讓更多小伙伴認識你,還能領取20積分哦,立即完善>

3天內不再提示

基于 Infineon PSoC?62開發板的信號處理前端虛擬示波器-工具集

RTThread物聯網操作系統 ? 來源:未知 ? 2023-08-16 20:15 ? 次閱讀
一、前言
本項目基于英飛凌PSoC6 RT-Thread開發板實現了信號處理前端-一個信號處理的工具集。包括虛擬示波器音頻采集分析,諧波分析,周期幅值相位分析,數字濾波,極值檢測,可上位機可視化和命令行人機交互,可以方便繼續擴展相關功能,繼續豐富工具集。wKgZomToRnOAesyNAAExUFSjJj4041.png視頻:https://www.bilibili.com/video/BV1PM4y147v1/代碼倉庫:https://gitee.com/qinyunti/infineon-psoc62.git 二.移植DSP算法
2.1添加代碼git clone https://github.com/ARM-software/CMSIS_5.gitCMSIS_5CMSISDSP下是相關文件,Source下是源碼wKgZomToRnOAKrdXAAA585KWFls161.jpg將DSP文件夾復制到自己的工程目錄中,只保留Include,PrivateInclude,Source三個文件夾wKgZomToRnOAcZdOAAAgt3v1Iw4368.jpgwKgZomToRnOAcZdOAAAgt3v1Iw4368.jpgSource下的每個子文件夾都是一類算法,里面的每個c都對應一個計算函數,并且有一個總文件包括其中所有的單個.c,比如BasicMathFunctions.c中wKgZomToRnOAHpLeAAByEZp5dEo698.jpg刪除這些總的.c,避免編譯重復刪除以下文件和所有的非.c和.h文件

		1BasicMathFunctions:BasicMathFunctions.c,BasicMathFunctionsF16.c  2BayesFunctions:BayesFunctions.c,BayesFunctionsF16.c  3CommonTables:CommonTables.c,CommonTablesF16.c  4ComplexMathFunctions:ComplexMathFunctions.c,ComplexMathFunctionsF16.c  5ControllerFunctions:ControllerFunctions.c  6DistanceFunctions:DistanceFunctions.c,DistanceFunctionsF16.c  7FastMathFunctions:FastMathFunctions.c,FastMathFunctionsF16.c  8FilteringFunctions:FilteringFunctions.c,FilteringFunctionsF16.c  9InterpolationFunctions:InterpolationFunctions.c,InterpolationFunctionsF16.c 10MatrixFunctions:MatrixFunctions.c,MatrixFunctionsF16.c 11QuaternionMathFunctions:QuaternionMathFunctions.c 12StatisticsFunctions:StatisticsFunctions.c,StatisticsFunctionsF16.c 13SupportFunctions:SupportFunctions.c,SupportFunctionsF16.c 14SVMFunctions:SVMFunctions.c,SVMFunctionsF16.c 15TransformFunctions:TransformFunctions.c,TransformFunctionsF16.c,arm_bitreversal2.S wKgZomToRnOAaA6vAACg8Vcl9Vs015.jpg工程設置添加相關頭文件包含路徑
		2.2測試復制CMSIS_5CMSISDSPExamplesARMarm_fft_bin_example下的arm_fft_bin_data.carm_fft_bin_example_f32.c到自己的工程目錄arm_fft_bin_example_f32.c下的int32_t main(void)改為int32_t ffttest_main(void)并添加#define SEMIHOSTING以使能printf打印,我們已經重定向實現了printf打印到串口。由于 arm_cfft_f32(&varInstCfftF32, testInput_f32_10khz, ifftFlag, doBitReverse);會修改testInput_f32_10khz的內容,所以添加一個緩存,以便能重復測試

		1float32_ttesttmp_f32_10khz[2048]; 2/*ProcessthedatathroughtheCFFT/CIFFTmodule*/ 3memcpy(testtmp_f32_10khz,testInput_f32_10khz,sizeof(testInput_f32_10khz)); 4arm_cfft_f32(&varInstCfftF32,testtmp_f32_10khz,ifftFlag,doBitReverse); 5/*ProcessthedatathroughtheComplexMagnitudeModulefor 6calculatingthemagnitudeateachbin*/ 7arm_cmplx_mag_f32(testtmp_f32_10khz,testOutput,fftSize); 在自己的main函數中申明并調用int32_t ffttest_main(void);ffttest_main();編譯運行可以看到串口打印SUCCESS說明測試OK。將輸入輸出數據打印

		1printf("SUCCESS ");  2for(inti=0;i 3{  4if(i2)  5{  6printf("/*%f,%f*/ ",testInput_f32_10khz[i],testOutput[i]);  7}  8else  9{ 10printf("/*%f,%f*/ ",testInput_f32_10khz[i],0.0); 11} 12} 使用serialstudio可視化顯示,可以看到計算結果FFT頻率明顯的峰值

		wKgZomToRnSAOJOBAACPjIRdtsg761.jpg

		
			三.音頻采集
3.1原理圖從原理圖看到有6路模擬輸入,分別對應P10.0~P10.5, VREF為模擬參考電壓。wKgZomToRnSAL2KUAABfT2wQacg834.jpgwKgZomToRnSANybeAACo7SSVwkM417.jpg使用的是MAX4466的MIC,接到ADC0,如下圖所示wKgZomToRnSAab26AADL_PiUxb4588.jpgwKgZomToRnSADbj2AABP72Xp4zw570.jpgwKgZomToRnSAX9qhAACnU2Fye04562.jpg ?3.2配置模擬采集引腳 wKgZomToRnSAIKiLAACQYqRTOsc532.jpg ?3.3代碼Adc.c

			1#include"cy_pdl.h"  2#include"cyhal.h"  3#include"cybsp.h"  4#include"cy_retarget_io.h"  5#defineVPLUS_CHANNEL_0(P10_0)  6/*Conversionfactor*/  7#defineMICRO_TO_MILLI_CONV_RATIO(1000u)  8/*Acquistiontimeinnanosecond*/  9#defineACQUISITION_TIME_NS(116680u)  10/*ADCScandelayinmillisecond*/  11#defineADC_SCAN_DELAY_MS(200u)  12/*******************************************************************************  13*EnumeratedTypes  14*******************************************************************************/  15/*ADCChannelconstants*/  16enumADC_CHANNELS  17{  18CHANNEL_0=0,  19NUM_CHANNELS  20}adc_channel;  21/*******************************************************************************  22*GlobalVariables  23*******************************************************************************/  24/*ADCObject*/  25cyhal_adc_tadc_obj;  26/*ADCChannel0Object*/  27cyhal_adc_channel_tadc_chan_0_obj;  28/*DefaultADCconfiguration*/  29constcyhal_adc_config_tadc_config={  30.continuous_scanning=false,//ContinuousScanningisdisabled  31.average_count=1,//Averagecountdisabled  32.vref=CYHAL_ADC_REF_VDDA,//VREFforSingleendedchannelsettoVDDA  33.vneg=CYHAL_ADC_VNEG_VSSA,//VNEGforSingleendedchannelsettoVSSA  34.resolution=12u,//12-bitresolution  35.ext_vref=NC,//Noconnection  36.bypass_pin=NC};//Noconnection  37/*Asynchronousreadcompleteflag,usedinEventHandler*/  38staticboolasync_read_complete=true;  39#defineNUM_SCAN(1000)  40#defineNUM_CHANNELS(1)  41/*Variabletostoreresultsfrommultiplechannelsduringasynchronousread*/  42int32_tresult_arr[NUM_CHANNELS*NUM_SCAN]={0};  43staticvoidadc_event_handler(void*arg,cyhal_adc_event_tevent)  44{  45if(0u!=(event&CYHAL_ADC_ASYNC_READ_COMPLETE))  46{  47/*Setasyncreadcompleteflagtotrue*/  48async_read_complete=true;  49}  50}  51intadc_init(void)  52{  53/*Variabletocapturereturnvalueoffunctions*/  54cy_rslt_tresult;  55/*InitializeADC.TheADCblockwhichcanconnecttothechannel0inputpinisselected*/  56result=cyhal_adc_init(&adc_obj,VPLUS_CHANNEL_0,NULL);  57if(result!=CY_RSLT_SUCCESS)  58{  59printf("ADCinitializationfailed.Error:%ld ",(longunsignedint)result);  60CY_ASSERT(0);  61}  62/*ADCchannelconfiguration*/  63constcyhal_adc_channel_config_tchannel_config={  64.enable_averaging=false,//Disableaveragingforchannel  65.min_acquisition_ns=ACQUISITION_TIME_NS,//Minimumacquisitiontimesetto1us  66.enabled=true};//SamplethischannelwhenADCperformsascan  67/*Initializeachannel0andconfigureittoscanthechannel0inputpininsingleendedmode.*/  68result=cyhal_adc_channel_init_diff(&adc_chan_0_obj,&adc_obj,VPLUS_CHANNEL_0,  69CYHAL_ADC_VNEG,&channel_config);  70if(result!=CY_RSLT_SUCCESS)  71{  72printf("ADCfirstchannelinitializationfailed.Error:%ld ",(longunsignedint)result);  73CY_ASSERT(0);  74}  75/*Registeracallbacktohandleasynchronousreadcompletion*/  76cyhal_adc_register_callback(&adc_obj,&adc_event_handler,result_arr);  77/*Subscribetotheasyncreadcompleteeventtoprocesstheresults*/  78cyhal_adc_enable_event(&adc_obj,CYHAL_ADC_ASYNC_READ_COMPLETE,CYHAL_ISR_PRIORITY_DEFAULT,true);  79printf("ADCisconfiguredinmultichannelconfiguration. ");  80printf("Channel0isconfiguredinsingleendedmode,connectedtothe ");  81printf("channel0inputpin.Provideinputvoltageatthechannel0inputpin ");  82return0;  83}  84intadc_samp(void)  85{  86/*Variabletocapturereturnvalueoffunctions*/  87cy_rslt_tresult;  88/*VariabletostoreADCconversionresultfromchannel0*/  89int32_tadc_result_0=0;  90/*Clearasyncreadcompleteflag*/  91async_read_complete=false;  92/*Initiateanasynchronousreadoperation.Theeventhandlerwillbecalled  93*whenitiscomplete.*/  94memset(result_arr,0,sizeof(result_arr));  95cyhal_gpio_write_internal(CYBSP_USER_LED,true);  96result=cyhal_adc_read_async_uv(&adc_obj,NUM_SCAN,result_arr);  97if(result!=CY_RSLT_SUCCESS)  98{  99printf("ADCasyncreadfailed.Error:%ld ",(longunsignedint)result); 100CY_ASSERT(0); 101} 102while(async_read_complete==false); 103cyhal_gpio_write_internal(CYBSP_USER_LED,false); 104/* 105*Readdatafromresultlist,inputvoltageintheresultlistisin 106*microvolts.Convertitmillivoltsandprintinputvoltage 107* 108*/ 109for(inti=0;i110{ 111adc_result_0=result_arr[i]/MICRO_TO_MILLI_CONV_RATIO; 112printf("/*%4ld*/ ",(longint)adc_result_0); 113} 114return0; 115}

		Adc.h

			1#ifndefADC_H 2#defineADC_H 3intadc_init(void); 4intadc_samp(void); 5#endif

		Main.c調用adc_init();adc_samp();
			3.4時鐘時鐘源是100Mhz,12分頻=8.33M,滿足1.8MHz~18MHz之間的要求默認是按照8M配置

		wKgZomToRnSABHRuAACJsoLRvJc093.jpg

		wKgZomToRnWAU7B7AABiJ-TBAgw773.jpg
			?3.5采樣時間采樣前后翻轉LED用示波器測量時間

			1intadc_samp(void)  2{  3/*Variabletocapturereturnvalueoffunctions*/  4cy_rslt_tresult;  5/*VariabletostoreADCconversionresultfromchannel0*/  6int32_tadc_result_0=0;  7/*Clearasyncreadcompleteflag*/  8async_read_complete=false;  9/*Initiateanasynchronousreadoperation.Theeventhandlerwillbecalled 10*whenitiscomplete.*/ 11memset(result_arr,0,sizeof(result_arr)); 12cyhal_gpio_write_internal(CYBSP_USER_LED,true); 13result=cyhal_adc_read_async_uv(&adc_obj,NUM_SCAN,result_arr); 14if(result!=CY_RSLT_SUCCESS) 15{ 16printf("ADCasyncreadfailed.Error:%ld ",(longunsignedint)result); 17CY_ASSERT(0); 18} 19while(async_read_complete==false); 20cyhal_gpio_write_internal(CYBSP_USER_LED,false); 21/* 22*Readdatafromresultlist,inputvoltageintheresultlistisin 23*microvolts.Convertitmillivoltsandprintinputvoltage 24* 25*/ 26for(inti=0;i27{ 28adc_result_0=result_arr[i]/MICRO_TO_MILLI_CONV_RATIO; 29printf("/*%4ld*/ ",(longint)adc_result_0); 30} 31return0; 32}

		采樣1000次,分別設置采樣時間為2uS和1uS對比。#define ACQUISITION_TIME_NS (2000u)10.28mS

		wKgZomToRnWAfAMaAABkvoWNWT0667.jpg#define ACQUISITION_TIME_NS (1000u)9.32mS

		wKgZomToRnWAL1UTAABskInWyZ8332.jpg10.28-9.32=0.96mS 1000次約1mS,1次剛好是1uS。而1000次除去采樣時間其他時間為8.32mS,即一次8.32uS。因為前面設置了時鐘為8.33MHz, 從前面時序一節可以看到,除去采樣時間,其他轉換時間等需要14個CLK,所以需要14/8.33uS=1.7uS. 剩余的8.32-1.7為數據搬運,軟件處理等時間。
			3.6 采樣值正確性1.545V和示波器采集為1.54V差不多是正確的,這里沒有高精度萬用表就不對測試精度了,只測試了正確性。wKgZomToRnWAL3ByAAB1jy7ZOdE552.jpg

		wKgZomToRnWASVRqAABdXRAQons955.jpg

		
			3.7音頻采集一次采集1000次然后串口打印,使用SerialStudio可視化顯示

			1intadc_samp(void)  2{  3/*Variabletocapturereturnvalueoffunctions*/  4cy_rslt_tresult;  5/*VariabletostoreADCconversionresultfromchannel0*/  6int32_tadc_result_0=0;  7/*Clearasyncreadcompleteflag*/  8async_read_complete=false;  9/*Initiateanasynchronousreadoperation.Theeventhandlerwillbecalled 10*whenitiscomplete.*/ 11memset(result_arr,0,sizeof(result_arr)); 12cyhal_gpio_write_internal(CYBSP_USER_LED,true); 13result=cyhal_adc_read_async_uv(&adc_obj,NUM_SCAN,result_arr); 14if(result!=CY_RSLT_SUCCESS) 15{ 16printf("ADCasyncreadfailed.Error:%ld ",(longunsignedint)result); 17CY_ASSERT(0); 18} 19while(async_read_complete==false); 20cyhal_gpio_write_internal(CYBSP_USER_LED,false); 21/* 22*Readdatafromresultlist,inputvoltageintheresultlistisin 23*microvolts.Convertitmillivoltsandprintinputvoltage 24* 25*/ 26for(inti=0;i27{ 28adc_result_0=result_arr[i]/MICRO_TO_MILLI_CONV_RATIO; 29printf("/*%4ld*/ ",(longint)adc_result_0); 30} 31return0; 32}

		
			四.信號處理前端
4.1 電能質量,諧波分析4.1.1添加命令行在電能檢測應用中,電能質量一項分析即諧波分析,諧波分量大,說明電能質量不好,基于本板信號處理前端也實現了該功能。shell_fun.h

			1voidFftFun(void*param); shell_fun.c

			1include"fft.h" shell_cmd_list中添加一行

			1{(constuint8_t*)"fft",FftFun,"fft"},/*打印幫助信息*/ 添加命令執行函數

			1voidFftFun(void*param) 2{ 3fft_main(); 4}
			4.1.2添加實現Fft.c

			1#include"arm_math.h"  2#include"arm_const_structs.h"  3#include  4#defineTEST_LENGTH_SAMPLES2048  5externfloat32_ttestInput_f32_10khz[TEST_LENGTH_SAMPLES];  6staticfloat32_ttestOutput[TEST_LENGTH_SAMPLES/2];  7staticuint32_tfftSize=1024;  8staticuint32_tifftFlag=0;  9staticuint32_tdoBitReverse=1; 10staticarm_cfft_instance_f32varInstCfftF32; 11staticinttestIndex=0; 12staticfloattesttmp_f32_10khz[2048]; 13staticint32_tadcbuffer[2048]; 14int32_tfft_main(void) 15{ 16arm_statusstatus; 17float32_tmaxValue; 18status=ARM_MATH_SUCCESS; 19status=arm_cfft_init_f32(&varInstCfftF32,fftSize); 20//memcpy(testtmp_f32_10khz,testInput_f32_10khz,sizeof(testInput_f32_10khz)); 21adc_samp(adcbuffer,2048); 22for(inti=0;i<2048;i++) 23{ 24testtmp_f32_10khz[i]=(float)adcbuffer[i]; 25} 26arm_cfft_f32(&varInstCfftF32,testtmp_f32_10khz,ifftFlag,doBitReverse); 27arm_cmplx_mag_f32(testtmp_f32_10khz,testOutput,fftSize); 28/*CalculatesmaxValueandreturnscorrespondingBINvalue*/ 29arm_max_f32(testOutput,fftSize,&maxValue,&testIndex); 30int32_tout=0; 31for(inti=0;i32{ 33if(i>TEST_LENGTH_SAMPLES/2) 34{ 35out=testOutput[i-TEST_LENGTH_SAMPLES/2]/1024; 36} 37else 38{ 39out=testOutput[i]/1024; 40} 41printf("/*%ld,%ld*/ ",adcbuffer[i],out); 42} 43} 44/**endlink*/Fft.h

			1#ifndefFFT_H 2#defineFFT_H 3intfft_main(void); 4#endif

		測試

		wKgZomToRnWARPWbAAA4rRk9BL8816.jpg

		
			4.2 周期(頻率),幅值,相位分析4.2.1 原理FFT變換結果,幅值最大的橫坐標對應信號頻率,縱坐標對應幅度。幅值最大的為out[m]=val;則信號頻率f0=(Fs/N)m ,信號幅值Vpp=val/(N/2)。N為FFT的點數,Fs為采樣頻率。相位Pha=atan2(a, b)弧度制,其中ab是輸出虛數結果的實部和虛部。
			4.2.2添加命令行shell_fun.h

			1voidFrqFun(void*param); shell_fun.cinclude "frq.h"shell_cmd_list中添加一行

			1{(constuint8_t*)"frt",FrqFun,"frq"}, 添加命令執行函數

			1voidFrqFun(void*param) 2{ 3Frq_main(); 4}

		
			4.2.3實現代碼Frq.c

			1#include"arm_math.h"  2#include"arm_const_structs.h"  3#include  4#defineTEST_LENGTH_SAMPLES2048  5#defineFS10000  6externfloat32_ttestInput_f32_10khz[TEST_LENGTH_SAMPLES];  7staticfloat32_ttestOutput[TEST_LENGTH_SAMPLES/2];  8staticuint32_tfftSize=1024;  9staticuint32_tifftFlag=0; 10staticuint32_tdoBitReverse=1; 11staticarm_cfft_instance_f32varInstCfftF32; 12staticinttestIndex=0; 13staticfloattesttmp_f32_10khz[2048]; 14staticint32_tadcbuffer[2048]; 15int32_tfrq_main(void) 16{ 17arm_statusstatus; 18float32_tmaxValue; 19status=ARM_MATH_SUCCESS; 20status=arm_cfft_init_f32(&varInstCfftF32,fftSize); 21//memcpy(testtmp_f32_10khz,testInput_f32_10khz,sizeof(testInput_f32_10khz)); 22adc_samp(adcbuffer,2048); 23for(inti=0;i<2048;i++) 24{ 25testtmp_f32_10khz[i]=(float)adcbuffer[i]; 26} 27arm_cfft_f32(&varInstCfftF32,testtmp_f32_10khz,ifftFlag,doBitReverse); 28arm_cmplx_mag_f32(testtmp_f32_10khz,testOutput,fftSize); 29/*CalculatesmaxValueandreturnscorrespondingBINvalue*/ 30arm_max_f32(testOutput,fftSize,&maxValue,&testIndex); 31floatfreq=(FS/TEST_LENGTH_SAMPLES)*testIndex; 32floatvpp=maxValue/(TEST_LENGTH_SAMPLES/2); 33floatpha=atan2(testOutput[2*testIndex],testOutput[2*testIndex+1]); 34printf("freq=%f,vpp=%f,pha=%f ",freq,vpp,pha); 35} 36/**endlink*/Frq.h

			1#ifndefFRQ_H 2#defineFRQ_H 3intfrq_main(void); 4#endif

		
			4.2.4測試輸入frq開始測試印如下wKgZomToRnWAN7ArAAB53lnTnBg135.jpg實時采集測試此時采集的是音頻背景聲,噪聲很小,所以頻率為0

		wKgZomToRnWAYhSBAAB_1jG6J70546.jpg

		
			4.3數字濾波信號前端4.3.1原理CMSIS-DSP提供直接I型IIR庫支持Q7,Q15,Q31和浮點四種數據類型。其中Q15和Q31提供了快速版本。直接I型IIR濾波器是基于二階Biquad級聯的方式來實現的。每個Biquad由一個二階的濾波器組成:y[n] = b0 x[n] + b1 x[n-1] + b2 x[n-2] + a1 y[n-1] + a2 * y[n-2]直接I型算法每個階段需要5個系數和4個狀態變量。wKgZomToRnaARkLfAABBdSBW9wY673.jpgy[n] = b0 x[n] + b1 x[n-1] + b2 x[n-2] - a1 y[n-1] - a2 * y[n-2]matlab使用上面的公式實現,在使用fdatool工具箱生成的a系數需要取反才能用于直接I型IIR濾波器的函數中。高階IIR濾波器的實現是采用二階Biquad級聯的方式來實現的。其中參數numStages就是用來做指定二階Biquad的個數。比如8階IIR濾波器就可以采用numStages=4個二階Biquad來實現。wKgZomToRnaAZPpKAAA8t1K_CcI553.jpg如果要實現9階IIR濾波器就需要將numStages=5,這時就需要其中一個Biquad配置成一階濾波器(也就是b2=0,a2=0)。
			4.3.2添加命令行shell_fun.hvoid IirFun(void* param);shell_fun.c

			1include"iir.h" shell_cmd_list中添加一行

			1{(constuint8_t*)"iir",IirFun,“iir"}, 

		添加命令執行函數 

			1voidIirFun(void*param) 2{ 3Iir_main(); 4}
			4.3.3實現代碼Iir.c

			1#include"arm_math.h"  2#include"arm_const_structs.h"  3#include  4#defineTEST_LENGTH_SAMPLES2048  5#defineFS10000  6externfloat32_ttestInput_f32_10khz[TEST_LENGTH_SAMPLES];  7staticfloat32_ttestOutput[TEST_LENGTH_SAMPLES];  8staticuint32_tfftSize=1024;  9staticuint32_tifftFlag=0; 10staticuint32_tdoBitReverse=1; 11staticarm_cfft_instance_f32varInstCfftF32; 12staticinttestIndex=0; 13staticfloattesttmp_f32_10khz[2048]; 14staticint32_tadcbuffer[2048]; 15#definenumStages2/*2階IIR濾波的個數*/ 16#defineBLOCK_SIZE128/*調用一次arm_biquad_cascade_df1_f32處理的采樣點個數*/ 17uint32_tblockSize=BLOCK_SIZE; 18uint32_tnumBlocks=TEST_LENGTH_SAMPLES/BLOCK_SIZE;/*需要調用arm_biquad_cascade_df1_f32的次數*/ 19staticfloat32_tIIRStateF32[4*numStages];/*狀態緩存*/ 20/*巴特沃斯低通濾波器系數80Hz*/ 21constfloat32_tIIRCoeffs32LP[5*numStages]={ 221.0f,2.0f,1.0f,1.479798894397216679763573665695730596781f, 23-0.688676953053861784503908438637154176831f, 241.0f,2.0f,1.0f,1.212812092620218384908525877108331769705f, 25-0.384004162286553540894828984164632856846f 26}; 27int32_tiir_main(void) 28{ 29uint32_ti; 30arm_biquad_casd_df1_inst_f32S; 31float32_tScaleValue; 32float32_t*inputF32,*outputF32; 33/*初始化輸入輸出緩存指針*/ 34//memcpy(testtmp_f32_10khz,testInput_f32_10khz,sizeof(testInput_f32_10khz)); 35#if1 36adc_samp(adcbuffer,2048); 37for(inti=0;i<2048;i++) 38{ 39testtmp_f32_10khz[i]=(float)adcbuffer[i]; 40} 41#endif 42inputF32=testtmp_f32_10khz; 43outputF32=testOutput; 44/*初始化*/ 45arm_biquad_cascade_df1_init_f32(&S,numStages,(float32_t*)&IIRCoeffs32LP[0], 46(float32_t*)&IIRStateF32[0]); 47/*實現IIR濾波,這里每次處理1個點*/ 48for(i=0;i49{ 50arm_biquad_cascade_df1_f32(&S,inputF32+(i*blockSize),outputF32+(i*blockSize), 51blockSize); 52} 53/*放縮系數*/ 54ScaleValue=0.012f*0.42f; 55/*打印濾波后結果*/ 56for(i=0;i57{ 58printf("/*%f,%f*/ ",testtmp_f32_10khz[i],testOutput[i]*ScaleValue); 59} 60} 61/**endlink*/Iir.h

			1#ifndefIIR_H 2#defineIIR_H 3intiir_main(void); 4#endif

		
			4.3.4測試輸入iir回車,查看波形見視頻以下可以看到濾波導致了滯后,黃色線有滯后

		wKgZomToRnaACxvXAABU0wp83P0068.jpg

		wKgZomToRnaAQ-HLAAA7J-2D5i0599.jpg以下是實時采集濾波wKgZomToRnaAXvGOAABCxtthyV8336.jpg4.4 極大值檢測在電力等行業,分析電壓極值,是一項重要的參數分析,可以分析電壓的波動;示波器中也有自動測量極值的功能更。本板作為信號處理前端也實現了該功能。
			4.4.1 算法算法來源于論文https://www.mdpi.com/1999-4893/5/4/588/htm核心代碼如下

			1voidampd(int32_t*data,int32_tlen)  2{  3introw_sum;  4for(intk=1;k2+1;k++)  5{  6row_sum=0;  7for(inti=k;i 8{  9if((data[i]>data[i-k])&&(data[i]>data[i+k])) 10{ 11row_sum-=1; 12} 13} 14arr_rowsum[k-1]=row_sum; 15} 16intmin_index=argmin(arr_rowsum,len/2+1); 17max_window_length=min_index; 18for(intk=1;k1;k++) 19{ 20for(inti=k;i21{ 22if((data[i]>data[i-k])&&(data[i]>data[i+k])) 23{ 24p_data[i]+=1; 25} 26} 27} 28for(intk=0;k29{ 30if(p_data[k]==max_window_length) 31{ 32/*極大值*/ 33} 34} 35}4.4.2 添加命令行

			1{(constuint8_t*)"max",MaxFun,"max"},/*打印幫助信息*/ 2voidMaxFun(void*param) 3{ 4max_test(); 5} 6voidMaxFun(void*param);測試代碼如下,串口命令行輸入命令max,開始采集ADC值,并計算極值,打印到PC串口通過seraistudio可視化顯示

			1intmax_test(void)  2{  3for(inti=0;i<10;i++)  4{  5memset(p_data,0,sizeof(p_data));  6//adc_samp(sim_data_buffer,1000);  7sim_data();  8ampd(sim_data_buffer,sizeof(sim_data_buffer)/sizeof(sim_data_buffer[0]));  9for(intk=0;k<sizeof(sim_data_buffer)/sizeof(sim_data_buffer[0]);k++) 10{ 11if(p_data[k]==max_window_length) 12{ 13/*極大值*/ 14printf("/*%ld,%ld*/ ",sim_data_buffer[k],sim_data_buffer[k]); 15} 16else 17{ 18printf("/*%ld,%d*/ ",sim_data_buffer[k],0); 19} 20cyhal_system_delay_ms(10); 21} 22} 23return0; 24}

		
			4.4.3 測試效果如下IN是原始數據,MAX是檢測到的極大值,如果檢測極小值將原始數據取反即可。

		wKgZomToRnaAcTyBAABkghq9CX4885.jpgwKgZomToRnaAR-RrAABCqborvGI036.jpg檢測語音,效果如下wKgZomToRnaAdQO-AABPblN0HsA150.jpgwKgZomToRneAOs23AAA91IodeKs578.jpgwKgZomToRneAI4J5AAA_Ym_aWb8611.jpg
		?五.總結
得益于開發板出色的處理性能,和外設性能,以及官方可視化的代碼配置工具,可以方便的搭建開發環境,實現外設采集信號比如ADC,移植DSP庫,實現各種算法。本Demo實現了諧波分析,周期幅值相位分析,數字濾波,極大值檢測等功能,是一個小的工具集,還可以繼續擴展,設計了人機交互命令行,方便實用和測試,具備一定實用價值。
原文:https://club.rt-thread.org/ask/article/889e7c861bc93340.html
———————End——————
點擊閱讀原文進入官網

?

?


原文標題:基于 Infineon PSoC?62開發板的信號處理前端虛擬示波器-工具集

文章出處:【微信公眾號:RTThread物聯網操作系統】歡迎添加關注!文章轉載請注明出處。


聲明:本文內容及配圖由入駐作者撰寫或者入駐合作網站授權轉載。文章觀點僅代表作者本人,不代表電子發燒友網立場。文章及其配圖僅供工程師學習之用,如有內容侵權或者其他違規問題,請聯系本站處理。 舉報投訴
  • RT-Thread
    +關注

    關注

    31

    文章

    1305

    瀏覽量

    40363

原文標題:基于 Infineon PSoC?62開發板的信號處理前端虛擬示波器-工具集

文章出處:【微信號:RTThread,微信公眾號:RTThread物聯網操作系統】歡迎添加關注!文章轉載請注明出處。

收藏 人收藏

    評論

    相關推薦

    正點原子fpga開發板不同型號

    正點原子作為國內領先的FPGA開發板供應商,其產品線覆蓋了從入門級到高端應用的各個領域。這些開發板不僅適用于學術研究,還廣泛應用于工業控制、通信、圖像處理等多個領域。 1. 入門級開發板
    的頭像 發表于 11-13 09:30 ?1214次閱讀

    ARM開發板與FPGA的結合應用

    大的功能和更高的靈活性,滿足各種復雜應用的需求。 二、ARM開發板與FPGA的優勢 高性能:ARM處理器具有高性能的特點,可以處理復雜的算法和任務。FPGA則可以實現高速的數字信號
    的頭像 發表于 11-05 11:42 ?697次閱讀

    linux開發板如何編譯curl

    準備開發環境 確保你的開發板已經安裝了編譯 curl 所需的工具,包括 gcc 、 make 和可能的其他庫。如果開發板是交叉編譯環境(即你在一個平臺上編譯用于另一個平臺的程序),你
    的頭像 發表于 08-30 15:33 ?775次閱讀

    linux開發板和單片機開發的區別

    、PIC等,處理能力和內存容量相對較低。 操作系統 Linux開發板通常使用Linux操作系統,具有豐富的軟件資源和開發工具。單片機開發則通常使用嵌入式操作系統或無操作系統,軟件資源和
    的頭像 發表于 08-30 15:30 ?1164次閱讀

    虛擬示波器還是數字示波器好用

    。 工作原理 虛擬示波器通過與計算機連接的數據采集卡(DAQ)或其他接口設備,獲取外部信號,然后通過軟件進行處理和分析,最后在計算機屏幕上顯示波形。 特點 成本效益 :
    的頭像 發表于 07-17 17:16 ?886次閱讀

    【PHYTEC AM62x開發板試用】ETHRNET以太網口調試過程

    現在將我調試PHYTEC AM62X以太網口過程發表一下。 通過認真看了開發板接口Interface章節內容,一步一步掌握了許多接囗操作方法,比如HDMI,UART,USB,I2C等,這里著重介紹
    發表于 06-25 11:45

    【PHYTEC AM62x開發板試用】開箱

    本帖最后由 如去如來 于 2024-6-27 12:28 編輯 收到PHYTEC AM62x開發板已經兩周,忙于其他事情還沒來得及開箱測試,今天來補上。 自從1986年德國PHYTEC公司創立
    發表于 06-18 23:27

    示波器噪聲的信號處理技術研究

    在電子測量領域,示波器作為一種重要的測量工具,廣泛應用于各種電子信號的觀測和分析中。然而,在實際應用中,示波器往往會受到各種噪聲的干擾,導致測量結果的準確性和可靠性受到影響。因此,對
    的頭像 發表于 05-17 17:17 ?644次閱讀

    沒搶到?米爾AM62x開發板限量5折又來了

    AM62x替代AM335x處理器上次米爾AM62x開發板的5折活動沒搶到,今天又來了!!!精彩不容錯過,限量30PCSSitara作為TI處理
    的頭像 發表于 05-17 08:01 ?551次閱讀
    沒搶到?米爾AM<b class='flag-5'>62</b>x<b class='flag-5'>開發板</b>限量5折又來了

    鴻蒙OpenHarmony【輕量系統 環境搭建】 (基于Hi3861開發板

    除上述[安裝庫和工具]和[安裝編譯工具]外,針對Hi3861開發板還需要安裝特定的編譯工具
    的頭像 發表于 04-26 17:47 ?1975次閱讀
    鴻蒙OpenHarmony【輕量系統 環境搭建】 (基于Hi3861<b class='flag-5'>開發板</b>)

    玩轉PSoC 6 RTT積木式開發套件,實現毫米波雷達等實用功能

    本期英飛凌手工課,將由來自英飛凌的工程師Jenson給大家帶來PSoC62withCAPSENSEevaluationkit(下稱PSoC6RTT開發板)和其積木式開發套件的全面介紹。
    的頭像 發表于 03-20 08:35 ?1065次閱讀
    玩轉<b class='flag-5'>PSoC</b> 6 RTT積木式<b class='flag-5'>開發</b>套件,實現毫米波雷達等實用功能

    fpga開發板是什么?fpga開發板有哪些?

    FPGA開發板是一種基于FPGA(現場可編程門陣列)技術的開發平臺,它允許工程師通過編程來定義和配置FPGA芯片上的邏輯電路,以實現各種數字電路和邏輯功能。FPGA開發板通常包括FPGA芯片、時鐘模塊、電源模塊、輸入輸出接口等組
    的頭像 發表于 03-14 18:20 ?2192次閱讀

    英飛凌PSoC? 62開發板 -LCD ILI9341 8080 DEMO

    英飛凌開發板
    rx_ted
    發布于 :2024年03月07日 21:08:39

    【英飛凌測評】英飛凌PSoC? 62開發板試用報告-LCD ILI9341 8080 DEMO

    # Infineon Psoc6-evaluationkit-062S2 說明 ## 簡介 本文檔為 `RT-Thread` 為 `PSoC6 CY8CKIT-062S2-43012` 開發板
    的頭像 發表于 03-07 21:07 ?776次閱讀
    【英飛凌測評】英飛凌<b class='flag-5'>PSoC</b>? <b class='flag-5'>62</b><b class='flag-5'>開發板</b>試用報告-LCD ILI9341 8080 DEMO

    如何在PSoC62 Evaluation kit開發板上運行RT-Thread呢?

    本文主要介紹如何在PSoC62 Evaluation kit開發板上運行RT-Thread,如何連接CYW43012 WiFi藍牙模組,以及如何進行WiFi聯網、網絡延遲、網絡帶寬測試。
    的頭像 發表于 02-27 10:49 ?1093次閱讀
    如何在<b class='flag-5'>PSoC62</b> Evaluation kit<b class='flag-5'>開發板</b>上運行RT-Thread呢?
    主站蜘蛛池模板: 欧美精品久久久久性色AV苍井 | 最近中文字幕完整版免费 | 18禁无遮遮挡羞漫画免费阅读 | 精品一成人岛国片在线观看 | 老师系列高H文 | 伦理片免费秋霞e | 欧美ZC0O人与善交的最新章节 | 国产精品久久久久久亚洲毛片 | 亚洲精品国偷拍电影自产在线 | 在线欧美精品一区二区三区 | 香蕉鱼视频观看在线视频下载 | 一本道综合久久免费 | yellow视频免费观看 | 99热这里有精品 | 免费毛片在线播放 | 国产真实强被迫伦姧女在线观看 | 羞羞影院男女爽爽影院尤物 | 国产精品久久久久久久久99热 | xx顶级欧美熟妞xxhd | 久久无码人妻中文国产 | 蜜芽丅v新网站在线观看 | 蜜臀AV人妻久久无码精品麻豆 | 亚洲国产在线播放在线 | 精油按摩日本 | 国产性色AV内射白浆肛交后入 | 97在线看视频福利免费 | 国产在线精品一区二区在线看 | 青草在线在线d青草在线 | 美妇教师双飞后菊 | 国产午夜伦鲁鲁 | 最近中文字幕MV免费高清在线 | 国产欧美日韩国产高清 | 国产乱人精品视频AV麻豆 | 99国产精品久久久久久久日本竹 | 夜夜穞狠狠穞 | ppypp日本欧美一区二区 | 国产免费阿v精品视频网址 国产免费69成人精品视频 | 99re久久热在线播放快 | 日日AV夜夜添久久奶无码 | 日本少妇无码精品12P | 久久中文字幕亚洲 |