电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

数字信号处理第三版用MATLAB上机实验

数字信号处理第三版用MATLAB上机实验_第1页
1/16
数字信号处理第三版用MATLAB上机实验_第2页
2/16
数字信号处理第三版用MATLAB上机实验_第3页
3/16
实验二:时域采样与频域采样 一、 时域采样 1.用 MATLAB 编程如下: %1 时域采样序列分析fs=1000 A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=1000;T=1/fs; n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn); subplot(3,2,1);stem(n,xn);xlabel('n,fs=1000Hz'); ylabel('xn');title('xn'); subplot(3,2,2);plot(n,abs(Xk));xlabel('k,fs=1000Hz'); title('|X(k)|'); %1 时域采样序列分析fs=200 A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=200;T=1/fs; n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn);subplot(3,2,3);stem(n,xn);xlabel('n,fs=200Hz'); ylabel('xn');title('xn'); subplot(3,2,4);plot(n,abs(Xk));xlabel('k,fs=200Hz'); title('|X(k)|'); %1 时域采样序列分析fs=500 A=444.128; a=222.144; w=222.144; ts=64*10^(-3); fs=500;T=1/fs; n=0:ts/T-1; xn=A*exp((-a)*n/fs).*sin(w*n/fs); Xk=fft(xn); subplot(3,2,5);stem(n,xn);xlabel('n,fs=500Hz'); ylabel('xn');title('xn'); subplot(3,2,6);plot(n,abs(Xk));xlabel('k,fs=500Hz'); title('|X(k)|'); 2.经调试结果如下图: 020406080-2000200n,fs=1000Hzxnxn02040608005001000k,fs=1000Hz|X(k)|051015-2000200n,fs=200Hzxnxn0510150100200k,fs=200Hz|X(k)|010203040-2000200n,fs=500Hzxnxn0102030400500k,fs=500Hz|X(k)| 实验结果说明:对时域信号采样频率必须大于等于模拟信号频率的两倍以上,才 能使采样信号的频谱不产生混叠.fs=200Hz 时,采样信号的频谱产生了混叠,fs=500Hz 和 fs=1000Hz 时,大于模拟信号频率的两倍以上,采样信号的频谱不产生混叠。 二、 频域采样 1.用MATLAB 编程如下: %频域采样定理验证 M=26;N=32;n=0:M; n1=0:13;x1=n1+1; n2=14:26;x2=27-n2; x=[x1,x2];Xk=fft(x,512); X32k=fft(x,32); k=0:511;w=(pi/512)*k; subplot(321);stem(n,x);xlabel('n'); ylabel('xn');axis([0,31,0,15]); subplot(322);plot(w,abs(Xk));xlabel('k'); ylabel('|X(k)|');axis([0,1,0,200]) X16k=X32k(1:2:N); x32n=ifft(X32k);x16n=ifft(X16k,16); k1=0:31;k2=0:15; subplot(323);stem(k1,abs(X32k));xlabel('k'); ylabel('X32k');axis([0,31,0,200]); subplot(325);stem(k2,abs(X16k));xlabel('k'); ylabel('|X(k)|');axis([0,15,0,200]) n=0:31; subp...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

数字信号处理第三版用MATLAB上机实验

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部