每帧噪声的能量的概率密度函数

close all;  
clear all;
clc;  
data1=textread('d.txt');
% data=data.^2
data1=data1(:,1:end-1);

data2=textread('c.txt');
% data=data.^2
data2=data2(:,1:end-1);
data2=data2'

car_num = size(data1,1);%帧数
day_num = size(data1,2);%采样点的个数
% % hour_num = 0.1;%？？？？
%  
for i=1:car_num%车序号  这个地方不用变，还是一帧
    %每辆车一条曲线
    count = 1;
    u = [];
    for j=1:day_num%天序号
        u(count,1) = i;%横坐标
         u(count,2) = data2(j);%纵坐标
         u(count,3) = data1(i,j);%空间中第三个点的值
           count = count + 1;
    end
 
    %绘制每辆车的曲线
    x=u(:,1);%车  %获取横轴
    y=u(:,2);%天  %获取纵坐标
    z=u(:,3);%时段 %获取z
    colorbar
    plot3(x,y,z);  %,'bo',x,y,z
    
    grid on;
    hold on;
%     axis([0, car_num+1,0,day_num,0,hour_num]);% 坐标轴的显示范围
%     set(gca, 'XTick', 0:car_num+1); % 车
%     set(gca, 'YTick', 0:day_num); % 天
%     set(gca, 'ZTick', 0:hour_num); % 时段
end
 
xlabel('帧数');
ylabel('能量范围');
zlabel('概率密度');
title('噪声能量概率密度图');