Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2012-02-13
帖子: 12
声望力: 14 ![]() |
![]()
我不知道是不是你不小心,Z的赋值式中有一个写错了,不该是a^2,该是a.^2下面是完整语句
Z=[ones(36,1),t,t.^2,l,l.^2,a,a.^2,b,b.^2,t.*l,t.*a,t.*b,t.*a.*b,l.*a,l.*b,l.*a.*b,a.*b,l.*t.*a.*b]; 另外补充一下你还可能犯的另外一个拟合错误,Z不是一个满秩阵,因此你可能做的结果并不是你想要的结果,我建议你用伪逆指令代替你的右除指令,比如说: c=pinv(Z)*pcrg; 这是一个良好的编程习惯,不会引起混淆。 另外告诉你求满秩方阵的逆矩阵的函数b=inv(A); |
![]() |
![]() |