自制示波器/频谱显示/DDS。
基于ucos+ARMcortexM4+FPGA。
带FIR，FFT,Kalman算法。
最高采样率50M。存储8K*2字节数据。每秒64次波形（VGA刷屏63次/秒，快了也显示不出）
fpga：数据采集，DDS波形产生，提供VGA时序驱动屏幕。
ARMcortexM4：外设接口控制，数据处理。
FPGA与CPU通过8位并行总线通信。CPU给FPGA内置RAM写 命令/数据
先上张图
频谱显示功能演示（FFT）对着数字信号处理教材自己写的，用递归和非递归两种方式实现，和MATLAB对照确认计算无误。
高通滤波器与低通滤波器功能演示（理想FIR滤波器+汉明窗）单独开任务计算FIR参数。每输入一次计算一次参数。
MATLAB仿真
卡尔曼滤波算法（kalman算法对于随机噪声为最优），一维卡尔曼算法，无转置矩阵。
MATLAB仿真，不同参数的kalman算法
任务调度部分代码。（#include fpga控制.h，数学算法.h，寄存器配置BSP.h）
#include <app.h>
#include "ucos_ii.h" //uC/OS-II系统函数头文件
#include "sysinit.h" //系统配置
#include "task_header.h" //任务相关宏 结构体
#include "fpga_data.h" //与FPGA通信部分
#include "bsp.h" //ARM驱动
#include "character.h" //字模
#include "dds_data.h" //dds数据
#include "mcu_math.h" //数学算法部分
//#include "arm_math.h"
/********临时定义变量******************/
void system_init_task(void *pdata){
//double cc;
pdata = pdata;
#if OS_CRITICAL_METHOD == 3 //关中断的方法为3
OS_CPU_SR cpu_sr;
#endif
OS_ENTER_CRITICAL();
OS_CPU_SysTickInit(CPU_frequency/1000); //初始化OS节拍无比重要
port_init(); //设置所有IO为低电平&&输出
//UART4_Init(115200);
init_LED();
key_init(); //将四个IO配备成中断
lcd_data_init(); //lcd 数据传输
PIT0_Init(CPU_frequency/64); //pit定时中断，VGA刷屏速率约为63次/秒
//PIT0_Init(5000000);
//初始化四个电位器AD
AD_init();
OS_EXIT_CRITICAL();
//计算sin相关参数，尼玛不要手动输入。
read_data_init();
LED1(0);
LED2(0);
LED3(0);
LED4(0);
//状态
trigger.state =0;
kalman.state = 0;
fft_wave.state =0;
freq_selc.state =0;
dds.state = 0;
trigger_lock.state =0;
//默认值
trigger.value =0;
kalman.value = 0;
fft_wave.value =0;
freq_selc.value =0;
dds.value = 0;
trigger_lock.value =0;
//光标
trigger.cursor =0;
kalman.cursor = 0;
fft_wave.cursor =0;
freq_selc.cursor =0;
dds.cursor = 0;
trigger_lock.cursor =0;
//双向循环链表
fft_wave.next_menu = &trigger;
trigger.last_menu = &fft_wave;
fft_wave.value_num = 2;
fft_wave.menu_num =1;
trigger.next_menu = &kalman;
kalman.last_menu = &trigger;
trigger.value_num = 3;//上升 下降 双向
trigger.menu_num =2;
kalman.next_menu = &freq_selc;
freq_selc.last_menu = &kalman;
kalman.value_num =2;
kalman.menu_num =3;
freq_selc.next_menu = &trigger_lock;
trigger_lock.last_menu = &freq_selc;
freq_selc.value_num = 3;//开启，关闭，确认输入
freq_selc.menu_num =4;
trigger_lock.next_menu = &dds;
dds.last_menu = &trigger_lock;
trigger_lock.value_num = 2;//锁定，解锁
trigger_lock.menu_num =5;
dds.next_menu = &fft_wave;
fft_wave.last_menu = &dds;
dds.value_num = 4;//方波，正弦，合成波（滤波演示），干扰波（卡尔曼演示）
dds.menu_num =6;
menu_p = &fft_wave;
menu_head = &fft_wave;
//key事件标志
KeyFlag = OSFlagCreate(0,&err);
Menu = OSFlagCreate(0,&err);
Data_Process = OSFlagCreate(0,&err);
//menu信号量，保护全局变量
MenuSem = OSSemCreate(1);
FirParamSem = OSSemCreate(1);
//这里换成互斥信号量，避免优先级反转标记一下。
//LCDdataSem = OSSemCreate(1);
LCDdataSem = OSMutexCreate(Mutex_prio,&err);
//发送参数的邮箱
ParaMbox = OSMboxCreate((void*)0);
FirMbox = OSMboxCreate((void*)0);
OSFlagPost(Menu,(OS_FLAGS)0x01,OS_FLAG_SET,NULL);
/*临时测试add与clk信号是否能用*/
OSTimeDly(5);
OSTaskDel(Init_task_prio);
}
void key_task(void *pdata)
{
pdata = pdata;
while(1)
{
//等待置1
OSFlagPend(KeyFlag,(OS_FLAGS)0x0f,OS_FLAG_WAIT_SET_OR,0,&err);
OSTimeDly(40);
if(!((GPIOD_PDIR>>8)&0x01)) //确认键
{
led_2 = !led_2;
while(!((GPIOD_PDIR>>8)&0x01))
OSTimeDly(1);
OSTimeDly(40);
LED2(led_2);
key_menu(1);
OSFlagPost(KeyFlag,(OS_FLAGS)0x01,OS_FLAG_CLR,NULL);
OSFlagPost(Menu,(OS_FLAGS)0x01,OS_FLAG_SET,NULL);
OSFlagPost(FirFlag,(OS_FLAGS)0x01,OS_FLAG_SET,NULL);
}
else if(!((GPIOD_PDIR>>10)&0x01)) //退出键
{
led_2 = !led_2;
while(!((GPIOD_PDIR>>10)&0x01))
OSTimeDly(1);
OSTimeDly(40);
LED2(led_2);
key_menu(2);
OSFlagPost(KeyFlag,(OS_FLAGS)0x02,OS_FLAG_CLR,NULL);
OSFlagPost(Menu,(OS_FLAGS)0x02,OS_FLAG_SET,NULL);
}
else if(!((GPIOD_PDIR>>12)&0x01)) //左
{
led_2 = !led_2;
while(!((GPIOD_PDIR>>12)&0x01))
OSTimeDly(1);
OSTimeDly(20);
LED2(led_2);
key_menu(3);
OSFlagPost(KeyFlag,(OS_FLAGS)0x04,OS_FLAG_CLR,NULL);
OSFlagPost(Menu,(OS_FLAGS)0x04,OS_FLAG_SET,NULL);
}
else if(!((GPIOD_PDIR>>14)&0x01)) //右
{
led_2 = !led_2;
while(!((GPIOD_PDIR>>14)&0x01))
OSTimeDly(1);
OSTimeDly(20);
LED2(led_2);
key_menu(4);
OSFlagPost(KeyFlag,(OS_FLAGS)0x08,OS_FLAG_CLR,NULL);
OSFlagPost(Menu,(OS_FLAGS)0x08,OS_FLAG_SET,NULL);
}
else
OSFlagPost(KeyFlag,(OS_FLAGS)0x0f,OS_FLAG_CLR,NULL);//清按键数据
}
}
void uart4_task(void *pdata)
{
pdata = pdata;
while(1)
{
/* LED1(0);
Uart4_SendByte(97);
OSTimeDly(100);
LED1(1);
OSTimeDly(100);*/
}
}
void menu_task(void * pdata)
{
//uint8 data;
#if OS_CRITICAL_METHOD == 3 //关中断的方法为3
OS_CPU_SR cpu_sr;
#endif
pdata = pdata;
//cpu_sr = cpu_sr;
while(1)
{
OSFlagPend(Menu,(OS_FLAGS)0x0f,OS_FLAG_WAIT_SET_OR+OS_FLAG_CONSUME,0,&err); //所有按键皆触发该事件
//OSSemPend(LCDdataSem,0,&err);
OSMutexPend(LCDdataSem,0,&err);
show_menu();
//OSSemPost(LCDdataSem);
OSMutexPost(LCDdataSem);
OSFlagPost(Menu,(OS_FLAGS)0x01,OS_FLAG_CLR,NULL);
}
}
void data_process_task(void * pdata)
{
#if OS_CRITICAL_METHOD == 3 //关中断的方法为3
OS_CPU_SR cpu_sr;
#endif
uint16 i,j,trigger_add,add;
uint8 data;
int8 tem;
uint8 wave_data[(LcdLength+FirLength-1)];
uint16 ReadLength;
struct complex_num fft_result[FFT_CNT];
uint8 fft_out[FFT_CNT];
int32 FirTem;
//uint8 wocao[4] = {0x01,0x04,0x02,0x03};
//uint8 sample_para,voltage_para,phase_para,show_para;
uint8 Para[4]; //sample show phase vlotage
//uint8 show_num,;
pdata = pdata;
while(1)
{
OSFlagPend(Data_Process,(OS_FLAGS)0x01,OS_FLAG_WAIT_SET_AND+OS_FLAG_CONSUME,0,&err);
//cc = arm_sin_f32(1.6);
if(!trigger_lock.value) //如果触发
{
//获取AD信号，分频
Para[0] = ADC_Getdata(sample_ctrl);
Para[0] = sample_algorithms(Para[0]);
debug2 = Para[0];
Set_fre_div(Para[0]);//1为2分频，2为四分频……
//开始采集
read_start(1);
//寻找触发源
trigger_add = trigger_seek();
}
else
{
}
Para[1] = (ADC_Getdata(show_ctrl)>>4)+1;
debug = Para[1];
Para[2] = ADC_Getdata(phase_ctrl);
Para[2] =(Para[2]>>2)<<2;
debug1 = Para[2];
OSMboxPost(ParaMbox,&Para[0]);
//phase_para = 0;
//debug = phase_para;
//发送消息
add = trigger_add+Para[2]*Para[1];
//读取更多数据以完成FIR滤波
if(freq_selc.value)
ReadLength = (LcdLength+FirLength-1);
else
ReadLength = LcdLength;
for(i=0;i<ReadLength;i++)
{
//OS_ENTER_CRITICAL();
add += Para[1];
wave_data[i] = read_data(add);
//OS_EXIT_CRITICAL();
}
read_start(0);
if(freq_selc.value==2)
{
for(i=0;i<LcdLength;i++)
{
FirTem=0;
for(j=0;j<FirLength;j++)
FirTem += wave_data[i+j]*FirParam[j];
wave_data[i]=FirTem>>10;
}
}
else if(freq_selc.value==1)
{
for(i=0;i<LcdLength;i++)
{
FirTem=0;
for(j=0;j<FirLength;j++)
FirTem += wave_data[i+j]*FirParam[j];
wave_data[i]=(FirTem>>10)-100;
}
}
else;
if(kalman.value)
KalmanFilter(&wave_data[0],KalmanInitX,KalmanInitP);
else;
if(fft_wave.value==1)
{
OS_ENTER_CRITICAL();
fft(&wave_data[0],&fft_result[0],FFT_LEN,FFT_CNT);
absfft(&fft_result[0],&fft_out[0],FFT_CNT);
OS_EXIT_CRITICAL();
//cc[0] = fft_out[0];
//cc[1] = fft_out[1];
//cc[2] = fft_out[2];
//cc[3] = fft_out[3];
OSMutexPend(LCDdataSem,0,&err);
//命令切换
data = (menu_p->state)?(0x80+clr_data_cnt):clr_data_cnt;
write_command(data);
//写波形数据
data = (menu_p->state)?(0x80+write_data_command):write_data_command;
write_command(data);
for(i=0;i<(FFT_CNT);i++)
{
tem = (fft_out[i+1]-fft_out[i])/3;
write_data(fft_out[i]);
write_data((fft_out[i]+tem));
write_data((fft_out[i]+tem+tem));
}
for(i=0;i<96;i++)
{
write_data(0x00);
}
OSMutexPost(LCDdataSem);
}
else
{
OSMutexPend(LCDdataSem,0,&err);
//命令切换
data = (menu_p->state)?(0x80+clr_data_cnt):clr_data_cnt;
write_command(data);
//写波形数据
data = (menu_p->state)?(0x80+write_data_command):write_data_command;
write_command(data);
for(i=0;i<480;i++)
{
write_data(wave_data[i]);
}
OSMutexPost(LCDdataSem);
}
OSFlagPost(Data_Process,(OS_FLAGS)0x01,OS_FLAG_CLR,NULL);
}
}
void dds_ctrl_task(void *pdata)
{
pdata = pdata;
while(1)
{
OSFlagPend(Menu,(OS_FLAGS)0x01,OS_FLAG_WAIT_SET_AND+OS_FLAG_CONSUME,0,&err); //确认键触发该事件
//OSFlagPend(Menu,(OS_FLAGS)0x01,OS_FLAG_WAIT_SET_AND,0,&err); //尝试清除/不清除标志
if((menu_p ==&dds)&&(!menu_p->state)) //在DDS选项下按了确认 并且状态非下拉
{
OSMutexPend(LCDdataSem,0,&err);
dds_data(dds.value);
OSMutexPost(LCDdataSem);
}
else
{
//OSFlagPost(Menu,(OS_FLAGS)0x01,OS_FLAG_CLR,NULL);
//标志不清除，留给menu事件用。
}
}
}
void FirPara_task(void *pdata)
{
pdata = pdata;
uint8 * ReceivePara;
uint32 FirFre;
while(1)
{
OSFlagPend(Menu,(OS_FLAGS)0x01,OS_FLAG_WAIT_SET_AND+OS_FLAG_CONSUME,0,&err); //确认键触发该事件
if((menu_p ==&freq_selc)&&(!menu_p->state))
//计算fir参数
{
ReceivePara = OSMboxPend(ParaMbox,0,&err);
if(freq_selc.value==1)
{
FirFre = ReceivePara[1]*3*(ReceivePara[0]-1)/25;
FirPara(255,FirFre,DataLength,&FirParam[0],FirLength);
}
else if(freq_selc.value==2)
{
FirFre = ReceivePara[1]*12*(ReceivePara[0]-1)/25;
FirPara(FirFre,0,DataLength,&FirParam[0],FirLength);
}
else;
}
OSTimeDly(10);
}
//参数计算 参数 采样频率 分频 DDS频率
//FirPara(uint8 H,uint8 L, uint16 N,float * fir_c,FirLength)
}
void Paradata_task(void *pdata)
{
uint8 * ReceivePara;
pdata =pdata;
while(1)
{
ReceivePara = OSMboxPend(ParaMbox,0,&err);
OSMutexPend(LCDdataSem,0,&err);
show_para(ReceivePara);
OSMutexPost(LCDdataSem);
OSTimeDly(3);
}
}