數(shù)字濾波算法.doc

打開方式: Microsoft Word 2003
資料大小: 40.5 KB
下載次數(shù): 27
上傳會(huì)員: admin
上傳時(shí)間: 2012-05-25
本站網(wǎng)址: www.www.9936d.com
查看附件所在的主題
下載附件

10種軟件濾波方法:

10種軟件濾波方法
1、限幅濾波法（又稱程序判斷濾波法）
A、方法：
      根據(jù)經(jīng)驗(yàn)判斷，確定兩次采樣允許的最大偏差值（設(shè)為A）
      每次檢測(cè)到新值時(shí)判斷：
      如果本次值與上次值之差<=A,則本次值有效
      如果本次值與上次值之差>A,則本次值無效,放棄本次值,用上次值代替本次值
B、優(yōu)點(diǎn)：
      能有效克服因偶然因素引起的脈沖干擾
C、缺點(diǎn)
      無法抑制那種周期性的干擾
      平滑度差

2、中位值濾波法
A、方法：
      連續(xù)采樣N次（N取奇數(shù)）
      把N次采樣值按大小排列
      取中間值為本次有效值
B、優(yōu)點(diǎn)：
      能有效克服因偶然因素引起的波動(dòng)干擾
      對(duì)溫度、液位的變化緩慢的被測(cè)參數(shù)有良好的濾波效果
C、缺點(diǎn)：
      對(duì)流量、速度等快速變化的參數(shù)不宜
3、算術(shù)平均濾波法
A、方法：
      連續(xù)取N個(gè)采樣值進(jìn)行算術(shù)平均運(yùn)算
      N值較大時(shí)：信號(hào)平滑度較高，但靈敏度較低
      N值較小時(shí)：信號(hào)平滑度較低，但靈敏度較高
      N值的選取：一般流量，N=12；壓力：N=4
B、優(yōu)點(diǎn)：
      適用于對(duì)一般具有隨機(jī)干擾的信號(hào)進(jìn)行濾波
      這樣信號(hào)的特點(diǎn)是有一個(gè)平均值，信號(hào)在某一數(shù)值范圍附近上下波動(dòng)
C、缺點(diǎn)：
      對(duì)于測(cè)量速度較慢或要求數(shù)據(jù)計(jì)算速度較快的實(shí)時(shí)控制不適用
      比較浪費(fèi)RAM

4、遞推平均濾波法（又稱滑動(dòng)平均濾波法）
A、方法：
      把連續(xù)取N個(gè)采樣值看成一個(gè)隊(duì)列
      隊(duì)列的長(zhǎng)度固定為N
      每次采樣到一個(gè)新數(shù)據(jù)放入隊(duì)尾,并扔掉原來隊(duì)首的一次數(shù)據(jù).(先進(jìn)先出原則)
      把隊(duì)列中的N個(gè)數(shù)據(jù)進(jìn)行算術(shù)平均運(yùn)算,就可獲得新的濾波結(jié)果
      N值的選取：流量，N=12；壓力：N=4；液面，N=4~12；溫度，N=1~4
B、優(yōu)點(diǎn)：
      對(duì)周期性干擾有良好的抑制作用，平滑度高
      適用于高頻振蕩的系統(tǒng)
C、缺點(diǎn)：
      靈敏度低
      對(duì)偶然出現(xiàn)的脈沖性干擾的抑制作用較差
      不易消除由于脈沖干擾所引起的采樣值偏差
      不適用于脈沖干擾比較嚴(yán)重的場(chǎng)合
      比較浪費(fèi)RAM

5、中位值平均濾波法（又稱防脈沖干擾平均濾波法）
A、方法：
      相當(dāng)于“中位值濾波法”+“算術(shù)平均濾波法”
      連續(xù)采樣N個(gè)數(shù)據(jù)，去掉一個(gè)最大值和一個(gè)最小值
      然后計(jì)算N-2個(gè)數(shù)據(jù)的算術(shù)平均值
      N值的選取：3~14
B、優(yōu)點(diǎn)：
      融合了兩種濾波法的優(yōu)點(diǎn)
      對(duì)于偶然出現(xiàn)的脈沖性干擾，可消除由于脈沖干擾所引起的采樣值偏差
C、缺點(diǎn)：
      測(cè)量速度較慢，和算術(shù)平均濾波法一樣
      比較浪費(fèi)RAM

