1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
|
clear all
close all
clc
y = load('DATA3/ISD/isdmax0010.m');
x = load('DATA3/ParametresSignal/Sa_f20010.m');
start_intervalx = min(x);
end_intervalx = max(x);
start_intervaly = min(y);
end_intervaly = max(y);
intervals = 50;
interval_lengthx = (end_intervalx - start_intervalx)/intervals;
interval_lengthy = (end_intervaly - start_intervaly)/intervals;
for i=1:intervals
leftx = start_intervalx + interval_lengthx*(i-1);
rightx = start_intervalx + interval_lengthx*i;
indx = x >= leftx & x < rightx;
lefty = start_intervaly + interval_lengthy*(i-1);
righty = start_intervaly + interval_lengthy*i;
indy = y >= lefty & y < righty;
m_xi = mean(x(indy));
m_yi = mean(y(indx));
std_xi = std(x(indy));
std_yi = std(y(indx));
plot(mid,m_yi,'k*')
plot(mid,m_yi+std_yi,'ro')
plot(mid,m_yi-std_yi,'ro')
line( [m_xi,m_xi],[m_yi-std_yi,m_yi+std_yi], 'Color','g' );
line( [m_xi-std_xi,m_xi+std_xi],[m_yi,m_yi], 'Color','c' ); |
Partager