function [Rs,Ps]=move_corr(X,Y,a)
%%%--- this function is for moving correlation
%%%--- X,Y is the imput data, a is the length of the moving window
%%%--- Rs is the spwarman correlation coeffcient, Ps is the p_vaue
read_segy_file_legacy
Function reads a SEG-Y file and outputs a
seismic structure
Last updated: April 6, 2007: fix bug in reading
of user-specified header
function [R,k,b] = msc(A)
% 多元散射校正
% 输入待处理矩阵，通过多元散射校正，求得校正后的矩阵
%% 获得矩阵行列数
[m,n] = size(A);
%% 求平均光谱
M = mean(A,2);
%% 利用最小二乘法求每一列的斜率k和截距b
for i = 1:n
a = polyfit(M,A(:,i),1);
if i == 1
k = a(1);
b = a(2);
...
How will future generations refer to our times? Will it be known as one
of space exploration, genetics, atomic energy or computing? Possibly,
but I think it is more likely to be ‘The age of communications’. Not since
printed books and newspapers were first introduced has there been
such an explosi ...
Over the past ten years there has been a revolution in the devel-
opment and acceptance of mobile products. In that period, cel-
lular telephony and consumer electronics have moved from the
realm of science fiction to everyday reality. Much of that revolu-
tion is unremarkable – we use wireless, in ...
When thinking about mobile radio engineers there is a tendency to
assume that the engineering function relates solely to the technical
aspects of the network, such as the equipment design or the network
design. That is certainly a key part of the role of a mobile radio engineer.
However,increasingly ...
When thinking about mobile radio engineers there is a tendency to
assume that the engineering function relates solely to the technical
aspects of the network, such as the equipment design or the network
design. That is certainly a key part of the role of a mobile radio engineer.
However,increasingly ...
function y=lagr(x0,y0,x)
%x0,y0为节点
%x是插值点
n=length(x0);
m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end
