第1个回答 2010-02-10
penalty=1e6;
A=[16,2,3,13;5,11,10,8;9,7,6,12;4,14,15,10];
b=[81;89;89;117];
x=A\b
%原始解是x=[1;2;3;4];
%下面限定x3=1
k=3; %代表x3
xk=1;%代表x3=1;
temp=zeros(1,length(b));
temp(k)=penalty;
%构造新的矩阵方程,当然也可以直接在A,b里修改
AA=A;
bb=b;
AA(k,:)=temp;
bb(k)=xk*penalty;
xx=AA\bb