matalab习题答案.docVIP

Collect(s,x)关于表达式s按照x合并同类项 syms x y; d=(x+y)*(x-y); bb=collect(d,x) Expand(s)对表达式s展开 syms x y; d=(x+y)*(x-y); bb=expand(d) Factor(s)因式分解 c=factor(‘x^4-y^4) [a,b]=numden(s)提取分子a和分母b syms x y; u=y/x+2*x/y; [a,b]=numden(u) Simple(s)化简 syms x; u=(x-2)*(x^2+2*x+5+6); a=simple(u) Compose(s)复合函数 syms x,y; f=2*x; g=y^2; bb=compose(f,g) Finverse(f,g)反函数f(g) syms x,y; y=2*x+3; bb=finverse(x,y) Subs(s,new,old)表达式中old被new取代 syms a b x; s=a*sin(b); u=cos(x); bb=subs(s,u,a) Poly(s)矩阵符号s的特征多项式 A=sym(‘[1 2 ; 5 9]’) V=ploy(A) Limit(s,x,x0,’right’)求s函数关于极限x趋于x0的右极限 Syms x; Limit(‘(2x^2+x)/x^2’,x,4,right) Diff(s,x,n)对s求关于x的n阶导数 syms y x; y=sin(x); diff(y,x,2) Int(s,x,a,b)被积函数s关于x的从a到b的积分 syms x y; y=sin(x); b=int(y,x,0,(1/2)*pi) Symsum(s,x,a,b)无穷级数的累加和，从第a项到第b项。 syms i s; s=2*i; b=symsum(s,i,1,5) Taylor(s,n,x,x0) syms x y; y=x^9; t=taylor(y,5,x,0) Solve(s1,s2)代数等式s1,s2的方程组的解 syms x y; [x y]=solve(y=x,x+y=1) Det(A)行列式 A=[1 2; 3 4] ; F=det(A) Inv(A)逆矩阵 A=[15 25 ;25 15] T=inv(A) [a,b]=Egi(A)特征值和特征向量 syms x y; a=[25 35 15;98 21 22 ;32 65 15]; [x y]=eig(a) Rref(A)行最简式 a=[1 2; 3 6]; c=rref(a) Trace(A)迹-主对角线之和 A=[15 16 17;22 33 25;68 2 55] Tr=trace(A) Rank(A)秩 A=[24 31 25;52 21 322 ;25 32 12;24 25 26]; R=rank(A) Cond(A)条件数（cond(A,inf)=norm(A,inf)*norm(inv(A),inf)） A=[12 25;36 25] X=cond(A) Norm(A)范数 A=[1 2 ;3 4] N=norm(a) [a,b]=Size(A)结构 syms x y; a=[25 35 15;98 21 22 ;32 65 15]; [x y]=size(a) Length(A)向量的维数 a=[1 2 3 4]; b=length(a) P=Polyfit(x,y,n),p1=polyval(p,x0)多项式拟合的系数；多项式在x0的值 [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub,x0)线性规划 求P条件数和讨论解的性态的MATLAB程序 function Acp=zpjxpb(A) Acw = cond (A, inf);Ac1= cond (A,1); Ac2= cond (A,2);Acf = cond (A,fro);dA=det(A); if (Acw50)(dA0.1) disp(请注意：AX=b是病态的，A的∞条件数,1条件数,2条件数,条件数和A的行列式的值依次如下：) Acp=[Acw Ac1 Ac2 Acf dA]; else disp( AX=b是良态的，A的∞条件数,1条件数,2条件数,条件数和A的行列式的值依次如下：) Acp=[Acw Ac1 Ac2 Acf dA]; end 程序：Acp =zpjxpb(hilb(7)); Acp,det(hilb(7)) 请问：这个矩阵A是否属于病态？