6、限幅平均濾波法
A、方法：
      相當(dāng)于“限幅濾波法”+“遞推平均濾波法”
      每次采樣到的新數(shù)據(jù)先進(jìn)行限幅處理，
      再送入隊(duì)列進(jìn)行遞推平均濾波處理
B、優(yōu)點(diǎn)：
      融合了兩種濾波法的優(yōu)點(diǎn)
      對(duì)于偶然出現(xiàn)的脈沖性干擾，可消除由于脈沖干擾所引起的采樣值偏差
C、缺點(diǎn)：
      比較浪費(fèi)RAM
7、一階滯后濾波法
A、方法：
      取a=0~1
      本次濾波結(jié)果=（1-a）*本次采樣值+a*上次濾波結(jié)果
B、優(yōu)點(diǎn)：
      對(duì)周期性干擾具有良好的抑制作用
      適用于波動(dòng)頻率較高的場(chǎng)合
C、缺點(diǎn)：
      相位滯后，靈敏度低
      滯后程度取決于a值大小
      不能消除濾波頻率高于采樣頻率的1/2的干擾信號(hào)

8、加權(quán)遞推平均濾波法
A、方法：
      是對(duì)遞推平均濾波法的改進(jìn)，即不同時(shí)刻的數(shù)據(jù)加以不同的權(quán)
      通常是，越接近現(xiàn)時(shí)刻的數(shù)據(jù)，權(quán)取得越大。
      給予新采樣值的權(quán)系數(shù)越大，則靈敏度越高，但信號(hào)平滑度越低
B、優(yōu)點(diǎn)：
      適用于有較大純滯后時(shí)間常數(shù)的對(duì)象
      和采樣周期較短的系統(tǒng)
C、缺點(diǎn)：
      對(duì)于純滯后時(shí)間常數(shù)較小，采樣周期較長(zhǎng)，變化緩慢的信號(hào)
      不能迅速反應(yīng)系統(tǒng)當(dāng)前所受干擾的嚴(yán)重程度，濾波效果差
9、消抖濾波法
A、方法：
      設(shè)置一個(gè)濾波計(jì)數(shù)器
      將每次采樣值與當(dāng)前有效值比較：
      如果采樣值＝當(dāng)前有效值，則計(jì)數(shù)器清零
      如果采樣值<>當(dāng)前有效值，則計(jì)數(shù)器+1，并判斷計(jì)數(shù)器是否>=上限N(溢出)
         如果計(jì)數(shù)器溢出,則將本次值替換當(dāng)前有效值,并清計(jì)數(shù)器
B、優(yōu)點(diǎn)：
      對(duì)于變化緩慢的被測(cè)參數(shù)有較好的濾波效果,
      可避免在臨界值附近控制器的反復(fù)開/關(guān)跳動(dòng)或顯示器上數(shù)值抖動(dòng)
C、缺點(diǎn)：
      對(duì)于快速變化的參數(shù)不宜
      如果在計(jì)數(shù)器溢出的那一次采樣到的值恰好是干擾值,則會(huì)將干擾值當(dāng)作有效值導(dǎo)入系統(tǒng)
10、限幅消抖濾波法
A、方法：
      相當(dāng)于“限幅濾波法”+“消抖濾波法”
      先限幅,后消抖
B、優(yōu)點(diǎn)：
      繼承了“限幅”和“消抖”的優(yōu)點(diǎn)
      改進(jìn)了“消抖濾波法”中的某些缺陷,避免將干擾值導(dǎo)入系統(tǒng)
C、缺點(diǎn)：
      對(duì)于快速變化的參數(shù)不宜
