MATLAB课堂习题2 云大.docVIP

课堂习题2 一、求方程组的解，用三种方法。 solve（方程组）、左除、左乘 (1)solve syms x y z eq1=sym(3*x+4*y-2*z-10); eq2=sym(45*x+5*y+4*z-23); eq3=sym(6*x+2*y-3*z-4); S=solve(eq1,eq2,eq3,x,y,z); disp([x,y,z]) disp([S.x,S.y,S.z]) xyz [ 58/447, 417/149, 118/149] (2)左除 A=sym([3 4 -2;45 5 4;6 2 -3]); b=sym([10;23;4]); X1=A\b X1 = 58/447 417/149 118/149 二、计算表达式与表达式的和、差、积、商。 简单 创建符号变量、符号表达式 clear syms x; f=x^3+sin(x)+3; g=x^2-cos(x)+5; a=simplify(f+g) b= simplify(f-g) c= simplify(f*g) d= simplify(f/g) 三、试用图解法回答：P236-8 （1）方程组有多少个实数解？ Ezplot函数，help课本有答案（大概6组解） Hold on 调整坐标轴范围 （2）求出离最近、且满足该方程组的一个近似解。 课本题 交点 （1）ezplot(y/(1+x^2+y^2)-0.1,[-2*pi,2*pi,-pi/3,3.5*pi])? hold?on? ezplot(sin(x+cos(y)),[-2*pi,2*pi,-pi/3,3.5*pi])?? （2）syms?x?y;? s=solve(y/(1+x^2+y^2)-0.1,sin(x+cos(y))) 可看到6个交点，及方程组有6个实数解 hold off grid on [xx,yy]=ginput1) 要解最接近x=0,y=0的解，首先将‘myfun8’另存为一个文件 function F=myfun8(x,y) F=[y/(1+x^2+y^2)-0.1;sin(x+cos(y))] end zoom on xy= ginput(1) f=fsolve(@myfun8,xy) [X,Y] = ginput(1) f=fsolve(@myfun8,[X,Y]) clear all syms x y s=solve(y/(1+x^2+y^2)-0.1,sin(x+cos(y))) 四、在区域，绘制的如图p5-6的三维（透视）网格曲面。（不得使用ezmesh） 了解ezmesh x=-3:0.1:3;y=x;[X,Y]=meshgrid(x,y); Z=4.*X.*exp(-X.^2-Y.^2); mesh(X,Y,Z) hidden off % colormap(cool), % shading interp, syms x y z % z=4.*x.*exp(-x.^2-y.^2); z=4*x*exp(-x^2-y^2); ezmesh(z,[-3,3]) hidden off