Matlab习题

^2-1.5*x).*(x+y>1);p=p+b*exp(-y.^2-6*x.^2).*(x+y>-1).*(x+y<=1);p=p+a*exp(-0.75*y.^2-3.75*x.^2+1.5*x).*(x+y<=-1);mesh(x,y,p)10.已知连续时间Lyapunov方程为AX+XA’=?C其中A=??????????087654321,C=??????????????????165622562452252.试通过lookfor和help的帮助用MATLAB求解。%page41,ex10lookforlyapunovhelplyap>>A=[123;456;780];C=[2-5-22;-5-24-56;-22-56-16];>>X=lyap(A,C)X=1.0000-1.0000-0.0000-1.00002.00001.0000-0.00001.00007.0000>>rank(A),rank([A,b])ans=3ans=3%可见方程组唯一解>>x=A\bx=-0.4706-0.29410(3)41321511112????????????????????????????xx%Exercise2(3)>>A=[41;32;1-5];b=[1;1;1];>>rank(A),rank([A,b])ans=2ans=3%可见方程组无解>>x=A\bx=0.3311-0.1219%最小二乘近似解(4)2111121111211231234???????????????????????????????????xxxx%Exercise2(4)>>a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[123]';%注意b的写法>>rank(a),rank([a,b])ans=3ans=3%rank(a)==rank([a,b])<4说明有无穷多解>>a\b