三.实验程序及结果
3.3产生信号
a,
n=0:1000;
f0=1000;
fsim=80000;
t=n./fsim;
x=cos(2*pi*f0*n/fsim);
subplot(211),stem(n,x);
title('SAMPLES'),ylabel('x(t)'),xlabel('n');
subplot(212),plot(t,x);
title('CONTINIOUS'),ylabel('x(t)'),xlabel('t')
c.
fmagplot函数m文件:
function fmagplot(xa,dt)
L=length(xa);
Nfft=round(2.^round(log2(5*L)));
Xa=fft(xa,Nfft);
range=0:(Nfft/4);
ff=range/Nfft/dt;
plot(ff/1000,abs(Xa(1:length(range))))
title('CONT-TIME FOURIER TRANSFORM(MAG)')
xlabel('FREQUENCY(kHz)'),grid;
Pause
fmagplot(x,1/fsim)
3.4 A/D转换
a.
fsim=80000;
fs=8000;
k=fsim/fs;
n=0:1000;
l=length(n);
x=cos(2*pi*1000*n/fsim);
y=x(1:k:l);
L=length(y);
n1=0:L-1;
stem(n1,y)
title('AD'),ylabel('x[n]'),xlabel('n')
b.
[H,W]=dtft(y,1000);
plot(W,abs(H));
title('DTFT AD'),ylabel('X(jw)'),xlabel('w')
3.5设计一个重建滤波器
B.
fsim=80000;
f0=1000;
fs=8000;
y2=zeros(1,L);
for i=0:L1-1;
y2(i*delta+1)=y(i+1);
end
fcut=2*(fs/2)/fsim;
[b,a]=cheby2(9,60,fcut);
x1=delta*filter(b,a,y2);
[c,d]=freqz(b,a,1000,'whole');
d(501:1000)=d(501:1000)-2*pi;
d=fftshift(d);
d=d/2/pi/1000*fsim;
c=fftshift(c);
plot(d,abs(c));ylabel('|H(j\omega)|');xlabel('FREQUENCY(kHz)');
3.6 D/A转换:
a.
y2=zeros(1,L);
for i=0:100;
y2(i*10+1)=y(i+1);
end
yy=filter(b,a,y2);
t=[0:0.001:1];
plot(t,yy)
ylabel('xr(t)'),xlabel('t')
b.
fmagplot(yy,0.0000125)
3.7测试混叠
function [x,x1,y,y1,y2]=testall(f0,fs,fsim)
n=0:1000;
fsim=80000;
dt=1/fsim;
t=n./fsim;
x=cos(2*pi*f0*n/fsim);
delta=round(fsim/fs);
L=length(x);
y1=zeros(1,L);
y1(1:delta:L)=x(1:delta:L);
y=x(1:delta:L);
L1=length(y);
n1=0:L1-1;
Nfft=round(2.^round(log2(5*L)));
Xa=fft(x,Nfft);
range=0:(Nfft/4);
ff=range/Nfft/dt;
[H,W]=dtft(y,L);
w1=W/2/pi/1000*fs;
y2=zeros(1,L);
for i=0:L1-1;
y2(i*delta+1)=y(i+1);
end
fcut=2*(fs/2)/fsim;
[b,a]=cheby2(9,60,fcut);
x1=delta*filter(b,a,y2);
M=length(x1);
Nfft1=round(2.^round(log2(5*M)));
X=fft(x1,Nfft1);
range1=0:(Nfft1/4);
ff1=range1/Nfft1/dt;
subplot(223);plot(t,x1);grid;title('Xr(t)');xlabel('t');
subplot(224);plot(ff1/1000,abs(X(1:length(range1))));title('|Xr(j\omega)|');xlabel('kHz');
subplot(221);plot(t,x);grid;title('X(t)');xlabel('t');
subplot(222);plot(ff/1000,abs(Xa(1:length(range))));title('|X(j\omega)|');xlabel('kHz');
a,
b,
f0=6kHz
f0=7kHz
f0=9kHz
f0=10kHz
f0=15kHz
本文档为【数字信号处理实验报告三】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。