完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
function [x,z]=fun5(A,B,C,N)
[m,n]=size(A); flag=1; z=0; while flag if C(:)>=0 flag=0; x=zeros(1,n); for i=1:m x(N(i))=B(i); end z=z; else for i=1:n if C(i)<0&A(1:m,i)<=0 disp('问题不可解'); flag=0; break; end end if flag temp=0; for i=1:n if C(i) inb=i; end end s=zeros(1,m); for i=1:m if A(i,inb)>0 s(i)=B(i)/A(i,inb); end end M=1000; for i=1:m if s(i)>0&s(i) outb=i+1; end end A(outb,:)=A(outb,:)/A(outb,inb); for i=1:m if i~=outb A(i,:)=A(i,:)-A(outb,:)*A(i,inb); B(i)=B(i)-B(outb)*A(i,inb); end end for i=1:n C(i)=C(i)-A(outb,:)*C(inb); %老提示这里维数不匹配,但这里C矩阵是1行n列,A矩阵这里是下标outb对应的行,也是一行n列 end z=z-B(outb)*C(inb); for i=1:m if i==outb N(i)=inb; end end end end end |
|
相关推荐
1个回答
|
|
帮顶!!!等高手过来!!!!!!
|
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
200圆!求助大佬给一份VSG并网和离网模式的simulink仿真
1215 浏览 0 评论
MATLAB(3)--矩阵的引用(sub2ind、ind2sub、reshape函数使用)
2034 浏览 0 评论
2377 浏览 0 评论
4068 浏览 1 评论
18222 浏览 0 评论
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-11-1 06:55 , Processed in 0.593380 second(s), Total 74, Slave 58 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号