10種軟件濾波方法的示例程序
OurWay 發(fā)表于 2005-9-2 22:24:00
10種軟件濾波方法的示例程序(JKRL)
假定從8位AD中讀取數(shù)據(jù)（如果是更高位的AD可定義數(shù)據(jù)類型為int),子程序?yàn)?/font>get_ad();
1、限副濾波
/*  A值可根據(jù)實(shí)際情況調(diào)整
value為有效值，new_value為當(dāng)前采樣值
濾波程序返回有效的實(shí)際值  */
#define A 10
char value;
char filter()
{
char  new_value;
new_value = get_ad();
if ( ( new_value - value > A ) || ( value - new_value > A )
   return value;
return new_value;

}
2、中位值濾波法
/*  N值可根據(jù)實(shí)際情況調(diào)整
排序采用冒泡法*/
#define N  11
char filter()
{
char value_buf[N];
char count,i,j,temp;
for ( count=0;count<N;count++)
{
   value_buf[count] = get_ad();
   delay();
}
for (j=0;j<N-1;j++)
{
   for (i=0;i<N-j;i++)
   {
      if ( value_buf>value_buf[i+1] )
      {
         temp = value_buf;
         value_buf = value_buf[i+1];
         value_buf[i+1] = temp;
      }
   }
}
return value_buf[(N-1)/2];
}
3、算術(shù)平均濾波法
#define N 12
char filter()
{
int  sum = 0;
for ( count=0;count<N;count++)
{
   sum + = get_ad();
   delay();
}
return (char)(sum/N);
}
4、遞推平均濾波法（又稱滑動(dòng)平均濾波法）
#define N 12
char value_buf[N];
char i=0;
char filter()
{
char count;
int  sum=0;
value_buf[i++] = get_ad();
if ( i == N ) i = 0;
for ( count=0;count<N,count++)
   sum = value_buf[count];
return (char)(sum/N);
}
5、中位值平均濾波法（又稱防脈沖干擾平均濾波法）
#define N 12
char filter()
{
char count,i,j;
char value_buf[N];
int  sum=0;
for  (count=0;count<N;count++)
{
   value_buf[count] = get_ad();
   delay();
}
for (j=0;j<N-1;j++)
{
   for (i=0;i<N-j;i++)
   {
      if ( value_buf>value_buf[i+1] )
      {
         temp = value_buf;
         value_buf = value_buf[i+1];
         value_buf[i+1] = temp;
      }
   }
}
for(count=1;count<N-1;count++)
   sum += value[count];
return (char)(sum/(N-2));
}
6、限幅平均濾波法
/*
*/
略參考子程序1、3
7、一階滯后濾波法
/* 為加快程序處理速度假定基數(shù)為100，a=0~100 */
#define a 50
char value;
char filter()
{
char  new_value;
new_value = get_ad();
return (100-a)*value + a*new_value;
}
8、加權(quán)遞推平均濾波法
/* coe數(shù)組為加權(quán)系數(shù)表，存在程序存儲(chǔ)區(qū)。*/
#define N 12
char code coe[N] = {1,2,3,4,5,6,7,8,9,10,11,12};
char code sum_coe = 1+2+3+4+5+6+7+8+9+10+11+12;
char filter()
{
char count;
char value_buf[N];
int  sum=0;
for (count=0,count<N;count++)
{
   value_buf[count] = get_ad();
   delay();
}
for (count=0,count<N;count++)
   sum += value_buf[count]*coe[count];
return (char)(sum/sum_coe);
}
9、消抖濾波法
#define N 12
char filter()
{
char count=0;
char new_value;
new_value = get_ad();
while (value !=new_value);
{
   count++;
   if (count>=N) return new_value;
   delay();
   new_value = get_ad();
}
return value;
}
10、限幅消抖濾波法
/*
*/
略參考子程序1、9
11、IIR 數(shù)字濾波器
A. 方法：
確定信號(hào)帶寬，濾之。
Y(n) = a1*Y(n-1) + a2*Y(n-2) + ... + ak*Y(n-k) + b0*X(n) + b1*X(n-1) + b2*X(n-2) + ... + bk*X(n-k)
B. 優(yōu)點(diǎn)：高通，低通，帶通，帶阻任意。設(shè)計(jì)簡(jiǎn)單(用matlab）
C. 缺點(diǎn)：運(yùn)算量大。

IIR濾波例子
int  BandpassFilter4(int InputAD4)
{
int  ReturnValue;
int  ii;
RESLO=0;
RESHI=0;
MACS=*PdelIn;
OP2=1068; //FilterCoeff4[4];
MACS=*(PdelIn+1);
OP2=8; //FilterCoeff4[3];
MACS=*(PdelIn+2);
OP2=-2001;//FilterCoeff4[2];
MACS=*(PdelIn+3);
OP2=8; //FilterCoeff4[1];
MACS=InputAD4;
OP2=1068; //FilterCoeff4[0];
MACS=*PdelOu;
OP2=-7190;//FilterCoeff4[8];
MACS=*(PdelOu+1);
OP2=-1973; //FilterCoeff4[7];
MACS=*(PdelOu+2);
OP2=-19578;//FilterCoeff4[6];
MACS=*(PdelOu+3);
OP2=-3047; //FilterCoeff4[5];
*p=RESLO;
*(p+1)=RESHI;
mytestmul<<=2;
ReturnValue=*(p+1);
for  (ii=0;ii<3;ii++)
{
   DelayInput[ii]=DelayInput[ii+1];
   DelayOutput[ii]=DelayOutput[ii+1];
   }
   DelayInput[3]=InputAD4;
   DelayOutput[3]=ReturnValue;

//  if (ReturnValue<0)
//  {
//  ReturnValue=-ReturnValue;
//  }
return ReturnValue;
}
　
//------------------------------------------------------------------------------
防脈沖干擾移動(dòng)平均值法數(shù)字濾波器的C語(yǔ)言算法及其實(shí)現(xiàn)

   在許多的數(shù)據(jù)采集系統(tǒng)中，現(xiàn)場(chǎng)的強(qiáng)電設(shè)備較多，不可避免
  地會(huì)產(chǎn)生尖脈沖干擾,這種干擾一般持續(xù)時(shí)間短，峰值大，對(duì)這樣
  的數(shù)據(jù)進(jìn)行數(shù)字濾波處理時(shí)，僅僅采用算術(shù)平均或移動(dòng)平均濾波
  時(shí)，盡管對(duì)脈沖干擾進(jìn)行了1/n的處理，但，其剩余值仍然較大。
  這種場(chǎng)合最好的策略是：將被認(rèn)為是受干擾的信號(hào)數(shù)據(jù)去掉，這
  就是防脈沖干擾平均值濾波法的原理。
防脈沖干擾平均值濾波法的算法是：對(duì)連續(xù)的n個(gè)數(shù)據(jù)進(jìn)行排序，
  去掉其中最大和最小的2個(gè)數(shù)據(jù)，將剩余數(shù)據(jù)示平均值。
在一般8051單片機(jī)的應(yīng)用中為了加快數(shù)據(jù)處理速度，n可以取值6。
  而對(duì)于具有較快速度的處理器，則n值可以適當(dāng)取大一些。但最好是
  n=2^k+2，  k為整數(shù),因?yàn)檫@樣在求平均值average=SUM/(n-2)=SUM/2^k時(shí)，
  可以寫成average=SUM>>k，用移位的方法，可以加快處理速度。

上述算法顯然還存在一個(gè)不足之處，就是每采集一個(gè)數(shù)據(jù)就要進(jìn)行一次排序，
  這樣會(huì)大量占用系統(tǒng)寶貴的時(shí)間。這可以通過存儲(chǔ)當(dāng)前數(shù)據(jù)中的最大值和最小值
  來改進(jìn)。具體做法是：
系統(tǒng)中用兩個(gè)變量來存儲(chǔ)當(dāng)前n個(gè)數(shù)據(jù)的最大值和最小值在這個(gè)數(shù)組中的
  偏移量(也就是數(shù)組下標(biāo)，存儲(chǔ)數(shù)組下標(biāo)而直接不存儲(chǔ)數(shù)據(jù)本身是因?yàn)椋涸谝话愕?/font>
  系統(tǒng)中，n不會(huì)超無符號(hào)短整形的表示范圍，因此用一個(gè)char形變量就可以存儲(chǔ)了
  而如果直接存儲(chǔ)數(shù)據(jù)本身，則許多情況下要用int形變量，甚至更長(zhǎng)的類型)。這樣
  只要在當(dāng)前輸入的數(shù)據(jù)將要覆蓋的數(shù)據(jù)正好是當(dāng)前的最大值或最小值時(shí)才在下個(gè)數(shù)
  組中查找最大值或最小值，而其他情況下則只要將輸入的數(shù)據(jù)與最大值和最小值比較
  就可以修改下最大值和最小值了，而且不用進(jìn)行數(shù)據(jù)排序。
  這個(gè)算法很簡(jiǎn)單，下面是對(duì)應(yīng)的C語(yǔ)言代碼實(shí)現(xiàn)，可以很方便的應(yīng)用的具體的51單片機(jī)，
  或其他處理器上，只須做少量的修改。

