Inaccuracy in solving simultaneous equations using matrix

2 次查看（过去 30 天）
Susmit Kumar Mishra
编辑： Walter Roberson 2020-6-3
So i was solving a system of n linear equations. My coefficient matrix is a tridiagonal one.
However as i was decreasing the values inside matrix or increasing n the error was increasing rapidly.
clear
n=1000;
B=(1:n);
B=B';
A=full(gallery('tridiag',n,0.341,0.232,0.741));
x=A\B;
c=A*x-B;
error=0;
for i=1:n
error=error+abs(c(i,1));
end
error
%error = 2.174626266011847e+155
Here the system is in the form Ax=B
ideally c should contain only zero.
Can anyone a suggest a method so that i can decrease the net error.
NOTE: I also tried the Thomas Algorithm even that gave an error of similar order .
  3 个评论
显示 1更早的评论
Susmit Kumar Mishra
Susmit Kumar Mishra 2020-6-3
How can one solve system of equations involving variable names in an array using symbolic toolbox.
Can you just give an example by solving above set of equations in symbolic toolbox and display the corresponding error.
Walter Roberson
Walter Roberson 2020-6-3
编辑：Walter Roberson 2020-6-3
n = 1000;
B = (1:n).';
A = sym( full(gallery('tridiag',n,0.341,0.232,0.741)) ); %11 seconds
x = A\B; %not fast!! 43 seconds
c = A*x-B;
error = sum(abs(c));
disp(error)

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Bjorn Gustavsson
Bjorn Gustavsson 2020-6-3
If you run this code-snippet:
N = round(logspace(1,3,7));
for i1 = 1:numel(N),
n = N(i1);
B=(1:n)';
A=full(gallery('tridiag',n,0.341,0.232,0.741));
[U,S,V] = svd(A);
ph(i1) = plot(diag(S),'.-','linewidth',2,'markersize',15,'color',rand(1,3));
title(n)
drawnow
end
% Pause
set(gca,'xscale','log')
% Pause
set(gca,'yscale','log')
You will see that the smallest eigenvalue is way smaller than the next smallest. Such matrices A are illconditioned and are a bit more problematic to solve. For more information and better tools to handle such problems have a look at regtools.
HTH

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

产品


版本

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by