小波分析如何将数值加密，有阴影

whmooncool

clear;
load zhunian.txt
sst = zhunian;
%------------------------------------------------------ Computation
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance) ;
n = length(sst);
dt =1/1;
time = [0:length(sst)-1]*dt + 1960.0 ;  
xlim = [1960,2001];  
pad = 1;      
dj = 1/4;   
s0 = 2*dt;   
j1 =fix((log(n*dt/s0)/log(2))/dj);  
lag1 = 0.32;  
mother = 'Morlet';
% Wavelet transform:
[wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);
power = (abs(wave)).^2 ;      
realpart=real(wave);
modulus=abs(wave);  
phase=atan2(imag(wave),real(wave));
%------------------------------------------------------ Plotting
% Plot wavelet coefficients of time series
subplot('position',[0.10 0.30 0.80 0.60])
levels = [1,2,4,8,16] ;
%level1 = [0,2,4,8];
v=[0,1,2,4];
Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period))));
[c,h]=contour(time,log2(period),realpart,log2(levels),'-');  
%[c,h]=contour(time,log2(period),realpart,level1,'k-');
clabel(c,h,v);
xlabel('Time (year)')
ylabel('Period (years)')
%title('a) Morlet Wavelet Real Part of the Temperature Anomaly Series')
set(gca,'XLim',xlim(:))
set(gca,'YLim',log2([min(period),max(period)]), ...
    'YDir','default', ...
'YTick',log2(Yticks(:)), ...
'YTickLabel',Yticks)
hold on
levels = [0.0625,0.125,0.25,0.5,1] ;
%level1=[-8,-4,-2,0];
v=[-2,-4];
[c,h]=contour(time,log2(period),realpart,log2(levels),':');
%[c,h]=contour(time,log2(period),realpart,level1,'k:');
clabel(c,h,v);
hold off

我现在想把图（无阴影）的图修改成有阴影的，无阴影的图值只有0,1,2，怎样能改成每隔0.2显示一条等值线？请高手们指导一下

edwardli

matlab 的图像属性设置窗口就可以

whmooncool

edwardli 发表于 2012-2-29 12:00

                               
登录/注册后可看大图

matlab 的图像属性设置窗口就可以

可不可以问下，那个matlab 的图像属性设置窗口在哪呢？谢谢

whmooncool

edwardli 发表于 2012-2-29 12:00

                               
登录/注册后可看大图

matlab 的图像属性设置窗口就可以

可不可以具体说下操作的地方，我对matlab还不熟悉。谢谢~

edwardli

whmooncool 发表于 2012-3-1 09:48

                               
登录/注册后可看大图

可不可以具体说下操作的地方，我对matlab还不熟悉。谢谢~

生成图片的窗口 菜单栏第二行最后一个按钮

edwardli

edwardli 发表于 2012-3-2 11:17

                               
登录/注册后可看大图

生成图片的窗口 菜单栏第二行最后一个按钮

好几种办法吧。。。好像程序中的counter语句就是画等值线的，你可以看看帮助。另外你的变量窗口会产生好多变量，选中realpart这个，在变量窗口的右上方，有一绘图按钮。一眼明了

whmooncool

edwardli 发表于 2012-3-2 11:33

                               
登录/注册后可看大图

好几种办法吧。。。好像程序中的counter语句就是画等值线的，你可以看看帮助。另外你的变量窗口会产生好多 ...

谢谢啦，我再试试~

liqix

楼主的这个程序与论坛中的其他程序有什么区别吗？

向毛主席保证

请问楼主lag1=0.32是如何选出的？这个值有什么特殊选择要求吗？