D=[1,1;2,4];
 B=[10,32]';
 u0=[0,0]';
 alpha=10;
 beta=10;
 deltat=0.0001;
 maxN=100000;
 V=u0;
 phi=u0;
 deltaut=0.0001;
 deltau=10;
 for k=1:maxN
     V=beta*u0;
     phi=alpha*(D*V-B);
     u1=(1-deltaut)*u0-deltaut*D'*phi;
     tmp=u1-u0;
     deltau=tmp'*tmp;
     if deltau<=eps
         break;
     end
     u0=
 
u1;
 end
 disp(k)
 disp(V)

---

相關推薦