没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
1、分别用短时 Fourier,Gabor 变换分析下列信号,要求提供程序,图形结果
并对它们的结果进行对比分析。采样频率 FS=1920HZ,采样长度 N=512.
Matlab 程序如下:
fs=1920;%采样频率
N=512; %采样长度
t=0:1/fs:(N-1)/fs; %时间序列
x1=(1+0.2*sin(2*pi*15*t)).*(cos(2*pi*30*t)+0.5*sin(2*pi*15*t))
+sin(2*pi*120*t);%信号
figure(1)
plot(t,x1);%画想(t)的图像
y1=fft(x1,N); %对信号进行快速 Fourier 变换
mag1=abs(y1);%求变换后的幅值
k=0:N-1;
f1=k*fs/N;
figure(2)
grid on
stem(f1,mag1);%绘制 N 点 DFI 的幅频特性图
xlabel('f1');
ylabel('幅值’);
axis([0,256,0,2*max(abs(y1))]);%x,y 的范围
grid on
figure(3)
h=window(321,'hamming');
sig=x1;
tfrstft(sig',1:512,512,h);%短时 Fourier 变换
xlabel('时间(秒)');
ylabel('频率(Hz)');
figure(4)
q=16;
h=window(211,'gauss');
h=h/norm(h);
tfrgabor(x1',128,q,h);%Gabor 变换
xlabel('时间(秒)');
ylabel('频率(Hz)');
1.1 信号的图形
图 1-1 信号时域波形图
1.2 信号 N 点的DFI幅频特性图
图 1-2 信号的幅频特性图
对 信 号 进 行 分 析 , 信 号 共 有 5 个 频 率 分 别 是 0
HZ,15HZ,30HZ,45HZ,120HZ,用火柴棍状表示出来。
1.3 短时 Fourier 变换图
图 1-3 短时傅里叶变换图
1.4 Gabor 变换图
图 1-4 Gabor 变换图
通过上面两图可以看出就显示两个频率,分别是 30HZ 和 120HZ,15HZ
比较模糊,而 0HZ 和 45HZ 的信号淹没了,经过分析原因可能一是信号的强
度不一样,显示的清晰度也不一样;二是采样频率过大显示的比较拥挤。
有以上两图知:短时 Fourier 变换和 Gabor 变换均能显示在特定时刻该信
号的频率,与 Fourier 变换相比具有定位的功能。但短时 Fourier 窗函数宽度
的选择对时间和频率分辨率的影响比较大,不能使时间分辨率和频率分辨率都
能提高。
2、分别用 wigner-ville 分布,伪 wigner-ville 分布,平滑伪 wigner-ville 分
布和 Cohen 分布分析下列信号:
其中 =0.25, =5s, 。要求提供图形结果,并对它们的结
果进行对比分析。
Matlab 程序如下:
a=0.25;
t0=5;
fs=10;
w0=1.57;
n=128;
t=0:1/fs:(n-1)/fs;
x=exp(-a*(t+t0).^2+w0*t*j)+exp(-a*(t-t0).^2+w0*t*j);
figure(1)
plot(abs(x));
x=x.';
set(gca,'xlim',[0,n])
set(gca,'xtick',[0:n/4:n])
剩余15页未读,继续阅读
资源评论
Dr.ZJH
- 粉丝: 5
- 资源: 41
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功