function x=super_matrix(A,b)
%edit by lancelot
%2009/09/10-- 21:00
%最小二乘法解超定方程组,最后使用列主元法解方程组
[m,n]=size(A);
AT=zeros(n,m);
for k=1:n
AT(k,:)=A(:,k);
end
AS=AT*A;
Ab=AT*b;
digits(31);
x=vpa(liezhuyuan(AS,Ab));
function x=liezhuyuan(A,b) %主函数
% A=[1,1,0,3;2,1,-1,1;3,-1,-1,3;-1,2,3,-1];b=[4;1;-1;4];liezhuyuan(A,b);
[m,n]=size(A);
nb=n+1;
Ab=[A b];
for i=1:m-1
[value,p]=max(abs(Ab(i:n,i)));
ip=p+i-1;
if ip~=i
Ab([i,ip],:)=Ab([ip,i],:);
end
value=Ab(i,i);
for k=i+1:m
Ab(k,i:nb)=Ab(k,i:nb)-(Ab(k,i)/value)*Ab(i,i:nb);
end
end
x=zeros(n,1);
x(n)=Ab(n,nb)/Ab(n,n);
for i=n-1:-1:1
x(i)=(Ab(i,nb)-Ab(i,i+1:n)*x(i+1:n,1))/Ab(i,i);
end
matlab,最下二乘法求解方程,如果返回计算结果的误差呢?比如是超定的方程
答案:2 悬赏:0 手机版
解决时间 2021-12-26 21:00
- 提问者网友:你挡着我发光了
- 2021-12-25 22:58
最佳答案
- 五星知识达人网友:低血压的长颈鹿
- 2021-12-26 00:31
这个不是开卷的嘛,问问老师然后记下来啊
全部回答
- 1楼网友:思契十里
- 2021-12-26 01:02
polyfit(x, y, n)
对描述n阶多项式y=f(x)的数据进行最小二乘曲线拟合
但是,你的数据只有一列,怎么拟合?
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
推荐资讯