4.2 下面三种表示方法有什么不同的含义? (1)f=3*x^2+5*x+2 (2)f='3*x^2+5*x+2' (3)x=sym('x') f=3*x^2+5*x+2 (1)f=3*x^2+5*x+2
表示在给定x时,将3*x^2+5*x+2的数值运算结果赋值给变量f,如果没有给定x则指示错误信息。 (2)f='3*x^2+5*x+2'
表示将字符串'3*x^2+5*x+2'赋值给字符变量f,没有任何计算含义,因此也不对字符串中的内容做任何分析。 (3)x=sym('x') f=3*x^2+5*x+2
表示x是一个符号变量,因此算式f=3*x^2+5*x+2就具有了符号函数的意义,f也自然成为符号变量了。
4.3 用符号函数法求解方程at+b*t+c=0。 >> r=solve('a*t^2+b*t+c=0','t') r =
[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]
4.4 用符号计算验证三角等式:
sin(?1)cos(?2)-cos(?1)sin(?2) =sin(?1-?2) >> syms phi1 phi2;
>> y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2)) y =
sin(phi1-phi2)
4.5 求矩阵A??2
?a11?a21a12?的行列式值、逆和特征根。 a22?? >> syms a11 a12 a21 a22; >> A=[a11,a12;a21,a22]
>> AD=det(A) % 行列式 >> AI=inv(A) % 逆 >> AE=eig(A) % 特征值 A = [ a11, a12] [ a21, a22] AD =
a11*a22-a12*a21 AI =
[ -a22/(-a11*a22+a12*a21), a12/(-a11*a22+a12*a21)] [ a21/(-a11*a22+a12*a21), -a11/(-a11*a22+a12*a21)] AE =
[ 1/2*a11+1/2*a22+1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)] [ 1/2*a11+1/2*a22-1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)]
4324.6 因式分解:x?5x?5x?5x?6
>> syms x;
>> f=x^4-5*x^3+5*x^2+5*x-6; >> factor(f) ans =
(x-1)*(x-2)*(x-3)*(x+1)
?a4.7 f???ax?e1?x?,用符号微分求df/dx。
?log(x)sin(x)?x2 >> syms a x;
>> f=[a, x^2, 1/x; exp(a*x), log(x), sin(x)]; >> df=diff(f) df =
[ 0, 2*x, -1/x^2] [ a*exp(a*x), 1/x, cos(x)]
??ax2?by?c?04.8 求代数方程组?关于x,y的解。
??x?y?0 >> S=solve('a*x^2+b*y+c=0','b*x+c=0','x','y'); >> disp('S.x=') , disp(S.x) >> disp('S.y=') , disp(S.y) S.x= -c/b S.y=
-c*(a*c+b^2)/b^3
4.9 符号函数绘图法绘制函数x=sin(3t)cos(t),y=sin(3t)sin(t)的图形,t的变化范围为[0,2?]。 >> syms t
>> ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),[0,pi])
4.10 绘制极坐标下sin(3*t)*cos(t)的图形。
>> syms t
>> ezpolar(sin(3*t)*cos(t))
第5章 基本图形处理功能
5.1 绘制曲线y?x3?x?1,x的取值范围为[-5,5]。 >> x=-5:0.2:5; >> y=x.^3+x+1; >> plot(x,y)
5.2 有一组测量数据满足y?e-at,t的变化范围为0~10,用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。
>> t=0:0.5:10; >> y1=exp(-0.1*t); >> y2=exp(-0.2*t); >> y3=exp(-0.5*t);
>> plot(t,y1,'-ob',t,y2,':*r',t,y3,'-.^g')
5.3 在5.1题结果图中添加标题y?e-at,并用箭头线标识出各曲线a的取值。 >> title('\\ity\\rm=e^{-\\itat}')
>> title('\\ity\\rm=e^{-\\itat}','FontSize',12)
>> text(t(6),y1(6),'\\leftarrow\\ita\\rm=0.1','FontSize',11) >> text(t(6),y2(6),'\\leftarrow\\ita\\rm=0.2','FontSize',11) >> text(t(6),y3(6),'\\leftarrow\\ita\\rm=0.5','FontSize',11)
5.4 在5.1题结果图中添加标题y?e-at和图例框。 >> title('\\ity\\rm=e^{-\\itat}','FontSize',12) >> legend('a=0.1','a=0.2','a=0.5')
5.5表中列出了4个观测点的6次测量数据,将数据绘制成为分组形式和堆叠形式的条形图。
观测点1 观测点2 观测点3 观测点4 >> bar(y)
第1次 3 6 9 6 第2次 6 7 7 4 第3次 7 3 2 3 第4次 4 2 5 2 第5次 2 4 8 7 第6次 8 7 4 4 >> y=[3 6 9 6;6 7 7 4;7 3 2 3;4 2 5 2;2 4 8 7;8 7 4 4];
>> bar(y,’stack’)
5.6 x= [66 49 71 56 38],绘制饼图,并将第五个切块分离出来。 >> x=[66 49 71 56 38]; >> L=[0 0 0 0 1]; >> pie(x,L)