tikhonov(U,s,V,b,lambda,x_0) Error
3 views (last 30 days)
Show older comments
Daniel Vasilaky
on 5 Jul 2013
Commented: ANA ROSHAN
on 29 Aug 2019
I cant get tikhonov(U,s,V,b,lambda,x_0) to run even for a simple example given in Christian Hansen Package. Here it is including the error:
EDU>> clear
EDU>> A= [
0.15 0.1
0.16 0.1
2.02 1.3
];
b =[
0.27
0.25
3.33];
x=[
1;1
];
b=A*x;
db=[.01; -.032; .01];
b1=b+db;
x2=pinv(A)*b1;
x_0=[0;0];
lambda=0.001;
[U,s,V]=svd(A);
[x_lambda,rho,eta]=tikhonov(U,s,V,b,lambda,x_0);
[x_lambda,rho,eta]=tikhonov(U,s,V,b1,lambda,x_0);
Index exceeds matrix dimensions.
Error in tikhonov (line 66)
if (nargin==6), omega = V\x_0; omega = omega(1:p); end
%******************************************************************
When I try it on
A=hilb(12)
x=[1;1;1;1;1;1;1;1;1;1;1;1
after computong
b=A*x
and
U,s,V
x_0=[0;0;0;0;0;0;0;0;0;0;0;0]
I get NaN for output.
Not sure what I'm doing wrong. Can someone give a simple example that runs?
Thank you
Accepted Answer
More Answers (1)
Yishi
on 21 Jan 2015
That works, I did not know either, but thanks for posting this comment
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!