Vẽ nhiều đường như một bản đồ nhiệt

Question

Tôi có một số lượng lớn (~ 1000) tệp từ trình ghi dữ liệu mà tôi đang cố gắng xử lý.

Nếu tôi muốn vẽ đồ thị xu hướng từ một đơn một trong những tập tin log tôi có thể làm điều đó bằng

plot(timevalues,datavalues) 

Tôi muốn để có thể xem tất cả các dòng cùng lúc trong một cách tương tự làm thế nào một dao động có một chế độ "kiên trì".

tôi có lẽ có thể rai với nhau một cái gì đó có sử dụng biểu đồ nhưng tôi hy vọng có giải pháp tồn tại từ trước hoặc thanh lịch hơn cho vấn đề này.

Answer 1

Bạn có thể thực hiện chính xác những gì bạn đang đề xuất, ví dụ: vẽ bản đồ nhiệt của tín hiệu.

Hãy xem xét những điều sau: Tôi sẽ xây dựng một tín hiệu thử nghiệm (ngoài sóng sin có biên độ khác nhau), sau đó tôi sẽ vẽ sơ đồ nhiệt qua hist3 và imagesc.

Ý tưởng là xây dựng một tín hiệu phụ trợ chỉ là vị trí của tất cả lịch sử thời gian của bạn (cả trong số x và y), sau đó trích xuất số liệu thống kê bivariate cơ bản.

% # Test signals 
xx = 0 : .01 : 2* pi; 
center = 1; 
eps_ = .2; 
amps = linspace(center - eps_ , center + eps_ , 100); 

% # the auxiliary signal will be stored in the following variables 
yy = []; 
xx_f = []; 

for A = amps 
    xx_f = [xx_f,xx]; 
    yy = [yy A*sin(xx)]; 
end 

% # final heat map 
colormap(hot) 
[N,C] = hist3([xx_f' yy'],[100 100]); 
imagesc(C{1},C{2},N') 

Bạn có thể sử dụng cũng jet colormap thay vì hot colormap cho dễ đọc. Trong biên độ sau đây là gaussian thay vì homogeneus.

Answer 2

đây là một "nguyên thủy" giải pháp mà chỉ được sử dụng hist:

%# generate some fake data 

x=-8:0.01:8; 
y=10*sinc(x); 
yy=bsxfun(@plus,y,0.1*randn(numel(x),1000)'); 
yy(randi(1000,1,200),:)= 5-randi(10)+ circshift(yy(randi(1000,1,200),:),[1 randi(numel(x),1,200)]); 

%# get plot limit parameters 

plot(x,yy) 
yl=get(gca,'Ylim'); 
xl=get(gca,'Xlim'); 
close all; 


%# set 2-d histogram ranges 

ybins=100; 
xbins=numel(x); 
yrange=linspace(yl(1),yl(2),ybins); 
xrange=linspace(xl(1),xl(2),xbins); 

%# prealocate 

m=zeros(numel(yrange),numel(xrange)); 

% build 2d hist 
for n=1:numel(x) 
    ind=hist(yy(:,n),yrange); 
    m(:,n)=m(:,n)+ind(:); 
end 

imagesc(xrange,yrange,m) 
set(gca,'Ydir','normal') 

Answer 3

Tại sao bạn không chuẩn hóa dữ liệu và sau đó thêm tất cả các dòng với nhau? Sau đó bạn có thể vẽ sơ đồ nhiệt từ datafile đơn.