tikhonov(U,s,V,b,lambda,x_0) Error
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Daniel Vasilaky
le 5 Juil 2013
Commenté : ANA ROSHAN
le 29 Août 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
2 commentaires
Réponse acceptée
Plus de réponses (1)
Yishi
le 21 Jan 2015
That works, I did not know either, but thanks for posting this comment
0 commentaires
Voir également
Catégories
En savoir plus sur Mathematics and Optimization dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!