matlab教程R答案(奇数题).docxVIP

习题1 1 与1500相同的是1.5e3 3 1+eps 5 clear的作用是清除matlab工作空间中保存的变量，clf的作用是清除图形窗，clc的作用是清除指令窗中显示内容。 7 前两个指令能实现目的。 习题2 1说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”对象？ 3/7+0.1, sym(3/7+0.1), vpa(sym(3/7+0.1)) 双精度 符号 符号 符号 3现有方程和，且在这两个方程中，都有的前提。请写出实现以下要求的程序： （1）只求的正实根，而不出现其他形式的根； （2）求的全部根。 （3）在解题结束后，请显示出对符号变量所作的限定性假设。 （提示：正确使用限定假设；清除多余的限定性假设） syms a x positive solve(x^3+(4-a)*x^2+(3-4*a)*x-3*a) syms x clear solve(x^2-a*x+a^2) 5求符号矩阵的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 A=sym([a11 a12 a13;a21 a22 a23;a31 a32 a33]) DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) 9求出的具有64位有效数字的积分值。（提示：int, vpa, ezplot） syms x y=exp(-abs(x))*abs(sin(x)) si=vpa(int(y,-10*pi,1.7*pi),64) 11在区间，画出曲线，并计算。（提示：int, subs） syms t x; f=sin(t)/t; y=int(f,t,0,x) y5=subs(y,x,sym(4.5)) ezplot(y,[0,2*pi]) 13有序列，，（在此，），求这两个序列的卷积。（提示：symsum, subs） syms a b k n x=a^k; h=b^k; w=symsum(subs(h,k,n)*subs(x,k,k-n),n,0,k) %据定义 y1=simple(w) 15求的Fourier变换。（提示：注意限定） syms A t w a=sym(a,positive); f=A*exp(-a*abs(t)); y=fourier(f,t,w) F=simple(y) 17求的Laplace反变换。 syms s t F=(s+3)/(s^3+3*s^2+6*s+4); f=simple(ilaplace(F,s,t)) 19求的Z变换表达式。 syms lambda k T z; f_k=k*exp(-lambda*k*T); F_z=simple(ztrans(f_k,k,z)) 23求微分方程的通解，并绘制任意常数为1时，如图p2-3所示的解曲线图形。（提示：通解中任意常数的替代；构造能完整反映所有解的统一表达式，然后绘图。） syms y x y=dsolve(0.2*y*Dy+0.25*x=0,x) yy=subs(y,C3,1) yy(1)^2==yy(2)^2 syms Y fxy=Y^2-yy(1)^2 clf ezplot(fxy,[-2,2,-2,2]) axis square grid on 25求边值问题的解。 [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g(0)=1) 习题3 1在MATLAB中，先运行指令A=magic(3), B=[1,2,1;3,4,3;5,6,7], C=reshape(1:6,3,2)生成阵列，，，然后根据运行结果回答以下问题： （1）计算A*B, B*A，这两个乘积相同吗？ （2）计算A\B, B/A，左除、右除结果相同吗？ （3）计算B( : ,[1,2]).*C和C.*B( : , [1,2])，这两个乘积相同吗？ （4）计算A\A和A.\A，这两个计算结果相同吗？ （5）计算A\eye(3)和inv(A)，这两个计算结果相同吗？ （提示：根据对计算结果的目测回答问题） （1）不相同，A*B= [41 56 53;53 68 67;41 56 45],B*A=[18 20 22;48 50 52;86 98 86] (2)不相同 (3)相同 (4)不相同 (5)相同 3在MATLAB中，先运行rng default, A=randn(3,3)+1j*randn(3,3)创建一个的复数阵列A，然后根据计算结果回答以下问题： （1）计算A，A., conj(A), conj(A.), 这四个计算结果相同吗？ （2）计算A*A, A.*A., A.*conj(A), 这三个计算结果相同吗？ （提示：根据对计算结果的目测回答