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) syms x1 x2 y1 y2; x1=0.4012;y1=0.8916;x2=1.5894;y2=0.8185; vpa(x1^2+y1^2-3*x1*y1^2) ans = -0.0008904940159999430449033752665855 vpa(x2^2+y2^2-3*x2*y2^2) ans = 0.0017159735499996386920429358724505 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、试用图解法求解下面的一元和二元方程，并验证得出的结果。 解：（1） ezplot(exp(-(x+1)^2+pi/2)*sin(5*x+2),[-4 4]) hold on, line([-4,4],[0,0]) （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中方程的根，并对得出的结果进行检验。 解：（1） 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 （2） syms x; y1=solve((x^2+y^2+x*y)*exp(-x^2-y^2-x*y)=0,y) y1 = - x/2 + (3^(1/2)*x*i)/2 - x/2 - (3^(1/2)*x*i)/2 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,[1;1;1;1]) Warning: Gradient must be provided for trust-region algorithm; using line-search algorithm instead. In fminunc at 367 Local minimum found. Optimization completed because the size of the gradient is less than the default value of the function tolerance. stopping criteria details x = 10.5464 111.2315 6.7823 -111.5041 试用图解法求解下面的非线性规划问题，并用数值求解算法验证结果。 解： [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) min(z(:)) ans = 5 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];