matlab课堂习题2云大.pdfVIP

课堂习题 2 3x  4y  2z  10   一、求方程组 45x  5y  4z  23 的解，用三种方法。   6x  2y  3z  4 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 二、计算表达式x 3  sin(x)  3 与表达式x 2  cos(x)  5 的和、差、积、商。 简单 创建符号变量、符号表达式 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  y   0.1 （1）方程组 (1 x 2  y 2 ) 有多少个实数解？  sin(x  cos(y ))  0  Ezplot 函数，help 课本有答案（大概 6 组解） Hold on 调整坐标轴范围 （2 ）求出离x  0, y  0 最近、且满足该方程组的一个近似解。 课本题 交点 （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))) 2 2 四、在区域x , y [3,3] ，绘制z  4xex 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*