MATLAB基验三及求解.docVIP

实验三、符号运算的应用 1、简化。（比较simple和simplify指令的区别） clear;syms x f = (1/x^3+6/x^2+8)^(1/3) f = (1/x^3+6/x^2+8)^(1/3) simple(f) %simple不返回指定的变量时，给出所有简化的结果，最后结果是最简的形式 simplify: ((1+6*x+8*x^3)/x^3)^(1/3) radsimp: 1/x*(1+6*x+8*x^3)^(1/3) combine(trig): ((1+6*x+8*x^3)/x^3)^(1/3) factor: ((1+6*x+8*x^3)/x^3)^(1/3) expand: (1/x^3+6/x^2+8)^(1/3) combine: (1/x^3+6/x^2+8)^(1/3) convert(exp): (1/x^3+6/x^2+8)^(1/3) convert(sincos): (1/x^3+6/x^2+8)^(1/3) convert(tan): (1/x^3+6/x^2+8)^(1/3) collect(x): (1/x^3+6/x^2+8)^(1/3) mwcos2sin: (1/x^3+6/x^2+8)^(1/3) ans = (1/x^3+6/x^2+8)^(1/3) B = simple(f) %把最简的形式赋给变量B B = (1/x^3+6/x^2+8)^(1/3) C = simplify(f) %数学上的恒等式化简 C = ((1+6*x+8*x^3)/x^3)^(1/3) 2、 clear;syms x a y = sin(a*x) y = sin(a*x) A = diff(y) A = cos(a*x)*a B = diff(y,a) B = cos(a*x)*x C = diff(y,2) C = -sin(a*x)*a^2 3、计算不定积分、定积分、反常积分 ， ，? clear;syms x f1 = (x^2+1)/(x^2-2*x+2)^2 f1 = (x^2+1)/(x^2-2*x+2)^2 I =int(f1) I = 1/4*(2*x-6)/(x^2-2*x+2)+3/2*atan(x-1) f2 = cos(x)/(sin(x)+cos(x)) f2 = cos(x)/(sin(x)+cos(x)) J =int(f2,0,pi/2) J = 1/4*pi f3 = exp(-x^2) f3 = exp(-x^2) K = int(f3,0,inf) K = 1/2*pi^(1/2) K = vpa(K,4) K = .8865 4、求级数的和S, 以及前十项的部分和S1． clear;syms n S = symsum(1/n^2,1,inf) S = 1/6*pi^2 S1 = symsum(1/n^2,1,10) S1 = 1968329/1270080 s1 = vpa(S1,4) s1 = 1.550 5、求一元二次方程的根． clear;syms a b c x f = a*x^2+b*x+c f = a*x^2+b*x+c X = solve(f) X = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) 6、求微分方程的通解． dsolve(Dy = x,x) ans = 1/2*x^2+C1 7、? 求微分方程的特解． y=dsolve(D2y=Dy+x,y(0)=1,Dy(0)=0,x) y = -1/2*x^2+exp(x)-x 8、求微分方程组的通解． [x,y]=dsolve(Dx =y+x,Dy = 2*x) x = -1/2*C1*exp(-t)+C2*exp(2*t) y = C1*exp(-t)+C2*exp(2*t) 9、求下列联立方程的解 （用数值和符号运算两种方法） 法一：数值方法 A = [3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; b = [4;-3;9;-8]; x = A\b x = -1.4841 -0.6816 0.5337 -1.2429 法二：符号方法 clear;syms x y z w f1 = 3*x+4*y-7*z-12*w-4 f2 = 5*x-7*y+4*z+2*w+3 f3 = x+8