赌场送彩金35|临海赌场看场子
查看: 204|回复: 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 幸运飞艇官方开奖历史结果 境外彩澳洲幸运8 快乐十分天津走势图今天 在线每日稳定更新 黑龙江哈尔滨麻将 曾道点特玄机 时时彩缩水app苹果 免费老时时票软件 2018香港马会奖卷开奖结果 云南时时中奖规则