matlab实验符号计算答案,实验7 Matlab符号计算.doc

该实验旨在掌握Matlab的符号计算功能,包括定义符号对象、符号表达式运算、求极限、导数、定积分和不定积分。实验中通过具体例子展示了如何进行因式分解、表达式简化、矩阵运算、求逆矩阵、下三角阵计算,以及求解极限、导数和积分问题。
摘要由CSDN通过智能技术生成

实验7 Matlab符号计算

实验目的:

掌握定义符号对象的方法;

掌握符号表达式的运算法则以及符号矩阵运算。

掌握求符号函数极限及导数的方法。

掌握求符号函数定积分和不定积分的方法。

实验内容:

已知x=6,y=5,利用符号表达式求

分解因式。

(1);

(2)

化简表达式

(1);

(2)

已知

完成下列运算:

(1).

(2)B的逆矩阵并验证结果。

(3)包括B矩阵主对角线元素的下三角阵。

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

,求和。

用符号方法求下列积分。

(1)

(2)

(3)

详细实验内容:

1.已知x=6,y=5,利用符号表达式求

>> syms x y

>> x=6;

>> y=5;

>> z=(x+1)./(sqrt(x+3)-sqrt(y))

z =

9.16311896062463

2.分解因式。

(1);

(2)

(1)>> clear

>> syms x y

>> factor(x.^4-y.^4)

ans =

(x-y)*(x+y)*(x^2+y^2)

(2)>> clear

>> syms x

>> factor(125.*x.^6+75.*x.^4+15.*x.^2+1)

ans =

(5*x^2+1)^3

3.化简表达式

(1);(2)

(1)>> clear

>> syms bata1 bata2

>> simplify(sin(bata1).*cos(bata2)-cos(bata1).*sin(bata2))

ans =

sin(bata1-bata2)

(2)>> clear

syms x

simplify((4.*x.^2+8.*x+3)./(2.*x+1))

ans =

2*x+3

4.已知

完成下列运算:

(1).

(2)B的逆矩阵并验证结果。

(3)包括B矩阵主对角线元素的下三角阵。

(1)>> p1=sym([0,1,0;1,0,0;0,0,1]);

p2=sym([1,0,0;0,1,0;1,0,1]);

A=sym('[a,b,c;d,e,f;g,h,i]');

>> B=p1*p2*A

B =

[ d, e, f]

[ a, b, c]

[ a+g, b+h, c+i]

(2)>> simplify(inv(B))

ans =

[ (i*b-c*h)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), (-e*c-i*e+f*b+f*h)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), -(-e*c+f*b)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b)]

[ -(i*a-c*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), -(-d*c-i*d+f*a+f*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), (-d*c+f*a)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b)]

[ (a*h-b*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), (-d*b-d*h+e*a+e*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), -(-d*b+e*a)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b)]

>> simplify(B*inv(B))

ans =

[ 1, 0, 0]

[ 0, 1, 0]

[ 0, 0, 1]

(3)>> B([4,7,8])=[0,0,0]

B =

[ d, 0, 0]

[ a, b, 0]

[ a+g, b+h, c+i]

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

,求和。

(1)>> clear

>> syms x;

fun=(x.*(exp(sin(x))+1)-2.*(exp(tan(x))-1))./(sin(x)).^3;

limit(fun,x,0)

ans =

-1/2

(2)>> clear

syms x;

fun=(sqrt(x)-sqrt(acos(x)))./sqrt(x+1);

limit(fun,x,1)

ans =

1/2*2^(1/2)

(3)>> clear

syms x;

y=(1-cos(2.*x))./x;

>> dx1=diff(y)

dx1 =

2*sin(2*x)/x-(1-cos(2*x))/x^2

>> diff(y,2)

ans =

4*co

  • 4
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值