测试技术基础作业.docxVIP

2016年编程作业：假定一个包含某种故障、异常或其他特征信息的连续时间信号，将信号离散化，即令，进行采样，以DFT编程计算为基础，合理地选用一种谱分析方法，在频域中识别出预设的特征。要求作出简要论述、给出程序和频谱图。一、给定故障信号函数定义一个干扰f(t)= : 我们令A1=1；对于这个函数我们能够做离散傅里叶变换，对其频谱进行分析。二、DFT编程变换算法的理解离散傅里叶变换（Discrete Fourier Transform，缩写为DFT），是傅里叶变换在时域和频域上都呈离散的形式，将信号的时域采样变换为其DTFT的频域采样。在形式上，变换两端（时域和频域上）的序列是有限长的，而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT，也应当将其看作其周期延拓的变换。在实际应用中通常采用快速傅里叶变换计算DFT。由于是离散傅里叶变化，所以将原有的积分号用累加来表示，自变量的取值也变成了离散的形式。于是离散傅里叶变换的公式如下：F(k)==|F(k)|=三、设计思路1．通过以上公式，可以选择用C语言实现得到离散变化过后的点，然后再把得到的点进行作图便能够得到对应函数的频谱图。2．也可以利用MATLAB里面带有的算法来完成，由于MATLAB能够很方便的直接得到图形，所以选择用matlab实现。四、代码粘贴1.matlab 代码：A1=1;N=10;t=0:1/N:1-1/N;f0=50;y2=A1*sin(2*pi*f0*t);??%采样i=0:N-1;subplot(3,3,4);stem(i,y2);xlabel(n;title(A1*sin(2*pi*f0*t)函数采样?f0=50);Y2=fft(y2,N);mag=abs(Y2);subplot(3,3,5);stem(i,mag);xlabel(K);title(A1*sin(2*pi*f0*t)函数DFT?f0=50);ang=angle(Y2)*180/pi;subplot(3,3,6);plot(i,ang);xlabel(K);ylabel(角度);title(A1*sin(2*pi*f0*t)函数DFT相位图);2：C语言编写的代码：#define m 8floatA1=1;floatm=10;x1[0]=1;y1[0]=1;for(i=1;im;i++){x1[i]= A1*sin(2*pi*f0*1/i);y1[i]= A1*sin(2*pi*f0*1/i);}for(i=0;im;i++){x2[i]=0;y2[i]=0;arg=2.0*3.141592654*(double)i/(double)m;for(k=0;km;k++){cosarg=cos(k*arg);sinarg=sin(k*arg);x2[i]+=x1[k]*cosarg;y2[i]+=-y1[k]*sinarg;}}for(i=0;im;i++){x1[i]=x2[i];y1[i]=y2[i];f[i]=sqrt(x1[i]*x1[i]+y1[i]*y1[i]);printf(%f\n ,f[i]);}五、结果分析取N=10；取N=20时：幅频特性：由上面的图像可以看出，幅频特性的点有对称性，且都是相对于N/2这个点对称，并且波动有规律，靠近中线的点的值都是最大的。相频特性:此干扰信号的相频差并不会超过但是相频图的波动较大并且没有明显规律。