DSP C語言的IIR濾波器程序源碼

ID:209315 · 發(fā)表于 2018-8-8 16:32

分享個(gè)DSP C語言的IIR濾波器

單片機(jī)源程序如下:

#include <math.h>
#include <stdio.h>
#include "fdacoefs.h"
//*********************************************
//fs=1000000,fpass=10000,fstop=50000(-80db)
//butterworth iir-df1
//*********************************************
#define fs 1000000
#define f1 4000//測試正弦頻率1
#define f2 50000//測試正弦頻率2
#define f_rec 5000//測試方波頻率
#define Impulse 1
#define Sine 2
#define Rect 3
#define pi 3.1415926
#define STG_NUM (MWSPT_NSEC-1)/2
#define DAT_NUM 800
float IIR_out[DAT_NUM];
float input_step[STG_NUM+1][3];// n; n-1; n-2
void iir_filter_init(float *ptr,int num)
{
int i;
for(i=0;i<num;i++)
{
*ptr++=0;
}
}
float iir_filter(float in)
{
int i;
float temp;
input_step[0][0] = in;
for(i=0;i<STG_NUM;i++)
{
temp =NUM[i*2][0]*(input_step[i][0]*NUM[i*2+1][0]+input_step[i][1]*NUM[i*2+1][1]+input_step[i][2]*NUM[i*2+1][2])-
input_step[i+1][1]*DEN[i*2+1][1]-input_step[i+1][2]*DEN[i*2+1][2];
temp *= DEN[i*2+1][0];
input_step[i][2]=input_step[i][1];
input_step[i][1]=input_step[i][0];
input_step[i+1][0]=temp;
}
return temp*NUM[MWSPT_NSEC-1][0];
}
float x;
int main()
{
int i,Test_type;
printf("請選擇測試類型：1、沖擊測試 2、正弦測試 3、方波測試\n");
scanf("%d",&Test_type);
while(1){
iir_filter_init(input_step[0],sizeof(input_step)/sizeof(input_step[0][0]));
switch(Test_type)
{
case Rect://方波輸入
for(i=0;i<DAT_NUM;i++)
{
if(((2*f_rec*i/fs)%2)==0)
x=5000;
else
x=-5000;
IIR_out[i]=iir_filter(x);
}
break;
case Impulse://沖激輸入
IIR_out[0]=iir_filter(5000);
for(i=1;i<DAT_NUM;i++)
IIR_out[i]=iir_filter(0);
break;
case Sine://正弦輸入
for(i=0;i<DAT_NUM;i++)
{
//***************************************************************
x=500*(cos(2*pi*f1*i/fs)+cos(2*pi*f2*i/fs)); //正弦測試
IIR_out[i]=iir_filter(x);
}
break;
}
//請?jiān)趐rintf前設(shè)置斷點(diǎn)，以便停止運(yùn)行查看輸出結(jié)果否則，看不到輸出
printf("請選擇測試類型：1、沖擊測試 2、正弦測試 3、方波測試\n");
scanf("%d",&Test_type);
}
while(1);
}

復(fù)制代碼

所有資料51hei提供下載:

IIR.7z (23.26 KB, 下載次數(shù): 78)

ID:95821 · 發(fā)表于 2018-8-10 11:55

學(xué)習(xí)了,正是在找的

ID:209315 · 發(fā)表于 2018-8-10 16:14

alai318 發(fā)表于 2018-8-10 11:55
學(xué)習(xí)了,正是在找的

努力加油

ID:95821 · 發(fā)表于 2018-8-12 12:36

剛學(xué)了FFT基2算法，F(xiàn)FT還沒學(xué)透，學(xué)好再來學(xué)濾波方面的

ID:209315 · 發(fā)表于 2018-8-14 14:05

alai318 發(fā)表于 2018-8-12 12:36
剛學(xué)了FFT基2算法，F(xiàn)FT還沒學(xué)透，學(xué)好再來學(xué)濾波方面的

老哥厲害了

ID:216962 · 發(fā)表于 2018-8-14 14:16

謝謝分享，論壇有您更精彩！

ID:424128 · 發(fā)表于 2018-11-10 18:29

老哥，這個(gè)程序適合用于DSP C54x做IIR濾波器的實(shí)物嗎，就是IIR數(shù)字濾波器的MATLAB仿真及DSP實(shí)現(xiàn)這種畢業(yè)設(shè)計(jì)題，能做出實(shí)物不？

ID:392405 · 發(fā)表于 2019-7-27 10:49

感謝樓主分享

ID:276779 · 發(fā)表于 2019-10-22 17:59

看一下正需要呢

帳號		自動登錄	找回密碼
密碼			立即注冊

久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

DSP C語言的IIR濾波器程序源碼

評分