MATLAB基础实验三及求解_第1页
MATLAB基础实验三及求解_第2页
MATLAB基础实验三及求解_第3页
MATLAB基础实验三及求解_第4页
MATLAB基础实验三及求解_第5页
已阅读5页,还剩7页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、实验三、符号运算的应用1、简化"壮+夕+乎+8(比较simple和simplify指令的区别)clear;syms xf = (1/xa3+6/xa2+8)a(1/3)f =(l/xx3+6/x2+8)人(1/3)simple(f) %simple不返回指定的变量时,给出所有简化的结果,最 后结果是最简的形式simplify:(l+6*x+8*xa3)/x3)"(1/3)radsimp:1/x*(l+6*x+8*x3)(1/3)combine(trig):(l+6*x+8*x3)/x3)人(1/3)factor:(l + 6*x+8*xa3)/xa3)a(1/3)expan

2、d:(1/xx3+6/xa2+8)人(1/3)combine:(l/x3+6/xx2+8)人(1/3)convert(exp):(1/xx3+6/xa2+8)人(1/3)convert(sincos):(1/xx3+6/xa2+8)人(1/3)convert(tan):(1/x"3 + 6/xa2 + 8)人(1/3)collect(x):(l/x3+6/xx2+8)人(1/3)mwcos2sin:(1/xx3+6/xa2+8)a(1/3)ans =(l/x3+6/xx2+8)人(1/3)=simple(f) %把最简的形式赋给变量b =(1/xa3 + 6/xa2 + 8)人(1/

3、3)c = simplify(f) %数学上的恒等式化简c =(l+6*x+8*xa3)/xa3)a(1/3)2、尸sinor,求心空,叟c = © 、dx da dx2clear;syms x a y = sin(a*x)y =sin(a*x)a = diff(y)a = cos(a*x)*ab = diff(y3a)b = cos(a*x)*xc = diff(y,2)-sin(a*x)*ax29.j (x2 - 2x+2)23、计算不定积分.定积分、反常积分k訂;ujoat rk/ 2 cos兀.dx j =dxjo sinx+cosxclear;syms xf1 = (xa2

4、+1)/(xa2-2*x+2)a2fl =(xa2+1)/(xx2-2*x+2)x21/4*(2*x-6)/(xx2-2*x+2)+3/2*atan(x-l)f2 = cos(x)/(sin(x)+cos(x)f2 =cos(x)/(sin(x)+cos(x)j =int(f250,pi/2)l/4*pif3 = exp(-xa2)f3 =exp(-x2)k = int(f3,0,inf)k =l/2*pix(l/2)k = vpa(k,4).88654、求级数工飞的和s,以及前十项的部分和sl n= “clear;syms ns = symsum(1/na2,15inf)l/6*pix2s1

5、 = symsum(1/na2,1,10)si =1968329/1270080s1 = vpa(s1,4)si =1 5505、求一元二次方程or,+bx + c = 0的根.f = sym(a*xa2+b*x+c=0)clear;syms a b c xf = a*xa2+b*x+cf =a*x2+b*x+cx = solve l/2/a*(-b+(b'2-4*a*c)(1/2)1/2/a*(b2-4*a*c)入(1/2)6、求微分方程y =兀的通解.dsolve(vdy = xjxjans = l/2*x2+cl7、求微分方程+x(o) = l,y(o) = o的特解.y=dso

6、lve(,d2y=dy+x,;y(0)=1,dy(0)=07x,)-l/2*x2+exp(x)-xx - y + x8、求微分方程组.的通解u =2兀x3y=dsolve(,dx =y+x,dy = 2*xf)-l/2*cl*exp(-t)+c2*exp(2*t)cl*exp(-t)+c2*exp(2*t)9、求下列联立方程的解(用数值和符号运算两种方法)3x + 4y -7z -12w = 45x - 7y + 4z + 2w = -3vx + 8z 5w =9-6x + 5y - 2z + low = -8法一:数值方法a = 3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6

7、 5 -2 10; b = 4;-3;9;-8;x = ab14841-068160 5337-12429法二:符号方法clear;syms x y z wf1 = 3*x+4*y-7*z-12*w-4 f2 = 5*x-7*y+4*z+2*w+3 f3 = x+8*z-5*w-9f4 = -6*x+5*y-2*z+10*w+8fl =3*x+4*y-7*z-12*w-4f2 =5*x-7*y+4*z+2*w+3f3 =x+8*z-5*w-9f4 =-6*x+5*y-2*z+10*w+8w x y z=solve(f 1 ,f2,f3,f4)w =-2932/2359-3501/2359-16

8、08/2359 z =1259/2359x= vpa(x;y;z;w,4)x =14846816 5337-1.24310、求极限m = limt ln(x + /z)-ln(x)l = lim力tohclear;syms x h nf1 = (log(x+h) log(x) / hfl =(log(x+h)-log(x)/hl = limit(f1,h,0)l =1/xf2 = (1 - x/n)an(1-x/n)人!1m = iimit(f2,n,inf)exp(-x)11、设a、b定义如下,试上机输出factor(a)和factor(b)的结果,并指 出那个结果才是12345678901

9、234567890的因式分解,为件么? a=symc12345678901234567890*); % 直接创建符号常量 b=sym(12345678901234567890); %将数值常量转换为符号常量, 有误差x1=factor(a)xl =(2)* (3)x2* (5)* (101)* (3803)* (3607)* (27961)* (3541)x2=factor(b)x2 = (2)11* (7)*(7324703)*(117570121)isequal(x1,x2)ans12、求函数/(x) = x4-3x + 8的函数值,其中x = l:10 (用数值和符号运算三种方法)符号,

10、置换f = sym(,xa4 -3*x + 8)xa4 -3*x + 8 fx = subs(f/x1:10)fx =columns1through561880252618columns6through1012862388408065429978f = symfxa4 -3*x + 81)f =xa4 -3*x + 8p = sym2poly(f)p =100-38fx = polyval(p,1:10)fx6188025261812862388408065429978逾亘:飙位,多项式求值p=10 0 -3 8p =10 0-38fx = polyval(p,1:10)fx =618802526181286238840806542997813、把(s +1)3(3, + 4屛+ 7)展成多项式的形式。(用数值和符号运算两种方法)毅仏$项式相乘n = conv(conv(1 151 1),conv(1 1,3 4 0 7)n =313212225217f = poly2sym(n/s,)f =3*s6+13*sx5+21*sa4+22*sa3+25*sa2+21*s+

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论