爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: ccnasyq

[源程序] MATLAB 小波分析源程序(可用) 简要注释

  [复制链接]

新浪微博达人勋

发表于 2019-11-4 08:47:59 | 显示全部楼层
新手学习,感谢楼主,特别棒!!!加油哈哈
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-5 13:38:59 | 显示全部楼层
看看你的代码,什么样子的?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-6 09:12:29 | 显示全部楼层
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2019-11-6 16:12:27 | 显示全部楼层
{:5_250:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2019-11-6 16:42:48 | 显示全部楼层
前来学习,谢谢!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-7 11:46:49 | 显示全部楼层
查看本帖隐藏内容
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-11 10:24:41 | 显示全部楼层
谢谢楼主分享,楼主真是好人啊啊啊啊啊
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-12 17:45:22 | 显示全部楼层
好东西,学习一下
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-13 09:11:18 | 显示全部楼层
% This is for power spectrum analysis.
clear;
load QR_jp0302;
for i=1:12:1200  
     jpn(i)=sum(QR_jp0302(i:i+11))/12;
%     jpn(i)=yjp(i);
end
   for j=1:100
       x(j)=jpn(12*(j-1)+1);
   end
y = fft(x);  %Fourier transform, default length = N
N = length(y);
y(1) = [];
power = abs(y(1:N/2)).^2;  
nyquist = 1/2;               
freq = (1:N/2)/(N/2)*nyquist;
period = 1./freq;
[mp,index] = max(power);
period(index);
[pxx,k]=psd(x,100,100,hanning(100),50,0.95,'none');
pxxgj=pxx.*((7.815+9.488)/2/3.5); psum=0;
%pxxgj=pxx.*((6.251+7.779)/2/3.5); psum=0; % alfa=0.10, r=3.5 freedum,
for i=1:length(pxx);psum=psum+pxxgj(i);end;
pxxc=psum/length(pxx); % 5% of Avg. PS line!
%[pzq,f] = periodogram(x,hanning(50),50,50-1);

subplot('position',[0.1 0.73 0.65 0.25])
plot(period,power,'r'),%axis([0 20 0 2e2])
xlabel('Period(/Years)')
ylabel('Power')
%title('Periodogram')
   
subplot('position',[0.1 0.40 0.65 0.25])%Plot Periodgram
plot(freq,power,'b') % ,grid on
%plot([0 freq],pzq,'b') % ,grid on
ylabel('Power')
xlabel('frequence')
hold on
%plot(freq,0.05*sum(power),'k-.')
hold off

%Plot PSDgram
subplot('position',[0.1 0.07 0.65 0.25])
plot(k,pxx,'m')%,axis([0 25 0 1]) % ,grid on
ylabel('Power Spectrum')
xlabel('Wave Number')
hold on
plot(xlim,pxxc+[0,0],'k:')
%plot(k,pxx(index+1)/(index+1),'r--')
hold off
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-11-13 21:44:12 | 显示全部楼层
谢谢楼主凤翔
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表