赌场送彩金35|临海赌场看场子
查看: 112|回复: 0|关注: 0
打印 上一主题 下一主题

[已答复] 快速离散S变换(FDST)处理电流信号得到电流谱能量

[复制链接]

新手

5 麦片

财?#25442;?#20998;


050


1

主题

1

帖子

0

最佳答案
跳转到指定楼层
1#
发表于 2019-5-30 10:19:37 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
对于时间序列x(t),以T秒为采样周期采样得到N个样本,离散S变换的表达式为:


老师给了我快速离散s变换的程序:
  1. function [S,N]=strans(h1)
  2. N=length(h1);
  3. H=[fft(h1) fft(h1) fft(h1)];
  4. a1=[0:(N/2)-1,-(N/2):-1];
  5. for n=1:N
  6.     warning off MATLAB:divideByZero
  7.     W=exp(-2*pi^2*a1.^2*0.7./((n-(N/2)-1).^2)); %the frequency index is not n.but make a modificatiom
  8.     S(n,:)=ifft(H(n+(N/2):n+(N/2)+N-1).*W); %the result of st first paragraph hz ,the second time .
  9. end
  10. S(N/2+1,:)=ifft(H(N+1:2*N).*[1 zeros(1,N-1)]);   
  11. return
复制代码

以及简单的应用示例:
  1. clear all
  2. close all
  3. N=256;
  4. x1=sin(2*pi*30*((1:N)/(N)))+0*randn(1,N);%+sin(2*pi*60*((1:N)/N))
  5. x1=[zeros(1,64) x1(1:128) zeros(1,64)];
  6. figure
  7. plot(abs(fft(x1)))
  8. [S,N]=strans(x1);
  9. y=sum(S,2);
  10. y1=real(ifft(fftshift(y)));
  11. figure
  12. pcolor(abs(S))
  13. shading interp
  14. figure
  15. plot(x1)
  16. figure
  17. plot(y1)
  18. hold off
复制代码

请问应该怎样通过老师给我的s变换程序得到图片中所说的信号谱能量呢



您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表
赌场送彩金35 韦斯卡对塞尔塔 乔治布莱顿 桑普多利亚vs热那亚直播 比利亚雷亚尔欧冠成绩 庆南fc是强队吗 威尼斯恋人 都灵vs恩波利预测 莱斯特城最不可思议夺冠 佩佩赫塔菲 吉达阿赫利0-0阿尔萨德足球直播