＃i nclude"stdio.h"
#define dtype unsigned int // 采集數(shù)據(jù)的數(shù)據(jù)類型
#define uint8 char
#define LEN  6 //移動(dòng)算術(shù)平均的個(gè)數(shù)+2=SHIFT<<2+2
#define SHIFT 2 //2^SHIFT
uint8 pdata; //移動(dòng)指針
uint8 pmax,pmin; //記錄數(shù)據(jù)表中最大值和最小值的位置,
   //在一般的數(shù)據(jù)采集系統(tǒng)中,數(shù)據(jù)的長(zhǎng)度>=8,
   //因此用指針記錄而不是直接記錄最大值和最小值
dtype datas[LEN];
dtype szlb(dtype _data)
{
/****************************/
/* 在調(diào)用此子程序前必須對(duì) */
/* pdata,datas[]數(shù)組,  */
/* pmax,pmin進(jìn)行初始化  */
/****************************/
uint8 i;
dtype average=0;  //清零，用來計(jì)算平均值
pdata=(pdata+1)%LEN; //指針下標(biāo)在0到LEN-1上滑動(dòng)
datas[pdata]=_data;  //采樣所得數(shù)據(jù)存入數(shù)據(jù)表中
for(i=0;i<LEN;i++)
      average+=datas; //求所有數(shù)據(jù)總和
/*******去除被認(rèn)為是脈沖的數(shù)據(jù)******/
if(_data>datas[pmax])
   pmax=pdata; //得到最大值的指針
else if(_data<datas[pmin])
   pmin=pdata; //得到最小值的指針
if(pdata==pmax) //如果當(dāng)前輸入值將存入當(dāng)前最大值的位置時(shí)
{    //由以上方法將不可行，必須從其他位置中查找極值
   for(i=0;i<LEN;i++)
      if(datas>datas[pmax])
      pmax=i;
}
else if(pdata==pmin)//如果當(dāng)前輸入值將存入當(dāng)前最大值的位置時(shí)
{    //由以上方法將不可行，必須從其他位置中查找極值
   for(i=0;i<LEN;i++)
      if(datas<datas[pmin])
         pmin=i;
}
average=average-datas[pmax]-datas[pmin];//減去脈沖
return (average>>SHIFT); //求算術(shù)平均值
}
/******以下是在VC＋＋6.0環(huán)境下運(yùn)行的測(cè)試程序**/
/***通過手動(dòng)輸入來模擬數(shù)據(jù)采集過程****/
void main()
{
uint8 i;
dtype _data;
pdata=0;
pmax=0;
pmin=0;
for(i=0;i<LEN;i++)
datas=0;
printf("數(shù)據(jù)：                         最大    最小\n");
while(1)
{
   scanf("%u",&_data);
   szlb(_data);
   for(i=0;i<LEN;i++)
   printf("%-3u  ",datas);
   printf(" %-3u    %-3u",datas[pmax],datas[pmin]);
   printf("\n");
}
}

