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

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

%平滑参数求解
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
b=x*inv(x'*x)*x'				    %计算平滑系数并在屏幕输出B矩阵