+ Trả lời chủ đề
Hiện kết quả từ 1 tới 2 của 2

Chủ đề: Phương pháp lặp đơn giải hệ phương trình tuyến tính X = BX + g

  1. #1
    Status : Administrator đang ẩn
    Tham gia ngày : Aug 2011
    Đến từ : Ninh Kiều - Cần Thơ
    Bài gửi : 1.450

    Phương pháp lặp đơn giải hệ phương trình tuyến tính X = BX + g

    Mã:
    Clear[n, B, X, g, X0, X1];
    ChuanB0[B_] := Max[Table[Sum[Abs[B[[i, j]]], {j, n}], {i, n}]];
    ChuanB1[B_] := Max[Table[Sum[Abs[B[[i, j]]], {i, n}], {j, n}]];
    B = Input["Nhap ma tran B"];
    g = Input["Nhap ma tran g"];
    Epxilon = Input["Nhap sai so gan dung"];
    X = Table[x, {i, n}];
    n = Length[g];
    Print["Giai he phuong trinh X = BX + g bang phuong phap lap"];
    Print[MatrixForm[X], " = ", MatrixForm[B], MatrixForm[X], "+",MatrixForm[g]];
    If[ChuanB0[B] < 1 && ChuanB1[B] < 1,
        		ϵ = 1;
        		X1 = Table[0, {i, n}];
        		While[ϵ > Epxilon,	
          			X0 = X1;
          			X1 = B.X0 + g;
          			Print["X1 = ", MatrixForm[X1]];
          			ϵ = Max[Abs[X1 - X0]];
          		];
        		Print["Nghiem cua he pttt: X = ", X1];
        	,
        		Print["Phuong phap lap khong hoi tu"];
        ];

  2. #2
    Status : Administrator đang ẩn
    Tham gia ngày : Aug 2011
    Đến từ : Ninh Kiều - Cần Thơ
    Bài gửi : 1.450
    Input:
    B = {{0, -0.06, 0.02}, {-0.03, 0, 0.05}, {-0.01, 0.02, 0}};
    g = {2, 3, 5};
    Epxilon = 0.0001;
    Output:
    Nghiem cua he pttt: X = {1.9092, 3.19496, 5.04481}

+ Trả lời chủ đề

Quyền viết bài

  • Bạn không thể gửi chủ đề mới
  • Bạn không thể gửi trả lời
  • Bạn không thể gửi file đính kèm
  • Bạn không thể sửa bài viết của mình
Trang Chủ Việc Làm Gia Sư Gia sư