《科技大学数字信号处理课程设计.doc》由会员分享，可在线阅读，更多相关《科技大学数字信号处理课程设计.doc（34页珍藏版）》请在三一文库上搜索。
1、数字信号处理课程设计第1章 信号的时域分析1.1 连续信号的时域分析用Matlab产生下列信号并与人工分析结果进行比较：(1) r (t)= t u (t) -1t10(2) x (t)=1+cos10 t -1t1(3) x (t)=(5e - t - 5e -3 t ) u (t) -1t5(4) x (t)=cos(2t) cos(20t) 0t5(5) x (t)=sin(t)/t -10t=0);plot(t,x);axis(-2,6,-0.1,1.1);（2）t=-1:0.001:1; x=1+cos(10*t);plot(t,x);ylabel(x(t);xlabel(t);（3
2、）t=0:0.001:5;x=t(t=0);plot(t,x);axis(-2,6,-0.1,1.1);t=0:0.1:10;m=(t=0);n=5*exp(-t)-5*exp(-3*t);x=n.*m;plot(t,x);（4）w0=2*pi; w1=20*pi;t=0:0.001:5; x=cos(w0*t).*cos(w1*t);plot(t,x);ylabel(x(t);xlabel(t);（5）t=-10:0.1:10;m=sin(t);x=m./t;plot(t,x);ylabel(x(t);xlabel(t);1.2 离散时间序列的时域分析及信号的运算1.使用Matlab产生下列
3、序列、作图并与理论值进行比较：（1）x(n)=2(n+n0) （2） x(n)=(0.9) n sin(0.25n)+cos(0.25n) n=-4:4;x=(0.9).n;y=sin(0.25*pi*n)+cos(0.25*pi*n);m=x.*y;stem(n,m);（3）已知LTI离散系统，x(n)=1 1 1,h(n)=0 1 2 3,求y(n)x=1,1,1;h=0,1,2,3;y=conv(x,h);subplot(3,1,1);stem(0:length(x)-1,x);ylabel(x);xlabel(Time index n);subplot(3,1,2);stem(0:le
4、ngth(h)-1,h);ylabel(h);xlabel(Time index n);subplot(3,1,3);stem(0:length(y)-1,y);ylabel(y=x*h);xlabel(Time index n)；（4）已知x(t)=e 2 tu (t), y(t)=e - tu (t) , 求：x(t) * y(t)t=-10:10;u=(t=0);m=exp(-2*t);n=exp(-1*t);x=m.*u;y=n.*u;h=conv(x,y);subplot(3,1,1);stem(0:length(x)-1,x);ylabel(x(n);xlabel(Time ind
5、ex n);subplot(3,1,2);stem(0:length(y)-1,y);ylabel(y(n);xlabel(Time index n);subplot(3,1,3);stem(0:length(h)-1,h);ylabel(h(n)=x(n)*y(n);xlabel(Time index n);（5）已知信号x(t)=(1+t/2)u(t+2)-u(t-2), 求x(t+2),x(t-2),x(-t),x(2t),-x(t)t=-10:10;m=(t=2);n=(t=-2);x=(1+(t./2).*(n-m);plot(t+2,x);t=-10:10;m=(t=2);n=(t
6、=-2);x=(1+(t./2).*(n-m);plot(t-2,x);t=-10:10;m=(t=2);n=(t=-2);x=(1+(t./2).*(n-m);plot(-t,x);t=-10:10;m=(t=2);n=(t=-2);x=(1+(t./2).*(n-m);plot(2*t,x);t=-10:10;m=(t=2);n=(t=-2);x=(1+(t./2).*(n-m);plot(t,-x);第2章 信号的频域分析 2.1 利用DFT分析连续信号频谱 1用fourier函数,理论上求下列连续时间信号的频谱。(1). 三角脉冲信号x1(t)= t=-2:0.1:2;x=tripul
7、s(t,2,0);plot(t,x);syms t wxt=sym(t+1)*Heaviside(t+1)-2*t*Heaviside(t)+(t-1)*Heaviside(t-1);Fw=fourier(xt,t,w);FFw=maple(convert,Fw,piecewise);FFP=abs(FFw);ezplot(FFP,-10*pi 10*pi)axis(-10*pi 10*pi 0 1)(2). 单边指数信号x2(t)=eu(t)N=256;fs=6;Ts=1/fs;t=(0:N-1)*Ts;x=exp(-t);y=fft(x);mag=Ts*abs(y);Ws=2*pi*fs;
8、w=(0:length(y)-1)*Ws/length(y);X=1./sqrt(w.2+1);w1=w(1:length(y)/2);plot(w1,mag(1:length(y)/2);xlabel(频率（弧度/秒）);ylabel(幅度谱);z=N= num2str(N) fs= num2str(fs) 的结果;title(exp(-t)的幅度谱);2用DFT计算下列信号的频谱：（1） T0=16; N=32;T=T0/N;t=0:T:T0;x=cos(pi/8)*t+pi/4);X=1/N*fft(x,N);f=1/T/N*(-N/2:(N/2-1);stem(f,abs(fftshi
9、ft(X);xlabel(Frequency(Hz);ylabel(magnitude);（2）T0=2.5; N=20;T=T0/N;t=0:T:T0;x=cos(0.8*pi)*t+cos(0.9*pi*t);X=1/N*fft(x,N);f=1/T/N*(-N/2:(N/2-1);stem(f,abs(fftshift(X);xlabel(Frequency(Hz);ylabel(magnitude);4产生一个淹没在噪声中的信号x(t)，例如由50Hz和120Hz的正弦信号以及一个零均值的随机噪声叠加而成。确定分析长度和取样速度，计算信号的频谱；计算其功率谱密度并作图，指出50Hz和1
10、20Hz的正弦成分以及噪声；详细列出检测信号的步骤和原理。 50赫兹的频率成分对应-50和50两个坐标点，120赫兹的频率成分对应-120和120两个坐标点的频率T0=1; N=241; T=T0/N; t=0:T:T0;x=sin(100*pi*t)+sin(240*pi*t)+ randn(size(t); ;Xm=fft(x,N)/N; f=(-(N-1)/2:(N-1)/2)/N/T; stem(f,abs(fftshift(Xm); xlabel(f (Hz);ylabel(magnitude); title(幅度谱);2.2 利用DFT分析离散序列频谱1 DFT计算序列的频谱；N=
11、51; n=-(N-1)/2:(N-1)/2; %若N为偶，n = -N / 2: (N / 2 - 1); x=(0.5).n.*(n=0); X=fft(x,N); omega = 2 * pi / N * n; subplot(2,1,1); stem(n,x); ylabel(x n); xlabel(Time n); subplot(2,1,2);stem(omega,real(fftshift(X);ylabel(Xk);xlabel(Frequency (rad);2.利用DFT计算序列的频谱；确定DFT计算的各参数（抽样间隔，截短长度，频谱分辨率等）； 答：抽样间隔T=0.01
12、s，截断长度Tp=3，N=300,取512fsam=100;Tp=3; N=512; T=1/fsam;t=0:T:Tp; x=exp(-2*t);X=T*fft(x,N);subplot(2,1,1);plot(t,x);xlabel(t);title(时域波形);w=(-N/2:N/2-1)*(2*pi/N)*fsam;y=1./(j*w+2);subplot(2,1,2);plot(w,abs(fftshift(X),w,abs(y),r-.);title(幅度谱);xlabel(w);legend(理论值,计算值,0); axis(-10,10,0,1.4);3 有限长序列， 0n31
13、，分别用N=32，N=60，N=120点DFT计算其频谱。要求：(1)确定DFT计算的各参数； (2)进行理论值与计算值比较，分析各信号频谱分析的计算精度； (3)详细列出利用DFT分析离散信号频谱的步骤； (4)写出实验原理。 N=32; n=0:N-1;x=cos(3*pi/8*n);X=1/N*fft(x,N);omega=2*pi/N*(n-N/2);subplot(2,1,1); stem(omega,abs(fftshift(X); axis(-pi,pi,0,1);ylabel(Magnitude); xlabel(Frequency (rad) );subplot(2,1,2)
14、; stem(omega,angle(fftshift(X); axis(-pi,pi,-4,4);ylabel(Phase); xlabel(Frequency (rad) );1. 既然可直接由Fourier变换的定义计算连续信号的傅里叶变换，为何利用DFT分析连续信号的频谱？ 答：因为有限长序列的DFT分析的是有限长序列，特别适合数字系统，且存在着快速算法，所以仍用DFT分析连续信号的频谱。2. 若信号持续时间无限，且无解析表达式，如何利用DFT分析其频谱？答：先用窗函数截取，将信号变为有限长的信号，然后利用DFT分析其频谱。3. 在利用DFT分析连续信号频谱时，会出现哪些误差？如何克服
15、或改善这些误差？ 答：混叠现象：减少时域抽样间隔可以克服这种误差。 频率泄漏：选择适合的窗函数。 栅栏现象：补零增加序列长度4. 在利用DFT分析连续信号频谱时，如何选择窗函数？ 答：加窗处理对频谱分析主要有两个方面的影响：频谱中出现多余的高频分量；谱线变成了具有一定宽度的谱峰。对于同一个窗函数，增加长度N虽然可以减少主瓣宽度，但是不减少旁瓣泄露。为了在两者之间取得平衡，应根据实际信号特性采用合适的窗函数。每种窗函数有其自身的特性，不同的窗函数适用于不同的应用。要选择正确的窗函数，必须先估计信号的频谱成份。如若信号中有许多远离被测频率的强干扰频率分量，应选择旁瓣衰减速度较快的窗函数；如果强干扰
16、频率分量紧邻被测频率时，应选择旁瓣峰值较小的窗函数；如果被测信号含有两个或两个以上的频率成份，应选用主瓣很窄的窗函数；如果是单一频率信号，且要求幅度精度较高，则推荐用宽主瓣的窗函数。对频带较宽或含有多个频率成份的信号则采用连续采样。 绝大多数应用采用汉宁（Hanning）窗即可得到满意的结果，因为它具有较好的频率分辨率和抑制频谱泄漏的能力。5. 讨论序列后补零对频谱分析结果的影响。 答：通过补零增加序列长度时，同时也减少了频谱分析的谱线间隔，从而计算出的频谱将会显示更多的细节。6. 窗函数对频谱分辨率有何影响？如何提高频谱分辨率？ 答：输入数据通过一个窗函数相当于原始数据的频谱与窗函数频谱的卷
17、积。窗函数的频谱由一个主瓣和几个旁瓣组成，主瓣以时域信号的每个频率成份为中心。旁瓣在主瓣的两侧以一定的间隔衰减至零。FFT 产生离散的频谱，出现在FFT 每个谱线的是在每个谱线上的连续卷积频谱。如果原始信号的频谱成份与FFT 中的谱线完全一致，这种情况下采样数据的长度为信号周期的整数倍，频谱中只有主瓣。没有出现旁瓣的原因是旁瓣正处在窗函数主瓣两侧采样频率间隔处的零分量点。如果时间序列的长度不是周期的整数倍，窗函数的连续频谱将偏离主瓣的中心，频率偏移量对应着信号频率和FFT 频率分辨率的差异，这个偏移导致了频谱中出现旁瓣，所以，窗函数的旁瓣特性直接影响着各频谱分量向相邻频谱的泄漏宽度。选取适当的
18、窗函数能够提高频谱分辨率。第3章 系统分析与设计3.1 连续系统分析2研究具有以下零极点连续系统： (1)1个极点s=-0.1，增益 k =1；a=1,-0.1;k=1;b=1*k;sys=tf(s);sys=tf(b,a);p=roots(a);z=roots(b);subplot(1,2,1);plot(real(p),imag(p),*);hold on;plot(real(z),imag(z),o);xlabel(t);title(零极点图);grid %运行程序发现不稳定，分析冲激响应 %冲激响应t=0:0.01:3;h=impulse(sys,t);subplot(1,2,2);p
19、lot(h);xlabel(t);title(冲激响应h(t);Grid(2)1个极点s =0，增益k =1； a=1,0;k=1;b=1*k;sys=tf(s);sys=tf(b,a);p=roots(a);z=roots(b);subplot(1,2,1);plot(real(p),imag(p),*);hold on;plot(real(z),imag(z),o);xlabel(t);title(零极点图);grid %运行程序发现不稳定，分析冲激响应 %冲激响应t=0:0.01:3;h=impulse(sys,t);subplot(1,2,2);plot(h);xlabel(t);ti
20、tle(冲激响应h(t);grid(3)2个共轭极点s=j5，增益 k =1 ；a=conv(1 5j,1 -5j);k=1;b=1*k;sys=tf(s);sys=tf(b,a);p=roots(a);z=roots(b);subplot(1,2,1);plot(real(p),imag(p),*);hold on;plot(real(z),imag(z),o);xlabel(t);title(零极点图);grid %运行程序发现不稳定，分析冲激响应 %冲激响应t=0:0.01:3;h=impulse(sys,t);subplot(1,2,2);plot(h);xlabel(t);title
21、(冲激响应h(t);grid(4)2个共轭极点s=-0 .1 j5，增益k =1； a=conv(1 -0.1+5j,1 -0.1-5j);k=1;b=1*k;sys=tf(s);sys=tf(b,a);p=roots(a);z=roots(b);subplot(1,2,1);plot(real(p),imag(p),*);hold on;plot(real(z),imag(z),o);xlabel(t);title(零极点图);grid %运行程序发现不稳定，分析冲激响应 %冲激响应t=0:0.01:3;h=impulse(sys,t);subplot(1,2,2);plot(h);xlab
22、el(t);title(冲激响应h(t);grid(5) 零点在s =0. 5，极点在s=-0 .1j5，增益k =1； a=conv(1 -0.1+5j,1 -0.1-5j);k=1;b=1 -0.5*k;sys=tf(s);sys=tf(b,a);p=roots(a);z=roots(b);subplot(2,2,1);plot(real(p),imag(p),*);hold on;plot(real(z),imag(z),o);xlabel(t);title(零极点图);grid %运行程序发现不稳定，分析冲激响应%冲激响应t=0:0.01:3;h=impulse(sys,t);subp
23、lot(2,2,2);plot(h);xlabel(t);title(冲激响应h(t);grid%幅频特性w=0:0.01:5;h=freqs(b,a,w);subplot(2,2,3);plot(w,abs(h);xlabel(w);title(幅频特性);grid(6)零点在s =0. 5，极点在s=0 .1j5，增益k =1； a=conv(1 0.1+5j,1 0.1-5j);k=1;b=1 -0.5*k;sys=tf(s);sys=tf(b,a);p=roots(a);z=roots(b);subplot(2,2,1);plot(real(p),imag(p),*);hold on;
24、plot(real(z),imag(z),o);xlabel(t);title(零极点图);grid %运行程序发现不稳定，分析冲激响应%冲激响应t=0:0.01:3;h=impulse(sys,t);subplot(2,2,2);plot(h);xlabel(t);title(冲激响应h(t);grid%幅频特性w=0:0.01:5;h=freqs(b,a,w);subplot(2,2,3);plot(w,abs(h);xlabel(w);title(幅频特性);grid3.2 离散系统分析1. 已知一个LTI系统的差分方程为： (1)初始条件y(-1)=1,y(-2)=2，输入x(n)=u
25、(n)，计算系统的零输入响应； (2)当下面三个信号分别通过系统，分别计算系统的响应： (3)指出这是一个什么特性的系统。 N=50;n=0:N-1;a=1,-1.143,0.4128;b=0.0675,0.1349,0.0675;x=zeros(N); zi=filtic(b,a,1,2);y=filter(b,a,x,zi);x1=cos(pi/10.*n);x2=cos(pi/5.*n);x3=cos(pi/10*7.*n);y1=filter(b,a,x1,zi);y2=filter(b,a,x2,zi);y3=filter(b,a,x3,zi);figure(1)subplot(4,
26、1,1);stem(n,y);xlabel(n);ylabel(y(n);title(零输入响应);subplot(4,1,2);stem(n,y1);xlabel(n);ylabel(y(n);title(x(n)=cos(pi/10.*n);subplot(4,1,3);stem(n,y2);xlabel(n);ylabel(y(n);title(x(n)=cos(pi/5.*n);subplot(4,1,4);stem(n,y3);xlabel(n);ylabel(y(n);title(x(n)=cos(pi/10*7.*n);figure(2)H,w=freqz(b,a); plot(
27、w,abs(H);xlabel(Frequency(rad);ylabel(Magnitude); title(Magnitude response);(1) 确定DFT计算的各参数（抽样间隔，截短长度，频谱分辨率等）； 答：抽样间隔T=0.01s，截断长度Tp=3，N=300,取512fsam=100;Tp=3; N=512; T=1/fsam;t=0:T:Tp; x=exp(-2*t);X=T*fft(x,N);subplot(2,1,1);plot(t,x);xlabel(t);title(时域波形);w=(-N/2:N/2-1)*(2*pi/N)*fsam;y=1./(j*w+2);s
28、ubplot(2,1,2);plot(w,abs(fftshift(X),w,abs(y),r-.);title(幅度谱);xlabel(w);legend(理论值,计算值,0); axis(-10,10,0,1.4);第4章 数字滤波器设计4.1 IIR数字滤波器设计2设计一个数字带通滤波器，使之满足以下指标：（1） 通带（2） =, =0.4（3） =0.8dB, 60dBomegac1=0.25*pi;omegac2=0.35*pi;omegar1=0.2*pi;omegar2=0.4*pi;Ap=0.8;Ar=60;T=5*pi*10(-7);wc1=(2/T)*tan(omegac1
29、/2);wc2=(2/T)*tan(omegac2/2);wr1=(2/T)*tan(omegar1/2);wr2=(2/T)*tan(omegar2/2);B=wc2-wc1;t=wc1*wc2;norm_wr1=(wr12)-(t)/(B*wr1);norm_wr2=(wr22)-(t)/(B*wr2);norm_wc1=(wc12)-(t)/(B*wc1);norm_wc2=(wc22)-(t)/(B*wc2);if(abs(norm_wr1)-abs(norm_wr2)0 norm_wr=abs(norm_wr2);else norm_wr=abs(norm_wr1);endnorm_
30、wc=1;N=buttord(norm_wc,norm_wr,Ap,Ar,s);b_LP,a_LP=butter(N,norm_wc,s);b_BP,a_BP=lp2bp(b_LP,a_LP,sqrt(t),B);Fs=1/T;b,a=bilinear(b_BP,a_BP,Fs);w=linspace(0,2*pi,500);h=freqz(b,a,w);plot(w,20*log10(abs(h);axis(0,2*pi,-300,50);grid on;xlabel(frequency(rad);ylabel(gain(dB);4.2 窗函数的特性分析1分析并绘出各种窗函数的时域特性；矩形
31、窗N=51;w=boxcar(N);Y=fft(w,256);subplot(2,1,1);stem(0:N-1,w);xlabel(w);ylabel(y);title(时域波形);subplot(2,1,2);Y0=abs(fftshift(Y);plot(-128:127,Y0);xlabel(W);ylabel(Y0);title(频谱波形);汉明窗N=51;k=0:N-1;w=0.54-0.46*cos(2*pi*k/(N-1);Y=fft(w,256);subplot(2,1,1);stem(0:N-1,w);xlabel(w);ylabel(y);title(时域波形);subp
32、lot(2,1,2);Y0=abs(fftshift(Y);plot(-128:127,Y0);xlabel(W);ylabel(Y0);title(频谱波形);汉宁窗N=51;k=0:N-1;w=1/2*(1-cos(2*pi*k/(N-1);Y=fft(w,256);subplot(2,1,1);stem(0:N-1,w);xlabel(w);ylabel(y);title(时域波形);subplot(2,1,2);Y0=abs(fftshift(Y);plot(-128:127,Y0);xlabel(W);ylabel(Y0);title(频谱波形);布拉克窗N=51;k=0:N-1;w
33、=0.42-0.5*cos(2*pi*k/(N-1)+0.08*cos(4*pi*k/(N-1);Y=fft(w,256);subplot(2,1,1);stem(0:N-1,w);xlabel(w);ylabel(y);title(时域波形);subplot(2,1,2);Y0=abs(fftshift(Y);plot(-128:127,Y0);xlabel(W);ylabel(Y0);title(频谱波形);三角形窗N=51;w=bartlett(N);Y=fft(w,256);subplot(2,1,1);stem(0:N-1,w);xlabel(w);ylabel(y);title(时
34、域波形);subplot(2,1,2);Y0=abs(fftshift(Y);plot(-128:127,Y0);xlabel(W);ylabel(Y0);title(频谱波形);凯撒窗N=51;beta=4;w=Kaiser(N,beta);Y=fft(w,256);subplot(2,1,1);stem(0:N-1,w);xlabel(w);ylabel(y);title(时域波形);subplot(2,1,2);Y0=abs(fftshift(Y);plot(-128:127,Y0);xlabel(W);ylabel(Y0);title(频谱波形); 4.3 FIR数字滤波器设计2. 分
35、别使用矩形窗、哈明窗、海宁窗设计一个N=10的FIR低通和高通滤波器， 截频为。 1)作出各滤波器的单位脉冲响应 2)作出各滤波器的幅频响应并比较各滤波器的通带纹波和阻带纹波。 3）若当输入为，计算各滤波器的输出并作出响应波形N=10;M=N-1;wc=pi/3;n=0:M;N1=2048;n1=0:N1-1;% LP% rectangular window b1=fir1(M,wc/pi,boxcar(N); H1,w=freqz(b1,wc,512); H1_db=20*log10(abs(H1); h1=ifft(H1,N1);% hamming window b2=fir1(M,wc/
36、pi,hamming(N); H2,w=freqz(b2,wc,512); H2_db=20*log10(abs(H2);h2=ifft(H2,N1);% hanning window b3=fir1(M,wc/pi,hanning(N); H3,w=freqz(b3,wc,512); H3_db=20*log10(abs(H3);h3=ifft(H3,N1);figure(1)c=plot(w,H1_db,w,H2_db,r,w,H3_db,g);figure(2)subplot(3,1,1);stem(n1,real(h1);axis(0 25 -0.1 0.2);subplot(3,1,
37、2);stem(n1,real(h2),r);axis(0 25 -0.1 0.2);subplot(3,1,3);stem(n1,real(h3),g);axis(0 25 -0.1 0.2);% HP% rectangular window b4=fir1(M,wc/pi,high,boxcar(N+1); H4,w=freqz(b4,wc,512); H4_db=20*log10(abs(H4); h4=ifft(H4,N1);% hamming window b5=fir1(M,wc/pi,high,hamming(N+1); H5,w=freqz(b5,wc,512); H5_db=
38、20*log10(abs(H5);h5=ifft(H5,N1);% hanning window b6=fir1(M,wc/pi,high,hanning(N+1); H6,w=freqz(b6,wc,512); H6_db=20*log10(abs(H6);h6=ifft(H6,N1);figure(3)c=plot(w,H4_db,w,H5_db,r,w,H6_db,g);figure(4)subplot(3,1,1);stem(n1,real(h4);axis(0 25 -0.2 0.2);subplot(3,1,2);stem(n1,real(h5),r);axis(0 25 -0.2
39、 0.2);subplot(3,1,3);stem(n1,real(h6),g);axis(0 25 -0.2 0.2);x=1+2*cos(pi/4*n1)+cos(pi/2*n1);y1=conv(x,h1);y2=conv(x,h2);y3=conv(x,h3);y4=conv(x,h4);y5=conv(x,h5);y6=conv(x,h6);n2=0:2*N1-2;figure(5)subplot(3,2,1);stem(n2,y1);axis(0 48 -0.2 1);subplot(3,2,2);stem(n2,y2);axis(0 48 -0.2 1);subplot(3,2,
40、3);stem(n2,y3);axis(0 48 -0.2 1);subplot(3,2,4);stem(n2,y4);axis(0 48 -1 2);subplot(3,2,5);stem(n2,y5);axis(0 48 -1 2);subplot(3,2,6);stem(n2,y6);axis(0 48 -1 2); 第5章 随机信号功率谱估计1某随机信号由两余弦信号与噪声构成：xk= cos (20k) +cos (40k) + s k，sk为均值为0，方差为1的高斯白噪声。 (1) 绘出此随机信号的时域波形； (2) 试分别用周期图法、平均周期图法和Welch法分析该序列的功率谱估计
41、。 Fs = 1000; % 抽样频率t = 0:1/Fs:1; % 抽样时间xn = cos(20*pi*t) + cos(40*pi*t) + randn(size(t); %粗略地估计xn的功率谱，做N=1024点FFT： Pxx = abs(fft(xn,1024).2/1001;subplot(3,3,1);plot(t,xn); xlabel(随机信号);grid on;subplot(3,3,2);plot(0:1023*Fs/1024,10*log10(Pxx);xlabel(利用公式); grid on;window=boxcar(1001);Pxx1,F1 = period
42、ogram(xn,window,1024,Fs);subplot(3,3,3);plot(F1, 10*log10(Pxx1); xlabel (利用函数periodogram); grid on;noverlap=500;Pxx2,F2 = psd(xn, 1024,Fs, window, noverlap);subplot(3,3,4);plot(F2, 10*log10(Pxx2); xlabel (利用函数psd); grid on;noverlap=500;Pxx3,F3 = pwelch(xn, window, noverlap, 1024,Fs);subplot(3,3,5);plot(F3, 10*log10(Pxx3); xlabel (利用函数pwelch);