算術(shù)平均濾波

所謂算術(shù)平均濾波就是把n個(gè)采樣值相加，然后取其算術(shù)平均值作為本次有效的采樣信號(hào)

PUSH  PSW ；現(xiàn)場(chǎng)保護(hù)

   PUSH  A

   MOV FLAG,#00H ；進(jìn)位位清零

   MOV R0,DATA  ；設(shè)置數(shù)據(jù)存儲(chǔ)區(qū)首址

   MOV R7,#08H ；設(shè)置采樣數(shù)據(jù)個(gè)數(shù)

   CLR A ；清累加器

LOOP: ADD A,@R0    ；兩數(shù)相加

   JNC    NEXT ；無進(jìn)位，轉(zhuǎn)NEXT

   INC    FLAG ；有進(jìn)位，進(jìn)位位加1

NEXT: INC    R0 ；數(shù)據(jù)指針加1

   DJNZ  R7, LOOP ；未加完，繼續(xù)加

   MOV R7,#03H ；設(shè)置循環(huán)次數(shù)

DIVIDE: MOV TEMP,A ；保存累加器中的內(nèi)容

   MOV A,FLAG ；累加結(jié)果除2

   CLR C

   RRC A

   MOV FLAG,A

   MOV A,TEMP

   RRC A

   DJNZ  R7,DIVIDE    ；未結(jié)束，繼續(xù)執(zhí)行

   MOV SAMP,A ；保存結(jié)果至SAMP中

   POP A ；恢復(fù)現(xiàn)場(chǎng)

   POP PSW

   RET

