研究生院matlab作业6.docVIP

MATLAB作业6 1、用图解的方式找到下面两个方程构成的联立方程的近似解。（注：在图上可用局部放大的方法精确读出交点值） ezplot(x^2+y^2-3*x*y^2); hold on ezplot(x^3-x^2=y^2-y) 2、在图形绘制语句中，若函数值为不定式NaN ，则相应的部分不绘制出来，试利用该规律绘制的表面图，并剪切下的部分。 [x,y]=meshgrid(-1:.1:1); z=sin(x.*y); ii=find(x.^2+y.^2=0.5^2); z(ii)=NaN; surf(x,y,z) 3、试用图解法求解下面的一元和二元方程，并验证得出的结果。 ezplot(exp(-(x+1)^2+pi/2)*sin(5*x+2)) ezsurf((x^2+y^2+x*y)*exp(-x^2-y^2-x*y)) [x,y]=meshgrid(-3:0.1:3); z=(0.1*x.^2+0.1*y.^2+x.*y).*exp(-x.^2-y.^2-x.*y); [C,h]=contour(x,y,z,[-0.1:0.05:0.1]); 4、用数值求解函数求解习题3中方程的根，并对得出的结果进行检验。 syms x; x1=solve(exp(-(x+1)^2+pi/2)*sin(5*x+2)) x1 = -2/5 subs(exp(-(x+1)^2+pi/2)*sin(5*x+2),x,x1) ans = 0 f=inline(exp(-(x+1).^2+pi/2).*sin(5*x+2),x); x2=fsolve(f,0) x2 = 0.22831852178755 subs(exp(-(x+1)^2+pi/2)*sin(5*x+2),x,x2) ans = 4.750949292642762e-008 x3=fsolve(f,-1) % 选择不同的初值可以得出其他的解 x3 = -1.02831853071796 subs(exp(-(x+1)^2+pi/2)*sin(5*x+2),x,x3) ans = -5.886413288211306e-016 syms x; y1=solve((x^2+y^2+x*y)*exp(-x^2-y^2-x*y)=0,y) y1 = (-1/2+1/2*i*3^(1/2))*x (-1/2-1/2*i*3^(1/2))*x y2=simple(subs((x^2+y^2+x*y)*exp(-x^2-y^2-x*y),y,y1)) y2 = 0 0 5、试求解下面的无约束最优化问题。 f=inline([100*(x(2)-x(1)^2)^2+(1-x(1))^2+,... 90*(x(4)-x(3)^2)+(1-x(3)^2)^2+,... 10.1*((x(2)-1)^2+(x(4)-1)^2)+,... 19.8*(x(2)-1)*(x(4)-1)],x); x=fminunc(f,ones(7,1)) x = 1.0e+002 * 0.10546446798713 1.11232066767234 0.06782323911149 -1.11504746457726 0.01000000000000 0.01000000000000 0.01000000000000 试用图解法求解下面的非线性规划问题，并用数值求解算法验证结果。 [x1,x2]=meshgrid(0:0.02:1,1:0.02:2); z=x1.^3+x2.^2+4*x1+4; ii=find(x1-x2+20); z(ii)=NaN; ii=find(-x1.^2+x2-10); z(ii)=NaN; ii=find(x10); z(ii)=NaN; ii=find(x20); z(ii)=NaN; surf(x1,x2,z) function [c,ce]=exc6f4(x) ce=[]; c=[x(1)^2-x(2)+1]; f_opt=inline(x(1)^3+x(2)^2+4*x(1)+4,x); A=[-1 1]; B=2; Aeq=[]; Beq=[]; xm=[0;0]; x=fmincon(f_opt,[0;1],A,B,Aeq,Beq,xm,[],exc6f4); 试求解此线性规划问题： f=[0,0,0,0,0,1,1]; Aeq=[1 1 1 1 0 0 0; -2 1 -1 0 0 -1 1; 0 3 1 0 1 0 1]; Beq=[4; 1; 9]; xm=[0;0;0;0;0;0;0]; A=[]; B=[]; x=linprog(f,A,B,Aeq,Beq,xm)