实验1matlab基础与符号计算.doc

2013年上学期MATLAB语言上机实验内容 姓名： 班级： 学号：08 实验1: matlab基础与符号计算 1、指令 clear all A= [2 2 ; 2 2]; A.^2 A^2 ans = 4 4 4 4 ans = 8 8 8 8 2、数组x=10:-2:2; y=1:5; x=10:-2:2; y=1:5; (x5)(y2) y==2|x~=2 ans = 0 0 0 1 1 ans = 1 1 1 1 0 （1）对于，求。（提示：理论结果为） clear syms k x f1=2/(2*k+1)*((x-1)/(x+1))^(2*k+1); s1=symsum(f1,k,0,inf) s1=simple(s1) s1 = (x-1)/(x+1)/((x^2-2*x+1)/(x^2+2*x+1))^(1/2)*log((1+((x^2-2*x+1)/(x^2+2*x+1))^(1/2))/(1-((x^2-2*x+1)/(x^2+2*x+1))^(1/2))) s1 = log(x) （2）求的符号解，并进而用该符号解求，，的准 确值。 （提示：注意subs的使用） clear syms x k f=x^k; s=symsum(f,k,0,inf) f1=simple(subs(f,x,-1/3)) f2=simple(subs(f,x,1/pi)) f3=simple(subs(f,x,3)) s = -1/(x-1) f1 = (-1/3)^k f2 = (5734161139222659/18014398509481984)^k f3 = 3^k （1）通过符号计算求的导数。（2）然后根据此结果，求和。 syms t y; y=abs(sin(t)); dy=diff(y) s1=limit(dy,t,0,left) s2=limit(dy,t,pi/2) dy = abs(1,sin(t))*cos(t) s1 = -1 s2 = 0 指令A= [2 1 ; 1 3]; A=[2 1;1 3]; A.^2-A A^2+A ans = 2 0 0 6 ans = 7 6 6 13 6、数组x=10:-2:2; y=5:-1:1; x=10:-2:2; y=5:-1:1; (x=4)(y=2) y==3|x~=3 ans = 1 1 1 1 0 ans = 1 1 1 1 1 7、计算二重积分。 clear syms x y f=x^2+y^2; f1=int(int(f,y,1,x^2),x,1,2) vf1=vpa(f1) f1 = 1006/105 vf1 = 9.5809523809523809523809523809524 8、在区间，画出曲线，并计算。 syms t x y=int(sin(t)/t,t,0,x) y1=subs(y,x,4.5) ezplot(y,[0,2*pi]) y = sinint(x) y1 = 1.6541 9、设系统的冲激响应为，求该系统在输入，作用下的输出。（提示：运用卷积进行计算） syms t tao ht=exp(-t); ut=cos(t); ub=subs(ut,t,tao)*subs(ht,t,t-tao); yt=simple(simple(int(ub,tao,0,t))) yt =1/2*cos(t)+1/2*sin(t)-1/2/exp(t) 10求的Fourier变换，并画出时的幅频谱。 clear all syms t A w; tao=sym(tao,positive); y=A*((1+t/tao)*(heaviside(t+tao)-heaviside(t))+(1-t/tao)*(heaviside(t)-heaviside(t-tao))); yw=fourier(y,t,w) s=simple(yw) yw2=subs(s,{A,tao},{2,2}); ezplot(abs(yw2)) grid yw = 4*A/tao/w^2*sin(1/2*tao*w)^2 s = 4*A/tao/w^2*sin(1/2*tao