限幅濾波程序

設(shè)在順序采樣時(shí)刻T1、T2、T3所采集的數(shù)據(jù)分別為Y1、Y2、Y3，則當(dāng)

∣Y2-Y1∣≤△y，則Y2作為采樣值；

∣Y2-Y1∣>△y，則保留Y2，但不作為采樣值，繼續(xù)采樣得Y3；

如果∣Y3-Y2∣≤△y，則Y3作為采樣值；

∣Y3-Y2∣>△y，則取作為采樣值。

程序

PUSH PSW    ；保護(hù)現(xiàn)場(chǎng)

   PUSH A

   CLR    C ；進(jìn)位標(biāo)志位清零

   MOV    DATA,DATA2

   MOV    A,DATA1

   SUBB A,DATA ；求Yn-1 -Yn

   JNC    COMPARE  ；如果Yn-1 -Yn≥0,轉(zhuǎn)COMPARE

   CPL    A ；如果Yn-1 -Yn<0, 求補(bǔ)

   INC    A

COMPARE:  CLR    C

   SUBB A,LIMIT    ；│Yn-Yn-1│和△y比較

   JC OVER ；如果│Yn-Yn-1│≤△y,DATA2→DATA

   MOV    DATA,DATA1    ；如果│Yn-Yn-1│>△y，DATA1→DATA

OVER: POP    A ；恢復(fù)現(xiàn)場(chǎng)

   POP    PSW

   RET    ；返回

如何選用濾波電容
  　濾波電容在開關(guān)電源中起著非常重要的作用，如何正確選擇濾波電容，尤其是輸出濾
波電容的選擇則是每個(gè)工程技術(shù)人員都十分關(guān)心的問題。50Hz工頻電路中使用的普通電解
電容器，其脈動(dòng)電壓頻率僅為100Hz，充放電時(shí)間是毫秒數(shù)量級(jí)。
為獲得更小的脈動(dòng)系數(shù)，所需的電容量高達(dá)數(shù)十萬μF，因此普通低頻鋁電解電容器的
目標(biāo)是以提高電容量為主，電容器的電容量、損耗角正切值以及漏電流是鑒別其優(yōu)劣的主
要參數(shù)。而開關(guān)電源中的輸出濾波電解電容器，其鋸齒波電壓頻率高達(dá)數(shù)十kHz，甚至是數(shù)
十MHz，這時(shí)電容量并不是其主要指標(biāo)，衡量高頻鋁電解電容優(yōu)劣的標(biāo)準(zhǔn)是“阻抗-頻率”
特性，要求在開關(guān)電源的工作頻率內(nèi)要有較低的等效阻抗，同時(shí)對(duì)于半導(dǎo)體器件工作時(shí)產(chǎn)
生的高頻尖峰信號(hào)具有良好的濾波作用。
普通的低頻電解電容器在10kHz左右便開始呈現(xiàn)感性，無法滿足開關(guān)電源的使用要求。
而開關(guān)電源專用的高頻鋁電解電容器有四個(gè)端子，正極鋁片的兩端分別引出作為電容器的
正極，負(fù)極鋁片的兩端也分別引出作為負(fù)極。電流從四端電容的一個(gè)正端流入，經(jīng)過電容
內(nèi)部，再?gòu)牧硪粋€(gè)正端流向負(fù)載；從負(fù)載返回的電流也從電容的一個(gè)負(fù)端流入，再?gòu)牧硪?/font>
個(gè)負(fù)端流向電源負(fù)端。由于四端電容具有良好的高頻特性，為減小電壓的脈動(dòng)分量以及抑
制開關(guān)尖峰噪聲提供了極為有利的手段。高頻鋁電解電容器還有多芯的形式，即將鋁箔分
成較短的若干段，用多引出片并聯(lián)連接以減小容抗中的阻抗成份。并且采用低電阻率的材
料作為引出端子，提高了電容器承受大電流的能力。

不錯(cuò)啊，，這個(gè)很好用

賬號(hào)		自動(dòng)登錄	找回密碼
密碼			注冊(cè)會(huì)員

板塊導(dǎo)航

熱門下載

數(shù)字濾波算法.doc