|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
下面是我自己编的用最小二乘法拟合曲线的程序
%用最小二乘法拟合曲线
function C=isploy(X,Y,N)
n=length(X);
A=zeros(N+1);
B=zeros(1,N+1);
for a=0:N
for b=0:N
A(b+1,a+1)=sum(X.^(a+b)) %构建系数矩阵
end
end
for c=0:N
B(c+1)=sum((X.^c).*Y)
end
B=B';
C=A\B;
运行后在工作空间会把每一次运行的A和B都显示出来,这是为什么啊?
下面是书上给的程序,没有问题
function C=isploy1(X,Y,M)
n=length(X);
B=zeros(1,M+1);
F=zeros(n,M+1);
for k=1:M+1
F(:,k)=X'.^(k-1);
end
A=F'*F;
B=F'*Y';
C=A\B;
C=flipud(C);
高手指教!先谢过了.
[ 本帖最后由 lxq 于 2006-10-10 23:19 编辑 ] |
|