多项式平滑公式系数程序(1).m

来自「本人利用MATLAB处理了并总结了工程上经常用到的“信号处理”的源代码」· M 代码 · 共 18 行

%平滑参数求解
m=input('m=');					%键盘输入m
k=input('k=');					%键盘输入k，k-1次多项式
n=2*m+1;
x=zeros(n,k);
for I=-m:m							%构造X矩阵行循环
	for J=0:k-1						%构造X矩阵列循环
		x(I+m+1,J+1)=I^J;		
    end
end
xt=x';							%求X的转置矩阵
xtx=xt*x;							%求XTX
xtx=inv(xtx);						%求(XTX)-1
xtt=xtx*xt;							%求(XTX)-1XT
b=x*xtt;							%求B=X(XTX)-1XT
b							        %屏幕输出B矩阵