代码搜索结果
找到约 10,000 项符合
0 的代码
gmnewxm0.m
function [Xm0]=gmnewXm0(x,T,N)
n=length(x);
m=N+n;
xT_1=gmT(T-1,x);
xT=gmT(T,x);
q=getab(xT_1,xT);
a=q(1);
b=q(2);
for i=1:m
xm1(i)=(x(1)-b/a)*exp(-a*(i-1))+b/a;
end
xm0=xm1;
Xm0=gm_T(
gmxm0.m
function [Xm0]=gmXm0(x)
n=length(x);
sum0=0;
x1=[];
for i=1:n
sum0=sum0+x(i);
x1=[x1 sum0];
end
for i=1:n-1
b(i)=-(x1(i)+x1(i+1))/2;
Yn(i)=x(i+1);
end
b;
B=[b;ones(size(b)