MATLAB实验《符号计算对象与符号微积分》

 1.已知x-6，y5，利用符号表达式求Z=(x+1)/√3+X-√y,提示:定义符号常数 x=sym(6)，y-sym(5)。

x=sym('6'),y=sym('5');
Z=(x+1)/(sqrt(3+x)-sqrt(y));

 

2.分解因式。

(1)x^4+y^4    (2) 5 135

syms x y;
t=sym('5135');
a=x^4-y^4;
factor(a)
factor(t) 

 

3.化简表达式。

(1) sin B1 cos B2-cos B2 sin B1   

syms B1 B2;
f= sin(B1)*cos(B2)-cos(B1)*sin(B2);
f1=simplify(f);
disp(f1);

 

（2）4x^2+8x+3 /(2x +1)

syms x;

f=(4*x^2+8*x+3)/(2*x+1);

f1=simplify(f);

disp(f1);

 

 4.已知

完成下列运算:
(1) B=P1·P2·A。
(2)B 的逆矩阵并验证结果。
(3)包括 B 矩阵主对角线元素的下三角阵
(4)B的行列式值。 

syms a b c d e f g h i
P1=[0,1,0;1,0,0;0,0,1];
P2=[1,0,0;0,1,0;1,0,1];
A=[a,b,c;d,e,f;g,h,i];
B=P1*P2*A
B1=inv(B)
C=B*B1
B2=tril(B)
B3=det(B)

 

5.用符号方法求下列极限或导数

syms x t a y z;
f1=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/sin(x)^3;  %(1)
limit(f1)
f2=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1);              %(2)
limit(f2,x,-1 ,'right')
y=(1-cos(2*x))/x;                                   %(3)
y1=diff(y)
y2=diff(y,2)
A=[a^x t^3;t*cos(x) log(x)];                        %(4)
Ax1=diff(A,x,1)
At2=diff(A,t,2)
Axt=diff(Ax1,t)
f=(x^2-2*x)*exp(-x^2-z^2-x*z);                      %(5)
Zx=-diff(f,x)/diff(f,z)
dfxz=diff(diff(f,x),z);
x=sym('0');
z=sym('1');
eval(dfxz)                                          %运算符号返回值 

 

6、 用符号方法求下列积分

x=sym('x');
f1=1/(1+x^4+x^8);               %(1)
f2=1/(asin(x))^2/sqrt(1-x^2);   %(2)
f3=(x^2+1)/(x^4+1);             %(3)
f4=exp(x)*(1+exp(x))^2;         %(4)
F1=int(f1)
F2=int(f2)
F3=int(f3,0,inf)
F4=int(f4,0,log(2))