simul.m
上传用户:zjhyt3
上传日期:2007-07-03
资源大小:89k
文件大小:2k
源码类别:

matlab例程

开发平台:

Matlab

  1. %  Script file: simul.m
  2. %
  3. %  Purpose: 
  4. %    This program solves a system of 8 linear equations in 8
  5. %    unknowns (a*x = b), using both full and sparse matrices.
  6. %
  7. %  Record of revisions:
  8. %      Date       Programmer          Description of change
  9. %      ====       ==========          =====================
  10. %    10/14/98    S. J. Chapman        Original code 
  11. %
  12. % Define variables:
  13. %   a            -- Coefficients of x (full matrix)
  14. %   as           -- Coefficients of x (sparse matrix)
  15. %   b            -- Constant coefficients (full matrix)
  16. %   bs           -- Constant coefficients (sparse matrix)
  17. %   x            -- Solution (full matrix)
  18. %   xs           -- Solution (sparse matrix)
  19. % Define coefficients of the equation a*x = b for
  20. % the full matrix solution.
  21. a = [  1.0  0.0  1.0  0.0  0.0  2.0  0.0 -1.0; ...
  22.        0.0  1.0  0.0  0.4  0.0  0.0  0.0  0.0; ...
  23.        0.5  0.0  2.0  0.0  0.0  0.0 -1.0  0.0; ...
  24.        0.0  0.0  0.0  2.0  0.0  1.0  0.0  0.0; ...
  25.        0.0  0.0  1.0  1.0  1.0  0.0  0.0  0.0; ...
  26.        0.0  0.0  0.0  1.0  0.0  1.0  0.0  0.0; ...
  27.        0.5  0.0  0.0  0.0  0.0  0.0  1.0  0.0; ...
  28.        0.0  1.0  0.0  0.0  0.0  0.0  0.0  1.0];
  29. b = [  3.0  2.0 -1.5  1.0 -2.0  1.0  1.0  1.0]'; 
  30. % Define coefficients of the equation a*x = b for
  31. % the sparse matrix solution.
  32. as = sparse(a);
  33. bs = sparse(b);
  34. % Solve the system both ways
  35. disp ('Full matrix solution:');
  36. x = ab
  37. disp ('Sparse matrix solution:');
  38. xs = asbs
  39. % Show workspace
  40. disp('Workspace contents after the solutions:')
  41. whos