发布时间 : 星期三 文章MATLAB符号运算练习题更新完毕开始阅读1b8f02cd5dbfc77da26925c52cc58bd63186933d
第五章 符合运算 练习题
3
2
1.求符号函数f=ax+by+cx+d分别对x,y进行三次微分；对 y进行定积分和不定积分，对y的定积分区间为（0，1）；对y趋向于1求极限。
>> syms x y
>> f=sym('a*x^3+b*y^2+c*x+d'); >> diff(f,x,3) ans = 6*a
>> diff(f,y,3) ans = 0
>> int(f,y,0,1) ans =
a*x^3 + c*x + b/3 + d >> int(f,y) ans =
(b*y^3)/3 + (a*x^3 + c*x + d)*y
>> limit(f,y,1) ans =
a*x^3 + c*x + b + d
2. 已知f=1/(1+x^2),g=sin(y),求复合函数f(g(y)).
>> syms x y
>> f=sym('1/(1+x^2)'); >> g=sym('sin(y)'); >> compose(f,g,x,y) ans =
1/(sin(y)^2 + 1)
?x2?2x?1?0?x?3z?4?3.求三元非线性方程组
?y*z??1的解。 ?>> f1=sym('x^2+2*x+1'); >> f2=sym('x+3*z-4'); >> f3=sym('y*z+1'); >> [x,y,z]=solve(f1,f2,f3)
x =
-1
y =
-3/5
z =
5/3
?dy?dx?z?cosx?4.解方程组?dz当y(0)=1,z(0)=0时，求微分方程组的解。
?dx?y?1??>> [y,z]=dsolve('Dy-z=cos(x), Dz+y=1','y(0)=1','z(0)=0','x') y =
sin(x)/2 + (x*cos(x))/2 + 1
z =
-(x*sin(x))/2
5.求级数1?122?132???1k2??和1+x+x2+…+xk+…的和。
>> syms k
>> symsum(1/(k^2),k,1,inf)
ans =
pi^2/6
>> syms x k
>> symsum(x^k,k,0,inf) ans =
piecewise([1 <= x, Inf], [abs(x) < 1, -1/(x - 1)])
??xdx6计算积分 21（1+x）
?>> f=sym('x^(1/2)/(1+x)^2'); >> int(f,1,inf) ans =
pi/4 + 1/2
?? xdz7计算积分 1 2 1+z
?>> syms x z
>> f=sym('x/(1+z^2)'); >> int(f,z,1,inf) ans =
(pi*x)/4
8求函数 1 ? 2 x 3 ? 3 1 ? 3 x ? x 2 的5阶泰勒级数展开式 x ?
>> syms x a
>>f=sym('(1-2*x+x^3)^(1/2)-(1-3*x+x^2)^(1/3)'); >>taylor(f,x,a,'order',5) ans =
(a^3 - 2*a + 1)^(1/2) - (a - x)^4*(((5*(3*a^2 - 2)*((3*(a^3 - 2*a + 1)^(1/2))/2 - (3*(3*a^2 - 2)*(3*a*(a^3 - 2*a + 1)^(1/2) - (3*a^2 - 2)^2/(4*(a^3 - 2*a + 1)^(1/2))))/(4*(a^3 - 2*a + 1))))/(6*(a^3 - 2*a + 1)) - (3*a^2 - 2)/(4*(a^3 - 2*a + 1)^(1/2)) + (3*a*(3*a*(a^3 - 2*a + 1)^(1/2) - (3*a^2 - 2)^2/(4*(a^3 - 2*a + 1)^(1/2))))/(2*(a^3 - 2*a + 1)))/(4*(a^3 - 2*a + 1)) + ((2*((2*(2*a - 3)^2)/(9*(a^2 - 3*a + 1)^(2/3)) - (2*(a^2 - 3*a + 1)^(1/3))/3))/(3*(a^2 - 3*a + 1)) + (8*(2*a - 3)*((2*a - 3)/(9*(a^2 - 3*a + 1)^(2/3)) - (5*((2*(2*a - 3)^2)/(9*(a^2 - 3*a + 1)^(2/3)) - (2*(a^2 - 3*a + 1)^(1/3))/3)*(2*a - 3))/(6*(a^2 - 3*a + 1))))/(9*(a^2 - 3*a + 1)))/(4*(a^2 - 3*a + 1))) - ((3*a^2 - 2)/(2*(a^3 - 2*a + 1)^(1/2)) - (2*a - 3)/(3*(a^2 - 3*a + 1)^(2/3)))*(a - x) - (a^2 - 3*a + 1)^(1/3) - (((2*a - 3)/(9*(a^2 - 3*a + 1)^(2/3)) - (5*((2*(2*a - 3)^2)/(9*(a^2 - 3*a + 1)^(2/3)) - (2*(a^2 - 3*a + 1)^(1/3))/3)*(2*a - 3))/(6*(a^2 - 3*a + 1)))/(3*(a^2 - 3*a + 1)) + ((3*(a^3 - 2*a + 1)^(1/2))/2 - (3*(3*a^2 - 2)*(3*a*(a^3 - 2*a + 1)^(1/2) - (3*a^2 - 2)^2/(4*(a^3 - 2*a + 1)^(1/2))))/(4*(a^3 - 2*a + 1)))/(3*(a^3 - 2*a + 1)))*(a - x)^3 + (a - x)^2*((3*a*(a^3 - 2*a + 1)^(1/2) - (3*a^2 - 2)^2/(4*(a^3 - 2*a + 1)^(1/2)))/(2*(a^3 - 2*a + 1)) + ((2*(2*a - 3)^2)/(9*(a^2 - 3*a + 1)^(2/3)) - (2*(a^2 - 3*a + 1)^(1/3))/3)/(2*(a^2 - 3*a + 1)))
xey9.计算下列函数的倒数：z?2，求zx’, zy’。
y>> syms x y
>> z=sym('x*exp(y)/(y^2)'); >> diff(z,x)
ans =
exp(y)/y^2
>> diff(z,y)
ans =