资源描述:
《MATLAB基础实验三及求解》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验三、符号运算的应用1、简化"壮+夕+乎+8(比较simple和simplify指令的区别)clear;symsxf=(1/xA3+6/xA2+8)A(1/3)f=(l/xx3+6/x^2+8)人(1/3)simple(f)%simple不返回指定的变量时,给出所有简化的结果,最后结果是最简的形式simplify:((l+6*x+8*xA3)/x^3)"(1/3)radsimp:1/x*(l+6*x+8*x^3)(1/3)combine(trig):((l+6*x+8*x^3)/x^3)人(1/3)facto
2、r:((l+6*x+8*xA3)/xA3)A(1/3)expand:(1/xx3+6/xa2+8)人(1/3)combine:(l/x^3+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/x^3+6/xx2+8)人(1/3)mwcos2sin:(1/xx3+6/xa2+8)a(1/3)ans=(l
3、/x^3+6/xx2+8)人(1/3)=simple(f)%把最简的形式赋给变量B=(1/xa3+6/xa2+8)人(1/3)c=simplify(f)%数学上的恒等式化简c=((l+6*x+8*xA3)/xA3)A(1/3)2、尸sinor,求心空,叟C=©、dxdadx2clear;symsxay=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、计算不定积
4、分.定积分、反常积分K訂;UJoATrK/2cos兀.dxJ=dxJosinx+cosxclear;symsxf1=(xA2+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(-x^2)K=int(f3,0,inf)K=l/2*p
5、ix(l/2)K=vpa(K,4).88654、求级数工飞的和S,以及前十项的部分和Sl・n=“clear;symsnS=symsum(1/nA2,15inf)l/6*pix2S1=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;symsabcxf=a*xA2+b*x+cf=a*x^2+b*x+cX=solve⑴l/2/a*(-b+(b'2-4*a*c
6、)(1/2))1/2/a*((b^2-4*a*c)入(1/2))6、求微分方程y=兀的通解.dsolve(vDy=xjxjans=l/2*x^2+Cl7、求微分方程+x(O)=l,y(O)=O的特解.y=dsolve(,D2y=Dy+x,;y(0)=1,Dy(0)=07x,)-l/2*x^2+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
7、、求下列联立方程的解(用数值和符号运算两种方法)3x+4y-7z-12w=45x-7y+4z+2w=-3Vx+8z・5w=9-6x+5y-2z+lOw=-8法一:数值方法A=[34-7-12;5-742;108-5;-65-210];b=[4;-3;9;-8];x=Ab■1・4841-0・68160・5337-1・2429法二:符号方法clear;symsxyzwf1=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*
8、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+8[wxyz]=solve(f1,f2,f3,f4)w=-2932/2359-3501/2359-1608/2359z=1259/2359X=vpa([x;y;z;w],4)x=■1・484■・6816・5337-1.24310、