1 bai roan cn cac bro jai quyet jup

[tranthanhquyen]

day la 1 bai toan thuoc ve chyen nhang cong trinh thuy loi bon em.Noi chung la em nho cac bro jup em cai vong lap While ay.ket qua thi dung roi nhung ma no ko chiu dung lai.
Nom na thi bai toan la the nay.Bai toan toi uu,Cac bac nhin thay cai ma tran K chua,Day,cu lap nhu the den thi nao cac phan tu o hang 1 cua K deu duong thi dung lai.

function thu
K=[1 -3 -5 0 0 0 0;0 1 0 1 0 0 4;0 0 1 0 1 0 6;0 3 2 0 0 1 18]
     disp('tim cot quay va bien thay the')
     m=3
     n=2
     while 1
         K(1,:)>0
        
       for i=2:n+1
     min(K(1,:));
     if K(1,i)==min(K(1,:))
         disp('cot quay la')
          cotquay=i
         disp('bien thay the la')
          bienthaythex=i-1
     end
  end
      if K(:,cotquay)<0
          disp('phuong an nghiem ban dau khong nam tren bien')
          disp('chon lai buoc 0')
         end
disp('tim hang quay va bien bi thay the')
     K4=K(:,cotquay)
     K5=K(:,m+n+2)
     teta=K5./K4
     kq=[];
for i=2:m+1
    if (teta(i,1)>0)&(teta(i,1)~=inf)
        kq = [kq, teta(i,:)];
    end
end
teta1=kq
tetamin=min(teta1)
for i=2:m+1
    if teta(i,1)==tetamin
        disp('hang quay la')
        hangquay=i
        disp('bien bi thay the la')
        bienbithaythex=i+1
    end
end

         K(hangquay,:)
         K(hangquay,cotquay)
        hesohangcobienmoi=K(hangquay,:)/K(hangquay,cotquay)
        if hangquay==m+1
            i=1:hangquay-1
              K(i,:)=K(i,:)-K(i,cotquay)*hesohangcobienmoi
              K(hangquay,:)=hesohangcobienmoi
        else
          for i=1:hangquay-1
              K(i,:)=K(i,:)-K(i,cotquay)*hesohangcobienmoi
              K(hangquay,:)=hesohangcobienmoi
          end
         
        for i=hangquay+1:m+1
              K(i,:)=K(i,:)-K(i,cotquay)*hesohangcobienmoi
              K(hangquay,:)=hesohangcobienmoi
        end
        end
       
   K
   if K(1,:)>0
   break
   end
  
        
         end
end

[thaithien]

Post lại đi tranthanhquyen.Phần code bạn nên đưa vào thẻ code

_____

đi chứ để như vậy không xem được.Nó ra mấy cái mặt cười rồi làm sao đoán được.

[ngohaibac]

Đã thêm thẻ code vào trong bài của bạn để chương trình rõ hơn. Đề nghị bạn gõ tiếng Việt (đã được tích hợp vào trong diễn đàn) để cho tiện việc đọc và tìm hiểu.

Tất cả mã lỗi, ... đều hiện lên cửa sổ command. Bạn cho chạy nếu thấy lỗi thì sẽ hiện ở đó cả.

Chúc bạn thành công.