滤波器实验报告.docVIP

数据处理算法实验报告 实验的目的: 熟悉掌握子程序、数组等C语言结构的使用方法。 熟悉各种数据处理办法。 实验用设备： PC系列微型计算机，操作系统使用DOS或Windows98。 实验要求： 完成如下数据处理算法,形成相应功能的函数: 算数平均值法 中位值算法 去掉极值的算数平均值算法 滑动平均滤波算法 加权平均滤波算法 数字低通滤波器 具有快速响应的数字滤波器 设计和计算： 从一文本中获取数据，并储存到数组中，传递给不同过程进行数据处理； 分别用不同的此程序实现不同的算法，并输出数据 2.1算数平均值法：【 void sspj(double *p,int n) 】 入口参数： *p : 数据数组指针 n : 数据个数 核心语句： 求平均数： for(i=0;in;i++) { temp+=*p; p++; } temp=temp/n; 作用： 将每一个数据求和，然后求平均数。 2.2中位值算法：【 void zwz(double *p,int n) 】 入口参数：（同上） 核心语句： 冒泡排序： for(i=n-1;i0;i--) for(j=0;ji;j++) { if(*(p+j)*(p+j+1)) { temp=*(p+j); *(p+j)=*(p+j+1); *(p+j+1)=temp; } } 中位值 = *(p+(n/2)+(n%2)-1) 作用： 对输入数据进行排序（从小到大），然后取排序位数最中间的数即为中位值。 2.3去掉极值的算数平均值算法：【void qjsspj(double *p,int n) 】 入口参数：（同上） 核心语句： 冒泡排序：（同上） 求算术平均： for(i=1;in-1;i++) { temp+=*(p+i); } temp=temp/(n-2); 作用： 对输入数据进行排序（从小到大），然后对数据第2位到第（n-1）位进行求和后取平均值。 2.4滑动平均滤波算法：【void hdpj(double *p,int n) 】 入口参数：（同上） 要求输入： it : 滑动平均的数据个数 核心语句： 滑动平均： for(i=it-1;in;i++) { temp=0; for(j=i-it+1;j=i;j++) { temp+=*(p+j); } printf(%d , %lf\n,i-it+1,temp/it); } 作用： 标记i从（it-1）到n滑动，每滑动一次对第（i-it+1）到第i个数据（共it个数据）求和并取平均数，达到滑动滤波的效果。 2.5加权平均滤波算法：【void jqpj(double *p,int n) 】 入口参数：（同上） 要求输入： it : 滑动加权平均的数据个数 a[i] : 每一个数据的权重（共it个数） 核心语句： 对权重进行处理，使其和为1： for(i=0;iit;i++) a[i]=a[i]/temp; 求加权平均： for(i=it-1;in;i++) { temp=0; for(j=i-it+1;j=i;j++) { temp+=(*(p+j))*(a[j+it-i-1]); } printf(%d , %lf\n,i-it+1,temp); } 作用： 在上一个“滑动平均滤波器”的基础上，对滑动的区间内的的每一个数据进行加权，然后求平均数，即为加权滑动平均滤波。 2.6数字低通滤波器：【void szdt(double *p,int n) 】 入口参数：（同上） 要求输入： high : 低通滤波器的数据截止上限 核心语句： 低通滤波器： for(i=0;in;i++) { if(*(p+i)-high0.00001) printf(%d , %lf\n,i,high); else printf(%d , %lf\n,i,*(p+i)); } 作用： 将每一个数据与截止上限进行比较，小于等于上限的直接输出，大于上线的数据进行截止，输出截止上限的数值，即为低通滤波。 2.7具有快速响应的数字滤波器