数字信号处理西电Word文档格式.docx
- 文档编号:20578621
- 上传时间:2023-01-24
- 格式:DOCX
- 页数:48
- 大小:1.04MB
数字信号处理西电Word文档格式.docx
《数字信号处理西电Word文档格式.docx》由会员分享,可在线阅读,更多相关《数字信号处理西电Word文档格式.docx(48页珍藏版)》请在冰豆网上搜索。
fs/2-fs/SPL;
plot(freqb,abs(sp))
频率'
频谱幅度'
set(gca,'
yy=abs(abs(sp)-3);
[aa,freqind]=min(yy);
(freqind-SPL/2)*fs/SPL
clc;
fs=1000;
实验三:
设
,编写MATLAB程序,计算:
(1)5点圆周卷积
;
(2)6点圆周卷积
(3)线性卷积
(4)画出的
和
时间轴对齐。
a=[1,2,2];
b=[1,2,3,4];
y1=cconv(a,b,5)
y2=cconv(a,b,6)
y3=conv(a,b)
figure
(1);
subplot(311)
stem(y1);
gridon
title('
五点圆周卷积y1(n)'
);
xlabel('
n'
),ylabel('
y1(n)'
axis([06015])
subplot(312)
stem(y2);
六点圆周卷积y2(n)'
y2(n)'
subplot(313)
stem(y3);
线性卷积y3(n)'
y3(n)'
实验四:
给定因果系统:
(1)求系统函数
并画出零极点示意图。
(2)画出系统的幅频特性
和相频特性
(3)求脉冲响应
并画序列图。
提示:
在MATLAB中,zplane(b,a)函数可画零极点图;
Freqz(b,a,N)可给出
范围内均匀间隔的
点频率响应的复振幅;
Impz(b,a,N)可求
的逆变换(即脉冲响应)。
clc
a=[1,0]
b=[1,-0.9]
figure
(1)
zplane(b,a);
零极点分布图'
w=[-3*pi:
0.01:
3*pi];
[h,phi]=freqz(b,a,w);
figure
(2);
subplot(2,1,1);
plot(w,abs(h));
gridon;
幅频特性'
f/Hz'
H(f)'
subplot(2,1,2);
plot(w,phi);
相频特性'
W(f)'
数字信号处理第二次实验
1.给定模拟信号
,对其进行采样,用DFT(FFT)进行信号频谱分析。
(1)确定最小采样频率和最小采样点数。
(2)若以
秒进行采样,至少需要取多少采样点?
(3)用DFT的点数
画出信号的
点DFT的幅度谱,讨论幅度谱结果。
(4)
分别为
,能否分辨出信号的所有频率分量。
(5)在(3)和(4)的条件下做补0FFT,分析结果。
(6)在不满足最小采样点数的情况下做补0DFT,观察是否可以分辨出两个频率分量。
程序如下:
clear
closeall
%
(1)确定最小采样频率和最小采样点数
w1=4*pi;
w2=8*pi;
f1=w1/(2*pi);
f2=w2/(2*pi);
disp('
最小采样频率:
'
fs1=2*max(f1,f2);
disp(fs1);
f=f2-f1;
最小采样点数:
N=ceil(fs1/f);
disp(N);
%
(2)t=0.01ns采样
T=0.01;
fs2=1/T;
以t=0.01ns采样,最少采样点数为:
N0=fs2/f;
disp(N0);
%(3)(4)N=50,100,64,60时的幅度谱
N1=50;
N2=100;
N3=64;
N4=60;
n1=0:
N1-1;
n2=0:
N2-1;
n3=0:
N3-1;
n4=0:
N4-1;
x1=2*cos(w1*n1*T)+5*cos(w2*n1*T);
x2=2*cos(w1*n2*T)+5*cos(w2*n2*T);
x3=2*cos(w1*n3*T)+5*cos(w2*n3*T);
x4=2*cos(w1*n4*T)+5*cos(w2*n4*T);
X1=abs(fft(x1,N1));
X2=abs(fft(x2,N2));
X3=abs(fft(x3,N3));
X4=abs(fft(x4,N4));
subplot(2,2,1);
stem(n1,X1,'
.'
N=50幅度谱'
ylabel('
X1'
subplot(2,2,2);
stem(n2,X2,'
N=100幅度谱'
X2'
subplot(2,2,3);
stem(n3,X3,'
N=64幅度谱'
X3'
subplot(2,2,4);
stem(n4,X4,'
N=60幅度谱'
X4'
%(5)补0DFT
N5=200;
n5=0:
N5-1;
X5=abs(fft(x1,N5));
X6=abs(fft(x2,N5));
X7=abs(fft(x3,N5));
X8=abs(fft(x4,N5));
figure
(2)
stem(n5,X5,'
补0后N=50幅度谱'
X5'
stem(n5,X6,'
补0后N=100幅度谱'
X6'
stem(n5,X7,'
补0后N=64幅度谱'
X7'
stem(n5,X8,'
补0后N=60幅度谱'
X8'
%(6)N=2时不满足最小采样点数
N6=2;
n9=0:
N6-1;
x9=2*cos(w1*n9*T)+5*cos(w2*n9*T);
X9=abs(fft(x9,N5));
figure(3)
stem(n5,X9,'
X9'
N=2时补0后的幅度谱'
运行结果:
2.设雷达发射线性调频信号
,采样率
,采样点数
回波信号
(1)画出
的频谱。
(2)利用DFT的时延性质产生
,比较直接在时域产生和在频域产生(再变换到时域)的结果是否相同。
(3)匹配滤波的结果是
,(“
”表示线性卷积)。
分别用直接线性卷积和DFT的卷积定理求解
比较二者结果,并记录两种方法的运行时间(用tic,toc指令)。
(4)画出
figure;
plot([-0.5:
1/(fft_num):
0.5-1/(fft_num)],...
fftshift(20*log10(abs(fft(ht,fft_num)))))%将线性调频信号转换到频域并将零频搬至频谱中央
axis([-0.50.51050])
归一化频率'
幅度/dB'
h(t)频谱'
第二问
时域构造回波信号,时延通过补零实现
s_shiyu=[zeros(1,shiyan1*fs),ht,zeros(1,N-shiyan1*fs)]+[zeros(1,shiyan2*fs),ht,zeros(1,N-shiyan2*fs)];
plot([0:
2*N-1],abs(s_shiyu))
axis([02*N-102.5])
距离单元'
时域法s(t)'
fftshift(20*log10(abs(fft(s_shiyu,fft_num)))))%将回波信号转换到频域并将零频搬至频谱中央
axis([-0.50.51060])
时域法s(t)频域'
频域构造回波信号,时延通过DFT时延性质产生
L=2*N;
P=fft(ht,L);
P_shiyan1=P.*exp(-j*2*pi*fs*[0:
L-1]/L*shiyan1);
%目标1频谱
P_shiyan2=P.*exp(-j*2*pi*fs*[0:
L-1]/L*shiyan2);
%目标2频谱
s_pinyu=ifft(P_shiyan1)+ifft(P_shiyan2);
2*N-1],abs(s_pinyu))
频域法s(t)'
fftshift(20*log10(abs(fft(s_pinyu,fft_num)))))
频域法s(t)频域'
2*N-1],abs(s_shiyu-s_pinyu))
axis([02*N-101])
时域法与频域法回波之差'
第三问
tic
y_shiyu=conv(s_pinyu,conj(fliplr(ht)));
%%%时域匹配滤波
toc
y_shiyu_quchu=y_shiyu(1,N:
end);
%%%%去暂态点
2*N-1],abs(y_shiyu_quchu))
时域匹配滤波'
y_pinyu=ifft(fft(s_pinyu).*conj(P));
频域匹配滤波
2*N-1],abs(y_pinyu))
频域匹配滤波'
第四问
figure
1/L:
0.5-1/L],fftshift(20*log10(abs(fft(y_pinyu)))))
y(t)频域'
)运行结果:
Elapsedtimeis0.235671seconds.
Elapsedtimeis0.001764seconds.
>
数字信号处理上机第三次实验
1.IR滤波器设计
(1)用matlab确定一个数字IIR低通滤波器所有四种类型的最低阶数。
指标如下:
40kHz的采样率,4kHz的通带边界频率,8kHz的阻带边界频率,0.5dB的通带波纹,40dB的最小阻带衰减。
并在同一张图中画出每种w。
(2)用matlab确定一个数字IIR高通滤波器所有四种类型的最低阶数。
3500Hz的采样率,1050Hz的通带边界频率,600Hz的阻带边界频率,1dB的通带波纹,50dB的最小阻带衰减。
并在同一张图中画出每种滤波器的频率响应。
(3)用matlab确定一个数字IIR带通滤波器所有四种类型的最低阶数。
7kHz的采样率,1.4kHz和2.1kHz的通带边界频率,1.05kHz和2.45kHz的阻带边界频率,0.4dB的通带波纹,50dB的最小阻带衰减。
(4)用matlab确定一个数字IIR带阻滤波器所有四种类型的最低阶数。
12kHz的采样率,2.1kHz和4.5kHz的通带边界频率,2.7kHz和3.9kHz的阻带边界频率,0.6dB的通带波纹,45dB的最小阻带衰减。
用到的函数:
butter,buttord,cheb2ord,chebl1,cheby2,
ellip,ellipord.
(1)clc
clearall
fc=40;
fp=4;
fs=8;
rp=0.5;
rs=40;
wp=2*pi*fp/fc;
ws=2*pi*fs/fc;
Inbuttord'
[n,wc]=buttord(wp,ws,rp,rs,'
s'
[b,a]=butter(n,wc,'
low'
w=0:
0.001:
6;
[h,w]=freqs(b,a,w);
h=20*log10(abs(h));
plot(w,h,'
r-'
Incheb1ord'
[n,wpo]=cheb1ord(wp,ws,rp,rs,'
[b,a]=cheby1(n,rp,wpo,'
holdon
b-'
Incheb2ord'
[n,wso]=cheb2ord(wp,ws,rp,rs,'
[b,a]=cheby2(n,rs,wso,'
k-'
Inellipord'
[n,wc]=ellipord(wp,ws,rp,rs,'
[b,a]=ellip(n,rp,rs,wc,'
m-'
滤波器的频率响应'
legend('
巴特沃斯'
切比雪夫Ⅰ型'
切比雪夫Ⅱ型'
椭圆'
w'
h'
巴特沃斯
n=
9
wc=
0.7533
切比雪夫Ⅰ型
5
wpo=
0.6283
切比雪夫Ⅱ型
wso=
1.2069
椭圆
4
(2)
clc
fc=3500;
fp=1050;
fs=600;
rp=1;
rs=50;
high'
axis([0,6,-100,0])
12
1.7402
7
1.8850
1.2049
(3)
fc=7;
fp=[1.4,2.1];
fs=[1.05,2.45];
rp=0.4;
[b,a]=butter(n,wc,'
[b,a]=cheby1(n,rp,wpo,'
[b,a]=cheby2(n,rs,wso,'
[b,a]=ellip(n,rp,rs,wc,'
axis([0,6,-100,20])
1.23051.9250
1.25661.8850
1.10502.1437
(4)
fc=12;
fp=[2.1,4.5];
fs=[2.7,3.9];
stop'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理