Các bác có thể giúp em cách tạo ma trận rồi mình hàm ED2(Pi,Gi,Pl) trên hàm M-file không? có nghĩa là ma trận Pdata,Gdata,Pldata thay đổi thì kết quả thay đổi mà khi em là thì nó bị lỗi. Em cám ơn các bác đã đọc bài này.
Code:
function ED2(Pdata,Gdata,Pldata) n=input('Nhap so to may phat= '); Pi=zeros(n,3); Gi=zeros(n,3); Pl=zeros(n,1); VT=zeros(n+1,n+1); VP=zeros(n+1,1); PT=zeros(n,n+1); KQ=zeros(n+1,1); Pload=sum(Pl(:,1).*(Pi(:,2)).^2); A=[polyder(Gi(1,:)),zeros(1,2)]; PT(1,1)=A(1,1); for k=2:n A=[polyder(Gi(k,:)),zeros(1,2)]; PT(k,k)=A(1,1); end tam=ones(n,1)-(2.*Pl.*Pi(:,2)); t=[zeros(n,n),tam]; VT=[(PT-t);[ones(1,n),0]]; VP=[(Gi(:,2));-(sum(Pi(:,2))+Pload)]; KQ=VT\-VP; disp(KQ)
Comment