MATLAB复习.docVIP

符号计算 试写出求三阶方程x3-44.5=0正实根的程序。注意：只要正实根。 x=solve(x^3-44.5=0); for k=1:length(x) xx(k) = isreal(x(k)); end x_real=x(xx) x_real = 3.5436705307503771775536687157232 试写出求三阶方程x3-27=0的正实根的程序。注意：只要正实根。 x=solve(x^3-27=0); for k=1:length(x) xx(k) = isreal(x(k)); end x_real=x(xx) x_real = 3 试求二阶方程x2-ax+ a2=0在a0时的根。(提示：关于符号变量的假设要注意) 试求二阶方程x3+4x+5=0的负实根。（提示：关于符号变量的假设要注意） x=solve(x^3+4*x+5=0); for k=1:length(x) xx(k) = isreal(x(k)); end x_real=x(xx) x_real = -1 通过符号计算求y(t)=|sin t|的导数，然后根据此结果，求和。 （p64） syms t f=abs(sin(t)); f1=diff(f) limit(f1,t,0,left) limit(f1,t,pi/2) f1 = abs(1,sin(t))*cos(t) ans = -1 ans = 0 通过符号计算求y(t)=|cos t|的导数，然后根据此结果，求和。syms t f=abs(cos(t)); f1=diff(f) limit(f1,t,0,left) limit(f1,t,pi/2) f1 = -abs(1,cos(t))*sin(t) ans = 0 ans = NaN 计算二重积分。（p69） syms x y f=x^2+y^2; f1=int(int(f,y,1,x^2),x,1,2) f1 = 1006/105 计算二重积分。 syms x y f=x+2*y^2; f1=int(int(f,y,x,3),x,0,5) f1 = -55/3 数值计算 求函数在区间[-6,0]中的最小值点。（提示：作图观察） x1=-6; x2=0; yx=inline((cos(2*t)).^2.*exp(0.06*t.^2)-1.5.*t.*sin(2*t)+0.7.*abs(t+3)) [xn0,fval]=fminbnd(yx,x1,x2) t=x1:0.1:x2; plot(t,yx(t)),hold on ,plot(xn0,fval,r*) yx = Inline function: yx(t) = (cos(2*t)).^2.*exp(0.06*t.^2)-1.5.*t.*sin(2*t)+0.7.*abs(t+3) xn0 = -3.9451 fval = -5.2489 求函数在区间[-6,0]中的最小值点。（提示：作图观察） x1=-5; x2=5; yx=inline((sin(5*t)).^2.*exp(0.06*t.^2)-1.5.*t.*cos(2*t)+1.8.*abs(t+0.5)) [xn0,fval]=fminbnd(yx,x1,x2) t=x1:0.1:x2; plot(t,yx(t)), hold on , plot(xn0,fval,r*) yx = Inline function: yx(t) = (sin(5*t)).^2.*exp(0.06*t.^2)-1.5.*t.*cos(2*t)+1.8.*abs(t+0.5) xn0 = -1.2850 fval = -0.1860 求矩阵Ax=b的解，A为3阶魔方阵，b=。 求矩阵Ax=b的解，A为4阶魔方阵，b=。 A=magic(4) b=(1:4) rref([A,b]) x=A\b A*x x=inv(A)*b A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15