数字信号实验三.docx
- 文档编号:26466232
- 上传时间:2023-06-19
- 格式:DOCX
- 页数:12
- 大小:856.87KB
数字信号实验三.docx
《数字信号实验三.docx》由会员分享,可在线阅读,更多相关《数字信号实验三.docx(12页珍藏版)》请在冰豆网上搜索。
数字信号实验三
数字信号处理实验三:
抽样
一、正弦信号混叠
(a)
>>n=[0:
80];
>>y=sin(2*pi*300/8000*n);
>>stem(n,y)
(b)plot(n,y)
(c)
M文件:
functionv=SIN(a,b,n,p)
nn=[0:
n];
v=sin(2*pi*a/b*nn);
subplot(410+p),stem(nn,v)
图一程序:
>>v=SIN(100,8000,80,1);
>>v=SIN(225,8000,80,2);
>>v=SIN(350,8000,80,3);
>>v=SIN(475,8000,80,4);
图二程序:
>>v=SIN(7525,8000,80,1);
>>v=SIN(7650,8000,80,2);
>>v=SIN(7775,8000,80,3);
>>v=SIN(7900,8000,80,4);
图三程序:
>>v=SIN(32100,8000,80,1);
>>v=SIN(32225,8000,80,2);
>>v=SIN(32350,8000,80,3);
>>v=SIN(32475,8000,80,4);
二、产生信号
T=12.5msf0=80Hz
(a)采样
M文件:
functionv=SIN(a,b,n)
nn=[0:
n-1];
v=cos(2*pi*a/b*nn);
stem(nn*1/80000,v)
v=SIN(80,80000,1000);
(b)
functionv=SIN(a,b,n)
nn=[0:
n-1];
v=cos(2*pi*a/b*nn);
plot(nn*1/80000,v)
v=SIN(80,80000,1000);
(c)
M文件:
functionfmapgplot(xa,dt)
L=length(xa);
Nifft=round(2.^round(log2(5*L)));
Xa=fft(xa,Nifft);
range=0:
(Nifft/4);
ff=range/Nifft/dt;
plot(ff/1000,abs(Xa(1+range)))
title('CONT-TIMEFOURIERTRANSFORM(MAG)')
xlabel('FREQUENCY(kHz)'),grid
pause
fmapgplot(v,5/80000);
>>v=SIN(8000,80000,1000);
>>fmapgplot(v,5/80000);
三、A/D转换
a.绘出离散时间信号
(1)l=20,n=500
>>fsim=80000;
>>fs=4000;
>>a=fsim/fs;
>>nn=0:
500;
>>l=length(nn);
>>x=cos(2*pi*500*nn/fsim);
>>y=x(1:
a:
l);
>>L=length(y);
>>n1=0:
L-1;
>>stem(n1,y)
>>title('A/D'),ylabel('x[n]'),xlabel('n')
(2)l=20,n=1500
>>fsim=80000;
>>fs=4000;
>>a=fsim/fs;
>>stem(n1,y)
>>title('A/D'),ylabel('x[n]'),xlabel('n')
>>nn=0:
1500;
>>l=length(nn);
>>x=cos(2*pi*1500*nn/fsim);
>>y=x(1:
a:
l);
>>L=length(y);
>>n1=0:
L-1;
>>stem(n1,y)
>>title('A/D'),ylabel('x[n]'),xlabel('n')
b.计算离散DTFT
(对第二个图)
>>[H,W]=dtft(y,1500);
>>plot(W,abs(H));
>>title('DTFTAD'),ylabel('X(jw)'),xlabel('w')
四、设计一个重建滤波器
fsim=80000;
f0=1000;
fs=8000;
y2=zeros(1,L);
fori=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)');
五、D/A转换
a.重建得到平滑输出
>>fcut=2*(fs/2)/fsim;
[b,a]=cheby2(9,60,fcut);
>>y2=zeros(1,L);
fori=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.绘出傅里叶变化
fmapgplot(yy,0.0000125)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 实验