数字信号处理实验指导书-程序

《数字信号处理实验指导书-程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、10数字信号处理实验指导书数字信号处理实验指导书课程组扬州大学信息工程学院2013年10月扬州大学信息工程学院10数字信号处理实验指导书2信号采样与重建-matlab仿真2.4实验内容实验内容一MATLAB参考程序和仿真内容%*******************************************************************%%f—余弦信号的频率%M—基2FFT幂次数N=2^M为采样点数，这样取值是为了便于作基2的FFT分析%2.采样频率Fs%*****************

2、**************************************************%functionsamples(f,Fs,M)N=2^M;%fft点数=取样总点数Ts=1/Fs;%取样时间间隔T=N*Ts;%取样总时间=取样总点数*取样时间间隔n=0:N-1;t=n*Ts;Xn=cos(2*f*pi*t);subplot(2,1,1);stem(t,Xn);axis([0T1.1*min(Xn)1.1*max(Xn)]);xlabel('t-->');ylabel('Xn');Xk=abs

3、(fft(Xn,N));subplot(2,1,2);stem(n,Xk);axis([0N1.1*min(Xk)1.1*max(Xk)]);xlabel('frequency-->');ylabel('!Xk!');%*******************************************************************%实验内容二MATLAB参考程序和仿真内容%信号取样与重构%k=1临界取样01欠取样display('Pleaseinputthevalueof

4、k');k=input('k=');wm=1;%信号带宽Ts=k*pi/wm;%取样间隔ws=2*pi/Ts;%取样角频率扬州大学信息工程学院10数字信号处理实验指导书wc=ws/2;%为简便，滤波器截止频率wc=ws/2n=-10:10;%定义取样点的数量m=fix(length(n)/2);%单边取样周期数,fix函数的功能为取整nTs=n*Ts;%计算每个取样点dt=0.05;t=-m*Ts:dt:m*Ts;%Sa(t)信号取样f=sinc(nTs/pi);%计算Sa(t)在取样点处的函数值，得到取样信号

5、%Sa(t)信号重构fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));%Sa(t)包络线f2=sinc(t/pi);subplot(311);stem(nTs,f);%绘制取样信号holdonplot(t,f2,'r:');%绘制包络线xlabel('t');title('对sa(t)信号进行取样');axis([-m*Tsm*Ts-0.51.2]);holdoffsubplot(312);h1=plot(t,f

6、a);%绘制重构信号holdonfori=-m*Ts:Ts:m*Tsft=sinc(i/pi)*sinc(wc/pi*(t-i));h2=plot(t,ft,'m:');%绘制重构信号的各个分量holdonyy=sinc(i/pi);yy1=yy:yy;plot(i,yy1,'o');endxlabel('t');title('重构信号');axis([-m*Tsm*Ts-0.51.2]);legend([h1,h2],'重构信号','分信号');holdoffsubplot(313);%计算原信号和重构信号之

7、间的误差error=abs(fa-f2);plot(t,error);axis([-m*Tsm*Tsmin(error)max(error)+0.1*max(error)]);xlabel('t');ylabel('error(t)');%End扬州大学信息工程学院10数字信号处理实验指导书4离散傅里叶变换及其快速算法—matlab仿真4.4实验内容MATLAB参考程序和仿真内容%******************************************************************

8、*%function[x]=ffts(mode,M)Nfft=2^M;x=zeros(1,Nfft);%定义一个长度为Nfft的一维全0数组ifmode==1forn=0:Nfft-1x(n+1)=sin(2*pi*n/Nfft);endend%定义一个长度为Nfft的单周期正弦序列ifmode==2forn=0:Nfft-1x(n+1)=sin(4*pi*n/Nfft);endend