1 / 14
文档名称:

现代信号处理例的题目及matlab代码实现.doc

格式:doc   大小:166KB   页数:14页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

现代信号处理例的题目及matlab代码实现.doc

上传人:511709291 2022/1/17 文件大小:166 KB

下载得到文件列表

现代信号处理例的题目及matlab代码实现.doc

文档介绍

文档介绍:word
word
1 / 14
word
《现代信号处理》期末考核作业
1 MATLAB仿真均值为0,方差为1的白噪声信号,信号长度N=1024,并用周期图法分别求500、1000和1500次实现的平均功率谱密度,画图。
n(*pi*n+2*pi*rand)+sin(*pi*n+2*pi*rand)+randn(1,N);
figure(2);
subplot(2,1,1);
plot(n,xn);
ylabel('幅值(V)');
xlabel('时间(s)');
title('原始信号');
Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);
f = 0:length(Pxx)-1;
subplot(212);plot(f,Pxx);
xlabel('频率/Hz');ylabel('功率谱/dB');
title('周期图N=64');
grid on;
%*******************N=64*********************
N = 128;
Nfft = 128;
word
word
5 / 14
word
n = 0:N-1;
xn = sin(*pi*n+2*pi*rand)+sin(*pi*n+2*pi*rand)+randn(1,N);
figure(3);
subplot(2,1,1);
plot(n,xn);
ylabel('幅值(V)');
xlabel('时间(s)');
title('原始信号');
Pxx = 10*log10(abs(fft(xn,Nfft).^2)/N);
f = 0:length(Pxx)-1;
subplot(212);plot(f,Pxx);
xlabel('频率/Hz');ylabel('功率谱/dB');
title('周期图N=128');
grid on;
%*******************N=64*********************
N = 1024;
Nfft =1024;
n = 0:N-1;
xn = sin(*pi*n+2*pi*rand)+sin(*pi*n+2*pi*rand)+randn(1,N);
figure(4);
subplot(2,1,1);
plot(n,xn);
ylabel('幅值(V)');
xlabel('时间(s)');
title('原始信号');
Pxx = 10*log10(abs(fft(xn,Nfft).^2)/N);
f = 0:length(Pxx)-1;
subplot(212);plot(f,Pxx);
xlabel('频率/Hz');ylabel('功率谱/dB');
title('周期图N=1024');
word
word
6 / 14
word
grid on;
实验结果如下:
数据长度N=16时的谱估计结果
数据长度N=64时的谱估计结果
word
word
7 / 14
word
数据长度N=128时的谱估计结果
数据长度N=1024时的谱估计结果
word
word
8 / 14
word
3计算并画图描绘如下函数离散时间傅立叶变换的幅度〔幅频〕特性:
n=0,1…31
n=其它

n=0,1…31
n=其它

频率X围[-2π,2π]
解答:
程序代码如下:
clear all;
close all;
N=32;
dt=1;%设置最大点数
n=[0:1:N-1];
k=[0:1:N-1];
word
word
9 / 14
word
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;%离散Fourier变换矩阵
t=n*dt;
%**************************xn信号
xn=exp(j*n*pi/5);
figure(1);
subplot(311);
plot(t,xn);
title('xn信号');
Xk=xn*WNnk;%对xn进展Fourier变换
magXk=abs(Xk);
phaXk=angle(Xk);
k=0:length(magXk)-1;
subplot(312);
plot(k/(N*dt),magXk*2/N);
title('xn信号的振幅谱')