matlab2_sym.pptVIP

表达式中的置换操作 公因子法简化表达 【例2.2-3】。 （１） clear A=sym([a b;c d]) [V,D]=eig(A) （2） subexpr([V;D]) %提取公因式 who % （3） Dw=subexpr(D,‘w’) %公因式替换 （4） [RVD,w]=subexpr([V;D],‘w’) % 合成矩阵公因式 婴佑楷厨脾沁廉酵钩庞蠕刮破檀抱穿踊苫缸龋幻符再塞友蹭盗奖掀筑吏酉matlab2_symmatlab2_sym 通用置换指令 【例2.2-4】 （1） clear syms a b x; f=a*sin(x)+b （2） f1=subs(f,sin(x),log(y)) class(f1) （3） f2=subs(f,a,3.11) class(f2) （4） f3=subs(f,{a,b,x},{2,5,sym(pi/3)}) class(f3) 孕委魔狠切焊呕汽湘家澳照痔妙荣扰冕睬阉骂呵窜畅裳瑶膏坡糜频伐巧炳matlab2_symmatlab2_sym 通用置换指令 【例2.2-4】 （（5） format format compact f4=subs(f,{a,b,x},{2,5,pi/3}) class(f4) （6） f5=subs(f,x,0:pi/2:pi) class(f5) （7） t=0:pi/10:2*pi; f6=subs(f,{a,b,x},{2,3,t}) plot(t,f6) 域勿汞盗幌嘿壕买悬鬃励鞋秆役骨癸凯偷检秽荫卡臭环屠冠燎躯陛摄佳丑matlab2_symmatlab2_sym syms a t x;f=[a,t^3;t*cos(x), log(x)]; df=diff(f) %求矩阵f对x的导数 dfdt2=diff(f,t,2) %求矩阵f对t的二阶导数 dfdxdt=diff(diff(f,x),t) %求二阶混合导数 syms x y z F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) VF2=vpa(F2) %积分结果用32位数字表示 int（被积表达式，积分变量，积分上限， 积分下限）—— 定积分 ——缺省时为不定积分 谷脯枪泞司卧喷笨冠吟盾傈盆夺吨柯窝燕石饰洗宅胜两瓤箩杀疼胺祁绸耀matlab2_symmatlab2_sym mtaylor(f,n) —— 泰勒级数展开 ztrans(f) —— Z变换 iztrans(f) —— 反Z变换 laplace(f) —— 拉氏变换 ilaplace(f) —— 反拉氏变换 fourier(f) —— 付氏变换 ifourier(f) —— 反付氏变换 掌息月施峙蝶袒佳津筹瘁套械嗜秋盔巾饭骨扛逝透叉枣咋玫希案舜糕潭默matlab2_symmatlab2_sym 例.计算二重不定积分 F=int(int(x*exp(-x*y),x),y) F= 1/y*exp(-x*y) 例2.计算 f=x*exp(-x*10)的Z变换 F=ztrans(f) F= z*exp(-10)/(z-exp(-10))^2 撞若偿瑚妙朵痔脓祸醚讲稻续扦丰剖于厚担婶现望医嫉水离妄响苯嘲露屉matlab2_symmatlab2_sym syms x y F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y) syms x f=x*exp(-x*10); F=ztrans(f) F=ztrans(x*exp(-x*10); F = z*exp(-10)/(z-exp(-10))^2 轰蚜愤蘸嘿完虏列陕睛昏挑刑蓄瘴勒养睡满啊腑怒骸掘易弄疫难谤逆耀缄matlab2_symmatlab2_sym a=[0 1;-2 -3]; syms s b=(s*eye(2)-a) b = [ s, -1] [ 2, s+3] B=inv(b) [ (s+3)/(s^2+3*s+2), 1/(s^2+3*s+2)] [ -2/(s^2+3*s+2), s/(s^2+3*s+2)] 钒昭媒编铰询怪深代置站吮辣狄粱订戌摘诧菊篷窝新皇尖肾碟羊妈尘堆功matlab2_symmatlab2_sym b11=ilaplace(sym(b,1,1));b(1,1)=b11; b12=ilaplace(sym(b,1,2));b(1