msxindexwriwavelet.m

用matlab实现股票分析中的WRI指标

clear all
nn=1075 %total plots 
mm=0
% mm=0  %delay plots of initial satate
delay=0
aa1=xlsread('e:\datafix\masangx\compdefs_temp.xls')
% aa1=xlsread('e:\课题\data\dafosi\strain-temperature307.xls')
% aa2=aa1(1:nn,16)
aa3=aa1(1:nn,6)


% aa2=medfilt1(aa2)
% aa3=medfilt1(aa3)

[C2,L2] = wavedec(aa3,5,'sym5')
aa5=wrcoef('a',C2,L2,'sym5',5)

%求残差
dd=aa3-aa5
zz3=dd.'

tt=1:nn-mm-delay

% macd(tt)=0
vdd1(tt)=0
vdd2(tt)=0
trma(tt)=0
% n为周期点数
n=144
m=6
k=5   % k为所取最大值或最小值进行平均的个数
kk=6   % 减kk是为消除指标的钝化，一般取10以下

for tt=(n+kk+1):m:nn-mm-delay
       sort1=sort(dd(tt-n-kk:tt-1-kk))
    
    vdd3(tt-n)=sum(sort1(n-k+1:n))/k %k个最大值平均
     vdd4(tt-n)=sum(sort1(1:k))/k  %k个最小值平均

    vdd2(tt)=100*(vdd3(tt-n)-dd(tt))/(vdd3(tt-n)-vdd4(tt-n)) %wri指标
   
end


%显示
subplot(2,1,1)
plot(dd,'r')
tt=n+1:nn
subplot(2,1,2)  %指标显示
plot(tt,vdd2(n+1:nn),'b')