实验6 Matlab数值计算(6)

实验6 Matlab数值计算

根。 6.

当x取矩阵A的每一元素时，求P(x)

的值。其中： 1

1.2 1.4

A

0.75

23.5

0

5

2.5

p1=[1,2,4,0,5]; p2=[1,2]; p3=[1,2,3]; p23=conv(p2,p3); disp('p(x)为：'); px=poly_sum(p1,p23) disp('p(x)的根为：'); r=roots(px)

x=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5]; disp('p(x)的值为：'); p=polyval(px,x) %%%%%

function p_out=poly_sum(p1, p2) % calculate the sum of two polynomials if length(p1)==length(p2) p_out=p1+p2;

elseif length(p1)>length(p2)

p2=[zeros(1,length(p1)-length(p2)), p2]; p_out=p1+p2; else

p1=[zeros(1,length(p2)-length(p1)), p1];

p_out=p1+p2; end

%%%%%%%%%%% >> p(x)为： px =

1 3 8 7 11 p(x)的根为： r =

-1.3840 + 1.8317i -1.3840 - 1.8317i -0.1160 + 1.4400i -0.1160 - 1.4400i p(x)的值为： p =

1.0e+003 *

0.0100 0.0382 0.0125 0.0223 0.0970 0.4122 0.0110 1.2460 0.1644

7. 求函数在指定点的数值导数。

（1）

f(x) sin2x cos2

x,x

,

64,3,

2f=inli

ne('(sin(x)).^2+(cos(x)).^2'); dx1=diff(f([pi/6,0.001+pi/6]))/0.001; dx2=diff(f([pi/4,0.001+pi/4]))/0.001; dx3=diff(f([pi/3,0.001+pi/3]))/0.001;