2.6458 2.8284 3.0000 >> cos(A) ans =
-0.9900 -0.9900 0.2837 -0.4161 -0.6536 0.9602 0.7539 -0.1455 -0.9111 >> pow2(A) ans =
8 8 32 4 16 64 128 256 512 >> 2.^A ans =
8 8 32 4 16 64 128 256 512 b.多项式的运算 >> a a =
1 2 3 4 >> b=[0,1] b =
0 1 加法:
>> c=a+[0 0 b] c =
1 2 3 5 乘法:
>> d=conv(a,b) d =
0 1 2 3 4 除法:
>> [div,rest]=deconv(d,a) div =
0 1 rest =
0 0 0 0 0 微分:
>> polyder(a) ans =
3 4 3 求根: >> roots(a) ans =
-1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i 求值:
>> a=[1,2,3,4] a =
1 2 3 4 >> b=[1,1;1,1] b =
1 1 1 1
>> polyvalm(a,4) % x=4时多项式的值 ans =
112 %与polyval(a,4)结果相同 >> polyval(a,b) %数组运算 ans =
10 10 10 10
>> polyvalm(a,b) %矩阵运算 ans =
15 11 11 15 多项式拟合: >> x=0:pi/20:pi/2; >> y=sin(x);
>> [p,s]=polyfit(x,y,5); >> x1=0:pi/30:pi*2; >> y1=sin(x1);
>>y2=p(1)*x1.^5+p(2)*x1.^4+p(3)*x1.^3+p(4)*x1.^2+p(5)*x1+p(6); >> plot(x1,y1,'b-',x1,y2,'r*')
>> legend('Original curve','Fitted curve') >> axis([0,7,-1.2,4])
上机习题(4)
( 1)在 MATLAB 6.0 的工作空间中用直接输入法建立如下两个矩阵,然后在矩阵编辑器中将矩阵A改为3行3列的矩阵,并将其保存。
?135?
??B?246 ?135???A???? ?789???246? ,
( 2)分别对(1)产生的两个矩阵 作加、减、乘和除(左除,右除)运算,同时运用数组运算法则进行运算,比较二者的计算结果有何异同。
( 3)利用矩阵生成函数建立一个对角线元素全部为 1 的 4 阶单位矩阵。 ( 4)利用矩阵生成函数建立一个 4×4 的随机矩阵。
b.由 sym 命令创建
>> f=sym('a*x^2+b*x+c') f =
a*x^2+b*x+c
>> f1=sym('a*x^2+b*x+c=0') f1 =
a*x^2+b*x+c=0
>> A=sym('[4+x, x^2, x;x^3, 5*x-3, x*a]') A =
[ 4+x, x^2, x] [ x^3, 5*x-3, x*a] c.由 syms 命令创建
>> syms x a b c >> f=a*x^2+b*x+c f =
a*x^2+b*x+c
>> syms x a
>> B=[4+x x^2 x;x^3 5*x-3 x*a] B =
[ 4+x, x^2, x] [ x^3, 5*x-3, x*a] 不能创建符号方程 b.符号矩阵的简化
因式分解 >> syms x
>> factor(x^9-1) ans =
(x-1)*(x^2+x+1)*(x^6+x^3+1)
符号矩阵展开 >> syms x y
>> expand((x+1)^3) ans =
x^3+3*x^2+3*x+1 >> expand(sin(x+y)) ans =
sin(x)*cos(y)+cos(x)*sin(y) 同类项合并
>> syms x y
>> collect(x^2*y+y*x-x^2-2*x) ans =
(y-1)*x^2+(y-2)*x
分式通分 >> syms x y
>> [n,d]=numden(x/y+y/x) n =
x^2+y^2 d = y*x
符号简化 >> syms x
>> simple(cos(x)^2-sin(x)^2) simplify: 2*cos(x)^2-1 radsimp:
cos(x)^2-sin(x)^2 combine(trig): cos(2*x) factor:
(cos(x)-sin(x))*(cos(x)+sin(x)) expand:
cos(x)^2-sin(x)^2 combine: cos(2*x)
上机练习(5)
43(1)设 f (x ) ? x5 ? x ( x ) ? x 4 ? x 2 ? ? 2 x ? x 2 ? 3 , gx ? 1对f(x)和g(x)
进行因式分解。
2 ,(2)运用 simplify 命令简化 x ? ? sin f ? cosx并与采用 simple 命令得出的结果相比较。
(3)计算下列极限 a. lim 2 b.
x?01?cosxxlimx?01?cos2xxsinx(4)建立一个符号表达式 y=sin((a+b)2*x+c),分别以变量 c 进行积分、以变量 x 从pi/2 到 pi 进行积分,
(5)求符号矩阵 F 对变量 x 的一阶微分、对变量 a 的二阶微分