永发信息网

matlab,最下二乘法求解方程,如果返回计算结果的误差呢?比如是超定的方程

答案:2  悬赏:0  手机版
解决时间 2021-12-26 21:00
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
最佳答案
这个不是开卷的嘛,问问老师然后记下来啊
全部回答
polyfit(x, y, n) 对描述n阶多项式y=f(x)的数据进行最小二乘曲线拟合 但是,你的数据只有一列,怎么拟合?
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
已知等腰三角形的一边长为4,它的其他两条边
哈里的别墅做的怎样?
请问,清哪位好心大大?兔υO?工資表,是一??物?
阅读下面一首词,然后回答问题。(8分)满江红
雷人是什么意思
把0.85、8.5%、0.805、8分之3这几个数从小到
The weather is very coid in winter (改为感
去日本用什么卡
科比退役最后一场比赛最后的演讲
读某区域等降水量图(单位:mm),回答问题。
怎么能自定义宏设置中间的按键延迟
台湾的人口和城市主要分布在A.东部地区B.南部
下列说法正确的是A.反应热是指反应过程中放出
农历吉日查询
要组成模42计数器,至少需要几个触发器
推荐资讯
凯迪仕指纹锁如何
求怀化到太原最佳的火车路线
一个女孩子说哇塞 我太有才了 什么意思
国内全案策划做的比较好的公司都有哪些呢?
分光计可以测量的物理量
请问13年锐志更换发动机电脑版需要电脑匹配吗
单选题一束平行单色光,通过双缝后,在屏上得
中国移动4G超级日租卡收费方式是什么
单选题在电磁学发展史上,提出分子电流假说的
单选题下列各句中加粗的成语使用正确的一句是
有一首歌歌词有“我的梦带我去。。。”是什么
龙泉派出所地址在哪?我要去那里办事
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?